test_digest.c

   1 /*
   2  * Test GOST 34.11 Digest operation
   3  *
   4  * Copyright (C) 2019-2020 Vitaly Chikunov <vt@altlinux.org>. All Rights Reserved.
   5  *
   6  * Contents licensed under the terms of the OpenSSL license
   7  * See https://www.openssl.org/source/license.html for details
   8  */
   9
  10 #include <openssl/engine.h>
  11 #include <openssl/opensslv.h>
  12 #include <openssl/evp.h>
  13 #include <openssl/rand.h>
  14 #include <openssl/err.h>
  15 #include <openssl/asn1.h>
  16 # include <openssl/hmac.h>
  17 #if OPENSSL_VERSION_MAJOR >= 3
  18 # include <openssl/core_names.h>
  19 #endif
  20 #include <openssl/obj_mac.h>
  21 #include <string.h>
  22 #include <stdlib.h>
  23 #if MIPSEL
  24 # include <sys/sysmips.h>
  25 #endif
  26 #ifndef EVP_MD_CTRL_SET_KEY
  27 # include "gost_lcl.h"
  28 #endif
  29
  30 /* Helpers to test OpenSSL API calls. */
  31 #define T(e) ({ if (!(e)) { \
  32                 ERR_print_errors_fp(stderr); \
  33                 OpenSSLDie(__FILE__, __LINE__, #e); \
  34             } \
  35         })
  36 #define TE(e) ({ if (!(e)) { \
  37                 ERR_print_errors_fp(stderr); \
  38                 fprintf(stderr, "Error at %s:%d %s\n", __FILE__, __LINE__, #e); \
  39                 return -1; \
  40             } \
  41         })
  42
  43 #define cRED    "\033[1;31m"
  44 #define cDRED   "\033[0;31m"
  45 #define cGREEN  "\033[1;32m"
  46 #define cDGREEN "\033[0;32m"
  47 #define cBLUE   "\033[1;34m"
  48 #define cDBLUE  "\033[0;34m"
  49 #define cMAGENT "\033[1;35m"
  50 #define cNORM   "\033[m"
  51 #define TEST_ASSERT(e) {if ((test = (e))) \
  52                  printf(cRED "  Test FAILED" cNORM "\n"); \
  53              else \
  54                  printf(cGREEN "  Test passed" cNORM "\n");}
  55
  56 /* To test older APIs. */
  57 #pragma GCC diagnostic ignored "-Wdeprecated-declarations"
  58
  59 /*
  60  * Test keys from both GOST R 34.12-2015 and GOST R 34.13-2015,
  61  * for 128-bit cipher (A.1).
  62  */
  63 static const char K[32] = {
  64     0x88,0x99,0xaa,0xbb,0xcc,0xdd,0xee,0xff,0x00,0x11,0x22,0x33,0x44,0x55,0x66,0x77,
  65     0xfe,0xdc,0xba,0x98,0x76,0x54,0x32,0x10,0x01,0x23,0x45,0x67,0x89,0xab,0xcd,0xef,
  66 };
  67
  68 /*
  69  * Key for 64-bit cipher (A.2).
  70  */
  71 static const char Km[32] = {
  72     0xff,0xee,0xdd,0xcc,0xbb,0xaa,0x99,0x88,0x77,0x66,0x55,0x44,0x33,0x22,0x11,0x00,
  73     0xf0,0xf1,0xf2,0xf3,0xf4,0xf5,0xf6,0xf7,0xf8,0xf9,0xfa,0xfb,0xfc,0xfd,0xfe,0xff,
  74 };
  75
  76 /*
  77  * Plaintext from GOST R 34.13-2015 A.1.
  78  * First 16 bytes is vector (a) from GOST R 34.12-2015 A.1.
  79  */
  80 static const char P[] = {
  81     0x11,0x22,0x33,0x44,0x55,0x66,0x77,0x00,0xff,0xee,0xdd,0xcc,0xbb,0xaa,0x99,0x88,
  82     0x00,0x11,0x22,0x33,0x44,0x55,0x66,0x77,0x88,0x99,0xaa,0xbb,0xcc,0xee,0xff,0x0a,
  83     0x11,0x22,0x33,0x44,0x55,0x66,0x77,0x88,0x99,0xaa,0xbb,0xcc,0xee,0xff,0x0a,0x00,
  84     0x22,0x33,0x44,0x55,0x66,0x77,0x88,0x99,0xaa,0xbb,0xcc,0xee,0xff,0x0a,0x00,0x11,
  85 };
  86
  87 /* Plaintext for 64-bit cipher (A.2) */
  88 static const char Pm[] = {
  89     0x92,0xde,0xf0,0x6b,0x3c,0x13,0x0a,0x59,0xdb,0x54,0xc7,0x04,0xf8,0x18,0x9d,0x20,
  90     0x4a,0x98,0xfb,0x2e,0x67,0xa8,0x02,0x4c,0x89,0x12,0x40,0x9b,0x17,0xb5,0x7e,0x41,
  91 };
  92
  93 /*
  94  * Kuznyechik OMAC1/CMAC test vector from GOST R 34.13-2015 А.1.6
  95  */
  96 static const char MAC_omac[] = { 0x33,0x6f,0x4d,0x29,0x60,0x59,0xfb,0xe3 };
  97
  98 /*
  99  * Magma OMAC1/CMAC test vector from GOST R 34.13-2015 А.2.6
 100  */
 101 static const char MAC_magma_omac[] = { 0x15,0x4e,0x72,0x10 };
 102
 103 /*
 104  * OMAC-ACPKM test vector from R 1323565.1.017-2018 A.4.1
 105  */
 106 static const char P_omac_acpkm1[] = {
 107     0x11,0x22,0x33,0x44,0x55,0x66,0x77,0x00,0xFF,0xEE,0xDD,0xCC,0xBB,0xAA,0x99,0x88,
 108     0x00,0x11,0x22,0x33,0x44,0x55,0x66,0x77,
 109 };
 110
 111 static const char MAC_omac_acpkm1[] = {
 112     0xB5,0x36,0x7F,0x47,0xB6,0x2B,0x99,0x5E,0xEB,0x2A,0x64,0x8C,0x58,0x43,0x14,0x5E,
 113 };
 114
 115 /*
 116  * OMAC-ACPKM test vector from R 1323565.1.017-2018 A.4.2
 117  */
 118 static const char P_omac_acpkm2[] = {
 119     0x11,0x22,0x33,0x44,0x55,0x66,0x77,0x00,0xFF,0xEE,0xDD,0xCC,0xBB,0xAA,0x99,0x88,
 120     0x00,0x11,0x22,0x33,0x44,0x55,0x66,0x77,0x88,0x99,0xAA,0xBB,0xCC,0xEE,0xFF,0x0A,
 121     0x11,0x22,0x33,0x44,0x55,0x66,0x77,0x88,0x99,0xAA,0xBB,0xCC,0xEE,0xFF,0x0A,0x00,
 122     0x22,0x33,0x44,0x55,0x66,0x77,0x88,0x99,0xAA,0xBB,0xCC,0xEE,0xFF,0x0A,0x00,0x11,
 123     0x33,0x44,0x55,0x66,0x77,0x88,0x99,0xAA,0xBB,0xCC,0xEE,0xFF,0x0A,0x00,0x11,0x22,
 124 };
 125
 126 static const char MAC_omac_acpkm2[] = {
 127     0xFB,0xB8,0xDC,0xEE,0x45,0xBE,0xA6,0x7C,0x35,0xF5,0x8C,0x57,0x00,0x89,0x8E,0x5D,
 128 };
 129
 130 /* Some other test vectors. */
 131 static const char etalon_M4[64] = { 0 };
 132
 133 static const char etalon_M5[] = {
 134     0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,
 135     0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x20,0x0a,
 136     0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,
 137     0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x20,0x0a,
 138 };
 139
 140 static const char etalon_M6[] = {
 141     0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,
 142     0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x20,0x0a,
 143     0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,
 144     0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x20,0x0a,
 145     0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,
 146     0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x20,0x0a,
 147     0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,
 148     0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x20,0x0a,
 149     0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,
 150     0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x20,0x0a,
 151     0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,
 152     0x37,0x38,0x39,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x30,0x20,0x0a,
 153 };
 154
 155 static const char etalon_carry[] = {
 156     0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,
 157     0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,
 158     0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,
 159     0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,0xee,
 160     0x16,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,
 161     0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,
 162     0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,
 163     0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x11,0x16,
 164 };
 165
 166 /* This is another carry test. */
 167 static const char ff[] = {
 168     0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,
 169     0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,
 170     0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,
 171     0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff
 172 };
 173
 174 struct hash_testvec {
 175     int nid;               /* OpenSSL algorithm numeric id. */
 176     const char *name;      /* Test name and source. */
 177     const char *plaintext; /* Input (of psize), NULL for synthetic test. */
 178     const char *digest;    /* Expected output (of EVP_MD_size or truncate). */
 179     const char *hmac;      /* Expected output for HMAC tests. */
 180     const char *key;       /* MAC key.*/
 181     int psize;             /* Input (plaintext) size. */
 182     int mdsize;            /* Compare to EVP_MD_size() if non-zero. */
 183     int truncate;          /* Truncated output (digest) size. */
 184     int key_size;          /* MAC key size. */
 185     int block_size;        /* Internal block size. */
 186     int acpkm;             /* The section size N (the number of bits that are
 187                               processed with one section key before this key is
 188                               transformed) (bytes) */
 189     int acpkm_t;           /* Master key (change) frequency T* (bytes) */
 190 };
 191
 192 static const struct hash_testvec testvecs[] = {
 193     { /* Test vectors from standards. */
 194         .nid = NID_id_GostR3411_2012_512,
 195         .name = "M1 from RFC 6986 (10.1.1) and GOST R 34.11-2012 (А.1.1)",
 196         .plaintext =
 197             "012345678901234567890123456789012345678901234567890123456789012",
 198         .psize = 63,
 199         .digest =
 200             "\x1b\x54\xd0\x1a\x4a\xf5\xb9\xd5\xcc\x3d\x86\xd6\x8d\x28\x54\x62"
 201             "\xb1\x9a\xbc\x24\x75\x22\x2f\x35\xc0\x85\x12\x2b\xe4\xba\x1f\xfa"
 202             "\x00\xad\x30\xf8\x76\x7b\x3a\x82\x38\x4c\x65\x74\xf0\x24\xc3\x11"
 203             "\xe2\xa4\x81\x33\x2b\x08\xef\x7f\x41\x79\x78\x91\xc1\x64\x6f\x48",
 204         .mdsize = 512 / 8,
 205         .block_size = 512 / 8,
 206     },
 207     {
 208         .nid = NID_id_GostR3411_2012_256,
 209         .name = "M1 from RFC 6986 (10.1.2) and GOST R 34.11-2012 (А.1.2)",
 210         .plaintext =
 211             "012345678901234567890123456789012345678901234567890123456789012",
 212         .psize = 63,
 213         .digest =
 214             "\x9d\x15\x1e\xef\xd8\x59\x0b\x89\xda\xa6\xba\x6c\xb7\x4a\xf9\x27"
 215             "\x5d\xd0\x51\x02\x6b\xb1\x49\xa4\x52\xfd\x84\xe5\xe5\x7b\x55\x00",
 216         .mdsize = 256 / 8,
 217         .block_size = 512 / 8,
 218     },
 219     {
 220         .nid = NID_id_GostR3411_2012_512,
 221         .name = "M2 from RFC 6986 (10.2.1) and GOST R 34.11-2012 (А.2.1)",
 222         .plaintext =
 223             "\xd1\xe5\x20\xe2\xe5\xf2\xf0\xe8\x2c\x20\xd1\xf2\xf0\xe8\xe1\xee"
 224             "\xe6\xe8\x20\xe2\xed\xf3\xf6\xe8\x2c\x20\xe2\xe5\xfe\xf2\xfa\x20"
 225             "\xf1\x20\xec\xee\xf0\xff\x20\xf1\xf2\xf0\xe5\xeb\xe0\xec\xe8\x20"
 226             "\xed\xe0\x20\xf5\xf0\xe0\xe1\xf0\xfb\xff\x20\xef\xeb\xfa\xea\xfb"
 227             "\x20\xc8\xe3\xee\xf0\xe5\xe2\xfb",
 228         .psize = 72,
 229         .digest =
 230             "\x1e\x88\xe6\x22\x26\xbf\xca\x6f\x99\x94\xf1\xf2\xd5\x15\x69\xe0"
 231             "\xda\xf8\x47\x5a\x3b\x0f\xe6\x1a\x53\x00\xee\xe4\x6d\x96\x13\x76"
 232             "\x03\x5f\xe8\x35\x49\xad\xa2\xb8\x62\x0f\xcd\x7c\x49\x6c\xe5\xb3"
 233             "\x3f\x0c\xb9\xdd\xdc\x2b\x64\x60\x14\x3b\x03\xda\xba\xc9\xfb\x28",
 234     },
 235     {
 236         .nid = NID_id_GostR3411_2012_256,
 237         .name = "M2 from RFC 6986 (10.2.2) and GOST R 34.11-2012 (А.2.2)",
 238         .plaintext =
 239             "\xd1\xe5\x20\xe2\xe5\xf2\xf0\xe8\x2c\x20\xd1\xf2\xf0\xe8\xe1\xee"
 240             "\xe6\xe8\x20\xe2\xed\xf3\xf6\xe8\x2c\x20\xe2\xe5\xfe\xf2\xfa\x20"
 241             "\xf1\x20\xec\xee\xf0\xff\x20\xf1\xf2\xf0\xe5\xeb\xe0\xec\xe8\x20"
 242             "\xed\xe0\x20\xf5\xf0\xe0\xe1\xf0\xfb\xff\x20\xef\xeb\xfa\xea\xfb"
 243             "\x20\xc8\xe3\xee\xf0\xe5\xe2\xfb",
 244         .psize = 72,
 245         .digest =
 246             "\x9d\xd2\xfe\x4e\x90\x40\x9e\x5d\xa8\x7f\x53\x97\x6d\x74\x05\xb0"
 247             "\xc0\xca\xc6\x28\xfc\x66\x9a\x74\x1d\x50\x06\x3c\x55\x7e\x8f\x50",
 248     },
 249     /* OMAC tests */
 250     {
 251         .nid = NID_grasshopper_mac,
 252         .name = "P from GOST R 34.13-2015 (А.1.6)",
 253         .plaintext = P,
 254         .psize = sizeof(P),
 255         .key = K,
 256         .key_size = sizeof(K),
 257         .digest = MAC_omac,
 258         .mdsize = 128 / 8,
 259         .truncate = sizeof(MAC_omac),
 260     },
 261     {
 262         .nid = NID_magma_mac,
 263         .name = "P from GOST R 34.13-2015 (А.2.6)",
 264         .plaintext = Pm,
 265         .psize = sizeof(Pm),
 266         .key = Km,
 267         .key_size = sizeof(Km),
 268         .digest = MAC_magma_omac,
 269         .mdsize = 64 / 8,
 270         .truncate = sizeof(MAC_magma_omac),
 271     },
 272     {
 273         .nid = NID_id_tc26_cipher_gostr3412_2015_kuznyechik_ctracpkm_omac,
 274         .name = "M from R 1323565.1.017-2018 (A.4.1)",
 275         .plaintext = P_omac_acpkm1,
 276         .psize = sizeof(P_omac_acpkm1),
 277         .key = K,
 278         .key_size = sizeof(K),
 279         .acpkm = 32,
 280         .acpkm_t = 768 / 8,
 281         .digest = MAC_omac_acpkm1,
 282         .mdsize = sizeof(MAC_omac_acpkm1),
 283     },
 284     {
 285         .nid = NID_id_tc26_cipher_gostr3412_2015_kuznyechik_ctracpkm_omac,
 286         .name = "M from R 1323565.1.017-2018 (A.4.2)",
 287         .plaintext = P_omac_acpkm2,
 288         .psize = sizeof(P_omac_acpkm2),
 289         .key = K,
 290         .key_size = sizeof(K),
 291         .acpkm = 32,
 292         .acpkm_t = 768 / 8,
 293         .digest = MAC_omac_acpkm2,
 294         .mdsize = sizeof(MAC_omac_acpkm2),
 295     },
 296     { /* HMAC tests. */
 297         .nid = NID_id_GostR3411_2012_512,
 298         .name = "HMAC_GOSTR3411_2012_512 from RFC 7836 (B) and R 50.1.113-2016 (A)",
 299         .plaintext =
 300             "\x01\x26\xbd\xb8\x78\x00\xaf\x21\x43\x41\x45\x65\x63\x78\x01\x00",
 301         .psize = 16,
 302         .key =
 303             "\x00\x01\x02\x03\x04\x05\x06\x07\x08\x09\x0a\x0b\x0c\x0d\x0e\x0f"
 304             "\x10\x11\x12\x13\x14\x15\x16\x17\x18\x19\x1a\x1b\x1c\x1d\x1e\x1f",
 305         .key_size = 32,
 306         .hmac =
 307             "\xa5\x9b\xab\x22\xec\xae\x19\xc6\x5f\xbd\xe6\xe5\xf4\xe9\xf5\xd8"
 308             "\x54\x9d\x31\xf0\x37\xf9\xdf\x9b\x90\x55\x00\xe1\x71\x92\x3a\x77"
 309             "\x3d\x5f\x15\x30\xf2\xed\x7e\x96\x4c\xb2\xee\xdc\x29\xe9\xad\x2f"
 310             "\x3a\xfe\x93\xb2\x81\x4f\x79\xf5\x00\x0f\xfc\x03\x66\xc2\x51\xe6",
 311     },
 312     {
 313         .nid = NID_id_GostR3411_2012_256,
 314         .name = "HMAC_GOSTR3411_2012_256 from RFC 7836 (B) and R 50.1.113-2016 (A)",
 315         .plaintext =
 316             "\x01\x26\xbd\xb8\x78\x00\xaf\x21\x43\x41\x45\x65\x63\x78\x01\x00",
 317         .psize = 16,
 318         .key =
 319             "\x00\x01\x02\x03\x04\x05\x06\x07\x08\x09\x0a\x0b\x0c\x0d\x0e\x0f"
 320             "\x10\x11\x12\x13\x14\x15\x16\x17\x18\x19\x1a\x1b\x1c\x1d\x1e\x1f",
 321         .key_size = 32,
 322         .hmac =
 323             "\xa1\xaa\x5f\x7d\xe4\x02\xd7\xb3\xd3\x23\xf2\x99\x1c\x8d\x45\x34"
 324             "\x01\x31\x37\x01\x0a\x83\x75\x4f\xd0\xaf\x6d\x7c\xd4\x92\x2e\xd9",
 325     },
 326     /* Other KATs. */
 327     {
 328         .nid = NID_id_GostR3411_2012_512,
 329         .name = "Zero length string (M3)",
 330         .plaintext = "",
 331         .psize = 0,
 332         .digest =
 333             "\x8e\x94\x5d\xa2\x09\xaa\x86\x9f\x04\x55\x92\x85\x29\xbc\xae\x46"
 334             "\x79\xe9\x87\x3a\xb7\x07\xb5\x53\x15\xf5\x6c\xeb\x98\xbe\xf0\xa7"
 335             "\x36\x2f\x71\x55\x28\x35\x6e\xe8\x3c\xda\x5f\x2a\xac\x4c\x6a\xd2"
 336             "\xba\x3a\x71\x5c\x1b\xcd\x81\xcb\x8e\x9f\x90\xbf\x4c\x1c\x1a\x8a",
 337         .mdsize = 512 / 8,
 338     },
 339     {
 340         .nid = NID_id_GostR3411_2012_256,
 341         .name = "Zero length string (M3)",
 342         .plaintext = "",
 343         .psize = 0,
 344         .digest =
 345             "\x3f\x53\x9a\x21\x3e\x97\xc8\x02\xcc\x22\x9d\x47\x4c\x6a\xa3\x2a"
 346             "\x82\x5a\x36\x0b\x2a\x93\x3a\x94\x9f\xd9\x25\x20\x8d\x9c\xe1\xbb",
 347         .mdsize = 256 / 8,
 348     },
 349     { /* M4 */
 350         .nid = NID_id_GostR3411_2012_512,
 351         .name = "64 bytes of zero (M4)",
 352         .plaintext = etalon_M4,
 353         .psize = sizeof(etalon_M4),
 354         .digest =
 355             "\xb0\xfd\x29\xac\x1b\x0d\xf4\x41\x76\x9f\xf3\xfd\xb8\xdc\x56\x4d"
 356             "\xf6\x77\x21\xd6\xac\x06\xfb\x28\xce\xff\xb7\xbb\xaa\x79\x48\xc6"
 357             "\xc0\x14\xac\x99\x92\x35\xb5\x8c\xb2\x6f\xb6\x0f\xb1\x12\xa1\x45"
 358             "\xd7\xb4\xad\xe9\xae\x56\x6b\xf2\x61\x14\x02\xc5\x52\xd2\x0d\xb7"
 359     },
 360     {
 361         .nid = NID_id_GostR3411_2012_256,
 362         .name = "64 bytes of zero (M4)",
 363         .plaintext = etalon_M4,
 364         .psize = sizeof(etalon_M4),
 365         .digest =
 366             "\xdf\x1f\xda\x9c\xe8\x31\x91\x39\x05\x37\x35\x80\x31\xdb\x2e\xca"
 367             "\xa6\xaa\x54\xcd\x0e\xda\x24\x1d\xc1\x07\x10\x5e\x13\x63\x6b\x95"
 368     },
 369     { /* M5 */
 370         .nid = NID_id_GostR3411_2012_512,
 371         .name = "64 bytes of (M5)",
 372         .plaintext = etalon_M5,
 373         .psize = sizeof(etalon_M5),
 374         .digest =
 375             "\x36\x3b\x44\x9e\xc8\x1a\xe4\x0b\x3a\x40\x7b\x12\x5c\x3b\x1c\x2b"
 376             "\x76\x8b\x50\x49\x6b\xcb\x5f\x69\x0b\x89\xe9\x00\x7b\x06\xe4\x08"
 377             "\x41\x82\xed\x45\xd4\x07\x2a\x67\xfe\xc9\xd3\x42\x1d\xab\x01\x3d"
 378             "\xa2\xaa\xbc\x1d\x65\x28\xe8\xe7\x75\xae\xc7\xb3\x45\x7a\xc6\x75"
 379     },
 380     {
 381         .nid = NID_id_GostR3411_2012_256,
 382         .name = "64 bytes of (M5)",
 383         .plaintext = etalon_M5,
 384         .psize = sizeof(etalon_M5),
 385         .digest =
 386             "\xf0\xa5\x57\xf6\xa0\x4a\x90\xab\x18\x11\xc1\xb6\xe9\xb0\x78\xe4"
 387             "\x16\x3b\x74\x03\x7c\x6c\xf5\x9f\x52\x44\x4a\x37\xf4\x8e\x11\xd8"
 388     },
 389     { /* M6 */
 390         .nid = NID_id_GostR3411_2012_512,
 391         .name = "(M6)",
 392         .plaintext = etalon_M6,
 393         .psize = sizeof(etalon_M6),
 394         .digest =
 395             "\x87\x81\xdf\xc8\x1d\x2d\xb6\xa4\x1d\x18\x57\xf3\x23\x0b\x3f\xfe"
 396             "\x2b\xda\x57\x42\x73\xea\x19\x47\x18\x9a\xaa\x54\x68\x47\x0d\xf1"
 397             "\xc4\xb3\x74\xb1\xa2\xb5\x6e\x59\xd1\x1d\x21\x3f\xea\x57\xe3\x51"
 398             "\x45\x43\xb0\xce\xd9\xb2\x0e\x55\x3a\xe6\x64\x25\xec\x90\x9c\xfd"
 399     },
 400     {
 401         .nid = NID_id_GostR3411_2012_256,
 402         .name = "(M6)",
 403         .plaintext = etalon_M6,
 404         .psize = sizeof(etalon_M6),
 405         .digest =
 406             "\x2f\x4f\x65\x1f\xe8\x8f\xea\x46\xec\x6f\x22\x23\x72\x8d\x8d\xff"
 407             "\x39\x68\x89\x35\x58\xef\x00\xa3\x10\xc2\x3e\x7d\x19\x23\xba\x0c"
 408     },
 409     { /* Carry */
 410         .nid = NID_id_GostR3411_2012_512,
 411         .name = "(carry)",
 412         .plaintext = etalon_carry,
 413         .psize = sizeof(etalon_carry),
 414         .digest =
 415             "\x8b\x06\xf4\x1e\x59\x90\x7d\x96\x36\xe8\x92\xca\xf5\x94\x2f\xcd"
 416             "\xfb\x71\xfa\x31\x16\x9a\x5e\x70\xf0\xed\xb8\x73\x66\x4d\xf4\x1c"
 417             "\x2c\xce\x6e\x06\xdc\x67\x55\xd1\x5a\x61\xcd\xeb\x92\xbd\x60\x7c"
 418             "\xc4\xaa\xca\x67\x32\xbf\x35\x68\xa2\x3a\x21\x0d\xd5\x20\xfd\x41"
 419     },
 420     {
 421         .nid = NID_id_GostR3411_2012_256,
 422         .name = "(carry)",
 423         .plaintext = etalon_carry,
 424         .psize = sizeof(etalon_carry),
 425         .digest =
 426             "\x81\xbb\x63\x2f\xa3\x1f\xcc\x38\xb4\xc3\x79\xa6\x62\xdb\xc5\x8b"
 427             "\x9b\xed\x83\xf5\x0d\x3a\x1b\x2c\xe7\x27\x1a\xb0\x2d\x25\xba\xbb"
 428     },
 429     { /* ff (Better carry test). */
 430         .nid = NID_id_GostR3411_2012_512,
 431         .name = "64 bytes of FF",
 432         .plaintext = ff,
 433         .psize = sizeof(ff),
 434         .digest =
 435             "\x41\x62\x9d\xe6\x77\xd7\xe8\x09\x0c\x3c\xd7\x0a\xff\xe3\x30\x0d"
 436             "\x1e\x1c\xfb\xa2\xdb\x97\x94\x5e\xc3\x7f\xeb\x4e\x13\x75\xbc\x02"
 437             "\xa5\x3f\x00\x37\x0b\x7d\x71\x5b\x07\xf3\x7f\x93\xca\xc8\x44\xef"
 438             "\xad\xbf\xd1\xb8\x5f\x9d\xda\xe3\xde\x96\x56\xc0\xe9\x5a\xff\xc7"
 439     },
 440     {
 441         .nid = NID_id_GostR3411_2012_256,
 442         .name = "64 bytes of FF",
 443         .plaintext = ff,
 444         .psize = sizeof(ff),
 445         .digest =
 446             "\x96\x4a\x5a\xb6\x02\x86\xf1\x06\x28\x87\x43\xe2\xfe\x1a\x42\x2d"
 447             "\x16\x08\x98\xca\x1b\xd5\x35\xe8\x31\xaa\x50\x0c\xfe\x34\xd7\xe8"
 448     },
 449     {
 450         .nid = NID_id_GostR3411_94,
 451         .name = "64 bytes of FF",
 452         .plaintext = ff,
 453         .psize = sizeof(ff),
 454         .digest =
 455             "\x58\x50\x4d\x26\xb3\x67\x7e\x75\x6b\xa3\xf4\xa9\xfd\x2f\x14\xb3"
 456             "\xba\x54\x57\x06\x6a\x4a\xa1\xd7\x00\x65\x9b\x90\xdc\xdd\xd3\xc6"
 457     },
 458     /* Synthetic tests. */
 459     {
 460         .nid = NID_id_GostR3411_2012_256,
 461         .name = "streebog256 synthetic test",
 462         .mdsize = 32,
 463         .block_size = 64,
 464         .digest =
 465             "\xa2\xf3\x6d\x9c\x42\xa1\x1e\xad\xe3\xc1\xfe\x99\xf9\x99\xc3\x84"
 466             "\xe7\x98\xae\x24\x50\x75\x73\xd7\xfc\x99\x81\xa0\x45\x85\x41\xf6"
 467     },
 468     {
 469         .nid = NID_id_GostR3411_2012_512,
 470         .name = "streebog512 synthetic test",
 471         .mdsize = 64,
 472         .block_size = 64,
 473         .digest =
 474             "\x1d\x14\x4d\xd8\xb8\x27\xfb\x55\x1a\x5a\x7d\x03\xbb\xdb\xfa\xcb"
 475             "\x43\x6b\x5b\xc5\x77\x59\xfd\x5f\xf2\x3b\x8e\xf9\xc4\xdd\x6f\x79"
 476             "\x45\xd8\x16\x59\x9e\xaa\xbc\xf2\xb1\x4f\xd0\xe4\xf6\xad\x46\x60"
 477             "\x90\x89\xf7\x2f\x93\xd8\x85\x0c\xb0\x43\xff\x5a\xb6\xe3\x69\xbd"
 478     },
 479     { 0 }
 480 };
 481
 482 static void hexdump(const void *ptr, size_t len)
 483 {
 484     const unsigned char *p = ptr;
 485     size_t i, j;
 486
 487     for (i = 0; i < len; i += j) {
 488         for (j = 0; j < 16 && i + j < len; j++)
 489             printf("%s%02x", j? "" : " ", p[i + j]);
 490     }
 491     printf("\n");
 492 }
 493
 494 static int do_hmac_old(int iter, const EVP_MD *type, const char *plaintext,
 495     const struct hash_testvec *t)
 496 {
 497     unsigned int len;
 498     unsigned char md[EVP_MAX_MD_SIZE];
 499     if (!iter)
 500         printf("[HMAC] ");
 501
 502     HMAC_CTX *ctx;
 503     T(ctx = HMAC_CTX_new());
 504     T(HMAC_Init_ex(ctx, t->key, t->key_size, type, NULL));
 505     T(HMAC_Update(ctx, (const unsigned char *)plaintext, t->psize));
 506     T(HMAC_Final(ctx, md, &len));
 507     HMAC_CTX_free(ctx);
 508
 509     if (t->mdsize)
 510         T(len == t->mdsize);
 511     if (memcmp(md, t->hmac, len) != 0) {
 512         printf(cRED "hmac mismatch" cNORM "\n");
 513         hexdump(t->hmac, len);
 514         hexdump(md, len);
 515         return 1;
 516     }
 517     return 0;
 518 }
 519
 520 #if OPENSSL_VERSION_MAJOR >= 3
 521 static int do_hmac_prov(int iter, const EVP_MD *type, const char *plaintext,
 522     const struct hash_testvec *t)
 523 {
 524     size_t len;
 525     unsigned char md[EVP_MAX_MD_SIZE];
 526     if (!iter)
 527         printf("[EVP_MAC] ");
 528
 529     EVP_MAC *hmac;
 530     T(hmac = EVP_MAC_fetch(NULL, "HMAC", NULL));
 531     EVP_MAC_CTX *ctx;
 532     T(ctx = EVP_MAC_CTX_new(hmac));
 533     OSSL_PARAM params[] = {
 534         OSSL_PARAM_utf8_string(OSSL_MAC_PARAM_DIGEST,
 535             (char *)EVP_MD_name(type), 0),
 536         OSSL_PARAM_END
 537     };
 538     T(EVP_MAC_init(ctx, (const unsigned char *)t->key, t->key_size, params));
 539     T(EVP_MAC_update(ctx, (unsigned char *)plaintext, t->psize));
 540     T(EVP_MAC_final(ctx, md, &len, EVP_MAX_MD_SIZE));
 541     EVP_MAC_CTX_free(ctx);
 542     EVP_MAC_free(hmac);
 543
 544     if (t->mdsize)
 545         T(len == t->mdsize);
 546     if (memcmp(md, t->hmac, len) != 0) {
 547         printf(cRED "hmac mismatch" cNORM "\n");
 548         hexdump(t->hmac, len);
 549         hexdump(md, len);
 550         return 1;
 551     }
 552     return 0;
 553 }
 554 #endif
 555
 556 static int do_hmac(int iter, const EVP_MD *type, const char *plaintext,
 557     const struct hash_testvec *t)
 558 {
 559     int ret;
 560
 561     /* Test old (deprecated) and (too) new APIs. */
 562     ret = do_hmac_old(iter, type, plaintext, t);
 563 #if OPENSSL_VERSION_MAJOR >= 3
 564     ret |= do_hmac_prov(iter, type, plaintext, t);
 565 #endif
 566
 567     return ret;
 568 }
 569
 570 /*
 571  * If we have OMAC1/CMAC test vector,
 572  * use CMAC provider to test it.
 573  */
 574 static int do_cmac_prov(int iter, const char *plaintext,
 575     const struct hash_testvec *t)
 576 {
 577 #if OPENSSL_VERSION_MAJOR >= 3
 578     char *mdname = NULL;
 579     /*
 580      * CMAC needs CBC.
 581      * Convert 'mac' digest to the underlying CBC cipher.
 582      */
 583     switch (t->nid) {
 584         case NID_grasshopper_mac:
 585             mdname = "kuznyechik-cbc";
 586             break;
 587         case NID_magma_mac:
 588             mdname = "magma-cbc";
 589             break;
 590         default:
 591             return 0;
 592     }
 593
 594     if (!iter)
 595         printf("[CMAC(%s)] ", mdname);
 596
 597     size_t len;
 598     unsigned char md[EVP_MAX_MD_SIZE];
 599     size_t mdsize = t->mdsize;
 600     if (t->truncate)
 601         mdsize = t->truncate;
 602
 603     EVP_MAC *cmac;
 604     T(cmac = EVP_MAC_fetch(NULL, "CMAC", NULL));
 605     EVP_MAC_CTX *ctx;
 606     T(ctx = EVP_MAC_CTX_new(cmac));
 607     OSSL_PARAM params[] = {
 608         OSSL_PARAM_utf8_string(OSSL_MAC_PARAM_CIPHER, mdname, 0),
 609         OSSL_PARAM_END
 610     };
 611     T(EVP_MAC_CTX_set_params(ctx, params));
 612     T(EVP_MAC_init(ctx, (const unsigned char *)t->key, t->key_size, params));
 613     T(EVP_MAC_update(ctx, (unsigned char *)plaintext, t->psize));
 614     T(EVP_MAC_final(ctx, md, &len, EVP_MAX_MD_SIZE));
 615     EVP_MAC_CTX_free(ctx);
 616     EVP_MAC_free(cmac);
 617
 618     /* CMAC provider will not respect outsize, and will output full block.
 619      * So, just compare until what we need. */
 620     T(mdsize <= len);
 621     if (memcmp(md, t->digest, mdsize) != 0) {
 622         printf(cRED "cmac mismatch" cNORM "\n");
 623         hexdump(t->digest, mdsize);
 624         hexdump(md, len);
 625         return 1;
 626     }
 627 #endif
 628     return 0;
 629 }
 630
 631 static int do_digest(int iter, const EVP_MD *type, const char *plaintext,
 632     const struct hash_testvec *t)
 633 {
 634     if (!iter)
 635         printf("[MD] ");
 636     if (t->mdsize)
 637         T(EVP_MD_size(type) == t->mdsize);
 638     size_t mdsize;
 639     if (t->truncate)
 640         mdsize = t->truncate;
 641     else
 642         mdsize = EVP_MD_size(type);
 643
 644     if (t->block_size)
 645         T(EVP_MD_block_size(type) == t->block_size);
 646     EVP_MD_CTX *ctx;
 647     T(ctx = EVP_MD_CTX_new());
 648     T(EVP_MD_CTX_init(ctx));
 649     T(EVP_DigestInit_ex(ctx, type, NULL));
 650     if (t->key)
 651         T(EVP_MD_CTX_ctrl(ctx, EVP_MD_CTRL_SET_KEY, t->key_size,
 652                 (void *)t->key));
 653     if (t->acpkm)
 654         T(EVP_MD_CTX_ctrl(ctx,
 655                 EVP_CTRL_KEY_MESH, t->acpkm,
 656                 t->acpkm_t? (void *)&t->acpkm_t : NULL));
 657     T(EVP_DigestUpdate(ctx, plaintext, t->psize));
 658
 659     unsigned int len;
 660     unsigned char md[EVP_MAX_MD_SIZE];
 661     if (EVP_MD_flags(EVP_MD_CTX_md(ctx)) & EVP_MD_FLAG_XOF) {
 662         T(EVP_DigestFinalXOF(ctx, md, mdsize));
 663         len = mdsize;
 664     } else {
 665         T(EVP_MD_CTX_size(ctx) == mdsize);
 666         T(EVP_DigestFinal_ex(ctx, md, &len));
 667     }
 668
 669     EVP_MD_CTX_free(ctx);
 670     T(len == mdsize);
 671     if (memcmp(md, t->digest, mdsize) != 0) {
 672         printf(cRED "digest mismatch" cNORM "\n");
 673         hexdump(t->digest, mdsize);
 674         hexdump(md, mdsize);
 675         return 1;
 676     }
 677
 678     return 0;
 679 }
 680
 681 static int do_test(const struct hash_testvec *tv)
 682 {
 683     int ret = 0;
 684
 685     const EVP_MD *type;
 686     T(type = EVP_get_digestbynid(tv->nid));
 687     const char *name = EVP_MD_name(type);
 688     printf(cBLUE "Test %s: %s: " cNORM, name, tv->name);
 689
 690     /* Test alignment problems. */
 691     int shifts = 32;
 692     int i;
 693     char *buf;
 694     T(buf = OPENSSL_malloc(tv->psize + shifts));
 695     for (i = 0; i < shifts; i++) {
 696         memcpy(buf + i, tv->plaintext, tv->psize);
 697         if (tv->hmac)
 698             ret |= do_hmac(i, type, buf + i, tv);
 699         else
 700             ret |= do_digest(i, type, buf + i, tv);
 701         /* Test CMAC provider. */
 702         ret |= do_cmac_prov(i, buf + i, tv);
 703         /* No need to continue loop on failure. */
 704         if (ret)
 705             break;
 706     }
 707     OPENSSL_free(buf);
 708
 709     if (!ret)
 710         printf(cGREEN "success" cNORM "\n");
 711     else
 712         printf(cRED "fail" cNORM "\n");
 713     return ret;
 714 }
 715
 716 #define SUPER_SIZE 256
 717 /*
 718  * For 256-byte buffer filled with 256 bytes from 0 to 255;
 719  * Digest them 256 times from the buffer end with lengths from 0 to 256,
 720  * and from beginning of the buffer with lengths from 0 to 256;
 721  * Each produced digest is digested again into final sum.
 722  */
 723 static int do_synthetic_once(const struct hash_testvec *tv, unsigned int shifts)
 724 {
 725     unsigned char *ibuf, *md;
 726     T(ibuf = OPENSSL_zalloc(SUPER_SIZE + shifts));
 727
 728     /* fill with pattern */
 729     unsigned int len;
 730     for (len = 0; len < SUPER_SIZE; len++)
 731             ibuf[shifts + len] = len & 0xff;
 732
 733     const EVP_MD *mdtype;
 734     T(mdtype = EVP_get_digestbynid(tv->nid));
 735     OPENSSL_assert(tv->nid == EVP_MD_type(mdtype));
 736     EVP_MD_CTX *ctx, *ctx2;
 737     T(ctx  = EVP_MD_CTX_new());
 738     T(ctx2 = EVP_MD_CTX_new());
 739     T(EVP_DigestInit(ctx2, mdtype));
 740     OPENSSL_assert(tv->nid == EVP_MD_CTX_type(ctx2));
 741     OPENSSL_assert(EVP_MD_block_size(mdtype) == tv->block_size);
 742     OPENSSL_assert(EVP_MD_CTX_size(ctx2) == tv->mdsize);
 743     OPENSSL_assert(EVP_MD_CTX_block_size(ctx2) == tv->block_size);
 744
 745     const unsigned int mdlen = EVP_MD_size(mdtype);
 746     OPENSSL_assert(mdlen == tv->mdsize);
 747     T(md = OPENSSL_zalloc(mdlen + shifts));
 748     md += shifts; /* test for output digest alignment problems */
 749
 750     /* digest cycles */
 751     for (len = 0; len < SUPER_SIZE; len++) {
 752         /* for each len digest len bytes from the end of buf */
 753         T(EVP_DigestInit(ctx, mdtype));
 754         T(EVP_DigestUpdate(ctx, ibuf + shifts + SUPER_SIZE - len, len));
 755         T(EVP_DigestFinal(ctx, md, NULL));
 756         T(EVP_DigestUpdate(ctx2, md, mdlen));
 757     }
 758
 759     for (len = 0; len < SUPER_SIZE; len++) {
 760         /* for each len digest len bytes from the beginning of buf */
 761         T(EVP_DigestInit(ctx, mdtype));
 762         T(EVP_DigestUpdate(ctx, ibuf + shifts, len));
 763         T(EVP_DigestFinal(ctx, md, NULL));
 764         T(EVP_DigestUpdate(ctx2, md, mdlen));
 765     }
 766
 767     OPENSSL_free(ibuf);
 768     EVP_MD_CTX_free(ctx);
 769
 770     T(EVP_DigestFinal(ctx2, md, &len));
 771     EVP_MD_CTX_free(ctx2);
 772
 773     if (len != mdlen) {
 774         printf(cRED "digest output len mismatch %u != %u (expected)" cNORM "\n",
 775             len, mdlen);
 776         goto err;
 777     }
 778
 779     if (memcmp(md, tv->digest, mdlen) != 0) {
 780         printf(cRED "digest mismatch" cNORM "\n");
 781
 782         unsigned int i;
 783         printf("  Expected value is: ");
 784         for (i = 0; i < mdlen; i++)
 785             printf("\\x%02x", md[i]);
 786         printf("\n");
 787         goto err;
 788     }
 789
 790     OPENSSL_free(md - shifts);
 791     return 0;
 792 err:
 793     OPENSSL_free(md - shifts);
 794     return 1;
 795 }
 796
 797 /* do different block sizes and different memory offsets */
 798 static int do_synthetic_test(const struct hash_testvec *tv)
 799 {
 800     int ret = 0;
 801
 802     printf(cBLUE "Test %s: " cNORM, tv->name);
 803     fflush(stdout);
 804
 805     unsigned int shifts;
 806     for (shifts = 0; shifts < 16 && !ret; shifts++)
 807         ret |= do_synthetic_once(tv, shifts);
 808
 809     if (!ret)
 810         printf(cGREEN "success" cNORM "\n");
 811     else
 812         printf(cRED "fail" cNORM "\n");
 813     return 0;
 814 }
 815
 816 int main(int argc, char **argv)
 817 {
 818     int ret = 0;
 819
 820 #if MIPSEL
 821     /* Trigger SIGBUS for unaligned access. */
 822     sysmips(MIPS_FIXADE, 0);
 823 #endif
 824     OPENSSL_add_all_algorithms_conf();
 825
 826     const struct hash_testvec *tv;
 827     for (tv = testvecs; tv->nid; tv++) {
 828         if (tv->plaintext)
 829             ret |= do_test(tv);
 830         else
 831             ret |= do_synthetic_test(tv);
 832     }
 833
 834     ENGINE *eng;
 835     ENGINE_DIGESTS_PTR fn_c;
 836     T(eng = ENGINE_by_id("gost"));
 837     T(fn_c = ENGINE_get_digests(eng));
 838     const int *nids;
 839     int n, k;
 840     n = fn_c(eng, NULL, &nids, 0);
 841     for (k = 0; k < n; ++k) {
 842         for (tv = testvecs; tv->nid; tv++)
 843             if (tv->nid == nids[k])
 844                 break;
 845         if (!tv->nid)
 846             printf(cMAGENT "Digest %s is untested!" cNORM "\n", OBJ_nid2sn(nids[k]));
 847     }
 848     ENGINE_free(eng);
 849
 850     if (ret)
 851         printf(cDRED "= Some tests FAILED!" cNORM "\n");
 852     else
 853         printf(cDGREEN "= All tests passed!" cNORM "\n");
 854     return ret;
 855 }