| /petsc/src/dm/dt/interface/ |
| H A D | dtaltv.c | 117 PetscReal prod; in PetscDTAltVApply() local
120 prod = permOdd ? -1. : 1.; in PetscDTAltVApply()
121 for (l = 0; l < k; l++) prod *= v[perm[l] * N + subset[l]]; in PetscDTAltVApply()
122 subsum += prod; in PetscDTAltVApply()
411 PetscReal prod; in PetscDTAltVPullback() local
415 prod = isOdd ? -ww[i] : ww[i]; in PetscDTAltVPullback()
416 for (l = 0; l < k; l++) prod *= L[subsetw[perm[l]] * N + subsetv[l]]; in PetscDTAltVPullback()
417 Lstarw[j] += prod; in PetscDTAltVPullback()
529 PetscReal prod; in PetscDTAltVPullbackMatrix() local
534 prod = isOdd ? -1. : 1.; in PetscDTAltVPullbackMatrix()
[all …]
|
| /petsc/src/tao/leastsquares/impls/pounders/ |
| H A D | gqt.c | 220 …PetscReal alpha, anorm, bnorm, parc, parf, parl, pars, par = *retpar, paru, prod, rxnorm, rznor… in gqt() local
323 PetscCallBLAS("BLASdot", prod = BLASdot_(&blasn, z, &blas1, x, &blas1) / delta); in gqt()
325 alpha = temp / (PetscAbs(prod) + PetscSqrtScalar(prod * prod + temp / delta)); in gqt()
326 if (prod >= 0) alpha = PetscAbs(alpha); in gqt()
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex13.c | 164 Mat prod; in test() local
231 … PetscCall(MatTransposeMatMult(mat_k_moment_T, AinvB, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &prod)); in test()
232 PetscCall(MatAXPY(prod, -1., mat_koszul, SAME_NONZERO_PATTERN)); in test()
233 PetscCall(MatNorm(prod, NORM_FROBENIUS, &frob_err)); in test()
236 PetscCall(MatDestroy(&prod)); in test()
|
| /petsc/src/mat/impls/aij/seq/seqviennacl/ |
| H A D | aijviennacl.cxx | 189 *ygpu = viennacl::linalg::prod(*viennaclstruct->compressed_mat, *xgpu); in MatMult_SeqAIJViennaCL()
191 *ygpu = viennacl::linalg::prod(*viennaclstruct->mat, *xgpu); in MatMult_SeqAIJViennaCL()
223 …if (a->compressedrow.use) *viennaclstruct->tempvec = viennacl::linalg::prod(*viennaclstruct->compr… in MatMultAdd_SeqAIJViennaCL()
224 else *viennaclstruct->tempvec = viennacl::linalg::prod(*viennaclstruct->mat, *xgpu); in MatMultAdd_SeqAIJViennaCL()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex7.c | 462 for (PetscInt d = 0, prod = 1; d < dim; prod *= extent[d], ++d) { in TestFreeField() local
463 const PetscInt idx = (v / prod) % extent[d]; in TestFreeField()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex11.h | 136 PetscReal prod = 0.0; in DotDIMReal() local
138 for (i = 0; i < DIM; i++) prod += x[i] * y[i]; in DotDIMReal()
139 return prod; in DotDIMReal()
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/ |
| H A D | dspacelagrange.c | 1055 Mat prod; in MatTensorAltV() local
1082 PetscCall(MatCreateSeqAIJ(PETSC_COMM_SELF, m, n, 0, nnz, &prod)); in MatTensorAltV()
1083 PetscCall(PetscObjectSetOptionsPrefix((PetscObject)prod, "altv_")); in MatTensorAltV()
1087 PetscCall(MatSetOption(prod, MAT_IGNORE_ZERO_ENTRIES, PETSC_FALSE)); in MatTensorAltV()
1165 …for (l = 0; l < Nk; l++) PetscCall(MatSetValue(prod, i, col * Nk + l, vals[l], INSERT_VALUES)); /*… in MatTensorAltV()
1176 PetscCall(MatAssemblyBegin(prod, MAT_FINAL_ASSEMBLY)); in MatTensorAltV()
1177 PetscCall(MatAssemblyEnd(prod, MAT_FINAL_ASSEMBLY)); in MatTensorAltV()
1178 *product = prod; in MatTensorAltV()
|
| /petsc/src/dm/impls/forest/p4est/ |
| H A D | pforest.h | 3635 PetscReal prod = 1.; in DMPforestMapCoordinates_Cell() local
3637 for (k = 0; k < P4EST_DIM; k++) prod *= (j & (1 << k)) ? eta[k] : (1. - eta[k]); in DMPforestMapCoordinates_Cell()
3638 for (k = 0; k < 3; k++) coordGuess[k] += prod * treeCoords[j][k]; in DMPforestMapCoordinates_Cell()
3653 PetscReal prod = 1.; in DMPforestMapCoordinates_Cell() local
3657 prod *= (k & (1 << l)) ? 1. : -1.; in DMPforestMapCoordinates_Cell()
3659 prod *= (k & (1 << l)) ? eta[l] : (1. - eta[l]); in DMPforestMapCoordinates_Cell()
3662 for (l = 0; l < 3; l++) dXdeta[l] += prod * treeCoords[k][l]; in DMPforestMapCoordinates_Cell()
3678 PetscReal prod = 1.; in DMPforestMapCoordinates_Cell() local
3680 for (l = 0; l < P4EST_DIM; l++) prod *= (k & (1 << l)) ? eta[l] : (1. - eta[l]); in DMPforestMapCoordinates_Cell()
3681 for (l = 0; l < 3; l++) coordGuess[l] += prod * treeCoords[k][l]; in DMPforestMapCoordinates_Cell()
|
| /petsc/share/petsc/matlab/ |
| H A D | laplacian.m | 167 maxeigs = prod(u);
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexcreate.c | 3151 PetscReal prod = 0.0; in DiffNormReal() local
3153 for (i = 0; i < dim; ++i) prod += PetscSqr(x[i] - y[i]); in DiffNormReal()
3154 return PetscSqrtReal(prod); in DiffNormReal()
3159 PetscReal prod = 0.0; in DotReal() local
3161 for (i = 0; i < dim; ++i) prod += x[i] * y[i]; in DotReal()
3162 return prod; in DotReal()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 471 (ignition/prod-fuel-cut 0.99)
2553 (dpm/subtet/max-turb-disp-face-normal-dot-prod -0.9)
8040 (coeff-ke-mu-prod 1.)
8077 (ke-prod-mod #f)
8134 (nut-prod-const 2.)
8135 (nut-prod-modify? #f)
|