| /petsc/src/ksp/ksp/interface/ |
| H A D | itres.c | 35 PetscErrorCode KSPInitialResidual(KSP ksp, Vec vsoln, Vec vt1, Vec vt2, Vec vres, Vec vb) in KSPInitialResidual() argument
43 PetscValidHeaderSpecific(vb, VEC_CLASSID, 6); in KSPInitialResidual()
49 PetscCall(VecCopy(vb, vt2)); in KSPInitialResidual()
60 PetscCall(VecCopy(vb, vt2)); in KSPInitialResidual()
62 PetscCall(PCDiagonalScaleLeft(ksp->pc, vb, vres)); in KSPInitialResidual()
64 PetscCall(KSP_PCApply(ksp, vb, vres)); in KSPInitialResidual()
67 PetscCall(PCApplySymmetricLeft(ksp->pc, vb, vres)); in KSPInitialResidual()
|
| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | eptorsion3.c | 261 PetscReal zero = 0.0, vb, vl, vr, vt, dvdx, dvdy, flin = 0.0, fquad = 0.0; in FormFunction() local
291 vb = zero; in FormFunction()
294 if (i < nx && j > 0) vb = x[k - nx]; in FormFunction()
298 dvdy = (v - vb) / hy; in FormFunction()
300 flin = flin - cdiv3 * (vb + vl + v); in FormFunction()
333 PetscReal vb, vl, vr, vt, dvdx, dvdy; in FormGradient() local
378 vb = zero; in FormGradient()
381 if (i < nx && j > 0) vb = x[k - nx]; in FormGradient()
385 dvdy = (v - vb) / hy; in FormGradient()
537 PetscReal v, vb, vl, vr, vt, hxhx, hyhy; in HessianProduct() local
[all …]
|
| H A D | eptorsion1.c | 245 PetscReal zero = 0.0, vb, vl, vr, vt, dvdx, dvdy, flin = 0.0, fquad = 0.0; in FormFunction() local
275 vb = zero; in FormFunction()
278 if (i < nx && j > 0) vb = x[k - nx]; in FormFunction()
282 dvdy = (v - vb) / hy; in FormFunction()
284 flin = flin - cdiv3 * (vb + vl + v); in FormFunction()
317 PetscReal vb, vl, vr, vt, dvdx, dvdy; in FormGradient() local
362 vb = zero; in FormGradient()
365 if (i < nx && j > 0) vb = x[k - nx]; in FormGradient()
369 dvdy = (v - vb) / hy; in FormGradient()
521 PetscReal v, vb, vl, vr, vt, hxhx, hyhy; in HessianProduct() local
[all …]
|
| H A D | eptorsion2.c | 205 PetscReal v, vb, vl, vr, vt, dvdx, dvdy; in FormFunctionGradient() local
279 vb = zero; in FormFunctionGradient()
282 if (i < mx && j > 0) vb = x[k - gxm]; in FormFunctionGradient()
286 dvdy = (v - vb) / hy; in FormFunctionGradient()
303 flin -= cdiv3 * (vb + vl + v); in FormFunctionGradient()
|
| H A D | eptorsion2f.F90 | 119 PetscReal v, vb, vl, vr, vt, dvdx
213 vb = zero
216 if (i < mx .and. j > 0) vb = lx_v(k + 1 - gxm)
220 dvdy = (v - vb)/hy
237 flin = flin - cdiv3*(vb + vl + v)
|
| /petsc/src/mat/tests/ |
| H A D | ex47.c | 10 PetscViewer va, vb, vc; in main() local
33 PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD, file, FILE_MODE_READ, &vb)); in main()
36 PetscCall(MatLoad(B, vb)); in main()
37 PetscCall(PetscViewerDestroy(&vb)); in main()
|
| /petsc/src/ts/tests/ |
| H A D | ex28.c | 216 static PetscReal ComputeCDF(PetscReal m, PetscReal n, PetscReal T, PetscReal va, PetscReal vb) in ComputeCDF() argument
220 PetscReal zb = alpha * vb; in ComputeCDF()
|
| /petsc/src/mat/graphops/coarsen/impls/hem/ |
| H A D | hem.c | 400 PetscReal va = ((Edge *)a)->weight, vb = ((Edge *)b)->weight; in gamg_hem_compare() local
401 return (va <= vb - MY_MEPS) ? 1 : (va > vb + MY_MEPS) ? -1 : 0; /* 0 for equal */ in gamg_hem_compare()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 63 const PetscScalar *va, *vb; in MatGetRowMaxAbs_MPIBAIJ() local
81 PetscCall(VecGetArrayRead(vB, &vb)); in MatGetRowMaxAbs_MPIBAIJ()
83 if (PetscAbsScalar(va[i]) < PetscAbsScalar(vb[i])) { in MatGetRowMaxAbs_MPIBAIJ()
84 vv[i] = vb[i]; in MatGetRowMaxAbs_MPIBAIJ()
88 …if (idx && PetscAbsScalar(va[i]) == PetscAbsScalar(vb[i]) && idxb[i] != -1 && idx[i] > bs * a->gar… in MatGetRowMaxAbs_MPIBAIJ()
93 PetscCall(VecRestoreArrayRead(vB, &vb)); in MatGetRowMaxAbs_MPIBAIJ()
|
| /petsc/src/mat/impls/dense/seq/cupm/ |
| H A D | matseqdensecupm.hpp | 1356 const auto vb = DeviceArrayWrite(dctx, B); in Copy() local
1365 …PetscCall(PetscCUPMMemcpy2DAsync(vb.data(), lda_b, va.data(), lda_a, m, n, cupmMemcpyDeviceToDevic… in Copy()
1367 … PetscCall(PetscCUPMMemcpyAsync(vb.data(), va.data(), m * n, cupmMemcpyDeviceToDevice, stream)); in Copy()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aij.c | 2250 const MatScalar *va, *vb; in MatIsTranspose_SeqAIJ() local
2261 PetscCall(MatSeqAIJGetArrayRead(B, &vb)); in MatIsTranspose_SeqAIJ()
2280 vr = vb[bptr[idc]]; in MatIsTranspose_SeqAIJ()
2294 PetscCall(MatSeqAIJRestoreArrayRead(B, &vb)); in MatIsTranspose_SeqAIJ()
2302 MatScalar *va, *vb; in MatIsHermitianTranspose_SeqAIJ() local
2317 vb = bij->a; in MatIsHermitianTranspose_SeqAIJ()
2332 vr = vb[bptr[idc]]; in MatIsHermitianTranspose_SeqAIJ()
|
| /petsc/src/mat/impls/dense/seq/ |
| H A D | dense.c | 2619 PetscScalar *vb; in MatCopy_SeqDense() local
2630 PetscCall(MatDenseGetArray(B, &vb)); in MatCopy_SeqDense()
2632 for (j = 0; j < n; j++) PetscCall(PetscArraycpy(vb + j * lda2, va + j * lda1, m)); in MatCopy_SeqDense()
2634 PetscCall(PetscArraycpy(vb, va, A->rmap->n * A->cmap->n)); in MatCopy_SeqDense()
2636 PetscCall(MatDenseRestoreArray(B, &vb)); in MatCopy_SeqDense()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiaij.c | 2163 const PetscScalar *va, *vb; in MatGetRowMaxAbs_MPIAIJ() local
2181 PetscCall(VecGetArrayRead(vB, &vb)); in MatGetRowMaxAbs_MPIAIJ()
2183 if (PetscAbsScalar(va[i]) < PetscAbsScalar(vb[i])) { in MatGetRowMaxAbs_MPIAIJ()
2184 vv[i] = vb[i]; in MatGetRowMaxAbs_MPIAIJ()
2188 …if (idx && PetscAbsScalar(va[i]) == PetscAbsScalar(vb[i]) && idxb[i] != -1 && idx[i] > a->garray[i… in MatGetRowMaxAbs_MPIAIJ()
2193 PetscCall(VecRestoreArrayRead(vB, &vb)); in MatGetRowMaxAbs_MPIAIJ()
|