/dports/math/py-spvcm/spvcm-0.3.0/spvcm/both_levels/generic/ |
H A D | sample.py | 24 logdet = splogdet(PsiRho) 29 return -.5*logdet -.5 * kernel + st.Log_Rho0(val) 45 logdet = splogdet(PsiLambda) 49 return -.5*logdet - .5*kernel + st.Log_Lambda0(val) 62 logdet = -splogdet(PsiLambdai) #negative because precision 66 return -.5 * logdet - .5 * kernel + st.Log_Lambda0(val) 79 logdet = -splogdet(PsiRhoi) 85 return -.5 * logdet - .5 * kernel + st.Log_Rho0(val)
|
/dports/science/qmcpack/qmcpack-3.11.0/src/QMCWaveFunctions/Fermion/ |
H A D | SPOSetProxyForMSD.cpp | 100 ValueMatrix_t& logdet, in evaluate() argument 110 logdet(i, p) = psiM(p, q); in evaluate() 120 ValueMatrix_t& logdet, in evaluate_notranspose() argument 130 logdet(p, i) = psiM(p, q); in evaluate_notranspose() 140 ValueMatrix_t& logdet, in evaluate() argument 150 logdet(i, p) = psiM(p, q); in evaluate() 160 ValueMatrix_t& logdet, in evaluate() argument 171 logdet(i, p) = psiM(p, q); in evaluate() 182 ValueMatrix_t& logdet, in evaluate_notranspose() argument 192 logdet(p, i) = psiM(p, q); in evaluate_notranspose() [all …]
|
H A D | SPOSetProxy.h | 62 ValueMatrix_t& logdet, 68 ValueMatrix_t& logdet, 74 ValueMatrix_t& logdet, 81 ValueMatrix_t& logdet, 89 ValueMatrix_t& logdet, 96 ValueMatrix_t& logdet,
|
H A D | SPOSetProxyForMSD.h | 70 ValueMatrix_t& logdet, 76 ValueMatrix_t& logdet, 82 ValueMatrix_t& logdet, 89 ValueMatrix_t& logdet, 97 ValueMatrix_t& logdet, 104 ValueMatrix_t& logdet,
|
H A D | SPOSetProxy.cpp | 66 ValueMatrix_t& logdet, in evaluate() argument 75 logdet(i, j) = psiM(j, i); in evaluate() 86 ValueMatrix_t& logdet, in evaluate_notranspose() argument 94 logdet(i, j) = psiM(i, j); in evaluate_notranspose() 104 ValueMatrix_t& logdet, in evaluate() argument 114 ValueMatrix_t& logdet, in evaluate() argument 125 ValueMatrix_t& logdet, in evaluate_notranspose() argument 135 ValueMatrix_t& logdet, in evaluate_notranspose() argument
|
/dports/math/cppad/CppAD-20210000.8/include/cppad/utility/ |
H A D | lu_solve.hpp | 271 Float &logdet ) in LuSolve() argument 315 logdet = Float(0); in LuSolve() 323 logdet = Float(0); in LuSolve() 329 { logdet += log( - pivot ); in LuSolve() 333 logdet += log( pivot ); in LuSolve()
|
/dports/math/SCIP/scip-7.0.3/src/cppad/utility/ |
H A D | lu_solve.hpp | 272 Float &logdet ) in LuSolve() argument 316 logdet = Float(0); in LuSolve() 324 logdet = Float(0); in LuSolve() 330 { logdet += log( - pivot ); in LuSolve() 333 else logdet += log( pivot ); in LuSolve()
|
/dports/math/tmv/tmv-0.75/src/ |
H A D | TMV_SymLDLDecompose.cpp | 278 logdet += TMV_LOG(absd); in NonBlockLDL_Decompose() 314 logdet += TMV_LOG(absd); in NonBlockLDL_Decompose() 426 ptrdiff_t* P, RT& logdet, T& signdet) in BlockLDL_Decompose() argument 772 ptrdiff_t* P, RT& logdet, T& signdet) in NonLapLDL_Decompose() argument 793 ptrdiff_t* P, RT& logdet, T& signdet) in LapLDL_Decompose() argument 993 logdet += TMV_LOG(absd); in LapLDL_Decompose() 1044 logdet += TMV_LOG(absd); in LapLDL_Decompose() 1056 ptrdiff_t* P, float& logdet, float& signdet) in LapLDL_Decompose() argument 1239 logdet += TMV_LOG(absd); in LapLDL_Decompose() 1290 logdet += TMV_LOG(absd); in LapLDL_Decompose() [all …]
|
H A D | TMV_SymCHD.cpp | 49 mutable RT logdet; member 64 zerodet(false), logdet(0), donedet(false) {} in HermCHDiv_Impl() 130 pimpl->logdet = DiagMatrixViewOf(pimpl->LLx.diag()).logDet(&s); in det() 131 pimpl->logdet *= RT(2); in det() 136 else return TMV_EXP(pimpl->logdet); in det() 144 pimpl->logdet = DiagMatrixViewOf(pimpl->LLx.diag()).logDet(&s); in logDet() 145 pimpl->logdet *= RT(2); in logDet() 153 return pimpl->logdet; in logDet()
|
H A D | TMV_SymBandCHD.cpp | 50 mutable RT logdet; member 71 zerodet(false), logdet(1), donedet(false) in HermBandCHDiv_Impl() 174 pimpl->logdet = DiagMatrixViewOf(pimpl->LLx.diag()).logDet(&s); in det() 175 if (pimpl->LLx.nlo() > 1) { pimpl->logdet *= RT(2); } in det() 180 else return TMV_EXP(pimpl->logdet); in det() 188 pimpl->logdet = DiagMatrixViewOf(pimpl->LLx.diag()).logDet(&s); in logDet() 189 if (pimpl->LLx.nlo() > 1) { pimpl->logdet *= RT(2); } in logDet() 197 return pimpl->logdet; in logDet()
|
/dports/science/qmcpack/qmcpack-3.11.0/src/QMCWaveFunctions/LCAO/ |
H A D | LCAOrbitalSetWithCorrection.cpp | 65 ValueMatrix_t& logdet, in evaluate_notranspose() argument 69 LCAOrbitalSet::evaluate_notranspose(P, first, last, logdet, dlogdet, d2logdet); in evaluate_notranspose() 72 cusp.add_vgl(P, iat, i, logdet, dlogdet, d2logdet); in evaluate_notranspose() 79 ValueMatrix_t& logdet, in evaluate_notranspose() argument 89 ValueMatrix_t& logdet, in evaluate_notranspose() argument
|
H A D | LCAOrbitalSet.h | 116 ValueMatrix_t& logdet, 123 ValueMatrix_t& logdet, 130 ValueMatrix_t& logdet, 235 ValueMatrix_t& logdet, 245 ValueMatrix_t& logdet, 258 ValueMatrix_t& logdet,
|
H A D | LCAOrbitalSetWithCorrection.h | 62 ValueMatrix_t& logdet, 69 ValueMatrix_t& logdet, 76 ValueMatrix_t& logdet,
|
H A D | LCAOrbitalSet.cpp | 391 ValueMatrix_t& logdet, in evaluate_vgl_impl() argument 395 const size_t output_size = logdet.cols(); in evaluate_vgl_impl() 396 std::copy_n(temp.data(0), output_size, logdet[i]); in evaluate_vgl_impl() 512 ValueMatrix_t& logdet, in evaluate_notranspose() argument 521 evaluate_vgl_impl(Temp, i, logdet, dlogdet, d2logdet); in evaluate_notranspose() 526 assert(logdet.cols() <= OrbitalSetSize); in evaluate_notranspose() 532 evaluate_vgl_impl(Tempv, i, logdet, dlogdet, d2logdet); in evaluate_notranspose() 540 ValueMatrix_t& logdet, in evaluate_notranspose() argument 554 assert(logdet.cols() <= OrbitalSetSize); in evaluate_notranspose() 568 ValueMatrix_t& logdet, in evaluate_notranspose() argument [all …]
|
/dports/math/py-spvcm/spvcm-0.3.0/spvcm/both_levels/se_se/ |
H A D | sample.py | 17 logdet = splogdet(PsiRhoi) 22 return .5*logdet -.5 * kernel + st.Log_Rho0(val) #since precision, no negative on ld 38 logdet = splogdet(PsiLambdai) 42 return .5*logdet - .5*kernel + st.Log_Lambda0(val)
|
/dports/math/dsdp/DSDP5.8/src/solver/ |
H A D | dualimpl.c | 261 int DSDPComputePotential(DSDP dsdp, DSDPVec y, double logdet, double *potential){ in DSDPComputePotential() argument 267 if (gap>0) dpotential=dsdp->rho*log(gap)-logdet; in DSDPComputePotential() 270 DSDPLogInfo(0,9,"Gap: %4.4e, Log Determinant: %4.4e, Log Gap: %4.4e\n",gap,logdet,log(gap)); in DSDPComputePotential() 287 int DSDPComputePotential2(DSDP dsdp, DSDPVec y, double mu, double logdet, double *potential){ in DSDPComputePotential2() argument 292 *potential=-(ddobj + mu*logdet)*dsdp->schurmu; in DSDPComputePotential2() 293 *potential=-(ddobj/mu + logdet)*dsdp->schurmu; in DSDPComputePotential2() 309 int DSDPSetY(DSDP dsdp, double beta, double logdet, DSDPVec ynew){ in DSDPSetY() argument 330 dsdp->logdet=logdet; in DSDPSetY() 331 info=DSDPComputePotential(dsdp,dsdp->y,dsdp->logdet,&dsdp->potential);DSDPCHKERR(info); in DSDPSetY()
|
H A D | dualalg.c | 27 double dstep,newpotential,logdet; in DSDPYStepLineSearch() local 42 info=DSDPComputeLogSDeterminant(dsdp,&logdet);DSDPCHKERR(info); in DSDPYStepLineSearch() 55 info=DSDPSetY(dsdp,dstep,logdet,dsdp->ytemp);DSDPCHKERR(info); in DSDPYStepLineSearch() 57 info=DSDPSetY(dsdp,0,dsdp->logdet,dsdp->y);DSDPCHKERR(info); in DSDPYStepLineSearch() 77 double dstep,newpotential,bdotdy,oldpotential,logdet; in DSDPYStepLineSearch2() local 93 info=DSDPComputeLogSDeterminant(dsdp,&logdet);DSDPCHKERR(info); in DSDPYStepLineSearch2() 107 info=DSDPSetY(dsdp,dstep,logdet,dsdp->ytemp);DSDPCHKERR(info); in DSDPYStepLineSearch2() 109 info=DSDPSetY(dsdp,0,dsdp->logdet,dsdp->y);DSDPCHKERR(info); in DSDPYStepLineSearch2() 339 info=DSDPComputeLogSDeterminant(dsdp,&dsdp->logdet);DSDPCHKERR(info); in DSDPResetY0() 340 info=DSDPSetY(dsdp,1.0,dsdp->logdet,dsdp->y);DSDPCHKERR(info); in DSDPResetY0() [all …]
|
/dports/math/casadi/casadi-3.5.5/external_packages/DSDP/src/solver/ |
H A D | dualimpl.c | 261 int DSDPComputePotential(DSDP dsdp, DSDPVec y, double logdet, double *potential){ in DSDPComputePotential() argument 267 if (gap>0) dpotential=dsdp->rho*log(gap)-logdet; in DSDPComputePotential() 270 DSDPLogInfo(0,9,"Gap: %4.4e, Log Determinant: %4.4e, Log Gap: %4.4e\n",gap,logdet,log(gap)); in DSDPComputePotential() 287 int DSDPComputePotential2(DSDP dsdp, DSDPVec y, double mu, double logdet, double *potential){ in DSDPComputePotential2() argument 292 *potential=-(ddobj + mu*logdet)*dsdp->schurmu; in DSDPComputePotential2() 293 *potential=-(ddobj/mu + logdet)*dsdp->schurmu; in DSDPComputePotential2() 309 int DSDPSetY(DSDP dsdp, double beta, double logdet, DSDPVec ynew){ in DSDPSetY() argument 330 dsdp->logdet=logdet; in DSDPSetY() 331 info=DSDPComputePotential(dsdp,dsdp->y,dsdp->logdet,&dsdp->potential);DSDPCHKERR(info); in DSDPSetY()
|
H A D | dualalg.c | 27 double dstep,newpotential,logdet; in DSDPYStepLineSearch() local 42 info=DSDPComputeLogSDeterminant(dsdp,&logdet);DSDPCHKERR(info); in DSDPYStepLineSearch() 55 info=DSDPSetY(dsdp,dstep,logdet,dsdp->ytemp);DSDPCHKERR(info); in DSDPYStepLineSearch() 57 info=DSDPSetY(dsdp,0,dsdp->logdet,dsdp->y);DSDPCHKERR(info); in DSDPYStepLineSearch() 77 double dstep,newpotential,bdotdy,oldpotential,logdet; in DSDPYStepLineSearch2() local 93 info=DSDPComputeLogSDeterminant(dsdp,&logdet);DSDPCHKERR(info); in DSDPYStepLineSearch2() 107 info=DSDPSetY(dsdp,dstep,logdet,dsdp->ytemp);DSDPCHKERR(info); in DSDPYStepLineSearch2() 109 info=DSDPSetY(dsdp,0,dsdp->logdet,dsdp->y);DSDPCHKERR(info); in DSDPYStepLineSearch2() 339 info=DSDPComputeLogSDeterminant(dsdp,&dsdp->logdet);DSDPCHKERR(info); in DSDPResetY0() 340 info=DSDPSetY(dsdp,1.0,dsdp->logdet,dsdp->y);DSDPCHKERR(info); in DSDPResetY0() [all …]
|
/dports/science/R-cran-eco/eco/src/ |
H A D | subroutines.c | 285 double logdet=0.0; in ddet() local 290 logdet += log(pdTemp[i][i]); in ddet() 294 return(2.0*logdet); in ddet() 296 return(exp(2.0*logdet)); in ddet() 305 double logdet=0.0; in ddet2D() local 310 logdet += log(pdTemp[i][i]); in ddet2D() 314 return(2.0*logdet); in ddet2D() 316 return(exp(2.0*logdet)); in ddet2D()
|
/dports/sysutils/istio/istio-1.6.7/vendor/github.com/openshift/api/vendor/gonum.org/v1/gonum/lapack/testlapack/ |
H A D | dlasq1.go | 77 var logdet float64 80 logdet += math.Log(math.Abs(v)) 82 if math.Exp(logdet) > 2 { 83 …("Incorrect singular value. n = %d, cas = %d, elem = %d, det = %v", n, cas, elem, math.Exp(logdet))
|
H A D | dsterf.go | 116 var logdet float64 119 logdet += math.Log(math.Abs(v)) 121 if math.Exp(logdet) > 2 { 122 t.Errorf("Incorrect singular value. n = %d, cas = %d, det = %v", n, cas, math.Exp(logdet))
|
/dports/math/cppad/CppAD-20210000.8/include/cppad/speed/ |
H A D | det_by_lu.hpp | 146 Scalar logdet; in operator ()() local 157 n_, m_, A_, B_, X_, logdet); in operator ()() 169 det = Scalar( signdet ) * exp( logdet ); in operator ()()
|
/dports/math/octave-forge-mvn/mvn/inst/ |
H A D | mvn_div_bc.m | 30 % 1/2*log(det(covsum)) - 1/4*m1.logdet - 1/4*m2.logdet; 33 1/2*log(det(covsum)) - 1/4*m1.logdet - 1/4*m2.logdet;
|
/dports/science/qmcpack/qmcpack-3.11.0/src/QMCWaveFunctions/BsplineFactory/ |
H A D | BsplineSet.h | 126 ValueMatrix_t& logdet, in evaluate_notranspose() argument 134 ValueVector_t v(logdet[i], OrbitalSetSize); in evaluate_notranspose() 193 ValueMatrix_t& logdet, in evaluate_notranspose() argument 199 ValueVector_t v(logdet[i], OrbitalSetSize); in evaluate_notranspose() 209 ValueMatrix_t& logdet, in evaluate_notranspose() argument 216 ValueVector_t v(logdet[i], OrbitalSetSize); in evaluate_notranspose()
|