| /petsc/src/mat/tests/ |
| H A D | ex74.c | 10 Mat sA, sB, sFactor, B, C; /* symmetric matrices */ in main() local
183 PetscCall(MatDuplicate(sA, MAT_COPY_VALUES, &sB)); in main()
184 PetscCall(MatEqual(sA, sB, &equal)); in main()
189 PetscCall(MatNorm(sB, NORM_FROBENIUS, &norm2)); in main()
193 PetscCall(MatNorm(sB, NORM_INFINITY, &norm2)); in main()
197 PetscCall(MatNorm(sB, NORM_1, &norm2)); in main()
203 PetscCall(MatGetInfo(sB, MAT_LOCAL, &minfo2)); in main()
211 PetscCall(MatGetSize(sB, &i, &j)); in main()
215 PetscCall(MatGetBlockSize(sB, &i)); in main()
232 PetscCall(MatDiagonalScale(sB, x, x)); in main()
[all …]
|
| H A D | ex75.c | 8 Mat A, sA, sB; in main() local
261 PetscCall(MatDuplicate(sA, MAT_COPY_VALUES, &sB)); in main()
262 PetscCall(MatEqual(sA, sB, &flg)); in main()
264 PetscCall(MatMultEqual(sA, sB, 5, &flg)); in main()
269 PetscCall(MatMultAddEqual(sA, sB, 5, &flg)); in main()
274 PetscCall(MatDestroy(&sB)); in main()
|
| /petsc/src/dm/dt/space/impls/tensor/ |
| H A D | spacetensor.c | 305 PetscReal *lpoints, *sB = NULL, *sD = NULL, *sH = NULL; in PetscSpaceEvaluate_Tensor() local
317 if (B || D || H) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc, MPIU_REAL, &sB)); in PetscSpaceEvaluate_Tensor()
343 PetscCall(PetscSpaceEvaluate(tens->tensspaces[s], npoints, lpoints, sB, sD, sH)); in PetscSpaceEvaluate_Tensor()
355 …for (PetscInt p = 0; p < npoints; p++) B[(pdim * p + i) * Nc + c] *= sB[(spdim * p + si) * sNc + s… in PetscSpaceEvaluate_Tensor()
378 …tscInt p = 0; p < npoints; p++) D[((pdim * p + i) * Nc + c) * Nv + der] *= sB[(spdim * p + si) * s… in PetscSpaceEvaluate_Tensor()
408 …p < npoints; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sB[(spdim * p + si) * s… in PetscSpaceEvaluate_Tensor()
425 if (B || D || H) PetscCall(DMRestoreWorkArray(dm, npoints * pdim * Nc, MPIU_REAL, &sB)); in PetscSpaceEvaluate_Tensor()
|
| /petsc/src/dm/dt/space/impls/sum/ |
| H A D | spacesum.c | 424 PetscReal *sB = NULL, *sD = NULL, *sH = NULL; in PetscSpaceEvaluate_Sum() local
436 if (B || D || H) PetscCall(DMGetWorkArray(dm, numelB, MPIU_REAL, &sB)); in PetscSpaceEvaluate_Sum()
453 … 0 || !sum->uniform) PetscCall(PetscSpaceEvaluate(sum->sumspaces[s], npoints, points, sB, sD, sH)); in PetscSpaceEvaluate_Sum()
468 if (B) B[BInd] = sB[sBInd]; in PetscSpaceEvaluate_Sum()
501 if (B || D || H) PetscCall(DMRestoreWorkArray(dm, numelB, MPIU_REAL, &sB)); in PetscSpaceEvaluate_Sum()
|
| /petsc/src/dm/field/impls/ds/ |
| H A D | dmfieldds.c | 499 PetscScalar *sB = (PetscScalar *)B; in DMFieldEvaluateFV_DS() local
505 for (j = 0; j < Nc; j++) sB[i * Nc + j] = 0.; in DMFieldEvaluateFV_DS()
508 …for (j = 0; j < Nc; j++) sB[i * Nc + j] += geom->detJ[i * Nq + k] * weights[k] * sqB[(i * Nq + k) … in DMFieldEvaluateFV_DS()
510 for (k = 0; k < Nc; k++) sB[i * Nc + k] /= vol; in DMFieldEvaluateFV_DS()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plextree.c | 84 PetscInt size, i, sA = -1, sB, sOrientB, sConeSize; in DMPlexReferenceTreeGetChildSymmetry_Default() local
102 …enceTreeGetChildSymmetry_Default(dm, parent, parentOrientA, 0, sA, parentOrientB, &sOrientB, &sB)); in DMPlexReferenceTreeGetChildSymmetry_Default()
105 PetscCall(DMPlexGetCone(dm, sB, &coneB)); in DMPlexReferenceTreeGetChildSymmetry_Default()
107 PetscCall(DMPlexGetConeOrientation(dm, sB, &oB)); in DMPlexReferenceTreeGetChildSymmetry_Default()
|
| /petsc/src/dm/impls/forest/p4est/ ! |
| H A D | pforest.h | 1943 PetscInt size, i, sA = -1, sB, sOrientB, sConeSize; in DMReferenceTreeGetChildSymmetry_pforest() local
1961 …enceTreeGetChildSymmetry_pforest(dm, parent, parentOrientA, 0, sA, parentOrientB, &sOrientB, &sB)); in DMReferenceTreeGetChildSymmetry_pforest()
1964 PetscCall(DMPlexGetCone(dm, sB, &coneB)); in DMReferenceTreeGetChildSymmetry_pforest()
1966 PetscCall(DMPlexGetConeOrientation(dm, sB, &oB)); in DMReferenceTreeGetChildSymmetry_pforest()
|