Home
last modified time | relevance | path

Searched refs:a_MSDptr (Results 1 – 12 of 12) sorted by relevance

/dports/math/cln/cln-1.3.6/src/base/digitseq/
H A Dcl_DS_sqrt.cc108 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 Dcl_DS_recip.cc23 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 Dcl_DS_recipsqrt.cc26 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 Dcl_DS_div.cc167 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 Dcl_DS.h2649 #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 Dcl_I_gcd.cc62 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 Dcl_I_xgcd.cc228 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 Dtimerecip2adic.cc24 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 Dtimediv2adic.cc27 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 Dintgcd.d596 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 Dintsqrt.d114 # 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 Dintdiv.d4 # 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]