/dports/math/libpgmath/flang-d07daf3/test/ncar_kernels/HOMME_grad_sphere/src/ |
H A D | kernel_gradient_sphere.F90 | 88 Dinv2(i,j,1,1,:) = Dinv(1,1,i,j,:) 89 Dinv2(i,j,1,2,:) = Dinv(1,2,i,j,:) 90 Dinv2(i,j,2,1,:) = Dinv(2,1,i,j,:) 91 Dinv2(i,j,2,2,:) = Dinv(2,2,i,j,:) 95 elem(ie)%dinv = Dinv(:,:,:,:,ie) 260 ds(i,j,1)=Dinv(1,1,i,j)*v1(i,j) + Dinv(2,1,i,j)*v2(i,j) 261 ds(i,j,2)=Dinv(1,2,i,j)*v1(i,j) + Dinv(2,2,i,j)*v2(i,j) 303 ds(i,j,1)=Dinv(1,1,i,j)*v1(i,j) + Dinv(2,1,i,j)*v2(i,j) 304 ds(i,j,2)=Dinv(1,2,i,j)*v1(i,j) + Dinv(2,2,i,j)*v2(i,j) 348 ds(i,j,1)=Dinv(i,j,1,1)*v1(i,j) + Dinv(i,j,2,1)*v2(i,j) [all …]
|
/dports/math/libpgmath/flang-d07daf3/test/ncar_kernels/HOMME_div_sphere/src/ |
H A D | kernel_divergence_sphere.F90 | 233 DinvTemp(:,:,1,1) = elem%Dinv(1,1,:,:) 234 DinvTemp(:,:,1,2) = elem%Dinv(1,2,:,:) 235 DinvTemp(:,:,2,1) = elem%Dinv(2,1,:,:) 236 DinvTemp(:,:,2,2) = elem%Dinv(2,2,:,:) 273 DinvTemp(:,:,1,1) = elem%Dinv(1,1,:,:) 274 DinvTemp(:,:,1,2) = elem%Dinv(1,2,:,:) 275 DinvTemp(:,:,2,1) = elem%Dinv(2,1,:,:) 334 … gv(i,j,1)=elem%metdet(i,j)*(elem%Dinv(1,1,i,j)*v(i,j,1) + elem%Dinv(1,2,i,j)*v(i,j,2)) 335 … gv(i,j,2)=elem%metdet(i,j)*(elem%Dinv(2,1,i,j)*v(i,j,1) + elem%Dinv(2,2,i,j)*v(i,j,2)) 387 … gv(i,j,1)=elem%metdet(i,j)*(elem%Dinv(1,1,i,j)*v(i,j,1) + elem%Dinv(1,2,i,j)*v(i,j,2)) [all …]
|
/dports/biology/ncbi-toolkit/ncbi/algo/blast/composition_adjustment/ |
H A D | optimize_target_freq.c | 398 free((*newton_system)->Dinv); in ReNewtonSystemFree() 399 (*newton_system)->Dinv = NULL; in ReNewtonSystemFree() 428 newton_system->Dinv = NULL; in ReNewtonSystemNew() 434 newton_system->Dinv = in ReNewtonSystemNew() 436 if (newton_system->Dinv == NULL) in ReNewtonSystemNew() 484 double * Dinv = newton_system->Dinv; in FactorReNewtonSystem() local 509 Dinv[i] = x[i] / (1 - eta); in FactorReNewtonSystem() 512 memcpy(Dinv, x, n * sizeof(double)); in FactorReNewtonSystem() 561 double * Dinv = newton_system->Dinv; in SolveReNewtonSystem() local 576 workspace[i] = x[i] * Dinv[i]; in SolveReNewtonSystem() [all …]
|
/dports/biology/ncbi-blast+/ncbi-blast-2.12.0+-src/c++/src/algo/blast/composition_adjustment/ |
H A D | optimize_target_freq.c | 398 free((*newton_system)->Dinv); in ReNewtonSystemFree() 399 (*newton_system)->Dinv = NULL; in ReNewtonSystemFree() 428 newton_system->Dinv = NULL; in ReNewtonSystemNew() 434 newton_system->Dinv = in ReNewtonSystemNew() 436 if (newton_system->Dinv == NULL) in ReNewtonSystemNew() 484 double * Dinv = newton_system->Dinv; in FactorReNewtonSystem() local 509 Dinv[i] = x[i] / (1 - eta); in FactorReNewtonSystem() 512 memcpy(Dinv, x, n * sizeof(double)); in FactorReNewtonSystem() 561 double * Dinv = newton_system->Dinv; in SolveReNewtonSystem() local 576 workspace[i] = x[i] * Dinv[i]; in SolveReNewtonSystem() [all …]
|
/dports/biology/ncbi-cxx-toolkit/ncbi_cxx--25_2_0/src/algo/blast/composition_adjustment/ |
H A D | optimize_target_freq.c | 398 free((*newton_system)->Dinv); in ReNewtonSystemFree() 399 (*newton_system)->Dinv = NULL; in ReNewtonSystemFree() 428 newton_system->Dinv = NULL; in ReNewtonSystemNew() 434 newton_system->Dinv = in ReNewtonSystemNew() 436 if (newton_system->Dinv == NULL) in ReNewtonSystemNew() 484 double * Dinv = newton_system->Dinv; in FactorReNewtonSystem() local 509 Dinv[i] = x[i] / (1 - eta); in FactorReNewtonSystem() 512 memcpy(Dinv, x, n * sizeof(double)); in FactorReNewtonSystem() 561 double * Dinv = newton_system->Dinv; in SolveReNewtonSystem() local 576 workspace[i] = x[i] * Dinv[i]; in SolveReNewtonSystem() [all …]
|
/dports/math/osqp/osqp-0.6.2/lin_sys/direct/qdldl/qdldl_sources/tests/ |
H A D | qdldl_tester.c | 75 QDLDL_float *Dinv; in ldl_factor_solve() local 103 Dinv = (QDLDL_float*)malloc(sizeof(QDLDL_float)*An); in ldl_factor_solve() 122 free(Lp); free(D); free(Dinv); free(etree); free(Lnz); in ldl_factor_solve() 136 factorStatus = QDLDL_factor(An,Ap,Ai,Ax,Lp,Li,Lx,D,Dinv,Lnz,etree,bwork,iwork,fwork); in ldl_factor_solve() 140 free(Lp); free(D); free(Dinv); free(etree); free(Lnz); in ldl_factor_solve() 148 QDLDL_solve(Ln,Lp,Li,Lx,Dinv,b); in ldl_factor_solve() 154 free(Lp); free(D); free(Dinv); free(etree); free(Lnz); in ldl_factor_solve()
|
/dports/math/py-qdldl/qdldl-0.1.5.post0/c/qdldl/tests/ |
H A D | qdldl_tester.c | 75 QDLDL_float *Dinv; in ldl_factor_solve() local 103 Dinv = (QDLDL_float*)malloc(sizeof(QDLDL_float)*An); in ldl_factor_solve() 122 free(Lp); free(D); free(Dinv); free(etree); free(Lnz); in ldl_factor_solve() 136 factorStatus = QDLDL_factor(An,Ap,Ai,Ax,Lp,Li,Lx,D,Dinv,Lnz,etree,bwork,iwork,fwork); in ldl_factor_solve() 140 free(Lp); free(D); free(Dinv); free(etree); free(Lnz); in ldl_factor_solve() 148 QDLDL_solve(Ln,Lp,Li,Lx,Dinv,b); in ldl_factor_solve() 154 free(Lp); free(D); free(Dinv); free(etree); free(Lnz); in ldl_factor_solve()
|
/dports/math/py-osqp/osqp-0.6.2.post0/osqp_sources/lin_sys/direct/qdldl/qdldl_sources/tests/ |
H A D | qdldl_tester.c | 75 QDLDL_float *Dinv; in ldl_factor_solve() local 103 Dinv = (QDLDL_float*)malloc(sizeof(QDLDL_float)*An); in ldl_factor_solve() 122 free(Lp); free(D); free(Dinv); free(etree); free(Lnz); in ldl_factor_solve() 136 factorStatus = QDLDL_factor(An,Ap,Ai,Ax,Lp,Li,Lx,D,Dinv,Lnz,etree,bwork,iwork,fwork); in ldl_factor_solve() 140 free(Lp); free(D); free(Dinv); free(etree); free(Lnz); in ldl_factor_solve() 148 QDLDL_solve(Ln,Lp,Li,Lx,Dinv,b); in ldl_factor_solve() 154 free(Lp); free(D); free(Dinv); free(etree); free(Lnz); in ldl_factor_solve()
|
/dports/devel/godot/godot-3.2.3-stable/thirdparty/opus/silk/float/ |
H A D | solve_LS_FLP.c | 44 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… 83 silk_float Dinv[ MAX_MATRIX_SIZE ]; /* inverse diagonal elements of D*/ in silk_solve_LDL_FLP() local 91 silk_LDL_FLP( A, M, &L[ 0 ][ 0 ], Dinv ); in silk_solve_LDL_FLP() 104 T[ i ] = T[ i ] * Dinv[ i ]; in silk_solve_LDL_FLP() 160 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… in silk_LDL_FLP() argument 190 Dinv[ j ] = ( silk_float )( 1.0f / temp ); in silk_LDL_FLP() 200 matrix_ptr( L, i, j, M ) = ( silk_float )( ( ptr1[ i ] - temp ) * Dinv[ j ] ); in silk_LDL_FLP()
|
/dports/devel/godot-tools/godot-3.2.3-stable/thirdparty/opus/silk/float/ |
H A D | solve_LS_FLP.c | 44 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… 83 silk_float Dinv[ MAX_MATRIX_SIZE ]; /* inverse diagonal elements of D*/ in silk_solve_LDL_FLP() local 91 silk_LDL_FLP( A, M, &L[ 0 ][ 0 ], Dinv ); in silk_solve_LDL_FLP() 104 T[ i ] = T[ i ] * Dinv[ i ]; in silk_solve_LDL_FLP() 160 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… in silk_LDL_FLP() argument 190 Dinv[ j ] = ( silk_float )( 1.0f / temp ); in silk_LDL_FLP() 200 matrix_ptr( L, i, j, M ) = ( silk_float )( ( ptr1[ i ] - temp ) * Dinv[ j ] ); in silk_LDL_FLP()
|
/dports/devel/godot2-tools/godot-2.1.6-stable/thirdparty/opus/silk/float/ |
H A D | solve_LS_FLP.c | 44 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… 83 silk_float Dinv[ MAX_MATRIX_SIZE ]; /* inverse diagonal elements of D*/ in silk_solve_LDL_FLP() local 91 silk_LDL_FLP( A, M, &L[ 0 ][ 0 ], Dinv ); in silk_solve_LDL_FLP() 104 T[ i ] = T[ i ] * Dinv[ i ]; in silk_solve_LDL_FLP() 160 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… in silk_LDL_FLP() argument 190 Dinv[ j ] = ( silk_float )( 1.0f / temp ); in silk_LDL_FLP() 200 matrix_ptr( L, i, j, M ) = ( silk_float )( ( ptr1[ i ] - temp ) * Dinv[ j ] ); in silk_LDL_FLP()
|
/dports/devel/godot2/godot-2.1.6-stable/thirdparty/opus/silk/float/ |
H A D | solve_LS_FLP.c | 44 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… 83 silk_float Dinv[ MAX_MATRIX_SIZE ]; /* inverse diagonal elements of D*/ in silk_solve_LDL_FLP() local 91 silk_LDL_FLP( A, M, &L[ 0 ][ 0 ], Dinv ); in silk_solve_LDL_FLP() 104 T[ i ] = T[ i ] * Dinv[ i ]; in silk_solve_LDL_FLP() 160 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… in silk_LDL_FLP() argument 190 Dinv[ j ] = ( silk_float )( 1.0f / temp ); in silk_LDL_FLP() 200 matrix_ptr( L, i, j, M ) = ( silk_float )( ( ptr1[ i ] - temp ) * Dinv[ j ] ); in silk_LDL_FLP()
|
/dports/www/firefox-legacy/firefox-52.8.0esr/media/libopus/silk/float/ |
H A D | solve_LS_FLP.c | 44 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… 83 silk_float Dinv[ MAX_MATRIX_SIZE ]; /* inverse diagonal elements of D*/ in silk_solve_LDL_FLP() local 91 silk_LDL_FLP( A, M, &L[ 0 ][ 0 ], Dinv ); in silk_solve_LDL_FLP() 104 T[ i ] = T[ i ] * Dinv[ i ]; in silk_solve_LDL_FLP() 160 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… in silk_LDL_FLP() argument 190 Dinv[ j ] = ( silk_float )( 1.0f / temp ); in silk_LDL_FLP() 200 matrix_ptr( L, i, j, M ) = ( silk_float )( ( ptr1[ i ] - temp ) * Dinv[ j ] ); in silk_LDL_FLP()
|
/dports/audio/sdl2_mixer/SDL2_mixer-2.0.4/external/opus-1.0.3/silk/float/ |
H A D | solve_LS_FLP.c | 44 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… 83 silk_float Dinv[ MAX_MATRIX_SIZE ]; /* inverse diagonal elements of D*/ in silk_solve_LDL_FLP() local 91 silk_LDL_FLP( A, M, &L[ 0 ][ 0 ], Dinv ); in silk_solve_LDL_FLP() 104 T[ i ] = T[ i ] * Dinv[ i ]; in silk_solve_LDL_FLP() 160 …silk_float *Dinv /* I/O Pointer to vector holding the inverse diagonal elements of… in silk_LDL_FLP() argument 190 Dinv[ j ] = ( silk_float )( 1.0f / temp ); in silk_LDL_FLP() 200 matrix_ptr( L, i, j, M ) = ( silk_float )( ( ptr1[ i ] - temp ) * Dinv[ j ] ); in silk_LDL_FLP()
|
/dports/math/osqp/osqp-0.6.2/lin_sys/direct/qdldl/qdldl_sources/julia/ |
H A D | QDLDL.jl | 15 Dinv::Array{Tf} field 26 (L, Dinv) = factor(triu(A)); 32 (L, Dinv) = factor(triu(A[p,p])); 57 (L, Dinv) = factor(Atr,logical) 91 QDLDL.Dinv,b) 123 Dinv = Array{Tv}(A.n) 128 Dinv[:] = 1 140 return L, Dinv 225 Dinv[1] = 1/D[1]; 343 Dinv[k]= 1/D[k] [all …]
|
/dports/math/fflas-ffpack/fflas-ffpack-2.4.3/fflas-ffpack/ffpack/ |
H A D | ffpack_fsytrf.inl | 40 typename Field::Element_ptr Dinvj = Dinv; 48 F.inv (Dinv[i*incDinv], *Ai); 53 typename Field::Element_ptr Dinvj = Dinv; 61 F.inv (Dinv[i*incDinv], *Ai); 72 typename Field::Element_ptr Ai = A, Dinvi = Dinv; 155 typename Field::Element_ptr bj = A+row, Dinvj=Dinv; 195 bj = A+row+i, Dinvj=Dinv; 697 FFLAS::fflas_delete(Dinv); 712 FFLAS::fflas_delete(Dinv); 729 // FFLAS::WriteMatrix(std::cerr<<"After fsytrf_nonunit Dinv = "<<std::endl,F,1,N,Dinv, 1); [all …]
|
/dports/math/osqp/osqp-0.6.2/lin_sys/direct/qdldl/qdldl_sources/examples/c/ |
H A D | example.c | 28 QDLDL_float *Dinv; in main() local 59 Dinv = (QDLDL_float*)malloc(sizeof(QDLDL_float)*An); in main() 84 QDLDL_factor(An,Ap,Ai,Ax,Lp,Li,Lx,D,Dinv,Lnz,etree,bwork,iwork,fwork); in main() 93 QDLDL_solve(Ln,Lp,Li,Lx,Dinv,x); in main() 122 print_arrayf(Dinv, An, "diag(D^{-1})"); in main() 139 free(Dinv); in main()
|
/dports/math/py-qdldl/qdldl-0.1.5.post0/c/qdldl/examples/ |
H A D | example.c | 28 QDLDL_float *Dinv; in main() local 59 Dinv = (QDLDL_float*)malloc(sizeof(QDLDL_float)*An); in main() 84 QDLDL_factor(An,Ap,Ai,Ax,Lp,Li,Lx,D,Dinv,Lnz,etree,bwork,iwork,fwork); in main() 93 QDLDL_solve(Ln,Lp,Li,Lx,Dinv,x); in main() 122 print_arrayf(Dinv, An, "diag(D^{-1})"); in main() 139 free(Dinv); in main()
|
/dports/math/py-osqp/osqp-0.6.2.post0/osqp_sources/lin_sys/direct/qdldl/qdldl_sources/examples/ |
H A D | example.c | 28 QDLDL_float *Dinv; in main() local 59 Dinv = (QDLDL_float*)malloc(sizeof(QDLDL_float)*An); in main() 84 QDLDL_factor(An,Ap,Ai,Ax,Lp,Li,Lx,D,Dinv,Lnz,etree,bwork,iwork,fwork); in main() 93 QDLDL_solve(Ln,Lp,Li,Lx,Dinv,x); in main() 122 print_arrayf(Dinv, An, "diag(D^{-1})"); in main() 139 free(Dinv); in main()
|
/dports/math/py-hdmedians/hdmedians-0.14.2/hdmedians/ |
H A D | geomedian.pyx | 104 Dinv[i] = 1. / Di 106 Dinv[i] = 0. 109 Dinvs = sum(Dinv) 112 W[i] = Dinv[i] / Dinvs 184 Dinv[i] = 1. / Di 186 Dinv[i] = 0. 188 Dinvs = sum(Dinv) 265 Dinv[i] = nan 268 Dinvs = nansum(Dinv) 346 Dinv[i] = nan [all …]
|
/dports/math/py-qdldl/qdldl-0.1.5.post0/cpp/ |
H A D | qdldl.cpp | 22 Dinv = new QDLDL_float[n]; in Solver() 64 D, Dinv, Lnz, in Solver() 86 QDLDL_solve(nx, Lp, Li, Lx, Dinv, work); in solve() 104 D, Dinv, Lnz, in update() 115 delete [] Dinv; in ~Solver()
|
/dports/biology/diamond/diamond-2.0.13/src/stats/ |
H A D | matrix_adjust_eigen.cpp | 64 Dinv(N,N), in ReNewtonSystem() 72 MatrixN Dinv; /**< The diagonal elements of the member 208 MatrixN& Dinv = newton_system.Dinv; in FactorReNewtonSystem() local 230 Dinv = x; in FactorReNewtonSystem() 231 Dinv /= 1 - eta; in FactorReNewtonSystem() 235 ScaledSymmetricProductA(W, Dinv); in FactorReNewtonSystem() 247 workspace.row(i) = Dinv.row(i).cwiseProduct(grad_re.segment(i * N, N)); in FactorReNewtonSystem() 250 workspace.row(i) = Dinv.row(i).cwiseProduct(grad_re.segment<N>(i * N)); in FactorReNewtonSystem() 267 const MatrixN& Dinv = newton_system.Dinv; in SolveReNewtonSystem() local 277 workspace = x.cwiseProduct(Dinv); in SolveReNewtonSystem() [all …]
|
H A D | matrix_adjust.cpp | 446 free((*newton_system)->Dinv); in ReNewtonSystemFree() 447 (*newton_system)->Dinv = NULL; in ReNewtonSystemFree() 476 newton_system->Dinv = NULL; in ReNewtonSystemNew() 482 newton_system->Dinv = in ReNewtonSystemNew() 484 if (newton_system->Dinv == NULL) in ReNewtonSystemNew() 532 double* Dinv = newton_system->Dinv; in FactorReNewtonSystem() local 557 Dinv[i] = x[i] / (1 - eta); in FactorReNewtonSystem() 561 memcpy(Dinv, x, n * sizeof(double)); in FactorReNewtonSystem() 611 double* Dinv = newton_system->Dinv; in SolveReNewtonSystem() local 626 workspace[i] = x[i] * Dinv[i]; in SolveReNewtonSystem() [all …]
|
/dports/math/armadillo/armadillo-10.7.1/include/armadillo_bits/ |
H A D | op_wishrnd_meat.hpp | 237 Mat<eT> Dinv; in apply_noalias_mode1() local 243 const bool chol_status = op_chol::apply_direct(Dinv, Tinv, 0); in apply_noalias_mode1() 247 return op_iwishrnd::apply_noalias_mode2(out, Dinv, df); in apply_noalias_mode1() 255 op_iwishrnd::apply_noalias_mode2(Mat<eT>& out, const Mat<eT>& Dinv, const eT df) in apply_noalias_mode2() argument 260 arma_debug_check( (Dinv.is_square() == false), "iwishrnd(): given matrix must be square sized" ); in apply_noalias_mode2() 262 if(Dinv.is_empty()) { out.reset(); return true; } in apply_noalias_mode2() 266 const bool wishrnd_status = op_wishrnd::apply_noalias_mode2(tmp, Dinv, df); in apply_noalias_mode2()
|
/dports/math/R-cran-RcppArmadillo/RcppArmadillo/inst/include/armadillo_bits/ |
H A D | op_wishrnd_meat.hpp | 237 Mat<eT> Dinv; in apply_noalias_mode1() local 243 const bool chol_status = op_chol::apply_direct(Dinv, Tinv, 0); in apply_noalias_mode1() 247 return op_iwishrnd::apply_noalias_mode2(out, Dinv, df); in apply_noalias_mode1() 255 op_iwishrnd::apply_noalias_mode2(Mat<eT>& out, const Mat<eT>& Dinv, const eT df) in apply_noalias_mode2() argument 260 arma_debug_check( (Dinv.is_square() == false), "iwishrnd(): given matrix must be square sized" ); in apply_noalias_mode2() 262 if(Dinv.is_empty()) { out.reset(); return true; } in apply_noalias_mode2() 266 const bool wishrnd_status = op_wishrnd::apply_noalias_mode2(tmp, Dinv, df); in apply_noalias_mode2()
|