| /dports/math/e-antic/e-antic-1.0.0-rc.13/libeantic/upstream/antic/ulong_extras/ |
| H A D | sqrtmod.c | 21 mp_limb_t pinv; in n_sqrtmod() local
50 pinv = n_preinvert_limb(p); in n_sqrtmod()
63 g = n_mulmod2_preinv(b, b, p, pinv); in n_sqrtmod()
68 g = n_powmod2_ui_preinv(2, (p - 1)/4, p, pinv); in n_sqrtmod()
69 return n_mulmod2_preinv(g, b, p, pinv); in n_sqrtmod()
80 b = n_powmod2_ui_preinv(a, p1, p, pinv); in n_sqrtmod()
87 g = n_powmod2_ui_preinv(k, p1, p, pinv); in n_sqrtmod()
88 res = n_powmod2_ui_preinv(a, (p1 + 1) / 2, p, pinv); in n_sqrtmod()
106 res = n_mulmod2_preinv(res, gpow, p, pinv); in n_sqrtmod()
107 g = n_mulmod2_preinv(gpow, gpow, p, pinv); in n_sqrtmod()
[all …]
|
| /dports/math/e-antic/flint2-ae7ec89/ulong_extras/ |
| H A D | sqrtmod.c | 21 mp_limb_t pinv; in n_sqrtmod() local
50 pinv = n_preinvert_limb(p); in n_sqrtmod()
63 g = n_mulmod2_preinv(b, b, p, pinv); in n_sqrtmod()
68 g = n_powmod2_ui_preinv(2, (p - 1)/4, p, pinv); in n_sqrtmod()
69 return n_mulmod2_preinv(g, b, p, pinv); in n_sqrtmod()
80 b = n_powmod2_ui_preinv(a, p1, p, pinv); in n_sqrtmod()
87 g = n_powmod2_ui_preinv(k, p1, p, pinv); in n_sqrtmod()
88 res = n_powmod2_ui_preinv(a, (p1 + 1) / 2, p, pinv); in n_sqrtmod()
106 res = n_mulmod2_preinv(res, gpow, p, pinv); in n_sqrtmod()
107 g = n_mulmod2_preinv(gpow, gpow, p, pinv); in n_sqrtmod()
[all …]
|
| /dports/math/flint2/flint-2.8.4/ulong_extras/ |
| H A D | sqrtmod.c | 21 mp_limb_t pinv; in n_sqrtmod() local
50 pinv = n_preinvert_limb(p); in n_sqrtmod()
63 g = n_mulmod2_preinv(b, b, p, pinv); in n_sqrtmod()
68 g = n_powmod2_ui_preinv(2, (p - 1)/4, p, pinv); in n_sqrtmod()
69 return n_mulmod2_preinv(g, b, p, pinv); in n_sqrtmod()
80 b = n_powmod2_ui_preinv(a, p1, p, pinv); in n_sqrtmod()
87 g = n_powmod2_ui_preinv(k, p1, p, pinv); in n_sqrtmod()
88 res = n_powmod2_ui_preinv(a, (p1 + 1) / 2, p, pinv); in n_sqrtmod()
106 res = n_mulmod2_preinv(res, gpow, p, pinv); in n_sqrtmod()
107 g = n_mulmod2_preinv(gpow, gpow, p, pinv); in n_sqrtmod()
[all …]
|
| /dports/math/py-qdldl/qdldl-0.1.5.post0/c/amd/src/ |
| H A D | perm.c | 32 void pinv(QDLDL_int const *p, QDLDL_int * pinv, QDLDL_int n) { in pinv() argument
33 for (QDLDL_int k = 0; k < n; k++) pinv[p[k]] = k; /* invert the permutation */ in pinv()
44 const QDLDL_int * pinv, in symperm() argument
51 j2 = pinv ? pinv[j] : j; /* column j of A is column j2 of C */ in symperm()
57 i2 = pinv ? pinv[i] : i; /* row i of A is row i2 of C */ in symperm()
64 j2 = pinv ? pinv[j] : j; /* column j of A is column j2 of C */ in symperm()
71 i2 = pinv ? pinv[i] : i; /* row i of A is row i2 in symperm()
|
| /dports/science/qmcpack/qmcpack-3.11.0/src/Numerics/ |
| H A D | determinant_operators.h | 72 inline void det_row_update(T* restrict pinv,
82 BLAS::gemv('T', m, m, c_ratio, pinv, m, tv, 1, const_traits<T>::zero(), temp, 1);
84 memcpy(rcopy, pinv + m * rowchanged, m * sizeof(T));
85 BLAS::ger(m, m, const_traits<T>::minus_one(), rcopy, 1, temp, 1, pinv, m);
89 inline void det_col_update(T* restrict pinv,
99 BLAS::gemv('N', m, m, c_ratio, pinv, m, tv, 1, T(), temp, 1);
101 BLAS::copy(m, pinv + colchanged, m, rcopy, 1);
102 BLAS::ger(m, m, -1.0, temp, 1, rcopy, 1, pinv, m);
140 inline T getRatioByColSubstitution(const T* restrict pinv, const T* restrict tc, int m, int colchan…
142 return BLAS::dot(m, pinv + colchanged, m, tc, 1);
[all …]
|
| /dports/games/naev/naev-0.7.0/lib/csparse/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const int *pinv, int values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/casadi/casadi-3.5.5/external_packages/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const int *pinv, int values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-csparse/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-klu/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-ldl/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-btf/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-amd/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-mongoose/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-rbio/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-colamd/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-camd/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-config/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-ccolamd/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-cholmod/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-umfpack/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-slip_lu/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-spqr/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-cxsparse/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/suitesparse-graphblas/SuiteSparse-5.10.1/CSparse/Source/ |
| H A D | cs_symperm.c | 3 cs *cs_symperm (const cs *A, const csi *pinv, csi values) in cs_symperm() argument
16 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
21 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
28 j2 = pinv ? pinv [j] : j ; /* column j of A is column j2 of C */ in cs_symperm()
33 i2 = pinv ? pinv [i] : i ; /* row i of A is row i2 of C */ in cs_symperm()
|
| /dports/math/igraph/igraph-0.9.5/vendor/cs/ |
| H A D | cs_lu.c | 9 CS_INT *Lp, *Li, *Up, *Ui, *pinv, *xi, *q, n, ipiv, k, top, p, i, col, lnz,unz; in cs_lu() local
19 N->pinv = pinv = cs_malloc (n, sizeof (CS_INT)) ; /* allocate result pinv */ in cs_lu()
20 if (!L || !U || !pinv) return (cs_ndone (N, NULL, xi, x, 0)) ; in cs_lu()
23 for (i = 0 ; i < n ; i++) pinv [i] = -1 ; /* no rows pivotal yet */ in cs_lu()
38 top = cs_spsolve (L, A, col, xi, x, pinv, 1) ; /* x = L\A(:,col) */ in cs_lu()
45 if (pinv [i] < 0) /* row i is not yet pivotal */ in cs_lu()
55 Ui [unz] = pinv [i] ; in cs_lu()
61 if (pinv [col] < 0 && CS_ABS (x [col]) >= a*tol) ipiv = col ; in cs_lu()
66 pinv [ipiv] = k ; /* ipiv is the kth pivot row */ in cs_lu()
72 if (pinv [i] < 0) /* x(i) is an entry in L(:,k) */ in cs_lu()
[all …]
|