/dports/graphics/blender/blender-2.91.0/intern/itasc/ |
H A D | WSDLSSolver.cpp | 85 e_scalar norm, mag, alpha, _qmax, Sinv, vmax, damp; in solve() local 90 Sinv = e_scalar(1.)/S; in solve() 121 M *= Sinv; in solve() 125 norm = fabs(Sinv*alpha*vmax); in solve() 127 damp = Sinv*alpha*_qmax/norm; in solve() 129 damp = Sinv*alpha; in solve()
|
/dports/science/qmcpack/qmcpack-3.11.0/src/Particle/LongRange/ |
H A D | LRBreakup.h | 158 if (Sinv[i] == 0.0) in DoBreakup() 170 coef *= Sinv[i]; in DoBreakup() 296 coef *= Sinv[i]; in DoGradBreakup() 538 if (Sinv[i] == 0.0) in DoBreakup() 548 coef *= Sinv[i]; in DoBreakup() 708 if (Sinv[i] == 0.0) in DoGradBreakup() 718 coef *= Sinv[i]; in DoGradBreakup() 878 if (Sinv[i] == 0.0) in DoStrainBreakup() 888 coef *= Sinv[i]; in DoStrainBreakup() 1143 if (Sinv[i] == 0.0) in DoAllBreakup() [all …]
|
/dports/science/openmx/openmx3.8/source/ |
H A D | TRAN_Calc_CurrentDensity.c | 13 dcomplex *Sinv) in TRAN_Calc_Sinv() argument 56 matrix1, &NUM_c, &zero, Sinv, &NUM_c); in TRAN_Calc_Sinv() 116 dcomplex *Sinv, in TRAN_Calc_CurrentDensity() argument 168 TRAN_MatProds5(NUM_c, "N", Sinv, VCC, GC, GammaR, GCdag, matrix1, matrix2); in TRAN_Calc_CurrentDensity() 169 if (minusk == 1) TRAN_MatProds5(NUM_c, "T", Sinv, VCC, GC, GammaR, GCdag, matrix1, matrix3); in TRAN_Calc_CurrentDensity() 176 TRAN_MatProds5(NUM_c, "N", GC, GammaR, GCdag, GammaL, Sinv, matrix1, matrix2); in TRAN_Calc_CurrentDensity() 184 TRAN_MatProds5(NUM_c, "N", GC, GammaL, GCdag, GammaR, Sinv, matrix1, matrix2); in TRAN_Calc_CurrentDensity() 256 dcomplex *Sinv, in TRAN_Calc_CurrentDensity_NC() argument 307 TRAN_MatProds5(NUM_c, "N", Sinv, VCC, GC, GammaR, GCdag, matrix1, matrix2); in TRAN_Calc_CurrentDensity_NC() 316 TRAN_MatProds5(NUM_c, "N", GC, GammaR, GCdag, GammaL, Sinv, matrix1, matrix2); in TRAN_Calc_CurrentDensity_NC() [all …]
|
/dports/science/dynare/dynare-4.6.4/matlab/partial_information/ |
H A D | PI_gensys.m | 61 Sinv=eye(FL_RANK); 63 Sinv(i,i)=1/S01(i,i); 65 F1=Sinv*U01'*a1*V01; 66 F2=Sinv*U01'*a1*V02; 67 F3=Sinv*U01'*a2*V01; 68 F4=Sinv*U01'*a2*V02; 69 F5=Sinv*U01'*PSI; 118 FF=F2; %=Sinv*U01'*a1*V02; 122 G31=-FF*G11+F4; % +Sinv*U01'*a2*V02; 123 G32=-FF*G12+F3; % +Sinv*U01'*a2*V01; [all …]
|
/dports/math/mfem/mfem-4.3/examples/ |
H A D | ex8.cpp | 168 BilinearForm *Sinv = new BilinearForm(test_space); in main() local 172 Sinv->AddDomainIntegrator(new InverseIntegrator(Sum)); in main() 173 Sinv->Assemble(); in main() 174 Sinv->Finalize(); in main() 184 SparseMatrix &matSinv = Sinv->SpMat(); in main() 278 delete Sinv; in main()
|
H A D | ex8p.cpp | 193 ParBilinearForm *Sinv = new ParBilinearForm(test_space); in main() local 197 Sinv->AddDomainIntegrator(new InverseIntegrator(Sum)); in main() 198 Sinv->Assemble(); in main() 199 Sinv->Finalize(); in main() 209 HypreParMatrix * matSinv = Sinv->ParallelAssemble(); delete Sinv; in main()
|
/dports/science/qmcpack/qmcpack-3.11.0/src/QMCHamiltonians/ |
H A D | ForceChiesaPBCAA.cpp | 57 Sinv.resize(N_basis, N_basis); in InitMatrix() 65 Sinv(k, j) = std::pow(Rcut, (m_exp + k + j + 3)) / static_cast<RealType>(m_exp + k + j + 3); in InitMatrix() 69 invert_matrix(Sinv, false); in InitMatrix() 71 MatrixOperators::product(Sinv, h.data(), c.data()); in InitMatrix() 241 tmp->Sinv.resize(N_basis, N_basis); in makeClone() 242 tmp->Sinv = Sinv; // terms in fitting polynomial in makeClone()
|
H A D | ForceCeperley.cpp | 65 Sinv.resize(N_basis, N_basis); in InitMatrix() 73 Sinv(k, j) = std::pow(Rcut, (m_exp + k + j + 3)) / static_cast<RealType>(m_exp + k + j + 3); in InitMatrix() 77 invert_matrix(Sinv, false); in InitMatrix() 79 MatrixOperators::product(Sinv, h.data(), c.data()); in InitMatrix()
|
/dports/math/py-statsmodels/statsmodels-0.13.1/statsmodels/sandbox/panel/ |
H A D | mixed.py | 81 def compute_P(self, Sinv): argument 88 self.P = self.W - np.dot(np.dot(t, Sinv), t.T) 130 def cov_random(self, D, Sinv=None): argument 142 if Sinv is not None: 143 self.compute_P(Sinv) 283 self.Sinv = L.pinv(S) 284 self.a = np.dot(self.Sinv, Y) 301 unit.compute_P(self.Sinv) 323 unit.compute_P(self.Sinv) 337 return self.Sinv [all …]
|
/dports/science/qwalk/mainline-1.0.1-300-g1b7e381/src/method/ |
H A D | Force_fitter.cpp | 38 Array2 <doublevar> Sinv(nexp, nexp); in setup() local 39 InvertMatrix(S, Sinv, nexp); in setup() 44 coeff(i)+=Sinv(i,j)*h(j); in setup()
|
H A D | Linear_optimization.cpp | 88 Array2 <doublevar> S(nparms+1,nparms+1),Sinv(nparms+1,nparms+1); in run() local 114 doublevar endiff= line_minimization(S,Sinv,H,alpha); in run() 165 doublevar find_directions(Array2 <doublevar> & S, Array2 <doublevar> & Sinv, in find_directions() argument 277 Array2 <doublevar> & Sinv, Array2 <doublevar> & H, Array1 <doublevar> & alpha) { in line_minimization() 287 doublevar psi_0_min=find_directions(S,Sinv,H,alpha_tmp,0.0,linear); in line_minimization() 288 doublevar psi_0_max=find_directions(S,Sinv,H,alpha_tmp,stabilmax,linear); in line_minimization() 309 guesspsi=find_directions(S,Sinv,H,alpha_tmp,guesstab,linear); in line_minimization() 328 doublevar prop_psi0=find_directions(S,Sinv,H,alphas(i),acc_stabils[i-1],linear); in line_minimization()
|
/dports/science/jdftx/jdftx-1.6.0/jdftx/fluid/ |
H A D | PCM_internal.cu | 83 void compute_kernel(int zBlock, const vector3<int> S, const vector3<> Sinv, const matrix3<> RTR, in compute_kernel() argument 86 compute_calc(i, iv, Sinv, RTR, nAtoms, x, nReps, reps, radius, shape, sigmaInv); in compute_kernel() 91 vector3<> Sinv(1./S[0], 1./S[1], 1./S[2]); in compute_gpu() local 93 …compute_kernel<<<glc.nBlocks,glc.nPerBlock>>>(zBlock, S, Sinv, RTR, nAtoms, x, nReps, reps, radius… in compute_gpu() 98 …void propagateGradient_kernel(int zBlock, const vector3<int> S, const vector3<> Sinv, const matrix… in propagateGradient_kernel() argument 102 …propagateGradient_calc(i, iv, Sinv, RTR, x, nReps, reps, radius, shape, E_shape, E_x, E_radius, si… in propagateGradient_kernel() 108 vector3<> Sinv(1./S[0], 1./S[1], 1./S[2]); in propagateGradient_gpu() local 110 …propagateGradient_kernel<<<glc.nBlocks,glc.nPerBlock>>>(zBlock, S, Sinv, RTR, x, nReps, reps, radi… in propagateGradient_gpu()
|
/dports/science/py-gpaw/gpaw-21.6.0/gpaw/test/parallel/ |
H A D | test_scalapack.py | 104 Sinv = dist.empty(dtype=dtype) 107 Sinv = dist.empty(dtype=dtype) 121 Glob2dist.redistribute(S0, Sinv, uplo='L') 134 scalapack_inverse(dist, Sinv, 'L') 140 Dist2glob.redistribute(Sinv, Sinv_test)
|
/dports/science/qmcpack/qmcpack-3.11.0/src/QMCHamiltonians/tests/ |
H A D | test_force.cpp | 103 REQUIRE(force2.Sinv.size() == force.Sinv.size()); in check_force_copy() 104 std::cout << force.Sinv << std::endl; in check_force_copy() 105 std::cout << force2.Sinv << std::endl; in check_force_copy() 106 for (int i = 0; i < force2.Sinv.rows(); i++) in check_force_copy() 108 for (int j = 0; j < force2.Sinv.cols(); j++) in check_force_copy() 111 REQUIRE(force2.Sinv(i, j) == Approx(force.Sinv(i, j))); in check_force_copy()
|
/dports/science/py-molmod/molmod-1.4.8/molmod/ |
H A D | unit_cells.py | 173 Sinv = np.zeros(S.shape, float) 176 Sinv[i] = 0.0 178 Sinv[i] = 1.0/S[i] 179 return np.dot(U*Sinv, Vt)*self.active
|
/dports/emulators/dolphin-emu/dolphin-3152428/Source/Core/Common/Crypto/ |
H A D | ec.cpp | 272 u8 Sinv[30]; in VerifySignature() local 274 bn_inv(Sinv, S, ec_N, 30); in VerifySignature() 279 bn_mul(w1, e, Sinv, ec_N, 30); in VerifySignature() 280 bn_mul(w2, R, Sinv, ec_N, 30); in VerifySignature()
|
/dports/math/py-statsmodels/statsmodels-0.13.1/statsmodels/base/ |
H A D | _constraints.py | 212 def transform_params_constraint(params, Sinv, R, q): argument 245 rsr = R.dot(Sinv).dot(R.T) 247 reduction = Sinv.dot(R.T).dot(np.linalg.solve(rsr, R.dot(params) - q))
|
/dports/graphics/geomview/geomview-1.9.5/src/lib/geometry/transform3/ |
H A D | tm3rotate.c | 177 Transform3 Sinv; in Tm3CarefulRotateTowardZ() local 204 Tm3Invert(S, Sinv); in Tm3CarefulRotateTowardZ() 206 Tm3Concat(T, Sinv, T); in Tm3CarefulRotateTowardZ()
|
/dports/science/minc2/minc-release-2.2.00/conversion/ana2mnc/ |
H A D | ana2mnc_xfm_reduce.pl | 688 my $Sinv = $M->shadow(); $Sinv->zero(); # Create and zero the inverse Scaling Matrix 689 $Sinv->assign(1, 1, $Sx); $Sx = 1/$Sx; 690 $Sinv->assign(2, 2, $Sy); $Sy = 1/$Sy; 691 $Sinv->assign(3, 3, $Sz); $Sz = 1/$Sz; 716 $M = $Sinv * $M;
|
/dports/games/kodi-addon-game.libretro.beetle-psx/game.libretro.beetle-psx-0.9.44.22-Matrix/depends/common/beetle-psx/deps/libkirk/ |
H A D | ec.c | 279 u8 Sinv[21]; in check_ecdsa() local 299 bn_mon_inv(Sinv, S, ec_N, 21); in check_ecdsa() 301 bn_mon_mul(w1, e, Sinv, ec_N, 21); in check_ecdsa() 303 bn_mon_mul(w2, R, Sinv, ec_N, 21); in check_ecdsa()
|
/dports/emulators/ppsspp/ppsspp-1.12.3/ext/libkirk/ |
H A D | ec.c | 312 u8 Sinv[21]; in check_ecdsa() local 332 bn_mon_inv(Sinv, S, ec_N, 21); in check_ecdsa() 334 bn_mon_mul(w1, e, Sinv, ec_N, 21); in check_ecdsa() 336 bn_mon_mul(w2, R, Sinv, ec_N, 21); in check_ecdsa()
|
/dports/emulators/ppsspp-qt5/ppsspp-1.12.3/ext/libkirk/ |
H A D | ec.c | 312 u8 Sinv[21]; in check_ecdsa() local 332 bn_mon_inv(Sinv, S, ec_N, 21); in check_ecdsa() 334 bn_mon_mul(w1, e, Sinv, ec_N, 21); in check_ecdsa() 336 bn_mon_mul(w2, R, Sinv, ec_N, 21); in check_ecdsa()
|
/dports/emulators/libretro-ppsspp/ppsspp-1.12.3/ext/libkirk/ |
H A D | ec.c | 312 u8 Sinv[21]; in check_ecdsa() local 332 bn_mon_inv(Sinv, S, ec_N, 21); in check_ecdsa() 334 bn_mon_mul(w1, e, Sinv, ec_N, 21); in check_ecdsa() 336 bn_mon_mul(w2, R, Sinv, ec_N, 21); in check_ecdsa()
|
/dports/games/libretro-beetle_psx/beetle-psx-libretro-3ec155d/deps/libkirk/ |
H A D | ec.c | 279 u8 Sinv[21]; in check_ecdsa() local 299 bn_mon_inv(Sinv, S, ec_N, 21); in check_ecdsa() 301 bn_mon_mul(w1, e, Sinv, ec_N, 21); in check_ecdsa() 303 bn_mon_mul(w2, R, Sinv, ec_N, 21); in check_ecdsa()
|
/dports/science/helfem/HelFEM-21461e9/src/diatomic/ |
H A D | completeness.cpp | 58 arma::mat Sinv(Sinvh*arma::trans(Sinvh)); in main() local 103 arma::mat Y(arma::diagvec(Plcao*Sinv*arma::trans(Plcao))); in main()
|