Home
last modified time | relevance | path

Searched refs:array (Results 1 – 25 of 341) sorted by relevance

12345678910>>...14

/petsc/src/mat/ftn-mod/
H A Dpetscmatmod.F9060 subroutine MatDenseGetArrayF901d(v, array, ierr) argument
61 PetscScalar, pointer :: array(:)
64 call MatDenseGetArray(v, array, ierr)
70 subroutine MatDenseRestoreArrayF901d(v, array, ierr) argument
71 PetscScalar, pointer :: array(:)
74 call MatDenseRestoreArray(v, array, ierr)
80 subroutine MatDenseGetArrayReadF901d(v, array, ierr) argument
81 PetscScalar, pointer :: array(:)
84 call MatDenseGetArrayRead(v, array, ierr)
90 subroutine MatDenseRestoreArrayReadF901d(v, array, ierr) argument
[all …]
H A Dpetscmat.h9028 subroutine MatDenseGetArray1d(A, array, ierr)
31 PetscScalar, pointer :: array(:)
35 subroutine MatDenseGetArray2d(A, array, ierr)
38 PetscScalar, pointer :: array(:, :)
45 subroutine MatDenseRestoreArray1d(A, array, ierr)
48 PetscScalar, pointer :: array(:)
52 subroutine MatDenseRestoreArray2d(A, array, ierr)
55 PetscScalar, pointer :: array(:, :)
62 subroutine MatDenseGetArrayRead1d(A, array, ierr)
65 PetscScalar, pointer :: array(:)
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex81a.c12 Mat array[Q * Q], A, a; in main() local
25 for (i = 0; i < Q * Q; ++i) array[i] = NULL; in main()
28 …(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main()
30 …TSC_COMM_WORLD, 2, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main()
32 …TSC_COMM_WORLD, 2, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main()
34 PetscCall(MatAssemblyBegin(array[(Q + 1) * i], MAT_FINAL_ASSEMBLY)); in main()
35 PetscCall(MatAssemblyEnd(array[(Q + 1) * i], MAT_FINAL_ASSEMBLY)); in main()
36 PetscCall(MatShift(array[(Q + 1) * i], 100 + i + 1)); in main()
38 PetscCall(MatDuplicate(array[(Q + 1) * i], MAT_COPY_VALUES, &a)); in main()
39 PetscCall(MatDestroy(array + (Q + 1) * i)); in main()
[all …]
H A Dex81.c11 Mat array[Q * Q], A, a; in main() local
33 for (i = 0; i < Q * Q; ++i) array[i] = NULL; in main()
36 …(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main()
38 …TSC_COMM_WORLD, 2, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main()
40 …TSC_COMM_WORLD, 2, PETSC_DECIDE, PETSC_DECIDE, size, size, 1, NULL, 0, NULL, array + (Q + 1) * i)); in main()
42 PetscCall(MatAssemblyBegin(array[(Q + 1) * i], MAT_FINAL_ASSEMBLY)); in main()
43 PetscCall(MatAssemblyEnd(array[(Q + 1) * i], MAT_FINAL_ASSEMBLY)); in main()
44 PetscCall(MatShift(array[(Q + 1) * i], 100 + i + 1)); in main()
46 PetscCall(MatDuplicate(array[(Q + 1) * i], MAT_COPY_VALUES, &a)); in main()
47 PetscCall(MatDestroy(array + (Q + 1) * i)); in main()
[all …]
/petsc/include/
H A Dpetscbt.h26 static inline PetscErrorCode PetscBTMemzero(PetscCount m, PetscBT array) in PetscBTMemzero() argument
28 return PetscArrayzero(array, PetscBTLength(m)); in PetscBTMemzero()
31 static inline PetscErrorCode PetscBTDestroy(PetscBT *array) in PetscBTDestroy() argument
33 return (*array) ? PetscFree(*array) : PETSC_SUCCESS; in PetscBTDestroy()
36 static inline PetscErrorCode PetscBTCreate(PetscCount m, PetscBT *array) in PetscBTCreate() argument
38 return PetscCalloc1(PetscBTLength(m), array); in PetscBTCreate()
46 static inline PetscByte PetscBTLookup(PetscBT array, PetscCount index) in PetscBTLookup() argument
48 return array[PetscBTIndex_Internal(index)] & PetscBTMask_Internal(index); in PetscBTLookup()
51 static inline PetscErrorCode PetscBTSet(PetscBT array, PetscCount index) in PetscBTSet() argument
54 array[PetscBTIndex_Internal(index)] |= PetscBTMask_Internal(index); in PetscBTSet()
[all …]
/petsc/src/sys/ftn-custom/
H A Df90_cwrap.c121 PetscErrorCode F90Array1dCreate(void *array, MPI_Datatype type, PetscInt start, PetscInt len, F90Ar… in F90Array1dCreate() argument
125 if (!len) array = PETSC_NULL_SCALAR_Fortran; in F90Array1dCreate()
126 f90array1dcreatescalar_(array, &start, &len, ptr PETSC_F90_2PTR_PARAM(ptrd)); in F90Array1dCreate()
128 if (!len) array = PETSC_NULL_REAL_Fortran; in F90Array1dCreate()
129 f90array1dcreatereal_(array, &start, &len, ptr PETSC_F90_2PTR_PARAM(ptrd)); in F90Array1dCreate()
131 if (!len) array = PETSC_NULL_INTEGER_Fortran; in F90Array1dCreate()
132 f90array1dcreateint_(array, &start, &len, ptr PETSC_F90_2PTR_PARAM(ptrd)); in F90Array1dCreate()
137 if (!len) array = PETSC_NULL_INTEGER_Fortran; in F90Array1dCreate()
138 f90array1dcreatempiint_(array, &start, &len, ptr PETSC_F90_2PTR_PARAM(ptrd)); in F90Array1dCreate()
140 f90array1dcreatefortranaddr_(array, &start, &len, ptr PETSC_F90_2PTR_PARAM(ptrd)); in F90Array1dCreate()
[all …]
/petsc/src/dm/ftn-mod/
H A Dpetscdmda.h9025 subroutine DMDAVecGetArray1(da, vec, array, ierr)
30 PetscScalar, pointer :: array(:)
33 subroutine DMDAVecGetArray2(da, vec, array, ierr)
38 PetscScalar, pointer :: array(:, :)
41 subroutine DMDAVecGetArray3(da, vec, array, ierr)
46 PetscScalar, pointer :: array(:, :, :)
49 subroutine DMDAVecGetArray4(da, vec, array, ierr)
54 PetscScalar, pointer :: array(:, :, :, :)
60 subroutine DMDAVecRestoreArray1(da, vec, array, ierr)
65 PetscScalar, pointer :: array(:)
[all …]
/petsc/lib/petsc/bin/
H A DPetscBinaryIO_tests.py19 array = np.array([1.1, 2.2, 3.3])
21 vec[...] = array
28 self.assertTrue(np.allclose(array, result))
32 array = np.array([1.1, 2.2, 3.3])
33 PetscBinaryIO().writeBinaryFile('test.dat', [array.view(Vec),])
41 self.assertTrue(np.allclose(array, vec[...]))
46 indices = np.array([3,4,5])
64 vals = np.array([1.1,2.1,3.1])
65 counts = np.array([0,2,3])
66 cols = np.array([0,1,1])
[all …]
/petsc/src/dm/impls/plex/
H A Dplexpoint.c62 PetscErrorCode DMPlexPointLocalRead(DM dm, PetscInt point, const PetscScalar *array, void *ptr) in DMPlexPointLocalRead() argument
68 PetscAssertPointer(array, 3); in DMPlexPointLocalRead()
71 *(const PetscScalar **)ptr = (start < end) ? array + start : NULL; in DMPlexPointLocalRead()
100 PetscErrorCode DMPlexPointLocalRef(DM dm, PetscInt point, PetscScalar *array, void *ptr) in DMPlexPointLocalRef() argument
106 PetscAssertPointer(array, 3); in DMPlexPointLocalRef()
109 *(PetscScalar **)ptr = (start < end) ? array + start : NULL; in DMPlexPointLocalRef()
166 …PlexPointLocalFieldRead(DM dm, PetscInt point, PetscInt field, const PetscScalar *array, void *ptr) in DMPlexPointLocalFieldRead() argument
172 PetscAssertPointer(array, 4); in DMPlexPointLocalFieldRead()
175 *(const PetscScalar **)ptr = array + start; in DMPlexPointLocalFieldRead()
197 PetscErrorCode DMPlexPointLocalFieldRef(DM dm, PetscInt point, PetscInt field, PetscScalar *array, … in DMPlexPointLocalFieldRef() argument
[all …]
/petsc/src/sys/ftn-src/
H A Df90_fwrap.F909 PetscScalar, target :: array(start:start + len1 - 1)
12 ptr => array
19 PetscReal, target :: array(start:start + len1 - 1)
22 ptr => array
29 PetscInt, target :: array(start:start + len1 - 1)
32 ptr => array
39 PetscMPIInt, target :: array(start:start + len1 - 1)
42 ptr => array
49 PetscFortranAddr, target :: array(start:start + len1 - 1)
52 ptr => array
[all …]
/petsc/src/dm/impls/da/
H A Ddagetarray.c43 PetscErrorCode DMDAVecGetArray(DM da, Vec vec, void *array) in DMDAVecGetArray() argument
50 PetscAssertPointer(array, 3); 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
99 PetscAssertPointer(array, 3); in DMDAVecRestoreArray()
116 PetscCall(VecRestoreArray1d(vec, gxm * dof, gxs * dof, (PetscScalar **)array)); in DMDAVecRestoreArray()
118 PetscCall(VecRestoreArray2d(vec, gym, gxm * dof, gys, gxs * dof, (PetscScalar ***)array)); in DMDAVecRestoreArray()
120 …tscCall(VecRestoreArray3d(vec, gzm, gym, gxm * dof, gzs, gys, gxs * dof, (PetscScalar ****)array)); in DMDAVecRestoreArray()
[all …]
/petsc/src/dm/tutorials/
H A Dswarm_ex2.c32 PetscReal *array; in ex2_1() local
33 PetscCall(DMSwarmGetField(dms, "viscosity", &bs, NULL, (void **)&array)); in ex2_1()
34 for (p = 0; p < nlocal; p++) array[p] = 11.1 + p * 0.1 + rank * 100.0; in ex2_1()
35 PetscCall(DMSwarmRestoreField(dms, "viscosity", &bs, NULL, (void **)&array)); in ex2_1()
39 PetscReal *array; in ex2_1() local
40 PetscCall(DMSwarmGetField(dms, "strain", &bs, NULL, (void **)&array)); in ex2_1()
42 array[bs * p + 0] = 2.0e-2 + p * 0.001 + rank * 1.0; in ex2_1()
43 array[bs * p + 1] = 2.0e-2 + p * 0.002 + rank * 1.0; in ex2_1()
44 array[bs * p + 2] = 2.0e-2 + p * 0.003 + rank * 1.0; in ex2_1()
46 PetscCall(DMSwarmRestoreField(dms, "strain", &bs, NULL, (void **)&array)); in ex2_1()
H A Dex13f90aux.F9068 subroutine petsc_to_local(da, vec, array, f, dof, stw) argument
72 PetscReal, pointer :: array(:, :, :, :)
77 PetscCall(DMDAVecGetArray(da, vec, array, ierr))
78 call transform_petsc_us(array, f, stw)
80 subroutine transform_petsc_us(array, f, stw) argument
83 PetscReal, intent(in), dimension(:, 1 - stw:, 1 - stw:, 1 - stw:) :: array
85 f(:, :, :, :) = array(:, :, :, :)
87 subroutine local_to_petsc(da, vec, array, f, dof, stw) argument
91 PetscReal, pointer :: array(:, :, :, :)
95 call transform_us_petsc(array, f, stw)
[all …]
/petsc/src/vec/vec/utils/matlab/
H A Dgcreatev.c11 const PetscScalar *array; in VecMatlabEnginePut_Default() local
15 PetscCall(VecGetArrayRead(vec, &array)); in VecMatlabEnginePut_Default()
22 PetscCall(PetscArraycpy(mxGetPr(mat), array, n)); in VecMatlabEnginePut_Default()
26 PetscCall(VecRestoreArrayRead(vec, &array)); in VecMatlabEnginePut_Default()
34 PetscScalar *array; in VecMatlabEngineGet_Default() local
38 PetscCall(VecGetArray(vec, &array)); in VecMatlabEngineGet_Default()
42 PetscCall(PetscArraycpy(array, mxGetPr(mat), n)); in VecMatlabEngineGet_Default()
43 PetscCall(VecRestoreArray(vec, &array)); in VecMatlabEngineGet_Default()
/petsc/src/sys/tests/
H A Dex41.c11 PetscInt n, off, array[4], na, nb, i, *marray, size; in main() local
66 PetscCall(PetscHSetIGetElems(ht, &off, array)); in main()
67 PetscCall(PetscSortInt(off, array)); in main()
69 PetscTestCheck(array[0] == 13); in main()
70 PetscTestCheck(array[1] == 42); in main()
71 PetscCall(PetscHSetIGetElems(ht, &off, array)); in main()
72 PetscCall(PetscSortInt(2, array + 2)); in main()
74 PetscTestCheck(array[0] == 13); in main()
75 PetscTestCheck(array[1] == 42); in main()
76 PetscTestCheck(array[0] == 13); in main()
[all …]
H A Dex14.c8 PetscScalar a, array[10]; in main() local
25 PetscCall(PetscOptionsScalarArray("-array", "Input a scalar array", "ex14.c", array, &n, NULL)); in main()
29 if (PetscImaginaryPart(array[i]) < 0.0) { in main()
30 …COMM_SELF, " %g - %gi\n", (double)PetscRealPart(array[i]), (double)PetscAbsReal(PetscImaginaryPart… in main()
32 …(PETSC_COMM_SELF, " %g + %gi\n", (double)PetscRealPart(array[i]), (double)PetscImaginaryPart(array… in main()
/petsc/src/vec/vec/tutorials/
H A Dex9.c24 PetscScalar value, *array, *tarray = 0; in main() local
88 PetscCall(VecGetArray(gxs, &array)); in main()
89 value = array[0]; in main()
90 PetscCall(VecRestoreArray(gxs, &array)); in main()
124 PetscCall(VecGetArray(lx, &array)); in main()
125 …ynchronizedPrintf(PETSC_COMM_WORLD, "%" PetscInt_FMT " %g\n", i, (double)PetscRealPart(array[i]))); in main()
126 PetscCall(VecRestoreArray(lx, &array)); in main()
134 PetscCall(VecGetArray(lx, &array)); in main()
135 for (i = 0; i < nghost; i++) array[nlocal + i] = rank ? (PetscScalar)4 : (PetscScalar)8; in main()
136 PetscCall(VecRestoreArray(lx, &array)); in main()
[all …]
/petsc/src/vec/vec/tests/
H A Dex27.c7 PetscScalar *array; in TestSetInf() local
13 PetscCall(VecGetArrayWrite(v, &array)); in TestSetInf()
14 PetscCall(PetscArrayzero(array, n)); in TestSetInf()
15 PetscCall(VecRestoreArrayWrite(v, &array)); in TestSetInf()
20 PetscCall(VecGetArrayRead(v, (const PetscScalar **)&array)); in TestSetInf()
22 const PetscScalar x = array[i]; in TestSetInf()
26 PetscCall(VecRestoreArrayRead(v, (const PetscScalar **)&array)); in TestSetInf()
H A Dex53.c7 const PetscScalar *array; in CheckVecShift() local
12 PetscCall(VecGetArrayRead(v, &array)); in CheckVecShift()
14 const PetscScalar actual = array[i], expected = array_copy[i]; in CheckVecShift()
18 PetscCall(VecRestoreArrayRead(v, &array)); in CheckVecShift()
26 const PetscScalar *array; in main() local
44 PetscCall(VecGetArrayRead(x, &array)); in main()
45 PetscCall(PetscArraycpy(array_copy, array, n)); in main()
46 PetscCall(VecRestoreArrayRead(x, &array)); in main()
/petsc/src/dm/impls/da/tests/
H A Dex1.c9 PetscReal *array; in DMSwarmPrint() local
22 PetscCall(DMSwarmGetField(sw, coordFields[0], &bs, NULL, (void **)&array)); in DMSwarmPrint()
26 const PetscReal th = PetscAtan2Real(array[2 * p + 1], array[2 * p]) / PETSC_PI; in DMSwarmPrint()
27 …const PetscReal r = PetscSqrtReal(array[2 * p + 1] * array[2 * p + 1] + array[2 * p] * array[2 * … in DMSwarmPrint()
28 ….3f*pi cellid=%" PetscInt_FMT "\n", rank, pidArray[p], (double)array[2 * p], (double)array[2 * p +… in DMSwarmPrint()
30 PetscCall(DMSwarmRestoreField(sw, coordFields[0], &bs, NULL, (void **)&array)); in DMSwarmPrint()
/petsc/src/mat/utils/
H A Dfreespace.c11 a->array = a->array_head; in PetscFreeSpaceGet()
63 PetscInt row, nnz, *bj, *array, total, bi_temp; in PetscFreeSpaceContiguous_LU() local
73 array = (*head)->array_head; in PetscFreeSpaceContiguous_LU()
89 PetscCall(PetscArraycpy(bj, array, nnzL)); in PetscFreeSpaceContiguous_LU()
100 PetscCall(PetscArraycpy(bj, array + nnzL + 1, nnzU)); in PetscFreeSpaceContiguous_LU()
101 array += nnz; in PetscFreeSpaceContiguous_LU()
139 PetscInt row, nnz, *uj, *array, total; in PetscFreeSpaceContiguous_Cholesky() local
146 array = (*head)->array_head; in PetscFreeSpaceContiguous_Cholesky()
153 PetscCall(PetscArraycpy(uj, array + 1, nnz)); in PetscFreeSpaceContiguous_Cholesky()
154 uj[nnz] = array[0]; /* diagonal */ in PetscFreeSpaceContiguous_Cholesky()
[all …]
/petsc/src/mat/tests/
H A Dex16.c8 const PetscScalar *array; in CheckValues() local
12 PetscCall(MatDenseGetArrayRead(A, &array)); in CheckValues()
20 PetscReal w = PetscRealPart(array[ii + jj * lda]); in CheckValues()
24 PetscCall(MatDenseRestoreArrayRead(A, &array)); in CheckValues()
35 PetscScalar *array; in main() local
89 PetscCall(PetscMalloc1((rend - rstart) * N, &array)); in main()
90 for (i = 0; i < (rend - rstart) * N; i++) array[i] = (PetscReal)1; in main()
91 PetscCall(MatDensePlaceArray(A, array)); in main()
98 PetscCall(PetscFree(array)); in main()
124 PetscCall(PetscMalloc1(m * N + 1, &array)); in main()
[all …]
/petsc/src/vec/vec/utils/
H A Dvsection.c25 PetscScalar *array; in PetscSectionVecView() local
41 PetscCall(VecGetArray(v, &array)); in PetscSectionVecView()
42 PetscCall(PetscSectionArrayView_ASCII_Internal(s->field[f], array, PETSC_SCALAR, viewer)); in PetscSectionVecView()
43 PetscCall(VecRestoreArray(v, &array)); in PetscSectionVecView()
47 PetscCall(VecGetArray(v, &array)); in PetscSectionVecView()
48 PetscCall(PetscSectionArrayView_ASCII_Internal(s, array, PETSC_SCALAR, viewer)); in PetscSectionVecView()
49 PetscCall(VecRestoreArray(v, &array)); in PetscSectionVecView()
104 PetscScalar *baseArray, *array; in VecSetValuesSection() local
117 array = &baseArray[s->atlasOff[p]]; in VecSetValuesSection()
124 for (i = 0; i < dim; ++i) array[i] = values[i]; in VecSetValuesSection()
[all …]
/petsc/include/petsc/private/cpp/
H A Darray.hpp27 using array_return_type = std::array<typename return_type_helper<D, T...>::type, sizeof...(T)>;
30 inline constexpr std::array<T, NL + NR> concat_array_impl(const std::array<T, NL> &l, const std::ar… in concat_array_impl()
44 PETSC_NODISCARD inline constexpr auto concat_array(const std::array<T, NL> &l, const std::array<T, …
/petsc/src/vec/is/ao/interface/
H A Dao.c289 PetscErrorCode AOPetscToApplicationPermuteInt(AO ao, PetscInt block, PetscInt array[]) in AOPetscToApplicationPermuteInt() argument
293 PetscAssertPointer(array, 3); in AOPetscToApplicationPermuteInt()
294 PetscUseTypeMethod(ao, petsctoapplicationpermuteint, block, array); in AOPetscToApplicationPermuteInt()
324 PetscErrorCode AOApplicationToPetscPermuteInt(AO ao, PetscInt block, PetscInt array[]) in AOApplicationToPetscPermuteInt() argument
328 PetscAssertPointer(array, 3); in AOApplicationToPetscPermuteInt()
329 PetscUseTypeMethod(ao, applicationtopetscpermuteint, block, array); in AOApplicationToPetscPermuteInt()
359 PetscErrorCode AOPetscToApplicationPermuteReal(AO ao, PetscInt block, PetscReal array[]) in AOPetscToApplicationPermuteReal() argument
363 PetscAssertPointer(array, 3); in AOPetscToApplicationPermuteReal()
364 PetscUseTypeMethod(ao, petsctoapplicationpermutereal, block, array); in AOPetscToApplicationPermuteReal()
394 PetscErrorCode AOApplicationToPetscPermuteReal(AO ao, PetscInt block, PetscReal array[]) in AOApplicationToPetscPermuteReal() argument
[all …]

12345678910>>...14