Home
last modified time | relevance | path

Searched refs:logdet (Results 1 – 25 of 190) sorted by relevance

12345678

/dports/math/py-spvcm/spvcm-0.3.0/spvcm/both_levels/generic/
H A Dsample.py24 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 DSPOSetProxyForMSD.cpp100 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 DSPOSetProxy.h62 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 DSPOSetProxyForMSD.h70 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 DSPOSetProxy.cpp66 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 Dlu_solve.hpp271 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 Dlu_solve.hpp272 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 DTMV_SymLDLDecompose.cpp278 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 DTMV_SymCHD.cpp49 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 DTMV_SymBandCHD.cpp50 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 DLCAOrbitalSetWithCorrection.cpp65 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 DLCAOrbitalSet.h116 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 DLCAOrbitalSetWithCorrection.h62 ValueMatrix_t& logdet,
69 ValueMatrix_t& logdet,
76 ValueMatrix_t& logdet,
H A DLCAOrbitalSet.cpp391 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 Dsample.py17 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 Ddualimpl.c261 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 Ddualalg.c27 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 Ddualimpl.c261 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 Ddualalg.c27 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 Dsubroutines.c285 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 Ddlasq1.go77 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 Ddsterf.go116 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 Ddet_by_lu.hpp146 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 Dmvn_div_bc.m30 % 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 DBsplineSet.h126 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()

12345678