xref: /freebsd/crypto/libecc/src/fp/fp_add.c (revision f0865ec9906d5a18fa2a3b61381f22ce16e606ad)
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