| /petsc/src/binding/petsc4py/test/ |
| H A D | test_vec.py | 17 self.vec = v 20 self.vec.destroy() 21 self.vec = None 25 self.vec.set(1) 26 vec = self.vec.duplicate() 27 self.assertFalse(self.vec.equal(vec)) 28 self.assertEqual(self.vec.sizes, vec.sizes) 29 del vec 32 self.vec.set(1) 33 vec = self.vec.duplicate() [all …]
|
| H A D | test_dmshell.py | 18 vec = PETSc.Vec().create(comm=self.COMM) 19 vec.setSizes((10, None)) 20 vec.setUp() 21 self.dm.setGlobalVector(vec) 23 self.assertEqual(vec.getSizes(), gvec.getSizes()) 24 self.assertEqual(vec.comm, gvec.comm) 39 vec = PETSc.Vec().create(comm=PETSc.COMM_SELF) 40 vec.setSizes((1 + 10 * self.COMM.rank, None)) 41 vec.setUp() 42 self.dm.setLocalVector(vec) [all …]
|
| /petsc/src/ts/utils/dmplexlandau/tutorials/output/ |
| H A D | ex1_re.out | 6 [0] <vec:seq> adaptToleranceFEM(): 0) Refine phase: Z1 7 [0] <vec:seq> adaptToleranceFEM(): 0) (debug) found Z1 cells 8 [0] <vec:seq> adaptToleranceFEM(): 0) (debug) found Z1 cells 9 [0] <vec:seq> adaptToleranceFEM(): 0) (debug) found Z1 cells 10 [0] <vec:seq> adaptToleranceFEM(): 0) (debug) found Z1 cells 11 [0] <vec:seq> adaptToleranceFEM(): 0) (debug) found Z1 cells 12 [0] <vec:seq> adaptToleranceFEM(): 0) (debug) found Z1 cells 13 [0] <vec:seq> adaptToleranceFEM(): 0) (debug) found Z1 cells 14 [0] <vec:seq> adaptToleranceFEM(): 0) (debug) found Z1 cells 15 [0] <vec:seq> adaptToleranceFEM(): 0) (debug) found Z1 cells [all …]
|
| /petsc/src/dm/impls/da/kokkos/ |
| H A D | dagetov.kokkos.cxx | 8 #define DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, d… argument 14 PetscCall(VecGetLocalSize(vec, &N)); \ 27 PetscErrorCode DMDAVecGetKokkosOffsetView_Private(DM da, Vec vec, PetscScalarKokkosOffsetView1DType… in DMDAVecGetKokkosOffsetView_Private() argument 34 PetscValidHeaderSpecific(vec, VEC_CLASSID, 2); in DMDAVecGetKokkosOffsetView_Private() 36 DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof); in DMDAVecGetKokkosOffsetView_Private() 38 if (overwrite) PetscCall(VecGetKokkosViewWrite(vec, &kv)); in DMDAVecGetKokkosOffsetView_Private() 39 else PetscCall(VecGetKokkosView(vec, &kv)); in DMDAVecGetKokkosOffsetView_Private() 46 PetscErrorCode DMDAVecRestoreKokkosOffsetView_Private(DM da, Vec vec, PetscScalarKokkosOffsetView1D… in DMDAVecRestoreKokkosOffsetView_Private() argument 52 PetscValidHeaderSpecific(vec, VEC_CLASSID, 2); in DMDAVecRestoreKokkosOffsetView_Private() 55 if (overwrite) PetscCall(VecRestoreKokkosViewWrite(vec, &kv)); in DMDAVecRestoreKokkosOffsetView_Private() [all …]
|
| /petsc/src/vec/vec/interface/ |
| H A D | vecreg.c | 52 PetscErrorCode VecSetType(Vec vec, VecType newType) in VecSetType() argument 66 PetscValidHeaderSpecific(vec, VEC_CLASSID, 1); in VecSetType() 68 PetscCall(VecGetType(vec, &curType)); in VecSetType() 72 PetscCall(PetscObjectTypeCompare((PetscObject)vec, newType, &match)); in VecSetType() 76 PetscCall(PetscObjectGetComm((PetscObject)vec, &comm)); in VecSetType() 83 if (size == 1) PetscCall(PetscObjectTypeCompare((PetscObject)vec, VECSEQ, &srcStandard)); in VecSetType() 84 else PetscCall(PetscObjectTypeCompare((PetscObject)vec, VECMPI, &srcStandard)); in VecSetType() 93 PetscCall(PetscObjectTypeCompareAny((PetscObject)vec, &match, seqType, mpiType, "")); in VecSetType() 108 if (size == 1) PetscCall(VecConvert_Seq_SeqCUDA_inplace(vec)); in VecSetType() 109 else PetscCall(VecConvert_MPI_MPICUDA_inplace(vec)); in VecSetType() [all …]
|
| H A D | veccreate.c | 42 PetscErrorCode VecCreate(MPI_Comm comm, Vec *vec) in VecCreate() argument 47 PetscAssertPointer(vec, 2); in VecCreate() 53 *vec = v; in VecCreate() 80 …CreateFromOptions(MPI_Comm comm, const char *prefix, PetscInt bs, PetscInt m, PetscInt n, Vec *vec) in VecCreateFromOptions() argument 83 PetscAssertPointer(vec, 6); in VecCreateFromOptions() 84 PetscCall(VecCreate(comm, vec)); in VecCreateFromOptions() 85 if (prefix) PetscCall(VecSetOptionsPrefix(*vec, prefix)); in VecCreateFromOptions() 86 PetscCall(VecSetBlockSize(*vec, bs)); in VecCreateFromOptions() 87 PetscCall(VecSetSizes(*vec, m, n)); in VecCreateFromOptions() 88 PetscCall(VecSetFromOptions(*vec)); in VecCreateFromOptions() [all …]
|
| H A D | vector.c | 39 PetscErrorCode VecStashGetInfo(Vec vec, PetscInt *nstash, PetscInt *reallocs, PetscInt *bnstash, Pe… in VecStashGetInfo() argument 42 PetscCall(VecStashGetInfo_Private(&vec->stash, nstash, reallocs)); in VecStashGetInfo() 43 PetscCall(VecStashGetInfo_Private(&vec->bstash, bnstash, breallocs)); in VecStashGetInfo() 117 PetscErrorCode VecAssemblyBegin(Vec vec) in VecAssemblyBegin() argument 120 PetscValidHeaderSpecific(vec, VEC_CLASSID, 1); in VecAssemblyBegin() 121 PetscValidType(vec, 1); in VecAssemblyBegin() 122 PetscCall(VecStashViewFromOptions(vec, NULL, "-vec_view_stash")); in VecAssemblyBegin() 123 PetscCall(PetscLogEventBegin(VEC_AssemblyBegin, vec, 0, 0, 0)); in VecAssemblyBegin() 124 PetscTryTypeMethod(vec, assemblybegin); in VecAssemblyBegin() 125 PetscCall(PetscLogEventEnd(VEC_AssemblyBegin, vec, 0, 0, 0)); in VecAssemblyBegin() [all …]
|
| /petsc/src/dm/ftn-mod/ |
| H A D | petscdmda.h90 | 25 subroutine DMDAVecGetArray1(da, vec, array, ierr) 29 Vec vec 33 subroutine DMDAVecGetArray2(da, vec, array, ierr) 37 Vec vec 41 subroutine DMDAVecGetArray3(da, vec, array, ierr) 45 Vec vec 49 subroutine DMDAVecGetArray4(da, vec, array, ierr) 53 Vec vec 60 subroutine DMDAVecRestoreArray1(da, vec, array, ierr) 64 Vec vec [all …]
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | petscvec.pxi | 223 cdef inline Vec ref_Vec(PetscVec vec): argument 225 ob.vec = vec 234 cdef Vec vec = type(self)() 235 CHKERR(VecDuplicate(self.vec, &vec.vec)) 236 CHKERR(VecCopy(self.vec, vec.vec)) 237 return vec 240 cdef Vec vec = <Vec> vec_pos(self) 241 CHKERR(VecScale(vec.vec, -1)) 242 return vec 245 cdef Vec vec = <Vec> vec_pos(self) [all …]
|
| H A D | Vec.pyx | 47 self.obj = <PetscObject*> &self.vec 48 self.vec = NULL 161 CHKERR(VecView(self.vec, vwr)) 173 CHKERR(VecDestroy(&self.vec)) 196 CHKERR(PetscCLEAR(self.obj)); self.vec = newvec 216 CHKERR(VecSetType(self.vec, cval)) 240 CHKERR(VecSetSizes(self.vec, n, N)) 242 CHKERR(VecSetBlockSize(self.vec, bs)) 280 CHKERR(PetscCLEAR(self.obj)); self.vec = newvec 316 CHKERR(PetscCLEAR(self.obj)); self.vec = newvec [all …]
|
| H A D | Scatter.pyx | 142 vec_from.vec, cisfrom, vec_to.vec, cisto, &newsct)) 206 def toAll(cls, Vec vec) -> tuple[Scatter, Vec]: argument 213 vec 218 The created scatter will have the same communicator of ``vec``. 230 vec.vec, &scatter.sct, &ovec.vec)) 234 def toZero(cls, Vec vec) -> tuple[Scatter, Vec]: argument 241 vec 246 The created scatter will have the same communicator of ``vec``. 258 vec.vec, &scatter.sct, &ovec.vec)) 282 CHKERR(VecScatterBegin(self.sct, vec_from.vec, vec_to.vec, [all …]
|
| /petsc/src/vec/vec/utils/tagger/tutorials/ |
| H A D | ex1.c | 5 static PetscErrorCode ISGetBlockGlobalIS(IS is, Vec vec, PetscInt bs, IS *isBlockGlobal) in ISGetBlockGlobalIS() argument 12 PetscCall(VecGetLayout(vec, &map)); in ISGetBlockGlobalIS() 19 …PetscCall(ISCreateBlock(PetscObjectComm((PetscObject)vec), bs, n, idxout, PETSC_OWN_POINTER, isBlo… in ISGetBlockGlobalIS() 25 Vec vec, tagged, untagged; in main() local 51 PetscCall(VecCreate(comm, &vec)); in main() 52 PetscCall(PetscObjectSetName((PetscObject)vec, "Vec to Tag")); in main() 53 PetscCall(VecSetBlockSize(vec, bs)); in main() 54 PetscCall(VecSetSizes(vec, PETSC_DECIDE, n)); in main() 55 PetscCall(VecSetUp(vec)); in main() 56 PetscCall(VecGetLocalSize(vec, &nloc)); in main() [all …]
|
| /petsc/src/dm/impls/da/ |
| H A D | dagetarray.c | 43 PetscErrorCode DMDAVecGetArray(DM da, Vec vec, void *array) in DMDAVecGetArray() argument 49 PetscValidHeaderSpecific(vec, VEC_CLASSID, 2); in DMDAVecGetArray() 56 PetscCall(VecGetLocalSize(vec, &N)); in DMDAVecGetArray() 67 PetscCall(VecGetArray1d(vec, gxm * dof, gxs * dof, (PetscScalar **)array)); in DMDAVecGetArray() 69 PetscCall(VecGetArray2d(vec, gym, gxm * dof, gys, gxs * dof, (PetscScalar ***)array)); in DMDAVecGetArray() 71 … PetscCall(VecGetArray3d(vec, gzm, gym, gxm * dof, gzs, gys, gxs * dof, (PetscScalar ****)array)); in DMDAVecGetArray() 92 PetscErrorCode DMDAVecRestoreArray(DM da, Vec vec, void *array) in DMDAVecRestoreArray() argument 98 PetscValidHeaderSpecific(vec, VEC_CLASSID, 2); in DMDAVecRestoreArray() 105 PetscCall(VecGetLocalSize(vec, &N)); in DMDAVecRestoreArray() 116 PetscCall(VecRestoreArray1d(vec, gxm * dof, gxs * dof, (PetscScalar **)array)); in DMDAVecRestoreArray() [all …]
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex21.c | 202 Vec vec; in main() local 271 PetscCall(VecCreate(comm, &vec)); in main() 274 PetscCall(VecSetSizes(vec, n, PETSC_DECIDE)); in main() 275 PetscCall(VecSetUp(vec)); in main() 277 PetscCall(DMGetGlobalVector(sdm, &vec)); in main() 279 PetscCall(PetscObjectSetName((PetscObject)vec, exampleVecName)); in main() 280 PetscCall(VecGetArrayWrite(vec, &array)); in main() 311 PetscCall(VecRestoreArrayWrite(vec, &array)); in main() 312 PetscCall(DMPlexGlobalVectorView(dm, viewer, sdm, vec)); in main() 314 PetscCall(VecDestroy(&vec)); in main() [all …]
|
| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex6.c | 8 Vec vec; in main() local 17 PetscCall(DMCreateGlobalVector(dm, &vec)); in main() 18 PetscCall(VecSet(vec, 1.234)); in main() 24 PetscCall(DMStagVecSplitToDMDA(dm, vec, DMSTAG_ELEMENT, -3, &da, &vecda)); in main() 33 PetscCall(DMStagVecSplitToDMDA(dm, vec, DMSTAG_ELEMENT, -5, &da, &vecda)); in main() 42 PetscCall(DMStagVecSplitToDMDA(dm, vec, DMSTAG_ELEMENT, -2, &da, &vecda)); in main() 51 PetscCall(DMStagVecSplitToDMDA(dm, vec, DMSTAG_FRONT_DOWN_LEFT, 2, &da, &vecda)); in main() 60 PetscCall(DMStagVecSplitToDMDA(dm, vec, DMSTAG_BACK_RIGHT, 1, &da, &vecda)); in main() 69 PetscCall(DMStagVecSplitToDMDA(dm, vec, DMSTAG_DOWN, 0, &da, &vecda)); in main() 74 PetscCall(VecDestroy(&vec)); in main()
|
| /petsc/src/vec/vec/utils/matlab/ |
| H A D | gcreatev.c | 10 Vec vec = (Vec)obj; in VecMatlabEnginePut_Default() local 15 PetscCall(VecGetArrayRead(vec, &array)); in VecMatlabEnginePut_Default() 16 PetscCall(VecGetLocalSize(vec, &n)); in VecMatlabEnginePut_Default() 26 PetscCall(VecRestoreArrayRead(vec, &array)); in VecMatlabEnginePut_Default() 33 Vec vec = (Vec)obj; in VecMatlabEngineGet_Default() local 38 PetscCall(VecGetArray(vec, &array)); in VecMatlabEngineGet_Default() 39 PetscCall(VecGetLocalSize(vec, &n)); in VecMatlabEngineGet_Default() 43 PetscCall(VecRestoreArray(vec, &array)); in VecMatlabEngineGet_Default()
|
| /petsc/src/vec/vec/tutorials/ |
| H A D | ex7.c | 24 Vec vec; in main() local 33 PetscCall(VecCreate(PETSC_COMM_WORLD, &vec)); in main() 34 PetscCall(VecSetSizes(vec, PETSC_DECIDE, m)); in main() 35 PetscCall(VecSetFromOptions(vec)); in main() 44 ex7f_(&vec, &fcomm); in main() 46 PetscCall(VecView(vec, PETSC_VIEWER_STDOUT_WORLD)); in main() 47 PetscCall(VecDestroy(&vec)); in main()
|
| H A D | ex6.c | 19 Vec vec; in main() local 58 PetscCall(VecCreate(PETSC_COMM_SELF, &vec)); in main() 59 PetscCall(VecSetSizes(vec, PETSC_DECIDE, m)); in main() 60 PetscCall(VecSetFromOptions(vec)); in main() 61 PetscCall(VecGetArray(vec, &avec)); in main() 71 PetscCall(VecRestoreArray(vec, &avec)); in main() 72 PetscCall(VecView(vec, PETSC_VIEWER_STDOUT_SELF)); in main() 75 PetscCall(VecDestroy(&vec)); in main()
|
| H A D | ex6f.F90 | 16 Vec vec 62 PetscCallA(VecCreate(PETSC_COMM_SELF, vec, ierr)) 63 PetscCallA(VecSetSizes(vec, PETSC_DECIDE, m, ierr)) 65 PetscCallA(VecSetFromOptions(vec, ierr)) 67 PetscCallA(VecGetArray(vec, avec, ierr)) 81 PetscCallA(VecRestoreArray(vec, avec, ierr)) 82 PetscCallA(VecView(vec, PETSC_VIEWER_STDOUT_SELF, ierr)) 86 PetscCallA(VecDestroy(vec, ierr))
|
| H A D | ex7f.F90 | 4 subroutine ex7f(vec, comm) argument 13 Vec vec 22 PetscCall(VecSet(vec, two, ierr)) 29 PetscCall(ex7c(vec, comm, ierr)) 38 PetscCall(VecView(vec, PETSC_VIEWER_STDOUT_WORLD, ierr))
|
| /petsc/interfaces/lua/src/ |
| H A D | luasrc.c | 7 Vec vec; in lua_VecCreate() local 9 ierr = VecCreate(PETSC_COMM_SELF, &vec); in lua_VecCreate() 10 lua_pushlightuserdata(L, vec); in lua_VecCreate() 17 Vec vec; in lua_VecSetSize() local 21 vec = (Vec)lua_touserdata(L, 1); in lua_VecSetSize() 23 ierr = VecSetSizes(vec, n, n); in lua_VecSetSize()
|
| /petsc/src/vec/vec/utils/ |
| H A D | vecio.c | 13 PetscErrorCode VecView_Binary(Vec vec, PetscViewer viewer) in VecView_Binary() argument 21 PetscCheckSameComm(vec, 1, viewer, 2); in VecView_Binary() 25 PetscCall(VecGetLayout(vec, &map)); in VecView_Binary() 34 PetscCall(VecGetArrayRead(vec, &xarray)); in VecView_Binary() 36 PetscCall(VecRestoreArrayRead(vec, &xarray)); in VecView_Binary() 45 PetscCallMPI(MPI_Comm_rank(PetscObjectComm((PetscObject)vec), &rank)); in VecView_Binary() 49 PetscCall(PetscObjectGetName((PetscObject)vec, &name)); in VecView_Binary() 57 PetscCall(PetscObjectGetOptionsPrefix((PetscObject)vec, &pre)); in VecView_Binary() 58 …PETSC_COMM_SELF, info, "-%svecload_block_size %" PetscInt_FMT "\n", pre ? pre : "", vec->map->bs)); in VecView_Binary() 63 static PetscErrorCode VecLoad_Binary(Vec vec, PetscViewer viewer) in VecLoad_Binary() argument [all …]
|
| /petsc/src/dm/impls/moab/ |
| H A D | dmmbvec.cxx | 38 … moab::Tag tag, const moab::Range *range, PetscBool is_global_vec, PetscBool destroy_tag, Vec *vec) in DMMoabCreateVector() argument 43 PetscCall(DMCreateVector_Moab_Private(dm, tag, range, is_global_vec, destroy_tag, vec)); in DMMoabCreateVector() 60 PetscErrorCode DMMoabGetVecTag(Vec vec, moab::Tag *tag) in DMMoabGetVecTag() argument 69 PetscCall(PetscObjectQuery((PetscObject)vec, "MOABData", (PetscObject *)&moabdata)); in DMMoabGetVecTag() 89 PetscErrorCode DMMoabGetVecRange(Vec vec, moab::Range *range) in DMMoabGetVecRange() argument 98 PetscCall(PetscObjectQuery((PetscObject)vec, "MOABData", (PetscObject *)&moabdata)); in DMMoabGetVecRange() 120 PetscErrorCode DMMoabVecGetArray(DM dm, Vec vec, void *array) in DMMoabVecGetArray() argument 133 PetscValidHeaderSpecific(vec, VEC_CLASSID, 2); in DMMoabVecGetArray() 138 PetscCall(PetscObjectQuery((PetscObject)vec, "MOABData", (PetscObject *)&moabdata)); in DMMoabVecGetArray() 147 PetscCall(DMGlobalToLocalBegin(dm, vec, INSERT_VALUES, vmoab->local)); in DMMoabVecGetArray() [all …]
|
| /petsc/src/mat/interface/ |
| H A D | matnull.c | 92 Vec vec[6]; in MatNullSpaceCreateRigidBody() local 110 PetscCall(VecCreate(PetscObjectComm((PetscObject)coords), &vec[0])); in MatNullSpaceCreateRigidBody() 111 PetscCall(VecSetSizes(vec[0], dim * n, dim * N)); in MatNullSpaceCreateRigidBody() 112 PetscCall(VecSetBlockSize(vec[0], dim)); in MatNullSpaceCreateRigidBody() 113 PetscCall(VecSetUp(vec[0])); in MatNullSpaceCreateRigidBody() 114 for (i = 1; i < nmodes; i++) PetscCall(VecDuplicate(vec[0], &vec[i])); in MatNullSpaceCreateRigidBody() 115 for (i = 0; i < nmodes; i++) PetscCall(VecGetArray(vec[i], &v[i])); in MatNullSpaceCreateRigidBody() 148 for (i = 0; i < nmodes; i++) PetscCall(VecRestoreArray(vec[i], &v[i])); in MatNullSpaceCreateRigidBody() 152 PetscCall(VecMDot(vec[i], i, vec, dots)); in MatNullSpaceCreateRigidBody() 154 PetscCall(VecMAXPY(vec[i], i, dots, vec)); in MatNullSpaceCreateRigidBody() [all …]
|
| /petsc/lib/petsc/bin/ |
| H A D | PetscBinaryIO_tests.py | 20 vec = PETSc.Vec().createSeq(3) 21 vec[...] = array 23 vec.view(viewer) 25 vec.destroy() 35 vec = PETSc.Vec().createSeq(3) 36 vec.set(0.) 38 vec.load(viewer) 41 self.assertTrue(np.allclose(array, vec[...])) 42 vec.destroy()
|