Home
last modified time | relevance | path

Searched refs:v2 (Results 1 – 25 of 75) sorted by relevance

123

/petsc/src/vec/is/sf/tests/
H A Dex1f.F9012 Vec v1, v2
28 PetscCallA(VecScatterCreateToAll(v1, toall, v2, ierr))
30 PetscCallA(VecScatterBegin(toall, v1, v2, INSERT_VALUES, SCATTER_FORWARD, ierr))
31 PetscCallA(VecScatterEnd(toall, v1, v2, INSERT_VALUES, SCATTER_FORWARD, ierr))
35 PetscCallA(VecDestroy(v2, ierr))
37 PetscCallA(VecScatterCreateToAll(v1, toall, v2, ierr))
38 PetscCallA(VecScatterBegin(toall, v1, v2, INSERT_VALUES, SCATTER_FORWARD, ierr))
39 PetscCallA(VecScatterEnd(toall, v1, v2, INSERT_VALUES, SCATTER_FORWARD, ierr))
42 PetscCallA(PetscObjectSetName(v2, 'v2', ierr))
43 PetscCallA(VecView(v2, PETSC_VIEWER_STDOUT_SELF, ierr))
[all …]
/petsc/src/mat/impls/aij/seq/
H A Dinode.c374 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()
415 v2 = v1 + n; in MatMult_SeqAIJ_Inode()
425 sum2 += v2[0] * tmp0 + v2[1] * tmp1; in MatMult_SeqAIJ_Inode()
426 v2 += 2; in MatMult_SeqAIJ_Inode()
431 sum2 += *v2++ * tmp0; in MatMult_SeqAIJ_Inode()
435 v1 = v2; /* Since the next block to be processed starts there*/ in MatMult_SeqAIJ_Inode()
442 v2 = v1 + n; in MatMult_SeqAIJ_Inode()
443 v3 = v2 + n; in MatMult_SeqAIJ_Inode()
453 sum2 += v2[0] * tmp0 + v2[1] * tmp1; in MatMult_SeqAIJ_Inode()
[all …]
/petsc/src/vec/vec/tests/
H A Dex30f.F9017 Vec v1, v2, v3
34 PetscCallA(VecCreateFromOptions(PETSC_COMM_WORLD, PETSC_NULL_CHARACTER, i1, nlocal, n, v2, ierr))
60 PetscCallA(VecScale(v2, num, ierr))
83 PetscCallA(VecScatterCreate(v1, fromis, v2, tois, scat1, ierr))
88 PetscCallA(VecScatterBegin(scat1, v1, v2, ADD_VALUES, SCATTER_FORWARD, ierr))
89 PetscCallA(VecScatterEnd(scat1, v1, v2, ADD_VALUES, SCATTER_FORWARD, ierr))
97 PetscCallA(PetscObjectSetName(v2, 'V2', ierr))
98 PetscCallA(VecView(v2, PETSC_VIEWER_STDOUT_WORLD, ierr))
108 PetscCallA(VecDestroy(v2, ierr))
/petsc/src/mat/tests/output/
H A Dex97_1.out3 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 Dex99_1.out3 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/src/vec/vec/utils/
H A Dprojection.c28 const PetscScalar *v1, *v2; in VecWhichEqual() local
41 v2 = v1; in VecWhichEqual()
44 PetscCall(VecGetArrayRead(Vec2, &v2)); in VecWhichEqual()
50 if (v1[i] == v2[i]) { in VecWhichEqual()
60 PetscCall(VecRestoreArrayRead(Vec2, &v2)); in VecWhichEqual()
94 const PetscScalar *v1, *v2; in VecWhichLessThan() local
107 v2 = v1; in VecWhichLessThan()
110 PetscCall(VecGetArrayRead(Vec2, &v2)); in VecWhichLessThan()
116 if (PetscRealPart(v1[i]) < PetscRealPart(v2[i])) { in VecWhichLessThan()
126 PetscCall(VecRestoreArrayRead(Vec2, &v2)); in VecWhichLessThan()
[all …]
/petsc/src/mat/tests/
H A Dex232.c11 PetscBool v2 = PETSC_FALSE; in main() local
15 PetscCall(PetscOptionsGetBool(NULL, NULL, "-v2", &v2, NULL)); in main()
19 if (!v2) { in main()
27 if (v2) { in main()
H A Dex96.c48 Vec x, v1, v2, v3, v4; in main() local
147 PetscCall(VecDuplicate(v1, &v2)); in main()
184 PetscCall(MatMult(A_tmp, v1, v2)); /* v2 = A*P*x */ in main()
186 PetscCall(VecAXPY(v1, none, v2)); in main()
188 PetscCall(VecNorm(v2, NORM_1, &norm_tmp1)); in main()
238 PetscCall(MatMult(A, v1, v2)); /* v2 = A*P*x */ in main()
240 PetscCall(MatMultTranspose(P, v2, v3)); /* v3 = Pt*A*P*x */ in main()
260 PetscCall(VecDestroy(&v2)); in main()
/petsc/src/mat/impls/submat/
H A Dsubmat.c73 static PetscErrorCode MatMultAdd_SubMatrix(Mat N, Vec v1, Vec v2, Vec v3) in MatMultAdd_SubMatrix() argument
81 if (v1 == v2) { in MatMultAdd_SubMatrix()
83 } else if (v2 == v3) { in MatMultAdd_SubMatrix()
85 PetscCall(VecScatterBegin(Na->lrestrict, v2, Na->lwork, INSERT_VALUES, SCATTER_REVERSE)); in MatMultAdd_SubMatrix()
86 PetscCall(VecScatterEnd(Na->lrestrict, v2, Na->lwork, INSERT_VALUES, SCATTER_REVERSE)); in MatMultAdd_SubMatrix()
94 PetscCall(VecScatterBegin(Na->lrestrict, v2, Na->lwork2, INSERT_VALUES, SCATTER_REVERSE)); in MatMultAdd_SubMatrix()
95 PetscCall(VecScatterEnd(Na->lrestrict, v2, Na->lwork2, INSERT_VALUES, SCATTER_REVERSE)); in MatMultAdd_SubMatrix()
117 static PetscErrorCode MatMultTransposeAdd_SubMatrix(Mat N, Vec v1, Vec v2, Vec v3) in MatMultTransposeAdd_SubMatrix() argument
125 if (v1 == v2) { in MatMultTransposeAdd_SubMatrix()
127 } else if (v2 == v3) { in MatMultTransposeAdd_SubMatrix()
[all …]
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A DScatter.pyx119 >>> v2 = PETSc.Vec().createWithArray([0, 0, 0])
120 >>> sct = PETSc.Scatter().create(v1,None,v2,None)
121 >>> sct.scatter(v1,v2) # v2[:] = v1[:]
122 >>> sct.scatter(v2,v1,mode=revmode) # v1[:] = v2[:]
126 >>> v2 = PETSc.Vec().createWithArray([0, 0])
128 >>> sct = PETSc.Scatter().create(v1,is1,v2,None)
129 >>> sct.scatter(v1,v2) # v2[:] = v1[3:0:-2]
130 >>> sct.scatter(v2,v1,mode=revmode) # v1[3:0:-2] = v2[:]
/petsc/src/vec/vec/tutorials/
H A Dex20f90.F9018 PetscReal norm, v, v1, v2, tol
171 v2 = abs(norm - 3.0*sqrt(nfloat))
172 if (v2 > -tol .and. v2 < tol) v2 = 0.0
173 if (rank == 0) write (6, 220) v, v1, v2
H A Dex1f90.F9017 PetscReal norm, v, v1, v2
169 v2 = abs(norm - 3.0*sqrt(nfloat))
170 if (v2 > -1.d-10 .and. v2 < 1.d-10) v2 = 0.0
171 if (rank == 0) write (6, 220) v, v1, v2
H A Dex1.c16 PetscReal norm, v, v1, v2, maxval; in main() local
150 v2 = norm - 3.0 * PetscSqrtReal((PetscReal)n); in main()
151 if (v2 > -PETSC_SMALL && v2 < PETSC_SMALL) v2 = 0.0; in main()
152 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "VecMAXPY %g %g %g \n", (double)v, (double)v1, (double)v2)… in main()
/petsc/src/binding/petsc4py/test/
H A Dtest_optdb.py131 for v1, v2 in zip(arrayv, av):
133 self.assertEqual(v1.item(), pyt(v2))
138 for v1, v2 in zip(arrayv, av):
140 self.assertEqual(v1.item(), pyt(v2))
145 for v1, v2 in zip(arrayv, defarrayv):
147 self.assertEqual(v1.item(), pyt(v2))
H A Dtest_vec.py85 v2 = v1.duplicate()
87 v2.set(2)
88 v1.swap(v2)
91 self.assertEqual(v2[idx], 1)
214 w2, v2 = w[e - 1], v[e - 1]
216 self.assertEqual(w2, v2)
320 v2 = PETSc.Vec().createWithArray(a, size=5, comm=PETSc.COMM_SELF)
324 self.assertTrue(v2.size == 5)
329 a2 = v2.getDict()['__array__']
338 v2 = PETSc.Vec().createWithArray(a, size=(5, None), comm=PETSc.COMM_WORLD)
[all …]
H A Dtest_nsp.py8 for v1, v2 in zip(seq1, seq2):
9 if abs(v1 - v2) > 1e-5:
/petsc/src/mat/tutorials/
H A Dex9.c15 Vec x, y, v, v2, z, z2; in main() local
43 PetscCall(VecDuplicate(v, &v2)); in main()
123 PetscCall(MatMult(B, x, v2)); in main()
126 PetscCall(VecAXPY(v2, -1.0, v)); in main()
127 PetscCall(VecNorm(v2, NORM_2, &rnorm)); in main()
149 PetscCall(VecDestroy(&v2)); in main()
/petsc/src/snes/tutorials/
H A Dex74ftube.h1 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…
H A Dex59.c37 PetscScalar v2; in main() local
96 v2 = a * xp + (b - a) * PetscPowScalar(xp, k); in main()
97 PetscCall(VecSetValues(x, 1, &i, &v2, INSERT_VALUES)); in main()
104 v2 = xx[i] * sperturb; in main()
105 PetscCall(VecSetValues(x, 1, &i, &v2, INSERT_VALUES)); in main()
H A Dex8.c42 const PetscScalar v2 = v1 / (1 + ux2 + uy2); in hessian_11_2d() local
44 jac[0] = v1 - v2 * ux2; in hessian_11_2d()
45 jac[1] = -v2 * uxy; in hessian_11_2d()
46 jac[2] = -v2 * uxy; in hessian_11_2d()
47 jac[3] = v1 - v2 * uy2; in hessian_11_2d()
/petsc/src/snes/tutorials/network/water/
H A Dwaterreaddata.c136 double v1, v2, v3, v4, v5, v6; in WaterReadData() local
194 ndata = sscanf(line, "%d %lf %lf %d", &id, &v1, &v2, &pattern); in WaterReadData()
199 junction->demand = (PetscScalar)v2; in WaterReadData()
228 …ndata = sscanf(line, "%d %lf %lf %lf %lf %lf %lf %d", &id, &v1, &v2, &v3, &v4, &v5, &v6, &… in WaterReadData()
233 tank->initlvl = (PetscScalar)v2; in WaterReadData()
249 …ndata = sscanf(line, "%d %d %d %lf %lf %lf %lf %s", &id, &node1, &node2, &v1, &v2, &v3, … in WaterReadData()
254 pipe->diam = (PetscScalar)v2; in WaterReadData()
291 ndata = sscanf(line, "%d %lf %lf", &icurve_id, &v1, &v2); in WaterReadData()
295 curve_y = (PetscScalar)v2; in WaterReadData()
/petsc/src/sys/tests/linter/output/
H A DtestReferences.patch4 -PetscErrorCode testOutOfLineReference(PetscViewer v, PetscViewer v2)
5 +static PetscErrorCode testOutOfLineReference(PetscViewer v, PetscViewer v2)
/petsc/src/ksp/ksp/utils/lmvm/tests/
H A Dlmvm_copy_test.c123 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/vec/is/sf/tests/output/
H A Dex1f_1.out1 Vec Object: v2 1 MPI process
/petsc/src/mat/impls/cdiagonal/
H A Dcdiagonal.c52 static PetscErrorCode MatMultAdd_ConstantDiagonal(Mat mat, Vec v1, Vec v2, Vec v3) in MatMultAdd_ConstantDiagonal() argument
57 if (v2 == v3) { 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()
73 PetscCall(VecAXPBYPCZ(v3, PetscConj(ctx->diag), 1.0, 0.0, v1, v2)); in MatMultHermitianTransposeAdd_ConstantDiagonal()

123