1*f0865ec9SKyle Evans /*
2*f0865ec9SKyle Evans * Copyright (C) 2017 - This file is part of libecc project
3*f0865ec9SKyle Evans *
4*f0865ec9SKyle Evans * Authors:
5*f0865ec9SKyle Evans * Ryad BENADJILA <ryadbenadjila@gmail.com>
6*f0865ec9SKyle Evans * Arnaud EBALARD <arnaud.ebalard@ssi.gouv.fr>
7*f0865ec9SKyle Evans * Jean-Pierre FLORI <jean-pierre.flori@ssi.gouv.fr>
8*f0865ec9SKyle Evans *
9*f0865ec9SKyle Evans * Contributors:
10*f0865ec9SKyle Evans * Nicolas VIVET <nicolas.vivet@ssi.gouv.fr>
11*f0865ec9SKyle Evans * Karim KHALFALLAH <karim.khalfallah@ssi.gouv.fr>
12*f0865ec9SKyle Evans *
13*f0865ec9SKyle Evans * This software is licensed under a dual BSD and GPL v2 license.
14*f0865ec9SKyle Evans * See LICENSE file at the root folder of the project.
15*f0865ec9SKyle Evans */
16*f0865ec9SKyle Evans #include <libecc/fp/fp_add.h>
17*f0865ec9SKyle Evans #include <libecc/nn/nn_add.h>
18*f0865ec9SKyle Evans
19*f0865ec9SKyle Evans /*
20*f0865ec9SKyle Evans * Compute out = in1 + in2 mod p. 'out' parameter must have been initialized
21*f0865ec9SKyle Evans * by the caller. Returns 0 on success, -1 on error.
22*f0865ec9SKyle Evans *
23*f0865ec9SKyle Evans * Aliasing is supported.
24*f0865ec9SKyle Evans */
fp_add(fp_t out,fp_src_t in1,fp_src_t in2)25*f0865ec9SKyle Evans int fp_add(fp_t out, fp_src_t in1, fp_src_t in2)
26*f0865ec9SKyle Evans {
27*f0865ec9SKyle Evans int ret, cmp;
28*f0865ec9SKyle Evans
29*f0865ec9SKyle Evans ret = fp_check_initialized(out); EG(ret, err);
30*f0865ec9SKyle Evans ret = fp_check_initialized(in1); EG(ret, err);
31*f0865ec9SKyle Evans ret = fp_check_initialized(in2); EG(ret, err);
32*f0865ec9SKyle Evans
33*f0865ec9SKyle Evans MUST_HAVE(((&(in1->ctx->p)) == (&(in2->ctx->p))), ret, err);
34*f0865ec9SKyle Evans MUST_HAVE(((&(in1->ctx->p)) == (&(out->ctx->p))), ret, err);
35*f0865ec9SKyle Evans FORCE_USED_VAR(cmp); /* silence warning when macro results in nothing */
36*f0865ec9SKyle Evans SHOULD_HAVE(!nn_cmp(&in1->fp_val, &(in1->ctx->p), &cmp) && (cmp < 0), ret, err);
37*f0865ec9SKyle Evans SHOULD_HAVE(!nn_cmp(&in2->fp_val, &(in2->ctx->p), &cmp) && (cmp < 0), ret, err);
38*f0865ec9SKyle Evans
39*f0865ec9SKyle Evans ret = nn_mod_add(&(out->fp_val), &(in1->fp_val),
40*f0865ec9SKyle Evans &(in2->fp_val), &(in1->ctx->p));
41*f0865ec9SKyle Evans
42*f0865ec9SKyle Evans err:
43*f0865ec9SKyle Evans return ret;
44*f0865ec9SKyle Evans }
45*f0865ec9SKyle Evans
46*f0865ec9SKyle Evans /*
47*f0865ec9SKyle Evans * Compute out = in + 1 mod p. 'out' parameter must have been initialized
48*f0865ec9SKyle Evans * by the caller. Returns 0 on success, -1 on error.
49*f0865ec9SKyle Evans *
50*f0865ec9SKyle Evans * Aliasing is supported.
51*f0865ec9SKyle Evans */
fp_inc(fp_t out,fp_src_t in)52*f0865ec9SKyle Evans int fp_inc(fp_t out, fp_src_t in)
53*f0865ec9SKyle Evans {
54*f0865ec9SKyle Evans int ret, cmp;
55*f0865ec9SKyle Evans
56*f0865ec9SKyle Evans ret = fp_check_initialized(in); EG(ret, err);
57*f0865ec9SKyle Evans ret = fp_check_initialized(out); EG(ret, err);
58*f0865ec9SKyle Evans
59*f0865ec9SKyle Evans MUST_HAVE(((&(in->ctx->p)) == (&(out->ctx->p))), ret, err);
60*f0865ec9SKyle Evans FORCE_USED_VAR(cmp); /* silence warning when macro results in nothing */
61*f0865ec9SKyle Evans SHOULD_HAVE(!nn_cmp(&in->fp_val, &(in->ctx->p), &cmp) && (cmp < 0), ret, err);
62*f0865ec9SKyle Evans
63*f0865ec9SKyle Evans ret = nn_mod_inc(&(out->fp_val), &(in->fp_val), &(in->ctx->p));
64*f0865ec9SKyle Evans
65*f0865ec9SKyle Evans err:
66*f0865ec9SKyle Evans return ret;
67*f0865ec9SKyle Evans }
68*f0865ec9SKyle Evans
69*f0865ec9SKyle Evans /*
70*f0865ec9SKyle Evans * Compute out = in1 - in2 mod p. 'out' parameter must have been initialized
71*f0865ec9SKyle Evans * by the caller. Returns 0 on success, -1 on error.
72*f0865ec9SKyle Evans *
73*f0865ec9SKyle Evans * Aliasing is supported.
74*f0865ec9SKyle Evans */
fp_sub(fp_t out,fp_src_t in1,fp_src_t in2)75*f0865ec9SKyle Evans int fp_sub(fp_t out, fp_src_t in1, fp_src_t in2)
76*f0865ec9SKyle Evans {
77*f0865ec9SKyle Evans int ret, cmp;
78*f0865ec9SKyle Evans
79*f0865ec9SKyle Evans ret = fp_check_initialized(out); EG(ret, err);
80*f0865ec9SKyle Evans ret = fp_check_initialized(in1); EG(ret, err);
81*f0865ec9SKyle Evans ret = fp_check_initialized(in2); EG(ret, err);
82*f0865ec9SKyle Evans
83*f0865ec9SKyle Evans MUST_HAVE(((&(in1->ctx->p)) == (&(in2->ctx->p))), ret, err);
84*f0865ec9SKyle Evans MUST_HAVE(((&(in1->ctx->p)) == (&(out->ctx->p))), ret, err);
85*f0865ec9SKyle Evans FORCE_USED_VAR(cmp); /* silence warning when macro results in nothing */
86*f0865ec9SKyle Evans SHOULD_HAVE(!nn_cmp(&in1->fp_val, &(in1->ctx->p), &cmp) && (cmp < 0), ret, err);
87*f0865ec9SKyle Evans SHOULD_HAVE(!nn_cmp(&in2->fp_val, &(in2->ctx->p), &cmp) && (cmp < 0), ret, err);
88*f0865ec9SKyle Evans
89*f0865ec9SKyle Evans ret = nn_mod_sub(&(out->fp_val), &(in1->fp_val),
90*f0865ec9SKyle Evans &(in2->fp_val), &(in1->ctx->p));
91*f0865ec9SKyle Evans
92*f0865ec9SKyle Evans err:
93*f0865ec9SKyle Evans return ret;
94*f0865ec9SKyle Evans }
95*f0865ec9SKyle Evans
96*f0865ec9SKyle Evans /*
97*f0865ec9SKyle Evans * Compute out = in - 1 mod p. 'out' parameter must have been initialized
98*f0865ec9SKyle Evans * by the caller. Returns 0 on success, -1 on error.
99*f0865ec9SKyle Evans *
100*f0865ec9SKyle Evans * Aliasing is supported.
101*f0865ec9SKyle Evans */
fp_dec(fp_t out,fp_src_t in)102*f0865ec9SKyle Evans int fp_dec(fp_t out, fp_src_t in)
103*f0865ec9SKyle Evans {
104*f0865ec9SKyle Evans int ret, cmp;
105*f0865ec9SKyle Evans
106*f0865ec9SKyle Evans ret = fp_check_initialized(out); EG(ret, err);
107*f0865ec9SKyle Evans ret = fp_check_initialized(in); EG(ret, err);
108*f0865ec9SKyle Evans
109*f0865ec9SKyle Evans MUST_HAVE(((&(in->ctx->p)) == (&(out->ctx->p))), ret, err);
110*f0865ec9SKyle Evans FORCE_USED_VAR(cmp); /* silence warning when macro results in nothing */
111*f0865ec9SKyle Evans SHOULD_HAVE(!nn_cmp(&in->fp_val, &(in->ctx->p), &cmp) && (cmp < 0), ret, err);
112*f0865ec9SKyle Evans
113*f0865ec9SKyle Evans ret = nn_mod_dec(&(out->fp_val), &(in->fp_val), &(in->ctx->p));
114*f0865ec9SKyle Evans
115*f0865ec9SKyle Evans err:
116*f0865ec9SKyle Evans return ret;
117*f0865ec9SKyle Evans }
118*f0865ec9SKyle Evans
119*f0865ec9SKyle Evans /*
120*f0865ec9SKyle Evans * Compute out = -in mod p = (p - in) mod p. 'out' parameter must have been
121*f0865ec9SKyle Evans * initialized by the caller. Returns 0 on success, -1 on error.
122*f0865ec9SKyle Evans *
123*f0865ec9SKyle Evans * Aliasing is supported.
124*f0865ec9SKyle Evans */
fp_neg(fp_t out,fp_src_t in)125*f0865ec9SKyle Evans int fp_neg(fp_t out, fp_src_t in)
126*f0865ec9SKyle Evans {
127*f0865ec9SKyle Evans int ret, cmp;
128*f0865ec9SKyle Evans
129*f0865ec9SKyle Evans ret = fp_check_initialized(in); EG(ret, err);
130*f0865ec9SKyle Evans ret = fp_check_initialized(out); EG(ret, err);
131*f0865ec9SKyle Evans
132*f0865ec9SKyle Evans MUST_HAVE(((&(in->ctx->p)) == (&(out->ctx->p))), ret, err);
133*f0865ec9SKyle Evans FORCE_USED_VAR(cmp); /* silence warning when macro results in nothing */
134*f0865ec9SKyle Evans SHOULD_HAVE(!nn_cmp(&in->fp_val, &(in->ctx->p), &cmp) && (cmp < 0), ret, err);
135*f0865ec9SKyle Evans
136*f0865ec9SKyle Evans ret = nn_sub(&(out->fp_val), &(in->ctx->p), &(in->fp_val));
137*f0865ec9SKyle Evans
138*f0865ec9SKyle Evans err:
139*f0865ec9SKyle Evans return ret;
140*f0865ec9SKyle Evans }
141