/dports/math/gmp/gmp-6.2.1/mpz/ |
H A D | lucmod.c | 62 mpz_set_si (Qk, Q); in mpz_lucas_mod() 65 mpz_set_ui (Qk, 1); /* U2 = 1 */ in mpz_lucas_mod() 80 mpz_mul (T1, Qk, Qk); /* U_{k+1}^2 */ in mpz_lucas_mod() 81 mpz_sub (Qk, V, Qk); /* |U_{k+1} - U_k| */ in mpz_lucas_mod() 82 mpz_mul (T2, Qk, Qk); /* |U_{k+1} - U_k|^2 */ in mpz_lucas_mod() 83 mpz_mul (Qk, V, V); /* U_k^2 */ in mpz_lucas_mod() 86 mpz_submul_ui (T1, Qk, Q); in mpz_lucas_mod() 98 mpz_tdiv_r (Qk, T1, n); in mpz_lucas_mod() 102 res = SIZ (Qk) == 0; in mpz_lucas_mod() 114 mpz_mul (T1, Qk, Qk); /* P^2 U_k^2 = U_k^2 */ in mpz_lucas_mod() [all …]
|
H A D | stronglucas.c | 62 mpz_stronglucas (mpz_srcptr x, mpz_ptr V, mpz_ptr Qk) in mpz_stronglucas() argument 159 if (mpz_lucas_mod (V, Qk, Q, b0, n, T1, T2) == 0) in mpz_stronglucas() 165 mpz_submul_ui (T2, Qk, 2); in mpz_stronglucas() 170 mpz_mul (T2, Qk, Qk); in mpz_stronglucas() 171 mpz_tdiv_r (Qk, T2, n); in mpz_stronglucas()
|
/dports/math/kalker/gmp-mpfr-sys-21966f4bfb56c87d407eb14c72e92ef4e55856aa/gmp-6.2.1-c/mpz/ |
H A D | lucmod.c | 62 mpz_set_si (Qk, Q); in mpz_lucas_mod() 65 mpz_set_ui (Qk, 1); /* U2 = 1 */ in mpz_lucas_mod() 80 mpz_mul (T1, Qk, Qk); /* U_{k+1}^2 */ in mpz_lucas_mod() 81 mpz_sub (Qk, V, Qk); /* |U_{k+1} - U_k| */ in mpz_lucas_mod() 82 mpz_mul (T2, Qk, Qk); /* |U_{k+1} - U_k|^2 */ in mpz_lucas_mod() 83 mpz_mul (Qk, V, V); /* U_k^2 */ in mpz_lucas_mod() 86 mpz_submul_ui (T1, Qk, Q); in mpz_lucas_mod() 98 mpz_tdiv_r (Qk, T1, n); in mpz_lucas_mod() 102 res = SIZ (Qk) == 0; in mpz_lucas_mod() 114 mpz_mul (T1, Qk, Qk); /* P^2 U_k^2 = U_k^2 */ in mpz_lucas_mod() [all …]
|
H A D | stronglucas.c | 62 mpz_stronglucas (mpz_srcptr x, mpz_ptr V, mpz_ptr Qk) in mpz_stronglucas() argument 159 if (mpz_lucas_mod (V, Qk, Q, b0, n, T1, T2) == 0) in mpz_stronglucas() 165 mpz_submul_ui (T2, Qk, 2); in mpz_stronglucas() 170 mpz_mul (T2, Qk, Qk); in mpz_stronglucas() 171 mpz_tdiv_r (Qk, T2, n); in mpz_stronglucas()
|
/dports/math/dune-pdelab/dune-pdelab-20c7085389d3eb4f8ca99e1bc60f7fa6036536c8/dune/pdelab/test/ |
H A D | testopbfem.cc | 108 … MassMatrixTest<double,1,1,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() 109 … MassMatrixTest<double,2,1,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() 110 … MassMatrixTest<double,3,1,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() 111 … MassMatrixTest<double,4,1,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() 112 … MassMatrixTest<double,5,1,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() 113 … MassMatrixTest<double,6,1,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() 114 … MassMatrixTest<double,7,1,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() 116 … MassMatrixTest<double,1,2,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() 117 … MassMatrixTest<double,2,2,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() 118 … MassMatrixTest<double,3,2,Dune::GeometryType::cube,Dune::GMPField<512>,BasisType::Qk >::compute(); in testmassmatrix() [all …]
|
/dports/math/libflame/libflame-5.2.0/old/lapack-test/play/ |
H A D | qr.c | 23 FLA_Obj A, Ak, T, Tk, D, Dk, A_copy, A_recovered, R, Q, Qk, W, x, y, z; in main() local 67 FLA_Obj_create( datatype, m, k, 0,0, &Qk ); FLA_Set_to_identity( Qk ); in main() 89 FLA_Apply_Q_UT_create_workspace_side( FLA_LEFT, Tk, Qk, &W ); in main() 91 Ak, Tk, W, Qk ); in main() 96 FLA_Set( FLA_ZERO, Qk ); in main() 97 FLA_Copy( Ak, Qk ); in main() 98 FLA_QR_UT_form_Q( Qk, Tk, Qk ); in main() 103 FLA_ONE, Qk, Qk, FLA_ZERO, Dk ); in main() 143 FLA_Obj_fshow( stdout, " - Qk - ", Qk, "% 6.4e", "------"); in main() 160 FLA_Obj_free( &Qk ); in main()
|
H A D | lq.c | 24 FLA_Obj A, Ak, T, Tk, D, Dk, A_copy, A_recovered, L, Q, Qk, W, x, y, z; in main() local 68 FLA_Obj_create( datatype, k, n, 0, 0, &Qk ); FLA_Set_to_identity( Qk ); in main() 94 FLA_Set( FLA_ZERO, Qk ); in main() 95 FLA_Copy( Ak, Qk ); in main() 96 FLA_LQ_UT_form_Q( Ak, Tk, Qk ); in main() 98 FLA_Apply_Q_UT_create_workspace( Tk, Qk, &W ); in main() 100 Ak, Tk, W, Qk ); in main() 106 FLA_ONE, Qk, Qk, FLA_ZERO, Dk ); in main() 156 FLA_Obj_fshow( stdout, " - Qk - ", Qk, "% 6.4e", "------"); in main() 172 FLA_Obj_free( &Qk ); in main()
|
/dports/math/py-sympy/sympy-1.9/sympy/ntheory/ |
H A D | primetest.py | 193 Qk = Q 212 if Qk == 1: 216 Qk = 1 225 Qk = -1 230 V = (V*V - 2*Qk) % n 231 Qk *= Qk 240 Qk *= Q 241 Qk %= n 242 return _int_tuple(U % n, V % n, Qk) 392 V = (V*V - 2*Qk) % n [all …]
|
/dports/misc/vxl/vxl-3.3.2/contrib/oul/ouel/ |
H A D | kalman_filter.cxx | 248 KalmanFilter::predict(const vnl_matrix<double> &Qk) in predict() argument 251 if ((Qk.rows()!=num_signal_dimensions)|| in predict() 252 (Qk.cols()!=num_signal_dimensions)) in predict() 263 P = A*P*A.transpose() + Qk; in predict() 282 KalmanFilter::predict(const vnl_matrix<double> &Qk, in predict() argument 286 if ((Qk.rows()!=num_signal_dimensions)|| in predict() 287 (Qk.cols()!=num_signal_dimensions)) in predict() 306 P = A*P*A.transpose() + Qk; in predict() 329 const vnl_matrix<double> &Qk in update_predict() argument 333 return predict(Qk); in update_predict()
|
H A D | kalman_filter.h | 101 vnl_matrix<double> predict(const vnl_matrix<double> &Qk); 107 vnl_matrix<double> predict(const vnl_matrix<double> &Qk, 113 const vnl_matrix<double> &Qk);
|
/dports/science/qmcpack/qmcpack-3.11.0/src/AFQMC/Hamiltonians/ |
H A D | rotateHamiltonian_Helper2.hpp | 48 MatQk const& Qk, in count_Qk_x_Rl() argument 54 assert(Qk.size(0) == Ta.size(0)); in count_Qk_x_Rl() 55 assert(Qk.size(1) == Rl.size(0)); in count_Qk_x_Rl() 68 nwork = std::min(int(Qk.size(0)), ncores); in count_Qk_x_Rl() 70 std::tie(ka0, kaN) = FairDivideBoundary(coreid, int(Qk.size(0)), nwork); in count_Qk_x_Rl() 78 ma::product(Qk, Rl(Rl.extension(0), {bl0, blN}), Ta(Ta.extension(0), {bl0, blN})); in count_Qk_x_Rl() 207 MatQk const& Qk, in Qk_x_Rl() argument 214 assert(Qk.size(0) == Ta.size(0)); in Qk_x_Rl() 215 assert(Qk.size(1) == Rl.size(0)); in Qk_x_Rl() 228 nwork = std::min(int(Qk.size(0)), ncores); in Qk_x_Rl() [all …]
|
H A D | rotateHamiltonian.hpp | 270 shmSpMatrix Qk({dummy_nrow, dummy_ncol}, shared_allocator<SPComplexType>{TG.Node()}); in rotateHijkl() local 272 std::fill_n(Qk.origin(), Qk.num_elements(), SPComplexType(0.0)); in rotateHijkl() 345 int n0_, n1_, sz_ = Qk.size(0); in rotateHijkl() 348 ma::transpose(Qk.sliced(n0_, n1_), Rl(Rl.extension(0), {n0_, n1_})); in rotateHijkl() 533 … std::copy(Qk.origin() + bi * maxnk * NEL0 * nvec, Qk.origin() + (bi * maxnk + nk) * NEL0 * nvec, in rotateHijkl() 654 … std::copy(Qk.origin() + bi * maxnk * NEL0 * nvec, Qk.origin() + (bi * maxnk + nk) * NEL0 * nvec, in rotateHijkl() 758 std::fill_n(Qk.origin(), Qk.num_elements(), SPComplexType(0.0)); in rotateHijkl_single_node() 792 int n0_, n1_, sz_ = Qk.size(0); in rotateHijkl_single_node() 795 ma::transpose(Qk.sliced(n0_, n1_), Rl(Rl.extension(0), {n0_, n1_})); in rotateHijkl_single_node() 859 Qk.sliced(NAEA * NMO + k0 * NAEB, NAEA * NMO + kN * NAEB), in rotateHijkl_single_node() [all …]
|
/dports/misc/vxl/vxl-3.3.2/contrib/rpl/rrel/ |
H A D | rrel_muse_table.cxx | 86 double pk, qk, Qk, pQk, Qk_prime, Qk_dprime, Qk_tprime, vrk; in calculate_standard_dev() local 92 Qk = expected_kth; // ak(k, N); // inverse cdf of absolute residuals in calculate_standard_dev() 95 pQk = std::exp( -0.5 * Qk*Qk) * std::sqrt(2.0 / vnl_math::pi); in calculate_standard_dev() 106 Qk_dprime = Qk/(pQk*pQk); in calculate_standard_dev() 109 Qk_tprime = ( 1.0 + 2.0 * Qk*Qk ) / (pQk*pQk*pQk); in calculate_standard_dev()
|
/dports/math/p5-Math-Prime-Util/Math-Prime-Util-0.73/ |
H A D | primality.c | 417 V = mulsubmod(V, V, addmod(Qk,Qk,n), n); in lucas_seq() 418 Qk = sqrmod(Qk, n); in lucas_seq() 425 Qk = mulmod(Qk, Qmod, n); in lucas_seq() 622 Qk = mont_sqrmod(Qk,n); in is_lucas_pseudoprime() 629 Qk = mont_mulmod(Qk, montQ, n); in is_lucas_pseudoprime() 645 Qk = mont_sqrmod(Qk,n); in is_lucas_pseudoprime() 658 Qk = mont_sqrmod(Qk,n); in is_lucas_pseudoprime() 693 V = mulsubmod(V, V, addmod(Qk,Qk,n), n); in is_lucas_pseudoprime() 694 Qk = sqrmod(Qk, n); in is_lucas_pseudoprime() 707 V = mulsubmod(V, V, addmod(Qk,Qk,n), n); in is_lucas_pseudoprime() [all …]
|
/dports/science/nwchem/nwchem-7b21660b82ebd85ef659f6fba7e1e73433b0bd0a/src/NWints/ecp/ |
H A D | ecp_contract.f | 4 subroutine ecp_contract (n_ab,n_c,m_count,Qk,coeff_c,Q) argument 22 double precision Qk(n_ab,n_c,m_count),Q(n_ab,m_count), local 33 Q(i_ab,m) = Q(i_ab,m)+Qk(i_ab,i_c,m)*coeff_c(i_c)
|
/dports/science/nwchem-data/nwchem-7.0.2-release/src/NWints/ecp/ |
H A D | ecp_contract.f | 4 subroutine ecp_contract (n_ab,n_c,m_count,Qk,coeff_c,Q) argument 22 double precision Qk(n_ab,n_c,m_count),Q(n_ab,m_count), 33 Q(i_ab,m) = Q(i_ab,m)+Qk(i_ab,i_c,m)*coeff_c(i_c)
|
/dports/graphics/gmt/gmt-6.3.0/doc/examples/ex40/ |
H A D | ex40.bat | 8 gmt spatial @GSHHS_h_Australia.txt -fg -Qk > centroid.txt 9 …gmt spatial @GSHHS_h_Australia.txt -fg -Qk | gawk "{printf \"Full area = %%.0f km@+2@+\n\", $3}" >… 17 …gmt spatial T100k.txt -fg -Qk | gawk "{printf \"Reduced area = %%.0f km@+2@+\n\", $3}" > area_T100… 29 …gmt spatial T500k.txt -fg -Qk | gawk "{printf \"Reduced area = %%.0f km@+2@+\n\", $3}" > area_T500…
|
H A D | ex40.sh | 13 gmt spatial @GSHHS_h_Australia.txt -fg -Qk > centroid.txt 22 …gmt spatial T100k.txt -fg -Qk | $AWK '{printf "Reduced area = %.0f km@+2@+\n", $3}' > area_T100k.t… 34 …gmt spatial T500k.txt -fg -Qk | $AWK '{printf "Reduced area = %.0f km@+2@+\n", $3}' > area_T500k.t…
|
/dports/math/p5-Math-Prime-Util-GMP/Math-Prime-Util-GMP-0.52/ |
H A D | primality.c | 286 mpz_t Qk, mpz_t t) in lucas_seq() argument 310 mpz_set_si(Qk, Q); in lucas_seq() 359 mpz_submul_ui(V, Qk, 2); in lucas_seq() 361 mpz_mul(Qk, Qk, Qk); /* Q2k = Qk^2 */ in lucas_seq() 376 mpz_mul_si(Qk, Qk, Q); in lucas_seq() 378 mpz_mod(Qk, Qk, n); in lucas_seq() 567 mpz_t v, k, V, U, Qk, t; in llr() local 743 mpz_t d, U, V, Qk, t; in _GMP_is_lucas_pseudoprime() local 792 mpz_submul_ui(V, Qk, 2); in _GMP_is_lucas_pseudoprime() 794 mpz_mulmod(Qk, Qk, Qk, n, t); in _GMP_is_lucas_pseudoprime() [all …]
|
/dports/print/scribus-devel/scribus-1.5.7/scribus/third_party/lib2geom/ |
H A D | sbasis.cpp | 441 SBasis Pk=Linear(1)-g,Qk=g,sg=Pk*Qk; in compose_inverse() local 443 Qk.truncate(order); in compose_inverse() 445 Qk.resize(order,Linear(0.)); in compose_inverse() 453 double q10 = Qk.at(k)[0];// p10*a + q10*b = r10 in compose_inverse() 454 double q01 = Qk.at(k)[1];// & in compose_inverse() 469 r=r-Pk*a-Qk*b; in compose_inverse() 472 Qk=Qk*sg; in compose_inverse() 474 Qk.truncate(order); in compose_inverse()
|
/dports/shells/zsh/zsh-5.8.1/Completion/Zsh/Type/ |
H A D | _command_names | 30 'builtins:builtin command:compadd -Qk builtins' 32 'aliases:alias:compadd -Qk aliases' 34 'reserved-words:reserved word:compadd -Qk reswords'
|
/dports/math/lib2geom/lib2geom-1.1/src/2geom/ |
H A D | sbasis.cpp | 626 SBasis Pk=Linear(1)-g,Qk=g,sg=Pk*Qk; in compose_inverse() local 628 Qk.truncate(order); in compose_inverse() 630 Qk.resize(order,Linear(0.)); in compose_inverse() 641 double q10 = Qk.at(k)[0];// p10*a + q10*b = r10 in compose_inverse() 642 double q01 = Qk.at(k)[1];// & in compose_inverse() 656 r=r-Pk*a-Qk*b; in compose_inverse() 659 Qk=Qk*sg; in compose_inverse() 662 Qk.resize(order,Linear(0.)); in compose_inverse()
|
/dports/math/mfem/mfem-4.3/tests/unit/fem/ |
H A D | test_fe.cpp | 97 REQUIRE( fe.Space() == (int) FunctionSpace::Qk ); 161 REQUIRE( fe.Space() == (int) FunctionSpace::Qk ); 193 REQUIRE( fe.Space() == (int) FunctionSpace::Qk ); 286 REQUIRE( fe.Space() == (int) FunctionSpace::Qk ); 348 REQUIRE( fe.Space() == (int) FunctionSpace::Qk ); 410 REQUIRE( fe.Space() == (int) FunctionSpace::Qk ); 472 REQUIRE( fe.Space() == (int) FunctionSpace::Qk ); 568 REQUIRE( fe.Space() == (int) FunctionSpace::Qk ); 632 REQUIRE( fe.Space() == (int) FunctionSpace::Qk ); 664 REQUIRE( fe.Space() == (int) FunctionSpace::Qk );
|
/dports/math/freefem++/FreeFem-sources-4.6/plugin/seq/ |
H A D | isoline.cpp | 620 R2 Qk[6]; in operator ( )() local 622 int np = IsoLineK(Pk, fk, Qk, i1, i2, eps); in operator ( )() 639 P.push_back(R2_I2(Qk[i])); in operator ( )() 652 cout << " +++ " << Qk[0] << " -> " << Qk[1] << " :: " << p[0] << " -> " << p[1] << endl; in operator ( )() 656 *fff << Qk[0] << "\n" in operator ( )() 657 << Qk[1] << "\n" in operator ( )() 658 << ((Qk[0] * 0.4 + Qk[1] * .6) + R2(Qk[0], Qk[1]).perp( ) * .4) << "\n\n"; in operator ( )() 680 R2 Qk[2]; in operator ( )() local 700 cout << " add : " << Qk[0] << ", " << i1[0] << ',' << i2[0] << " -> " << Qk[1] in operator ( )() 711 P.push_back(R2_I2(Qk[i])); in operator ( )() [all …]
|
/dports/graphics/inkscape/inkscape-1.1_2021-05-24_c4e8f9ed74/src/livarot/ |
H A D | PathSimplify.cpp | 238 g_free(data.Qk); in DoSimplify() 293 Q[0] += N13 (tk[i]) * Qk[i]; in FitCubic() 294 Q[1] += N23 (tk[i]) * Qk[i]; in FitCubic() 311 Q[0] += N13 (tk[i]) * Qk[i]; in FitCubic() 332 data.Qk = (double *) g_realloc(data.Qk, data.maxPt * sizeof(double)); in ExtendFit() 752 double *Qk; // les Qk in AttemptSimplify() local 837 g_free(Qk); in AttemptSimplify() 880 g_free(Qk); in AttemptSimplify() 1001 g_free(Qk); in AttemptSimplify() 1093 g_free(Qk); in AttemptSimplify() [all …]
|