Home
last modified time | relevance | path

Searched refs:vecs (Results 1 – 25 of 57) sorted by relevance

123

/petsc/src/ksp/ksp/utils/lmvm/ !
H A Dlmbasis.c28 basis->vecs = backing; in LMBasisCreate()
37 PetscValidLogicalCollectiveInt(basis->vecs, idx, 2); in LMBasisGetVec_Internal()
38 …PetscCheck(idx < basis->k, PetscObjectComm((PetscObject)basis->vecs), PETSC_ERR_ARG_OUTOFRANGE, "A… in LMBasisGetVec_Internal()
40 …PetscCheck(idx >= earliest, PetscObjectComm((PetscObject)basis->vecs), PETSC_ERR_ARG_OUTOFRANGE, "… in LMBasisGetVec_Internal()
44 PetscCall(MatDenseGetColumnVecRead(basis->vecs, idx % basis->m, single)); in LMBasisGetVec_Internal()
46 PetscCall(MatDenseGetColumnVecWrite(basis->vecs, idx % basis->m, single)); in LMBasisGetVec_Internal()
64 PetscCall(MatDenseRestoreColumnVecRead(basis->vecs, idx % basis->m, single)); in LMBasisRestoreVec()
66 PetscCall(MatDenseRestoreColumnVecWrite(basis->vecs, idx % basis->m, single)); in LMBasisRestoreVec()
118 …PetscCheck(basis->work_vecs_in_use == NULL, PetscObjectComm((PetscObject)basis->vecs), PETSC_ERR_A… in LMBasisDestroy()
130 …PetscCheck(basis->work_rows_in_use == NULL, PetscObjectComm((PetscObject)basis->vecs), PETSC_ERR_A… in LMBasisDestroy()
[all …]
H A Dlmproducts.c16 PetscValidHeaderSpecific(basis->vecs, MAT_CLASSID, 1); in LMProductsCreate()
17 …PetscCheck(block_type >= 0 && block_type < LMBLOCK_END, PetscObjectComm((PetscObject)basis->vecs),… in LMProductsCreate()
21 PetscCall(MatGetLocalSize(basis->vecs, NULL, &m_local)); in LMProductsCreate()
26 PetscCall(MatCreateVecs(basis->vecs, &(*dots)->diagonal_global, NULL)); in LMProductsCreate()
36 PetscCall(MatGetVecType(basis->vecs, &vec_type)); in LMProductsCreate()
37 …PetscCall(MatCreateDenseFromVecType(PetscObjectComm((PetscObject)basis->vecs), vec_type, m_local, … in LMProductsCreate()
96 MPI_Comm comm = PetscObjectComm((PetscObject)X->vecs); in LMProductsUpdate_Internal()
174 PetscCall(MatGetSize(X->vecs, &N, NULL)); in LMProductsUpdate_Internal()
H A Dlmbasis.h20 Mat vecs; // Dense matrix backing storage of vectors member
/petsc/src/mat/interface/ !
H A Dmatnull.c56 …rrorCode MatNullSpaceGetVecs(MatNullSpace sp, PetscBool *has_const, PetscInt *n, const Vec *vecs[]) in MatNullSpaceGetVecs() argument
62 if (vecs) *vecs = sp->vecs; in MatNullSpaceGetVecs()
196 for (i = 0; i < sp->n; i++) PetscCall(VecView(sp->vecs[i], viewer)); in MatNullSpaceView()
230 PetscErrorCode MatNullSpaceCreate(MPI_Comm comm, PetscBool has_cnst, PetscInt n, const Vec vecs[], … in MatNullSpaceCreate() argument
237 if (n) PetscAssertPointer(vecs, 4); in MatNullSpaceCreate()
238 for (i = 0; i < n; i++) PetscValidHeaderSpecific(vecs[i], VEC_CLASSID, 4); in MatNullSpaceCreate()
243 PetscCall(VecLockReadPush(vecs[i])); in MatNullSpaceCreate()
250 PetscCall(VecNorm(vecs[i], NORM_2, &norm)); in MatNullSpaceCreate()
251 …eal(norm - 1) <= PETSC_SQRT_MACHINE_EPSILON, PetscObjectComm((PetscObject)vecs[i]), PETSC_ERR_ARG_… in MatNullSpaceCreate()
256 PetscCall(VecSum(vecs[i], &sum)); in MatNullSpaceCreate()
[all …]
/petsc/src/vec/vec/tutorials/ !
H A Dex16f.F907 Vec, pointer, dimension(:) :: vecs
32 allocate (vecs(0:2))
34 vecs(0) = s
35 vecs(1) = r
49 PetscCallA(VecStrideGatherAll(v, vecs, INSERT_VALUES, ierr))
54 PetscCallA(VecStrideScatterAll(vecs, v, ADD_VALUES, ierr))
60 deallocate (vecs)
H A Dex16.c14 Vec v, s, r, vecs[2]; /* vectors */ in main() local
39 vecs[0] = s; in main()
40 vecs[1] = r; in main()
55 PetscCall(VecStrideGatherAll(v, vecs, INSERT_VALUES)); in main()
60 PetscCall(VecStrideScatterAll(vecs, v, ADD_VALUES)); in main()
/petsc/src/ksp/ksp/impls/gmres/pipefgmres/ !
H A Dpipefgmresimpl.h36 #define VEC_TEMP pipefgmres->vecs[0] /* work space */
37 #define VEC_TEMP_MATOP pipefgmres->vecs[1] /* work space */
38 #define VEC_Q pipefgmres->vecs[2] /* work space - Q pipelining var */
39 #define VEC_W pipefgmres->vecs[3] /* work space - W pipelining var */
41 #define VEC_VV(i) pipefgmres->vecs[VEC_OFFSET + i] /* use to access othog basis vectors. Note the o…
/petsc/src/mat/interface/ftn-custom/ !
H A Dzmatnullf.c13 …egetvecs_(MatNullSpace *sp, PetscBool *HAS_CNST, PetscInt *N, F90Array1d *vecs, PetscErrorCode *ie… in matnullspacegetvecs_() argument
24 *ierr = F90Array1dCreate(tvecs, MPIU_FORTRANADDR, 1, n, vecs PETSC_F90_2PTR_PARAM(ptrd)); in matnullspacegetvecs_()
27 …torevecs_(MatNullSpace *sp, PetscBool *HAS_CNST, PetscInt *N, F90Array1d *vecs, PetscErrorCode *ie… in matnullspacerestorevecs_() argument
29 *ierr = F90Array1dDestroy(vecs, MPIU_FORTRANADDR PETSC_F90_2PTR_PARAM(ptrd)); in matnullspacerestorevecs_()
/petsc/src/vec/vec/tests/ !
H A Dex28.c25 Vec x, y, vecs[40]; in main() local
135 PetscCall(VecCreate(PETSC_COMM_WORLD, vecs + i)); in main()
136 PetscCall(VecSetSizes(vecs[i], PETSC_DECIDE, n)); in main()
137 PetscCall(VecSetFromOptions(vecs[i])); in main()
139 PetscCall(VecSet(vecs[i], value)); in main()
141 for (i = 0; i < 39; i++) PetscCall(VecDotBegin(vecs[i], vecs[i + 1], results + i)); in main()
144 PetscCall(VecDotEnd(vecs[i], vecs[i + 1], results + i)); in main()
147 for (i = 0; i < 40; i++) PetscCall(VecDestroy(&vecs[i])); in main()
/petsc/src/ksp/ksp/impls/gmres/pgmres/ !
H A Dpgmresimpl.h20 #define VEC_TEMP pgmres->vecs[0] /* work space */
21 #define VEC_TEMP_MATOP pgmres->vecs[1] /* work space */
22 #define VEC_VV(i) pgmres->vecs[VEC_OFFSET + i] /* use to access othog basis vectors */
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/ !
H A Dpetscdmcomposite.pxi22 cdef PetscVec *vecs
36 self.vecs_mem = oarray_p(empty_p(self.nlocs), NULL, <void**>&self.vecs)
45 CHKERR(DMCompositeGetAccessArray(self.dm, self.gvec, self.nlocs, self.locs, self.vecs))
46 self.access = [ref_Vec(self.vecs[i]) for i from 0 <= i < n]
52 CHKERR(DMCompositeRestoreAccessArray(self.dm, self.gvec, self.nlocs, self.locs, self.vecs))
H A DVec.pyx954 vecs: Sequence[Vec],
963 vecs
976 vecs = list(vecs)
979 assert len(isets) == len(vecs)
983 cdef Py_ssize_t i, m = len(vecs)
989 for i from 0 <= i < m: cvecs[i] = (<Vec?>vecs[i]).vec
1902 def mDot(self, vecs: Sequence[Vec], out: ArrayScalar | None = None) -> ArrayScalar:
1909 vecs
1919 cdef PetscInt nv=<PetscInt>len(vecs), no=0
1925 v[i] = (<Vec?>(vecs[i])).vec
[all …]
/petsc/src/ksp/ksp/impls/gmres/fgmres/ !
H A Dfgmresimpl.h35 #define VEC_TEMP fgmres->vecs[0] /* work space */
36 #define VEC_TEMP_MATOP fgmres->vecs[1] /* work space */
38 fgmres->vecs[VEC_OFFSET + i] /* use to access
/petsc/src/vec/vec/interface/ftn-custom/ !
H A Dzvectorf90.c117 Vec *vecs; in vecdestroyvecs_() local
120 *ierr = F90Array1dAccess(ptr, MPIU_FORTRANADDR, (void **)&vecs PETSC_F90_2PTR_PARAM(ptrd)); in vecdestroyvecs_()
123 PETSC_FORTRAN_OBJECT_F_DESTROYED_TO_C_NULL(&vecs[i]); in vecdestroyvecs_()
124 *ierr = VecDestroy(&vecs[i]); in vecdestroyvecs_()
126 PETSC_FORTRAN_OBJECT_C_NULL_TO_F_DESTROYED(&vecs[i]); in vecdestroyvecs_()
130 *ierr = PetscFree(vecs); in vecdestroyvecs_()
/petsc/src/ksp/ksp/impls/gmres/lgmres/
H A Dlgmresimpl.h47 #define VEC_TEMP lgmres->vecs[0] /* work space */
48 #define VEC_TEMP_MATOP lgmres->vecs[1] /* work space */
49 #define VEC_VV(i) lgmres->vecs[VEC_OFFSET + i] /* use to access othog basis vectors */
/petsc/src/ksp/ksp/impls/gmres/
H A Dgmresimpl.h32 Vec *vecs; /* the work vectors */ \
86 #define VEC_TEMP gmres->vecs[0]
87 #define VEC_TEMP_MATOP gmres->vecs[1]
88 #define VEC_VV(i) gmres->vecs[VEC_OFFSET + i]
/petsc/src/ksp/ksp/impls/gmres/agmres/
H A Dagmresimpl.h42 #define VEC_TMP agmres->vecs[0]
43 #define VEC_TMP_MATOP agmres->vecs[1]
44 #define VEC_V(i) agmres->vecs[VEC_OFFSET + i]
/petsc/src/ksp/ksp/impls/gmres/dgmres/
H A Ddgmresimpl.h61 #define VEC_TEMP dgmres->vecs[0]
62 #define VEC_TEMP_MATOP dgmres->vecs[1]
63 #define VEC_VV(i) dgmres->vecs[VEC_OFFSET + i]
/petsc/src/dm/impls/composite/
H A Dpack.c247 …e DMCompositeGetAccessArray(DM dm, Vec pvec, PetscInt nwanted, const PetscInt wanted[], Vec vecs[]) in DMCompositeGetAccessArray() argument
279 vecs[wnum++] = v; in DMCompositeGetAccessArray()
309 …ompositeGetLocalAccessArray(DM dm, Vec pvec, PetscInt nwanted, const PetscInt wanted[], Vec vecs[]) in DMCompositeGetLocalAccessArray() argument
343 vecs[wnum++] = v; in DMCompositeGetLocalAccessArray()
417 …CompositeRestoreAccessArray(DM dm, Vec pvec, PetscInt nwanted, const PetscInt wanted[], Vec vecs[]) in DMCompositeRestoreAccessArray() argument
435 PetscCall(VecResetArray(vecs[wnum])); in DMCompositeRestoreAccessArray()
436 if (readonly) PetscCall(VecLockReadPop(vecs[wnum])); in DMCompositeRestoreAccessArray()
437 PetscCall(DMRestoreGlobalVector(link->dm, &vecs[wnum])); in DMCompositeRestoreAccessArray()
466 …ositeRestoreLocalAccessArray(DM dm, Vec pvec, PetscInt nwanted, const PetscInt wanted[], Vec *vecs) in DMCompositeRestoreLocalAccessArray() argument
484 PetscCall(VecResetArray(vecs[wnum])); in DMCompositeRestoreLocalAccessArray()
[all …]
/petsc/src/ts/impls/bdf/
H A Dbdf.c137 Vec *vecs = bdf->work; in TSBDF_VecLTE() local
146 PetscCall(VecMAXPY(lte, n + 1, alpha, vecs)); in TSBDF_VecLTE()
155 Vec *vecs = bdf->work + 1; in TSBDF_Extrapolate() local
162 PetscCall(VecMAXPY(X, n, alpha, vecs)); in TSBDF_Extrapolate()
171 Vec *vecs = bdf->work; in TSBDF_Interpolate() local
177 PetscCall(VecMAXPY(X, n, alpha, vecs)); in TSBDF_Interpolate()
190 Vec vecs[7]; in TSBDF_PreSolve() local
196 for (i = 1; i < n; i++) vecs[i] = bdf->transientvar ? bdf->tvwork[i] : bdf->work[i]; in TSBDF_PreSolve()
198 PetscCall(VecMAXPY(V0, n - 1, alpha + 1, vecs + 1)); in TSBDF_PreSolve()
/petsc/src/mat/ftn-mod/
H A Dpetscmat.h9011 subroutine MatNullSpaceRestoreVecs(sp, has_const, n, vecs, ierr)
17 Vec, pointer :: vecs(:)
/petsc/src/binding/petsc4py/src/petsc4py/lib/_pytypes/viewer/
H A Dpetscpyvista.py122 vecs = np.zeros((scalars[1].shape[0] // scalars[2], 3))
123 vecs[:, 0:2] = scalars[1].reshape(-1, scalars[2])
124 grid.point_data[scalars[0]] = vecs
/petsc/src/ksp/ksp/utils/lmvm/dense/
H A Ddenseqn.c150 Mat Sfull = lmvm->basis[LMBASIS_S]->vecs; in MatAllocate_LMVMDQN_Internal()
274 Mat Sfull = lmvm->basis[LMBASIS_S]->vecs; in MatUpdate_LMVMDQN()
275 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatUpdate_LMVMDQN()
671 Mat Sfull = lmvm->basis[LMBASIS_S]->vecs; in MatLMVMDBFGSUpdateMultData()
672 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatLMVMDBFGSUpdateMultData()
766 Mat Sfull = lmvm->basis[LMBASIS_S]->vecs; in MatSolve_LMVMDBFGS()
767 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatSolve_LMVMDBFGS()
855 Mat Sfull = lmvm->basis[LMBASIS_S]->vecs; in MatMult_LMVMDBFGS()
856 Mat Yfull = lmvm->basis[LMBASIS_Y]->vecs; in MatMult_LMVMDBFGS()
1096 Mat Sfull = lmvm->basis[LMBASIS_S]->vecs; in MatLMVMDDFPUpdateSolveData()
[all …]
/petsc/src/snes/tests/
H A Dex17.c465 Vec *vecs; in block_system() local
466 PetscCall(VecNestGetSubVecs(x, NULL, &vecs)); in block_system()
467 bv = vecs[0]; in block_system()
474 bv = vecs[1]; in block_system()
/petsc/src/ts/impls/implicit/alpha/
H A Dalpha1.c259 Vec vecs[3]; in TSEvaluateWLTE_Alpha() local
263 vecs[0] = th->X1; in TSEvaluateWLTE_Alpha()
264 vecs[1] = th->X0; in TSEvaluateWLTE_Alpha()
265 vecs[2] = th->vec_sol_prev; in TSEvaluateWLTE_Alpha()
267 PetscCall(VecMAXPY(Y, 3, scal, vecs)); in TSEvaluateWLTE_Alpha()

123