| /petsc/src/vec/vec/tests/ |
| H A D | ex30f.F90 | 17 Vec v1, v2, v3 35 PetscCallA(VecCreateSeq(PETSC_COMM_SELF, n, v3, ierr)) 61 PetscCallA(VecScale(v3, num, ierr)) 84 PetscCallA(VecScatterCreate(v1, fromis, v3, tois, scat2, ierr)) 91 PetscCallA(VecScatterBegin(scat2, v1, v3, ADD_VALUES, SCATTER_FORWARD, ierr)) 92 PetscCallA(VecScatterEnd(scat2, v1, v3, ADD_VALUES, SCATTER_FORWARD, ierr)) 101 PetscCallA(PetscObjectSetName(v3, 'V3', ierr)) 102 PetscCallA(VecView(v3, PETSC_VIEWER_STDOUT_SELF, ierr)) 109 PetscCallA(VecDestroy(v3, ierr))
|
| /petsc/src/mat/tests/output/ |
| H A D | ex99_1.out | 3 MatMultAdd v2==v3: norm difference < sqrt(eps_machine) 4 MatMultAdd v2!=v3: norm difference < sqrt(eps_machine) 5 MatMultTransposeAdd v2==v3: norm difference < sqrt(eps_machine) 6 MatMultTransposeAdd v2!=v3: norm difference < sqrt(eps_machine)
|
| H A D | ex97_1.out | 3 MatMultAdd v2==v3: norm difference < sqrt(eps_machine) 4 MatMultAdd v2!=v3: norm difference < sqrt(eps_machine) 5 MatMultTransposeAdd v2==v3: norm difference < sqrt(eps_machine) 6 MatMultTransposeAdd v2!=v3: norm difference < sqrt(eps_machine)
|
| /petsc/doc/install/ |
| H A D | download.md | 79 Releases (for example: 3.20.0, 3.21.0, etc. with corresponding Git tags v3.20.0, v3.21.0, etc): 84 Patch updates (for example: 3.21.1, 2.21.2, etc. with corresponding Git tags v3.21.1, v3.21.2, etc) 94 - May (features added since v3.21.0) main branch 95 - May (bug fixes since v3.21.1) release branch 96 - April end (bug fixes since v3.21.0) v3.21.1 tag and tarball 97 - March end (features added after v3.20.0) v3.21.0 tag and tarball 98 - March end (bug fixes since v3.20.5) v3.20.6 tag and tarball 100 - October end (bug fixes since v3.20.0) v3.20.1 tag and tarball 101 - September end (features added after v3.19.0) v3.20.0 tag and tarball
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | inode.c | 374 const MatScalar *v1, *v2, *v3, *v4, *v5; in MatMult_SeqAIJ_Inode() local 379 #pragma disjoint(*x, *y, *v1, *v2, *v3, *v4, *v5) in MatMult_SeqAIJ_Inode() 443 v3 = v2 + n; in MatMult_SeqAIJ_Inode() 455 sum3 += v3[0] * tmp0 + v3[1] * tmp1; in MatMult_SeqAIJ_Inode() 456 v3 += 2; in MatMult_SeqAIJ_Inode() 462 sum3 += *v3++ * tmp0; in MatMult_SeqAIJ_Inode() 467 v1 = v3; /* Since the next block to be processed starts there*/ in MatMult_SeqAIJ_Inode() 476 v3 = v2 + n; in MatMult_SeqAIJ_Inode() 477 v4 = v3 + n; in MatMult_SeqAIJ_Inode() 489 sum3 += v3[0] * tmp0 + v3[1] * tmp1; in MatMult_SeqAIJ_Inode() [all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex96.c | 48 Vec x, v1, v2, v3, v4; in main() local 229 PetscCall(VecCreate(PETSC_COMM_WORLD, &v3)); in main() 230 PetscCall(VecSetSizes(v3, n, PETSC_DECIDE)); in main() 231 PetscCall(VecSetFromOptions(v3)); in main() 232 PetscCall(VecDuplicate(v3, &v4)); in main() 240 PetscCall(MatMultTranspose(P, v2, v3)); /* v3 = Pt*A*P*x */ in main() 242 PetscCall(VecAXPY(v4, none, v3)); in main() 244 PetscCall(VecNorm(v3, NORM_1, &norm_tmp1)); in main() 251 PetscCall(VecDestroy(&v3)); in main()
|
| H A D | ex94.c | 26 Vec x, v1, v2, v3, v4; in main() local 394 PetscCall(VecCreate(PETSC_COMM_WORLD, &v3)); in main() 395 PetscCall(VecSetSizes(v3, n, PETSC_DECIDE)); in main() 396 PetscCall(VecSetFromOptions(v3)); in main() 397 PetscCall(VecDuplicate(v3, &v4)); in main() 405 PetscCall(MatMultTranspose(P, v2, v3)); /* v3 = Pt*A*P*x */ in main() 408 PetscCall(VecAXPY(v4, none, v3)); in main() 419 PetscCall(VecDestroy(&v3)); in main()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex74ftube.h | 1 double precision p1, p2, p3, p4, r1, r2, r3, r4 double precision e1, e2, e3, e4, v1, v2, v3, v4 dou… 6 …p2, p3, p4, r1, r2, r3, r4, e1, e2, e3, e4, v1, v2, v3, v4, ru1, ru2, ru3, ru4, et1, et2, et3, et4…
|
| /petsc/src/mat/impls/submat/ |
| H A D | submat.c | 73 static PetscErrorCode MatMultAdd_SubMatrix(Mat N, Vec v1, Vec v2, Vec v3) in MatMultAdd_SubMatrix() argument 83 } else if (v2 == v3) { in MatMultAdd_SubMatrix() 98 PetscCall(VecScatterBegin(Na->lrestrict, Na->lwork, v3, INSERT_VALUES, SCATTER_FORWARD)); in MatMultAdd_SubMatrix() 99 PetscCall(VecScatterEnd(Na->lrestrict, Na->lwork, v3, INSERT_VALUES, SCATTER_FORWARD)); in MatMultAdd_SubMatrix() 117 static PetscErrorCode MatMultTransposeAdd_SubMatrix(Mat N, Vec v1, Vec v2, Vec v3) in MatMultTransposeAdd_SubMatrix() argument 127 } else if (v2 == v3) { in MatMultTransposeAdd_SubMatrix() 142 PetscCall(VecScatterBegin(Na->rprolong, Na->rwork, v3, INSERT_VALUES, SCATTER_REVERSE)); in MatMultTransposeAdd_SubMatrix() 143 PetscCall(VecScatterEnd(Na->rprolong, Na->rwork, v3, INSERT_VALUES, SCATTER_REVERSE)); in MatMultTransposeAdd_SubMatrix()
|
| /petsc/src/ksp/ksp/utils/lmvm/tests/ |
| H A D | lmvm_copy_test.c | 123 Vec u, x, f, x_copy, f_copy, v1, v2, v3; in testMatLMVMCopy() local 154 PetscCall(testUnchangedBegin(rand, A_copy, &v1, &v2, &v3)); in testMatLMVMCopy() 157 PetscCall(testUnchangedEnd(rand, A_copy, &v1, &v2, &v3, &equal)); in testMatLMVMCopy() 160 PetscCall(testUnchangedBegin(rand, A, &v1, &v2, &v3)); in testMatLMVMCopy() 163 PetscCall(testUnchangedEnd(rand, A, &v1, &v2, &v3, &equal)); in testMatLMVMCopy()
|
| /petsc/src/mat/impls/cdiagonal/ |
| H A D | cdiagonal.c | 52 static PetscErrorCode MatMultAdd_ConstantDiagonal(Mat mat, Vec v1, Vec v2, Vec v3) in MatMultAdd_ConstantDiagonal() argument 57 if (v2 == v3) { in MatMultAdd_ConstantDiagonal() 58 PetscCall(VecAXPBY(v3, ctx->diag, 1.0, v1)); in MatMultAdd_ConstantDiagonal() 60 PetscCall(VecAXPBYPCZ(v3, ctx->diag, 1.0, 0.0, v1, v2)); in MatMultAdd_ConstantDiagonal() 65 static PetscErrorCode MatMultHermitianTransposeAdd_ConstantDiagonal(Mat mat, Vec v1, Vec v2, Vec v3) in MatMultHermitianTransposeAdd_ConstantDiagonal() argument 70 if (v2 == v3) { in MatMultHermitianTransposeAdd_ConstantDiagonal() 71 PetscCall(VecAXPBY(v3, PetscConj(ctx->diag), 1.0, v1)); in MatMultHermitianTransposeAdd_ConstantDiagonal() 73 PetscCall(VecAXPBYPCZ(v3, PetscConj(ctx->diag), 1.0, 0.0, v1, v2)); in MatMultHermitianTransposeAdd_ConstantDiagonal()
|
| /petsc/src/snes/tutorials/network/water/ |
| H A D | waterreaddata.c | 136 double v1, v2, v3, v4, v5, v6; in WaterReadData() local 228 …ndata = sscanf(line, "%d %lf %lf %lf %lf %lf %lf %d", &id, &v1, &v2, &v3, &v4, &v5, &v6, &… in WaterReadData() 234 tank->minlvl = (PetscScalar)v3; in WaterReadData() 249 …= sscanf(line, "%d %d %d %lf %lf %lf %lf %s", &id, &node1, &node2, &v1, &v2, &v3, &v4, pipe->stat); in WaterReadData() 255 pipe->roughness = (PetscScalar)v3; in WaterReadData()
|
| /petsc/src/dm/impls/moab/tests/output/ |
| H A D | ex2_1.out | 6 Field{2} = v3.
|
| /petsc/src/mat/utils/ |
| H A D | multequal.c | 440 Vec x, v1, v2, v3, v4, Cx, Bx; in MatProjMultEqual_Private() local 465 PetscCall(VecDuplicate(x, &v3)); in MatProjMultEqual_Private() 485 PetscCall(MatMult(B, v1, v3)); /* v3 = R*A*R^t*x */ in MatProjMultEqual_Private() 487 PetscCall(MatMultTranspose(B, v1, v3)); /* v3 = Bt*A*B*x */ in MatProjMultEqual_Private() 490 PetscCall(VecAXPY(v4, -1.0, v3)); in MatProjMultEqual_Private() 507 PetscCall(VecDestroy(&v3)); in MatProjMultEqual_Private()
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_vec.py | 321 v3 = PETSc.Vec().createWithArray(a, size=3, comm=PETSc.COMM_SELF) 325 self.assertTrue(v3.size == 3) 331 a3 = v3.getDict()['__array__'] 339 v3 = PETSc.Vec().createWithArray(a, size=(3, None), comm=PETSc.COMM_WORLD) 343 self.assertTrue(v3.local_size == 3) 349 a3 = v3.getDict()['__array__']
|
| /petsc/src/mat/impls/diagonal/ |
| H A D | diagonal.c | 89 static PetscErrorCode MatMultAdd_Diagonal(Mat mat, Vec v1, Vec v2, Vec v3) in MatMultAdd_Diagonal() argument 95 if (v2 != v3) { in MatMultAdd_Diagonal() 96 PetscCall(VecPointwiseMult(v3, ctx->diag, v1)); in MatMultAdd_Diagonal() 97 PetscCall(VecAXPY(v3, 1.0, v2)); in MatMultAdd_Diagonal() 100 PetscCall(VecDuplicate(v3, &w)); in MatMultAdd_Diagonal() 102 PetscCall(VecAXPY(v3, 1.0, w)); in MatMultAdd_Diagonal()
|
| /petsc/src/vec/vec/utils/ |
| H A D | projection.c | 1042 const PetscScalar *v1, *v2, *v3; in VecMedian() local 1094 PetscCall(VecGetArrayRead(Vec3, &v3)); in VecMedian() 1096 v3 = vmed; in VecMedian() 1099 …), PetscMin(PetscRealPart(v1[i]), PetscRealPart(v3[i]))), PetscMin(PetscRealPart(v2[i]), PetscReal… in VecMedian() 1104 if (VMedian != Vec3) PetscCall(VecRestoreArrayRead(Vec3, &v3)); in VecMedian()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact9.c | 886 __m256d a0, a1, a2, a3, a4, a5, w0, w1, w2, w3, s0, s1, s2, v0, v1, v2, v3; in MatSolve_SeqBAIJ_9_NaturalOrdering() local 1118 v3 = _mm256_set1_pd(ls[3]); in MatSolve_SeqBAIJ_9_NaturalOrdering() 1120 w0 = _mm256_fmadd_pd(a3, v3, w0); in MatSolve_SeqBAIJ_9_NaturalOrdering() 1122 w1 = _mm256_fmadd_pd(a4, v3, w1); in MatSolve_SeqBAIJ_9_NaturalOrdering() 1124 w2 = _mm256_fmadd_pd(a5, v3, w2); in MatSolve_SeqBAIJ_9_NaturalOrdering() 1154 v3 = _mm256_set1_pd(ls[7]); in MatSolve_SeqBAIJ_9_NaturalOrdering() 1156 w0 = _mm256_fmadd_pd(a3, v3, w0); in MatSolve_SeqBAIJ_9_NaturalOrdering() 1158 w1 = _mm256_fmadd_pd(a4, v3, w1); in MatSolve_SeqBAIJ_9_NaturalOrdering() 1160 w2 = _mm256_fmadd_pd(a5, v3, w2); in MatSolve_SeqBAIJ_9_NaturalOrdering()
|
| /petsc/src/mat/interface/ |
| H A D | matrix.c | 2835 PetscErrorCode MatMultAdd(Mat mat, Vec v1, Vec v2, Vec v3) in MatMultAdd() argument 2842 PetscValidHeaderSpecific(v3, VEC_CLASSID, 4); in MatMultAdd() 2849 …->rmap->n == v3->map->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Mat mat,Vec v3: local dim %" PetscIn… in MatMultAdd() 2851 …PetscCheck(v1 != v3, PetscObjectComm((PetscObject)mat), PETSC_ERR_ARG_IDN, "v1 and v3 must be diff… in MatMultAdd() 2854 PetscCall(PetscLogEventBegin(MAT_MultAdd, mat, v1, v2, v3)); in MatMultAdd() 2856 PetscUseTypeMethod(mat, multadd, v1, v2, v3); in MatMultAdd() 2858 PetscCall(PetscLogEventEnd(MAT_MultAdd, mat, v1, v2, v3)); in MatMultAdd() 2859 PetscCall(PetscObjectStateIncrease((PetscObject)v3)); in MatMultAdd() 2884 PetscErrorCode MatMultTransposeAdd(Mat mat, Vec v1, Vec v2, Vec v3) in MatMultTransposeAdd() argument 2893 PetscValidHeaderSpecific(v3, VEC_CLASSID, 4); in MatMultTransposeAdd() [all …]
|
| /petsc/doc/changes/ |
| H A D | 313.md | 130 accepts an optional threshold parameter (since v3.10) and 137 - Add SNESLineSearchGetType() in v3.13.1.
|
| H A D | 233.md | 159 - fftw interface is updated to use v3.2alpha2
|
| H A D | 323.md | 107 - v3.23.3: Add `TSGetIJacobianP()`
|
| /petsc/src/ksp/ksp/tests/benchmarkscatters/ |
| H A D | Baseline-MPICH-8 | 3 Using Petsc Development GIT revision: v3.6.2-1540-geddb79b GIT Date: 2015-11-11 14:39:31 -0600
|
| /petsc/src/snes/tutorials/output/ |
| H A D | ex19_tut_6.out | 10 Using Petsc Development GIT revision: v3.9.3-921-gfc2aa81 GIT Date: 2018-07-27 11:07:58 +0200
|
| /petsc/doc/manual/ |
| H A D | dmnetwork.md | 62 subnetwork 0: v0 — v1 — v2 — v3
|