/dports/math/cln/cln-1.3.6/src/base/digitseq/ |
H A D | cl_DS_sqrt.cc | 108 while ((a_len>0) && (mspref(a_MSDptr,0)==0)) { msshrink(a_MSDptr); a_len--; } in cl_UDS_sqrt() 138 #define a_MSDptr new_a_MSDptr in cl_UDS_sqrt() macro 152 cl_UDS_recipsqrt(a_MSDptr,2*n,c_MSDptr,n); in cl_UDS_sqrt() 175 if (!inc_loop_lsp(a_MSDptr mspop (n-1),n-1)) in cl_UDS_sqrt() 178 else if (test_loop_msp(a_MSDptr,n-1)) in cl_UDS_sqrt() 197 if (test_loop_msp(a_MSDptr mspop (n-1),n+1)) in cl_UDS_sqrt() 239 mspref(a_MSDptr,0) = highD(a_msdd); mspref(a_MSDptr,1) = lowD(a_msdd); in cl_UDS_sqrt() 244 mspref(a_MSDptr,0) = a_msd - x2hi; in cl_UDS_sqrt() 246 mspref(a_MSDptr,1) = a_2msd - x2lo; in cl_UDS_sqrt() 252 var uintD* a_mptr = a_MSDptr mspop 0; in cl_UDS_sqrt() [all …]
|
H A D | cl_DS_recip.cc | 23 extern void cl_UDS_recip (const uintD* a_MSDptr, uintC a_len, 81 void cl_UDS_recip (const uintD* a_MSDptr, uintC a_len, in cl_UDS_recip() argument 98 shiftrightcopy_loop_msp(a_MSDptr,x_MSDptr,a_len,1,0); in cl_UDS_recip() 102 shiftrightcopy_loop_msp(a_MSDptr,x_MSDptr,b_len,1,0) in cl_UDS_recip() 103 | ((mspref(a_MSDptr,b_len) & -bit(intDsize-3)) >> 1); in cl_UDS_recip() 106 { var uintD x1 = mspref(a_MSDptr,0); in cl_UDS_recip() 107 var uintD x2 = (a_len > 1 ? (mspref(a_MSDptr,1) & -bit(intDsize-3)) : 0); in cl_UDS_recip()
|
H A D | cl_DS_recipsqrt.cc | 26 extern void cl_UDS_recipsqrt (const uintD* a_MSDptr, uintC a_len, 74 void cl_UDS_recipsqrt (const uintD* a_MSDptr, uintC a_len, in cl_UDS_recipsqrt() argument 79 var const uintD* const x_MSDptr = a_MSDptr; in cl_UDS_recipsqrt()
|
H A D | cl_DS_div.cc | 167 void cl_UDS_divide (const uintD* a_MSDptr, uintC a_len, const uintD* a_LSDptr, in cl_UDS_divide() argument 172 while ((a_len>0) && (mspref(a_MSDptr,0)==0)) { msshrink(a_MSDptr); a_len--; } in cl_UDS_divide() 200 …{var uintD rest = divucopy_loop_msp(mspref(b_MSDptr,0),a_MSDptr,q_MSDptr,a_len); // Division durch… in cl_UDS_divide()
|
H A D | cl_DS.h | 2649 #define UDS_divide(a_MSDptr,a_len,a_LSDptr,b_MSDptr,b_len,b_LSDptr,q_,r_) \ argument 2653 cl_UDS_divide(a_MSDptr,_a_len,a_LSDptr,b_MSDptr,b_len,b_LSDptr,roomptr,q_,r_); 2654 extern void cl_UDS_divide (const uintD* a_MSDptr, uintC a_len, const uintD* a_LSDptr, 2668 #define UDS_sqrt(a_MSDptr,a_len,a_LSDptr,b_,squarep_zuweisung) \ argument 2672 squarep_zuweisung cl_UDS_sqrt(a_MSDptr,_a_len,a_LSDptr,b_); \ 2674 extern bool cl_UDS_sqrt (const uintD* a_MSDptr, uintC a_len, const uintD* a_LSDptr, DS* b_); 2679 extern void cl_UDS_recip (const uintD* a_MSDptr, uintC a_len, 2684 extern void cl_UDS_recipsqrt (const uintD* a_MSDptr, uintC a_len,
|
/dports/math/cln/cln-1.3.6/src/integer/gcd/ |
H A D | cl_I_gcd.cc | 62 var uintD* a_MSDptr; in gcd() local 105 swap(uintD*, a_MSDptr,b_MSDptr); in gcd() 112 while (mspref(a_MSDptr,0) == 0) { msshrink(a_MSDptr); a_len--; } // normalisieren in gcd() 306 var uintD* a_MSDptr; in gcd() local 342 swap(uintD*, a_MSDptr,b_MSDptr); in gcd() 355 { a0 = mspref(a_MSDptr,0); } in gcd() 357 { a0 = divu_loop_msp(b0,a_MSDptr,a_len); in gcd() 391 highlowDD(a_msd, mspref(a_MSDptr,1)); in gcd() 417 a_nsd = mspref(a_MSDptr,1); in gcd() 499 { mspref(a_MSDptr,0) -= 1; } in gcd() [all …]
|
H A D | cl_I_xgcd.cc | 228 var uintD* a_MSDptr; in xgcd() local 297 swap(uintD*, a_MSDptr,b_MSDptr); in xgcd() 317 a_msd = mspref(a_MSDptr,0); // führendes Digit von a in xgcd() 334 highlowDD(a_msd, (a_len==1 ? 0 : mspref(a_MSDptr,1))); in xgcd() 356 { a_msd |= mspref(a_MSDptr,1) >> a_msd_size; in xgcd() 362 { a_nsd = mspref(a_MSDptr,1); in xgcd() 388 a_nsd = mspref(a_MSDptr,1) << (intDsize-a_msd_size); in xgcd() 391 { a_nsd |= mspref(a_MSDptr,2) >> a_msd_size; in xgcd() 453 { mspref(a_MSDptr,0) -= 1; } in xgcd() 456 while (mspref(a_MSDptr,0)==0) { msshrink(a_MSDptr); a_len--; } in xgcd() [all …]
|
/dports/math/cln/cln-1.3.6/tests/ |
H A D | timerecip2adic.cc | 24 uintD* a_MSDptr; in main() local 28 num_stack_alloc(len,a_MSDptr=,a_LSDptr=); in main() 30 random_UDS(default_random_state,a_MSDptr,len); in main()
|
H A D | timediv2adic.cc | 27 uintD* a_MSDptr; in main() local 35 num_stack_alloc(a_len,a_MSDptr=,a_LSDptr=); in main() 39 random_UDS(default_random_state,a_MSDptr,a_len); in main()
|
/dports/lang/clisp/clisp-df3b9f6fdcff22832898e89a989eb499c0f842ed/src/ |
H A D | intgcd.d | 596 var uintD* a_MSDptr; in I_I_gcd_I() local 613 # Jetzt ist a = a_MSDptr/a_len/a_LSDptr, b = b_MSDptr/b_len/b_LSDptr, in I_I_gcd_I() 638 swap(uintD*, a_MSDptr,b_MSDptr); in I_I_gcd_I() 653 a_msd = a_MSDptr[0]; # führendes Digit von a in I_I_gcd_I() 671 highlowDD(a_msd, (a_len==1 ? 0 : a_MSDptr[1])); in I_I_gcd_I() 691 a_msd |= a_MSDptr[1] >> a_msd_size; in I_I_gcd_I() 726 *--a_MSDptr = 0; a_len++; in I_I_gcd_I() 759 a_MSDptr[0] -= 1; in I_I_gcd_I() 762 while (a_MSDptr[0]==0) { in I_I_gcd_I() 763 a_MSDptr++; a_len--; in I_I_gcd_I() [all …]
|
H A D | intsqrt.d | 114 # UDS_sqrt(a_MSDptr,a_len,a_LSDptr, &b, squarep=) 115 # > a_MSDptr/a_len/a_LSDptr: eine UDS 177 while ((a_len>0) && (a_MSDptr[0]==0)) { in UDS_sqrt_() 178 a_MSDptr++; a_len--; in UDS_sqrt_() 215 # Nun ist A = a_MSDptr/2n/.. in UDS_sqrt_() 256 a_MSDptr[0] = highD(a_msdd); a_MSDptr[1] = lowD(a_msdd); in UDS_sqrt_() 262 a_MSDptr[0] = a_msd - x2hi; in UDS_sqrt_() 264 a_MSDptr[0] -= 1; in UDS_sqrt_() 265 a_MSDptr[1] = a_2msd - x2lo; in UDS_sqrt_() 271 var uintD* a_mptr = &a_MSDptr[0]; in UDS_sqrt_() [all …]
|
H A D | intdiv.d | 4 # UDS_divide(a_MSDptr,a_len,a_LSDptr, b_MSDptr,b_len,b_LSDptr, &q,&r); 5 # Die UDS a = a_MSDptr/a_len/a_LSDptr (a>=0) wird durch 16 #define UDS_divide(a_MSDptr,a_len,a_LSDptr,b_MSDptr,b_len,b_LSDptr,q_,r_) \ argument 20 UDS_divide_(a_MSDptr,_a_len,a_LSDptr,b_MSDptr,b_len,b_LSDptr,roomptr,q_,r_); \ 22 local void UDS_divide_ (uintD* a_MSDptr, uintC a_len, uintD* a_LSDptr, 70 local void UDS_divide_ (uintD* a_MSDptr, uintC a_len, uintD* a_LSDptr, 75 # a normalisieren (a_MSDptr erhöhen, a_len erniedrigen): 76 while ((a_len>0) && (a_MSDptr[0]==0)) { 77 a_MSDptr++; a_len--; 108 var uintD rest = divucopy_loop_up(b_MSDptr[0],a_MSDptr,q_MSDptr,a_len); # Division durch b[0]
|