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 */ 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 */ 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 */ 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 */ 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 */ 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