Changeset 16965
- Timestamp:
- 11/27/13 12:06:34 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/shared/Bamg/BigPrimeNumber.cpp
r14949 r16965 7 7 /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/AGoodNumberPrimeWith)*/ 8 8 9 / /list of big prime numbers9 /*list of big prime numbers*/ 10 10 const long BigPrimeNumber[] ={ 567890359L, 11 11 567890431L, 567890437L, 567890461L, 567890471L, … … 13 13 567890591L, 567890599L, 567890621L, 567890629L , 0}; 14 14 15 / /initialize o and pi16 long o =0;17 long pi =BigPrimeNumber[1];15 /*initialize o and pi*/ 16 long o = 0; 17 long pi = BigPrimeNumber[1]; 18 18 19 / /loop until BigPrimeNumber[i]==0 (end of BigPrimeNumber)19 /*loop until BigPrimeNumber[i]==0 (end of BigPrimeNumber)*/ 20 20 for (int i=0; BigPrimeNumber[i]; i++){ 21 21 22 / /compute r, rest of the remainder of the division of BigPrimeNumber[i] by n22 /*compute r, remainder of the division of BigPrimeNumber[i] by n*/ 23 23 long r = BigPrimeNumber[i] % n; 24 24 25 25 /*compute oo = min ( r , n-r , |n - 2r|, |n-3r|)*/ 26 long oo = Min(Min(r,n-r),Min(Abs(n-2*r),Abs(n-3*r)));27 if (o < oo){28 o =oo;29 pi =BigPrimeNumber[i];26 long oo = Min(Min(r,n-r),Min(Abs(n-2*r),Abs(n-3*r))); 27 if(o < oo){ 28 o = oo; 29 pi = BigPrimeNumber[i]; 30 30 } 31 31 }
Note:
See TracChangeset
for help on using the changeset viewer.