math/mpreduce.h: Missing include files.
[u/mdw/catacomb] / math / ec-test.c
CommitLineData
bc985cef 1/* -*-c-*-
2 *
bc985cef 3 * Code for testing elliptic-curve stuff
4 *
5 * (c) 2004 Straylight/Edgeware
6 */
7
45c0fd36 8/*----- Licensing notice --------------------------------------------------*
bc985cef 9 *
10 * This file is part of Catacomb.
11 *
12 * Catacomb is free software; you can redistribute it and/or modify
13 * it under the terms of the GNU Library General Public License as
14 * published by the Free Software Foundation; either version 2 of the
15 * License, or (at your option) any later version.
45c0fd36 16 *
bc985cef 17 * Catacomb is distributed in the hope that it will be useful,
18 * but WITHOUT ANY WARRANTY; without even the implied warranty of
19 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20 * GNU Library General Public License for more details.
45c0fd36 21 *
bc985cef 22 * You should have received a copy of the GNU Library General Public
23 * License along with Catacomb; if not, write to the Free
24 * Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
25 * MA 02111-1307, USA.
26 */
27
bc985cef 28/*----- Header files ------------------------------------------------------*/
29
30#include <assert.h>
31#include <ctype.h>
32#include <stdio.h>
33#include <stdlib.h>
34#include <string.h>
35
36#include <mLib/alloc.h>
37#include <mLib/testrig.h>
38#include <mLib/sub.h>
39
40#include "ec.h"
41#include "ec-test.h"
42
43/*----- Cardboard cut-out elliptic curve ----------------------------------*/
44
45typedef struct ecctx {
46 ec_curve c;
47 unsigned long magic;
48 char *name;
49 ec_curve *real;
50} ecctx;
51
52#define MAGIC 0x3a1f0b07
53
54static void ecDESTROY(ec_curve *cc)
55{
56 ecctx *c = (ecctx *)cc;
57 xfree(c->name);
58 ec_destroycurve(c->real);
59 DESTROY(c);
60}
61
62#define UNOP(OP) \
63 static ec *ec##OP(ec_curve *cc, ec *d, const ec *p) { \
64 ecctx *c = (ecctx *)cc; \
65 return (EC_##OP(c->real, d, p)); \
66 }
67
68#define BINOP(OP) \
69 static ec *ec##OP(ec_curve *cc, ec *d, const ec *p, const ec *q) { \
70 ecctx *c = (ecctx *)cc; \
71 return (EC_##OP(c->real, d, p, q)); \
72 }
73
74UNOP(IN)
75UNOP(OUT)
76UNOP(FIX)
77UNOP(NEG)
78UNOP(DBL)
79BINOP(ADD)
80BINOP(SUB)
81
82#undef UNOP
83#undef BINOP
84
85static ec *ecFIND(ec_curve *cc, ec *d, mp *x)
86{
87 ecctx *c = (ecctx *)cc;
88 return (EC_FIND(c->real, d, x));
89}
90
91static int ecCHECK(ec_curve *cc, const ec *p)
92{
93 ecctx *c = (ecctx *)cc;
94 return (EC_CHECK(c->real, p));
95}
96
34e4f738 97static int ecSAMEP(ec_curve *cc, ec_curve *dd)
98{
99 ecctx *c = (ecctx *)cc, *d = (ecctx *)dd;
100 return (ec_samep(c->real, d->real));
101}
102
4e66da02 103static const ec_ops ecops = {
f94b972d 104 "cardboard",
34e4f738 105 ecDESTROY, ecSAMEP, ecIN, ecOUT, ecFIX,
bc985cef 106 ecFIND, ecNEG, ecADD, ecSUB, ecDBL, ecCHECK
107};
108
109static ec_curve *ec_cutout(ec_curve *real, const char *name)
110{
111 ecctx *c = CREATE(ecctx);
112 c->c.f = real->f;
113 c->c.ops = &ecops;
34e4f738 114 c->c.a = real->a;
115 c->c.b = real->b;
bc985cef 116 c->magic = MAGIC;
117 c->name = xstrdup(name);
118 c->real = real;
119 return (&c->c);
120}
121
122static const char *ec_name(ec_curve *cc)
123{
124 ecctx *c = (ecctx *)cc;
125 assert(c->magic == MAGIC);
126 return (c->name);
127}
128
432c4e18 129/*----- Test field types --------------------------------------------------*/
bc985cef 130
131static void ecvcvt(const char *buf, dstr *d)
132{
bc985cef 133 ec_curve *v;
432c4e18 134 qd_parse qd;
135
136 qd.p = buf;
137 qd.e = 0;
138 if ((v = ec_curveparse(&qd)) == 0) {
139 fprintf(stderr, "bad curve `%.*s|%s': %s\n",
bb77b1d1 140 (int)(qd.p - buf), buf, qd.p, qd.e);
432c4e18 141 exit(1);
bc985cef 142 }
bc985cef 143 dstr_ensure(d, sizeof(v));
144 *(ec_curve **)d->buf = ec_cutout(v, buf);
145 d->len += sizeof(v);
146}
147
148static void ecvdump(dstr *d, FILE *fp)
149{
150 ec_curve *v = *(ec_curve **)d->buf;
151 fprintf(fp, "%s", ec_name(v));
152}
153
4e66da02 154const test_type type_ecurve = { ecvcvt, ecvdump };
bc985cef 155
156static void eccvt(const char *p, dstr *d)
157{
158 ec *a;
432c4e18 159 qd_parse qd;
bc985cef 160
432c4e18 161 qd.p = p;
162 qd.e = 0;
bc985cef 163 dstr_ensure(d, sizeof(ec));
164 a = (ec *)d->buf;
165 d->len += sizeof(ec);
166 ec_create(a);
432c4e18 167 if (!ec_ptparse(&qd, a)) {
bb77b1d1
MW
168 fprintf(stderr, "bad point `%.*s|%s': %s\n",
169 (int)(qd.p - p), p, qd.p, qd.e);
432c4e18 170 exit(1);
171 }
bc985cef 172}
173
174static void ecdodump(ec *a, FILE *fp)
175{
176 if (EC_ATINF(a))
177 fputs("inf", fp);
178 else {
45c0fd36 179 fputs("0x", fp);
bc985cef 180 mp_writefile(a->x, fp, 16);
181 fputs(", 0x", fp);
182 mp_writefile(a->y, fp, 16);
183 }
184}
185
186static void ecdump(dstr *d, FILE *fp)
187{
188 ec *a = (ec *)d->buf;
189 ecdodump(a, fp);
190}
191
5685a696 192const test_type type_ec = { eccvt, ecdump };
bc985cef 193
194/*----- Testing elliptic curve functionality ------------------------------*/
195
196#ifdef TEST_RIG
197
198static void ecdestroy(ec_curve *c)
199{
200 field *f = c->f;
201 ec_destroycurve(c);
202 F_DESTROY(f);
203}
204
205#define UNOP(op) \
206 static int v##op(dstr v[]) \
207 { \
208 ec_curve *e = *(ec_curve **)v[0].buf; \
209 ec *a = (ec *)v[1].buf; \
210 ec *r = (ec *)v[2].buf; \
211 ec c = EC_INIT; \
212 int ok = 1; \
213 ec_##op(e, &c, a); \
214 if (!EC_EQ(r, &c)) { \
215 fprintf(stderr, #op "failed"); \
216 fprintf(stderr, "\ncurve = "); type_ecurve.dump(v, stderr); \
217 fprintf(stderr, "\n a = "); ecdodump(a, stderr); \
218 fprintf(stderr, "\n r = "); ecdodump(r, stderr); \
219 fprintf(stderr, "\n c = "); ecdodump(&c, stderr); \
220 fprintf(stderr, "\n"); \
221 ok = 0; \
222 } \
223 EC_DESTROY(a); EC_DESTROY(r); EC_DESTROY(&c); \
224 ecdestroy(e); \
225 return (ok); \
226 }
227
228#define BINOP(op) \
229 static int v##op(dstr v[]) \
230 { \
231 ec_curve *e = *(ec_curve **)v[0].buf; \
232 ec *a = (ec *)v[1].buf; \
233 ec *b = (ec *)v[2].buf; \
234 ec *r = (ec *)v[3].buf; \
235 ec c = EC_INIT; \
236 int ok = 1; \
237 ec_##op(e, &c, a, b); \
238 if (!EC_EQ(r, &c)) { \
239 fprintf(stderr, #op "failed"); \
240 fprintf(stderr, "\ncurve = "); type_ecurve.dump(v, stderr); \
241 fprintf(stderr, "\n a = "); ecdodump(a, stderr); \
242 fprintf(stderr, "\n b = "); ecdodump(b, stderr); \
243 fprintf(stderr, "\n r = "); ecdodump(r, stderr); \
244 fprintf(stderr, "\n c = "); ecdodump(&c, stderr); \
245 fprintf(stderr, "\n"); \
246 ok = 0; \
247 } \
248 EC_DESTROY(a); EC_DESTROY(b); EC_DESTROY(r); EC_DESTROY(&c); \
249 ecdestroy(e); \
250 return (ok); \
251 }
252
253UNOP(neg)
254UNOP(dbl)
255BINOP(add)
256BINOP(sub)
257
258static int vcheck(dstr v[])
259{
260 ec_curve *e = *(ec_curve **)v[0].buf;
261 ec *a = (ec *)v[1].buf;
262 int r = *(int *)v[2].buf;
263 int c;
264 int ok = 1;
265 c = ec_check(e, a);
266 if (r != c) {
267 fprintf(stderr, "check failed");
268 fprintf(stderr, "\ncurve = "); type_ecurve.dump(v, stderr);
45c0fd36
MW
269 fprintf(stderr, "\n a = "); ecdodump(a, stderr);
270 fprintf(stderr, "\n r = %d", r);
271 fprintf(stderr, "\n c = %d", c);
bc985cef 272 fprintf(stderr, "\n");
273 ok = 0;
274 }
275 EC_DESTROY(a);
276 ecdestroy(e);
277 return (ok);
278}
279
280static int vmul(dstr v[])
281{
282 ec_curve *e = *(ec_curve **)v[0].buf;
283 ec *a = (ec *)v[1].buf;
284 mp *n = *(mp **)v[2].buf;
285 ec *r = (ec *)v[3].buf;
286 ec c = EC_INIT;
287 int ok = 1;
288 ec_mul(e, &c, a, n);
289 if (!EC_EQ(r, &c)) {
290 fprintf(stderr, "mul failed");
291 fprintf(stderr, "\ncurve = "); type_ecurve.dump(v, stderr);
45c0fd36
MW
292 fprintf(stderr, "\n a = "); ecdodump(a, stderr);
293 fprintf(stderr, "\n n = "); mp_writefile(n, stderr, 10);
294 fprintf(stderr, "\n r = "); ecdodump(r, stderr);
295 fprintf(stderr, "\n c = "); ecdodump(&c, stderr);
bc985cef 296 fprintf(stderr, "\n");
297 ok = 0;
298 }
299 EC_DESTROY(a); EC_DESTROY(r); EC_DESTROY(&c); MP_DROP(n);
300 ecdestroy(e);
301 return (ok);
302}
303
304static int vfind(dstr v[])
305{
306 ec_curve *e = *(ec_curve **)v[0].buf;
307 mp *x = *(mp **)v[1].buf;
308 ec *r = (ec *)v[2].buf;
309 ec c = EC_INIT;
310 int ok = 1;
311 if (!ec_find(e, &c, x)) EC_SETINF(&c);
312 if (!EC_EQ(r, &c)) {
313 fprintf(stderr, "find failed");
314 fprintf(stderr, "\ncurve = "); type_ecurve.dump(v, stderr);
45c0fd36
MW
315 fprintf(stderr, "\n x = "); mp_writefile(x, stderr, 16);
316 fprintf(stderr, "\n r = "); ecdodump(r, stderr);
317 fprintf(stderr, "\n c = "); ecdodump(&c, stderr);
bc985cef 318 fprintf(stderr, "\n");
319 ok = 0;
320 }
321 MP_DROP(x); EC_DESTROY(r); EC_DESTROY(&c);
322 ecdestroy(e);
323 return (ok);
324}
325
326static test_chunk tests[] = {
327 { "neg", vneg, { &type_ecurve, &type_ec, &type_ec } },
328 { "dbl", vdbl, { &type_ecurve, &type_ec, &type_ec } },
329 { "add", vadd, { &type_ecurve, &type_ec, &type_ec, &type_ec } },
330 { "sub", vsub, { &type_ecurve, &type_ec, &type_ec, &type_ec } },
331 { "mul", vmul, { &type_ecurve, &type_ec, &type_mp, &type_ec } },
332 { "check", vcheck, { &type_ecurve, &type_ec, &type_int } },
333 { "find", vfind, { &type_ecurve, &type_mp, &type_ec } },
334 { 0, 0, { 0 } }
335};
336
337int main(int argc, char *argv[])
338{
339 sub_init();
0f00dc4c 340 test_run(argc, argv, tests, SRCDIR "/t/ec");
bc985cef 341 return (0);
342}
343
344#endif
345
346/*----- That's all, folks -------------------------------------------------*/