Home
last modified time | relevance | path

Searched refs:Qk (Results 1 – 25 of 793) sorted by relevance

12345678910>>...32

/dports/math/gmp/gmp-6.2.1/mpz/
H A Dlucmod.c62 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 Dstronglucas.c62 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 Dlucmod.c62 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 Dstronglucas.c62 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 Dtestopbfem.cc108 … 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 Dqr.c23 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 Dlq.c24 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 Dprimetest.py193 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 Dkalman_filter.cxx248 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 Dkalman_filter.h101 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 DrotateHamiltonian_Helper2.hpp48 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 DrotateHamiltonian.hpp270 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 Drrel_muse_table.cxx86 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 Dprimality.c417 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 Decp_contract.f4 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 Decp_contract.f4 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 Dex40.bat8 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 Dex40.sh13 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 Dprimality.c286 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 Dsbasis.cpp441 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_names30 '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 Dsbasis.cpp626 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 Dtest_fe.cpp97 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 Disoline.cpp620 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 DPathSimplify.cpp238 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 …]

12345678910>>...32