Home
last modified time | relevance | path

Searched refs:v1 (Results 1 – 25 of 84) sorted by relevance

1234

/petsc/src/vec/vec/tests/
H A Dex30f.F9017 Vec v1, v2, v3
33 …PetscCallA(VecCreateFromOptions(PETSC_COMM_WORLD, PETSC_NULL_CHARACTER, i1, nlocal2, n2, v1, ierr))
39 PetscCallA(VecSetValue(v1, row, num, INSERT_VALUES, ierr))
41 PetscCallA(VecSetValue(v1, row, num, INSERT_VALUES, ierr))
43 PetscCallA(VecSetValue(v1, row, num, INSERT_VALUES, ierr))
45 PetscCallA(VecSetValue(v1, row, num, INSERT_VALUES, ierr))
48 PetscCallA(VecSetValue(v1, row, num, INSERT_VALUES, ierr))
50 PetscCallA(VecSetValue(v1, row, num, INSERT_VALUES, ierr))
52 PetscCallA(VecSetValue(v1, row, num, INSERT_VALUES, ierr))
54 PetscCallA(VecSetValue(v1, row, num, INSERT_VALUES, ierr))
[all …]
H A Dbenchmark_veccreate.c12 PetscLogDouble v0, v1; in main() local
30 PetscCall(PetscTime(&v1)); in main()
32 …ll(PetscPrintf(PETSC_COMM_WORLD, "Iteration %" PetscInt_FMT ": Time= %g\n", i, (double)(v1 - v0))); in main()
/petsc/src/vec/is/sf/tests/
H A Dex1f.F9012 Vec v1, v2
20 …A(VecCreateFromOptions(PETSC_COMM_WORLD, PETSC_NULL_CHARACTER, i1, nlocal, PETSC_DECIDE, v1, ierr))
24 PetscCallA(VecSetValue(v1, row, num, INSERT_VALUES, ierr))
25 PetscCallA(VecAssemblyBegin(v1, ierr))
26 PetscCallA(VecAssemblyEnd(v1, ierr))
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))
37 PetscCallA(VecScatterCreateToAll(v1, toall, v2, ierr))
38 PetscCallA(VecScatterBegin(toall, v1, v2, INSERT_VALUES, SCATTER_FORWARD, ierr))
[all …]
/petsc/src/vec/vec/utils/
H A Dprojection.c28 const PetscScalar *v1, *v2; in VecWhichEqual() local
40 PetscCall(VecGetArrayRead(Vec1, &v1)); in VecWhichEqual()
41 v2 = v1; in VecWhichEqual()
43 PetscCall(VecGetArrayRead(Vec1, &v1)); in VecWhichEqual()
50 if (v1[i] == v2[i]) { in VecWhichEqual()
57 PetscCall(VecRestoreArrayRead(Vec1, &v1)); in VecWhichEqual()
59 PetscCall(VecRestoreArrayRead(Vec1, &v1)); in VecWhichEqual()
94 const PetscScalar *v1, *v2; in VecWhichLessThan() local
106 PetscCall(VecGetArrayRead(Vec1, &v1)); in VecWhichLessThan()
107 v2 = v1; in VecWhichLessThan()
[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()
387 v1 = &a->a[a->i[row]]; in MatMult_SeqAIJ_Inode()
389 …PetscPrefetchBlock(v1 + nsz * n, nsz * n, 0, PETSC_PREFETCH_HINT_NTA); /* Prefetch the values for … in MatMult_SeqAIJ_Inode()
402 sum1 += v1[0] * tmp0 + v1[1] * tmp1; in MatMult_SeqAIJ_Inode()
403 v1 += 2; in MatMult_SeqAIJ_Inode()
408 sum1 += *v1++ * tmp0; in MatMult_SeqAIJ_Inode()
415 v2 = v1 + n; in MatMult_SeqAIJ_Inode()
423 sum1 += v1[0] * tmp0 + v1[1] * tmp1; in MatMult_SeqAIJ_Inode()
424 v1 += 2; in MatMult_SeqAIJ_Inode()
[all …]
/petsc/src/binding/petsc4py/test/
H A Dtest_optdb.py131 for v1, v2 in zip(arrayv, av):
132 self.assertTrue(isinstance(v1.item(), pyt))
133 self.assertEqual(v1.item(), pyt(v2))
138 for v1, v2 in zip(arrayv, av):
139 self.assertTrue(isinstance(v1.item(), pyt))
140 self.assertEqual(v1.item(), pyt(v2))
145 for v1, v2 in zip(arrayv, defarrayv):
146 self.assertTrue(isinstance(v1.item(), pyt))
147 self.assertEqual(v1.item(), pyt(v2))
H A Dtest_vec.py84 v1 = self.vec
85 v2 = v1.duplicate()
86 v1.set(1)
88 v1.swap(v2)
90 self.assertEqual(v1[idx], 2)
213 w1, v1 = w[s], v[s]
215 self.assertEqual(w1, v1)
319 v1 = PETSc.Vec().createWithArray(a, comm=PETSc.COMM_SELF)
323 self.assertTrue(v1.size == 5)
327 a1 = v1.getDict()['__array__']
[all …]
/petsc/src/mat/tests/
H A Dex96.c48 Vec x, v1, v2, v3, v4; in main() local
143 PetscCall(VecCreate(PETSC_COMM_WORLD, &v1)); in main()
145 PetscCall(VecSetSizes(v1, m, PETSC_DECIDE)); in main()
146 PetscCall(VecSetFromOptions(v1)); in main()
147 PetscCall(VecDuplicate(v1, &v2)); in main()
183 PetscCall(MatMult(P, x, v1)); in main()
184 PetscCall(MatMult(A_tmp, v1, v2)); /* v2 = A*P*x */ in main()
185 PetscCall(MatMult(C, x, v1)); /* v1 = C*x */ in main()
186 PetscCall(VecAXPY(v1, none, v2)); in main()
187 PetscCall(VecNorm(v1, NORM_1, &norm_tmp)); in main()
[all …]
H A Dex94.c26 Vec x, v1, v2, v3, v4; in main() local
85 PetscCall(VecCreate(PETSC_COMM_WORLD, &v1)); in main()
87 PetscCall(VecSetSizes(v1, m, PETSC_DECIDE)); in main()
88 PetscCall(VecSetFromOptions(v1)); in main()
89 PetscCall(VecDuplicate(v1, &v2)); in main()
402 PetscCall(MatMult(P, x, v1)); in main()
403 PetscCall(MatMult(A, v1, v2)); /* v2 = A*P*x */ in main()
425 PetscCall(VecDestroy(&v1)); in main()
H A Dex18.c13 PetscScalar v, v0, v1, v2, a0 = 0.1, a, rhsval, *boundary_values, diag = 1.0; in main() local
70 v1 = -2.0 * a; in main()
75 v1 = -1.0 * a; in main()
80 PetscCall(MatSetValues(A, 1, &Ii, 1, &J, &v1, ADD_VALUES)); in main()
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A DScatter.pyx118 >>> v1 = PETSc.Vec().createWithArray([1, 2, 3])
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[:]
125 >>> v1 = PETSc.Vec().createWithArray([1, 2, 3, 4])
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 Dex1f90.F9017 PetscReal norm, v, v1, v2
166 v1 = abs(norm - 2.0*sqrt(nfloat))
167 if (v1 > -1.d-10 .and. v1 < 1.d-10) v1 = 0.0
171 if (rank == 0) write (6, 220) v, v1, v2
H A Dex20f90.F9018 PetscReal norm, v, v1, v2, tol
168 v1 = abs(norm - 2.0*sqrt(nfloat))
169 if (v1 > -tol .and. v1 < tol) v1 = 0.0
173 if (rank == 0) write (6, 220) v, v1, v2
H A Dex1.c16 PetscReal norm, v, v1, v2, maxval; in main() local
147 v1 = norm - 2.0 * PetscSqrtReal((PetscReal)n); in main()
148 if (v1 > -PETSC_SMALL && v1 < PETSC_SMALL) v1 = 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/snes/tutorials/output/
H A Dex28_glvis_composite_da_1d.out1 MFEM mesh v1.0
79 MFEM mesh v1.0
/petsc/src/ts/tutorials/output/
H A Dex12_glvis_da_2d_vect.out1 MFEM mesh v1.0
633 MFEM mesh v1.0
H A Dex12_glvis_da_2d_vect_ll.out1 MFEM mesh v1.0
633 MFEM mesh v1.0
/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()
198 junction->elev = (PetscScalar)v1; in WaterReadData()
212 ndata = sscanf(line, "%d %lf %d", &id, &v1, &pattern); in WaterReadData()
216 reservoir->head = (PetscScalar)v1; in WaterReadData()
228 …ndata = sscanf(line, "%d %lf %lf %lf %lf %lf %lf %d", &id, &v1, &v2, &v3, &v4, &v5, &v6, &… in WaterReadData()
232 tank->elev = (PetscScalar)v1; in WaterReadData()
249 …ndata = sscanf(line, "%d %d %d %lf %lf %lf %lf %s", &id, &node1, &node2, &v1, &v2, &v3, … in WaterReadData()
253 pipe->length = (PetscScalar)v1; in WaterReadData()
291 ndata = sscanf(line, "%d %lf %lf", &icurve_id, &v1, &v2); in WaterReadData()
[all …]
/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 Dex8.c41 const PetscScalar v1 = 1. / PetscSqrtReal(PetscAbsScalar(1 + ux2 + uy2)); in hessian_11_2d() local
42 const PetscScalar v2 = v1 / (1 + ux2 + uy2); in hessian_11_2d()
44 jac[0] = v1 - v2 * ux2; in hessian_11_2d()
47 jac[3] = v1 - v2 * uy2; in hessian_11_2d()
/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
79 PetscCall(VecScatterBegin(Na->rprolong, v1, Na->rwork, INSERT_VALUES, SCATTER_FORWARD)); in MatMultAdd_SubMatrix()
80 PetscCall(VecScatterEnd(Na->rprolong, v1, Na->rwork, INSERT_VALUES, SCATTER_FORWARD)); in MatMultAdd_SubMatrix()
81 if (v1 == v2) { in MatMultAdd_SubMatrix()
117 static PetscErrorCode MatMultTransposeAdd_SubMatrix(Mat N, Vec v1, Vec v2, Vec v3) in MatMultTransposeAdd_SubMatrix() argument
123 PetscCall(VecScatterBegin(Na->lrestrict, v1, Na->lwork, INSERT_VALUES, SCATTER_REVERSE)); in MatMultTransposeAdd_SubMatrix()
124 PetscCall(VecScatterEnd(Na->lrestrict, v1, Na->lwork, INSERT_VALUES, SCATTER_REVERSE)); in MatMultTransposeAdd_SubMatrix()
125 if (v1 == v2) { in MatMultTransposeAdd_SubMatrix()
/petsc/src/mat/utils/
H A Dmultequal.c440 Vec x, v1, v2, v3, v4, Cx, Bx; in MatProjMultEqual_Private() local
461 PetscCall(MatCreateVecs(B, &v1, &x)); in MatProjMultEqual_Private()
463 PetscCall(MatCreateVecs(B, &x, &v1)); in MatProjMultEqual_Private()
477 PetscCall(MatMultTranspose(B, x, v1)); in MatProjMultEqual_Private()
479 PetscCall(MatMult(B, x, v1)); in MatProjMultEqual_Private()
481 PetscCall(VecCopy(v1, Bx)); in MatProjMultEqual_Private()
483 PetscCall(VecCopy(v2, v1)); 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()
505 PetscCall(VecDestroy(&v1)); in MatProjMultEqual_Private()
/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/dm/impls/plex/tests/output/
H A Dex1_glvis_2d_quad_per_shift.out1 MFEM mesh v1.0
H A Dex1_glvis_2d_tet.out1 MFEM mesh v1.0

1234