index 45b417b..233e3f5 100644 (file)
*
*/

+
+
#include "cryptlib.h"
#include "bn_lcl.h"

@@ -261,7 +263,8 @@ BIGNUM *int_bn_mod_inverse(BIGNUM *in,
if (!BN_nnmod(B, B, A, ctx)) goto err;
}
sign = -1;
-       /* From  B = a mod |n|,  A = |n|  it follows that
+       /*-
+        * From  B = a mod |n|,  A = |n|  it follows that
*
*      0 <= B < A,
*     -sign*X*a  ==  B   (mod |n|),
@@ -278,7 +281,7 @@ BIGNUM *int_bn_mod_inverse(BIGNUM *in,

while (!BN_is_zero(B))
{
-                       /*
+                       /*-
*      0 < B < |n|,
*      0 < A <= |n|,
* (1) -sign*X*a  ==  B   (mod |n|),
@@ -325,7 +328,8 @@ BIGNUM *int_bn_mod_inverse(BIGNUM *in,
}

-                       /* We still have (1) and (2).
+                       /*-
+                        * We still have (1) and (2).
* Both  A  and  B  are odd.
* The following computations ensure that
*
@@ -361,7 +365,7 @@ BIGNUM *int_bn_mod_inverse(BIGNUM *in,
{
BIGNUM *tmp;

-                       /*
+                       /*-
*      0 < B < A,
* (*) -sign*X*a  ==  B   (mod |n|),
*      sign*Y*a  ==  A   (mod |n|)
@@ -408,7 +412,8 @@ BIGNUM *int_bn_mod_inverse(BIGNUM *in,
if (!BN_div(D,M,A,B,ctx)) goto err;
}

-                       /* Now
+                       /*-
+                        * Now
*      A = D*B + M;
* thus we have
* (**)  sign*Y*a  ==  D*B + M   (mod |n|).
@@ -421,7 +426,8 @@ BIGNUM *int_bn_mod_inverse(BIGNUM *in,
B=M;
/* ... so we have  0 <= B < A  again */

-                       /* Since the former  M  is now  B  and the former  B  is now  A,
+                       /*-
+                        * Since the former  M  is now  B  and the former  B  is now  A,
* (**) translates into
*       sign*Y*a  ==  D*A + B    (mod |n|),
* i.e.
@@ -474,7 +480,7 @@ BIGNUM *int_bn_mod_inverse(BIGNUM *in,
}
}

-       /*
+       /*-
* The while loop (Euclid's algorithm) ends when
*      A == gcd(a,n);
* we have
@@ -563,7 +569,8 @@ static BIGNUM *BN_mod_inverse_no_branch(BIGNUM *in,
if (!BN_nnmod(B, pB, A, ctx)) goto err;
}
sign = -1;
-       /* From  B = a mod |n|,  A = |n|  it follows that
+       /*-
+        * From  B = a mod |n|,  A = |n|  it follows that
*
*      0 <= B < A,
*     -sign*X*a  ==  B   (mod |n|),
@@ -574,7 +581,7 @@ static BIGNUM *BN_mod_inverse_no_branch(BIGNUM *in,
{
BIGNUM *tmp;

-               /*
+               /*-
*      0 < B < A,
* (*) -sign*X*a  ==  B   (mod |n|),
*      sign*Y*a  ==  A   (mod |n|)
@@ -589,7 +596,8 @@ static BIGNUM *BN_mod_inverse_no_branch(BIGNUM *in,
/* (D, M) := (A/B, A%B) ... */
if (!BN_div(D,M,pA,B,ctx)) goto err;

-               /* Now
+               /*-
+                * Now
*      A = D*B + M;
* thus we have
* (**)  sign*Y*a  ==  D*B + M   (mod |n|).
@@ -602,7 +610,8 @@ static BIGNUM *BN_mod_inverse_no_branch(BIGNUM *in,
B=M;
/* ... so we have  0 <= B < A  again */

-               /* Since the former  M  is now  B  and the former  B  is now  A,
+               /*-
+                * Since the former  M  is now  B  and the former  B  is now  A,
* (**) translates into
*       sign*Y*a  ==  D*A + B    (mod |n|),
* i.e.
@@ -630,7 +639,7 @@ static BIGNUM *BN_mod_inverse_no_branch(BIGNUM *in,
sign = -sign;
}

-       /*
+       /*-
* The while loop (Euclid's algorithm) ends when
*      A == gcd(a,n);
* we have