Home
last modified time | relevance | path

Searched refs:ia (Results 1 – 25 of 59) sorted by relevance

123

/petsc/src/mat/tests/
H A Dex98.c16 PetscInt *ia, *ja; in main() local
25 PetscCall(PetscMalloc1(5, &ia)); in main()
38 ia[0] = 0; in main()
39 ia[1] = 2; in main()
40 ia[2] = 5; in main()
41 ia[3] = 8; in main()
42 ia[4] = 10; in main()
58 ia[0] = 0; in main()
59 ia[1] = 3; in main()
60 ia[2] = 7; in main()
[all …]
H A Dex80.c17 PetscInt *ia, *ja, row; in main() local
28 PetscCall(PetscMalloc1(5, &ia)); in main()
41 ia[0] = 0; in main()
42 ia[1] = 2; in main()
43 ia[2] = 5; in main()
44 ia[3] = 8; in main()
45 ia[4] = 10; in main()
61 ia[0] = 0; in main()
62 ia[1] = 3; in main()
63 ia[2] = 7; in main()
[all …]
H A Dex82.c17 PetscInt *ia, *ja; in main() local
30 PetscCall(PetscMalloc1(5, &ia)); in main()
43 ia[0] = 0; in main()
44 ia[1] = 2; in main()
45 ia[2] = 5; in main()
46 ia[3] = 8; in main()
47 ia[4] = 10; in main()
63 ia[0] = 0; in main()
64 ia[1] = 3; in main()
65 ia[2] = 7; in main()
[all …]
H A Dex83.c17 PetscInt *ia, *ja, *indices_sc, isrows_localsize; in main() local
31 PetscCall(PetscMalloc1(5, &ia)); in main()
44 ia[0] = 0; in main()
45 ia[1] = 2; in main()
46 ia[2] = 5; in main()
47 ia[3] = 8; in main()
48 ia[4] = 10; in main()
65 ia[0] = 0; in main()
66 ia[1] = 3; in main()
67 ia[2] = 7; in main()
[all …]
H A Dex213.c16 PetscInt *ia, *ja, bs = 2; in main() local
35 PetscCall(PetscMalloc1(n + 1, &ia)); in main()
40 ia[0] = 0; in main()
42 ia[i] = ia[i - 1]; in main()
47 ja[ia[i]] = col; in main()
48 ia[i]++; in main()
54 ja[ia[i]] = col; in main()
55 ia[i]++; in main()
61 ja[ia[i]] = col; in main()
62 ia[i]++; in main()
[all …]
H A Dex79f.F9018 PetscInt, pointer :: ia(:), ja(:), icol(:)
38 PetscCallA(MatGetRowIJ(Ad, one, bb, bb, n, ia, ja, done, ierr))
41 write (7 + rank, *) 'row ', i + rstart, ' number nonzeros ', ia(i + 1) - ia(i)
42 do j = ia(i), ia(i + 1) - 1
46 PetscCallA(MatRestoreRowIJ(Ad, one, bb, bb, n, ia, ja, done, ierr))
51 PetscCallA(MatGetRowIJ(Ao, one, bb, bb, n, ia, ja, done, ierr))
54 write (7 + rank, *) 'row ', i + rstart, ' number nonzeros ', ia(i + 1) - ia(i)
55 do j = ia(i), ia(i + 1) - 1
60 PetscCallA(MatRestoreRowIJ(Ao, one, bb, bb, n, ia, ja, done, ierr))
H A Dmmio.c23 int *ia, *ja; in mm_read_unsymmetric_sparse() local
52 ia = (int *)malloc(nz * sizeof(int)); in mm_read_unsymmetric_sparse()
57 *I_ = ia; in mm_read_unsymmetric_sparse()
65 if (fscanf(f, "%d %d %lg\n", &ia[i], &ja[i], &val[i]) != 3) { in mm_read_unsymmetric_sparse()
69 ia[i]--; /* adjust from 1-based to 0-based */ in mm_read_unsymmetric_sparse()
205 int mm_read_mtx_crd_data(FILE *f, int M, int N, int nz, int ia[], int ja[], double val[], MM_typeco… in mm_read_mtx_crd_data() argument
210 …if (fscanf(f, "%d %d %lg %lg", &ia[i], &ja[i], &val[2 * i], &val[2 * i + 1]) != 4) return MM_PREMA… in mm_read_mtx_crd_data()
213 if (fscanf(f, "%d %d %lg\n", &ia[i], &ja[i], &val[i]) != 3) return MM_PREMATURE_EOF; in mm_read_mtx_crd_data()
219 if (fscanf(f, "%d %d", &ia[i], &ja[i]) != 2) return MM_PREMATURE_EOF; in mm_read_mtx_crd_data()
225 int mm_read_mtx_crd_entry(FILE *f, int *ia, int *ja, double *real, double *imag, MM_typecode matcod… in mm_read_mtx_crd_entry() argument
[all …]
H A Dmmloader.c8 PetscInt *ia, *ja; in MatCreateFromMTX() local
28 PetscCall(PetscMalloc4(nz, &ia, nz, &ja, nz, &val, M, &rownz)); in MatCreateFromMTX()
35 ninput = fscanf(file, "%d %d %lg\n", &ia[i], &ja[i], &val[i]); in MatCreateFromMTX()
37 ia[i]--; in MatCreateFromMTX()
40 rownz[ia[i]]++; in MatCreateFromMTX()
41 if (ja[i] != ia[i]) rownz[ja[i]]++; in MatCreateFromMTX()
44 else rownz[ia[i]]++; in MatCreateFromMTX()
68 …for (j = 0; j < nz; j++) PetscCall(MatSetValues(*A, 1, &ia[j], 1, &ja[j], &val[j], INSERT_VALUES)); in MatCreateFromMTX()
73 …for (j = 0; j < nz; j++) PetscCall(MatSetValues(*A, 1, &ja[j], 1, &ia[j], &val[j], INSERT_VALUES)); in MatCreateFromMTX()
78 PetscCall(MatSetValues(*A, 1, &ja[j], 1, &ia[j], &val[j], INSERT_VALUES)); in MatCreateFromMTX()
[all …]
H A Dex89.c15 const PetscInt *ia, *ja; in main() local
51 PetscCall(MatGetRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
54 for (i = 0; i < ia[nrows]; i++) array[i] = one; in main()
57 PetscCall(MatRestoreRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
61 PetscCall(MatGetRowIJ(AA, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
64 for (i = 0; i < ia[nrows]; i++) array[i] = one; in main()
67 PetscCall(MatRestoreRowIJ(AA, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
68 PetscCall(MatGetRowIJ(AB, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
71 for (i = 0; i < ia[nrows]; i++) array[i] = one; in main()
74 PetscCall(MatRestoreRowIJ(AB, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
H A Dex233.c7 PetscInt ia[3] = {0, 2, 4}; in main() local
24 for (i = 0; i < 3; i++) ia[i] = 0; in main()
39 PetscCall(MatMPIBAIJSetPreallocationCSR(ssbaij, 2, ia, ja, c)); in main()
40 PetscCall(MatMPISBAIJSetPreallocationCSR(ssbaij, 2, ia, ja, c)); in main()
H A Dex111.c48 const PetscInt *ia, *ja; in main() local
99 PetscCall(MatGetRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
102 for (i = 0; i < ia[nrows]; i++) array[i] = one; in main()
105 PetscCall(MatRestoreRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
109 PetscCall(MatGetRowIJ(AA, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
112 for (i = 0; i < ia[nrows]; i++) array[i] = one; in main()
115 PetscCall(MatRestoreRowIJ(AA, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
116 PetscCall(MatGetRowIJ(AB, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
119 for (i = 0; i < ia[nrows]; i++) array[i] = one; in main()
122 PetscCall(MatRestoreRowIJ(AB, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
/petsc/src/vec/is/ao/tests/
H A Dex5.c19 PetscInt *app_indices, *petsc_indices, *ia, *ia0; in main() local
57 PetscCall(PetscMalloc1(n_loc, &ia)); in main()
60 ia[0] = 0; in main()
61 ia[1] = -1; in main()
62 ia[2] = 1; in main()
63 ia[3] = 2; in main()
64 ia[4] = -1; in main()
65 ia[5] = 4; in main()
66 ia[6] = 5; in main()
67 ia[7] = 6; in main()
[all …]
/petsc/src/mat/impls/aij/seq/
H A Dij.c33 PetscInt *work, *ia, *ja, *j, i, nz, row, col; in MatToSymmetricIJ_SeqAIJ() local
37 PetscCall(PetscCalloc1(m + 1, &ia)); in MatToSymmetricIJ_SeqAIJ()
38 *iia = ia; in MatToSymmetricIJ_SeqAIJ()
42 ia[0] = shiftout; in MatToSymmetricIJ_SeqAIJ()
53 if (col != row) ia[row + 1]++; in MatToSymmetricIJ_SeqAIJ()
54 ia[col + 1]++; in MatToSymmetricIJ_SeqAIJ()
60 row = ia[i - 1]; in MatToSymmetricIJ_SeqAIJ()
61 ia[i] += row; in MatToSymmetricIJ_SeqAIJ()
66 nz = ia[m] + (!shiftin); in MatToSymmetricIJ_SeqAIJ()
/petsc/src/vec/is/ao/interface/
H A Dao.c134 PetscInt *ia; in AOPetscToApplicationIS() local
141 PetscCall(ISGetIndices(is, (const PetscInt **)&ia)); in AOPetscToApplicationIS()
143 PetscUseTypeMethod(ao, petsctoapplication, n, ia); in AOPetscToApplicationIS()
144 PetscCall(ISRestoreIndices(is, (const PetscInt **)&ia)); in AOPetscToApplicationIS()
177 PetscInt n, *ia; in AOApplicationToPetscIS() local
184 PetscCall(ISGetIndices(is, (const PetscInt **)&ia)); in AOApplicationToPetscIS()
186 PetscUseTypeMethod(ao, applicationtopetsc, n, ia); in AOApplicationToPetscIS()
187 PetscCall(ISRestoreIndices(is, (const PetscInt **)&ia)); in AOApplicationToPetscIS()
219 PetscErrorCode AOPetscToApplication(AO ao, PetscInt n, PetscInt ia[]) in AOPetscToApplication() argument
223 if (n) PetscAssertPointer(ia, 3); in AOPetscToApplication()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex83f.F9018 PetscInt, allocatable :: ia(:)
41 allocate (ia(n + 1), ja(nz), a(nz))
49 ia(1) = 0
50 ia(2) = 1
52 ia(i) = ia(i - 1) + 3
54 ia(n + 1) = ia(n) + 1
74 PetscCallA(MatCreateSeqAIJWithArrays(PETSC_COMM_SELF, n, n, ia, ja, a, J, ierr))
118 deallocate (ia, ja, a)
/petsc/src/mat/graphops/order/
H A Dspnd.c10 const PetscInt *ia, *ja; in MatGetOrdering_ND() local
15 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_ND()
18 PetscCall(MatGetRowIJ(B, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_ND()
22 PetscCall(SPARSEPACKgennd(&nrow, ia, ja, mask, perm, xls, ls)); in MatGetOrdering_ND()
24 PetscCall(MatRestoreRowIJ(B, 1, PETSC_TRUE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_ND()
27 PetscCall(MatRestoreRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_ND()
H A Dsprcm.c10 const PetscInt *ia, *ja; in MatGetOrdering_RCM() local
14 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_RCM()
18 PetscCall(SPARSEPACKgenrcm(&nrow, ia, ja, perm, mask, xls)); in MatGetOrdering_RCM()
19 PetscCall(MatRestoreRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_RCM()
H A Dsp1wd.c10 const PetscInt *ia, *ja; in MatGetOrdering_1WD() local
14 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_1WD()
18 PetscCall(SPARSEPACKgen1wd(&nrow, ia, ja, mask, &nblks, xblk, perm, xls, ls)); in MatGetOrdering_1WD()
19 PetscCall(MatRestoreRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_1WD()
H A Dwbm.c35 const PetscInt *ia, *ja; in MatGetOrdering_WBM() local
44 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_WBM()
46 nnz = ia[nrow]; in MatGetOrdering_WBM()
79 …PetscCallExternal(mc64ad_dist, &job, &nrow, &nnz, ia, ja, a, &num, perm, &liw, iw, &ldw, dw, icntl… in MatGetOrdering_WBM()
80 PetscCall(MatRestoreRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_WBM()
H A Dspqmd.c10 const PetscInt *ia, *ja; in MatGetOrdering_QMD() local
14 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_QMD()
21 …PetscCall(SPARSEPACKgenqmd(&nrow, ia, ja, perm, iperm, deg, marker, rchset, nbrhd, qsize, qlink, &… in MatGetOrdering_QMD()
22 PetscCall(MatRestoreRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_QMD()
/petsc/src/mat/graphops/order/metisnd/
H A Dmetisnd.c11 const PetscInt *ia, *ja; in MatGetOrdering_METISND() local
18 PetscCall(MatGetRowIJ(mat, 0, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_METISND()
21 PetscCall(MatGetRowIJ(B, 0, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_METISND()
45 PetscCall(PetscMalloc4(nrow + 1, &xadj, ia[nrow], &adjncy, nrow, &perm, nrow, &iperm)); in MatGetOrdering_METISND()
51 for (i = ia[j]; i < ia[j + 1]; i++) { in MatGetOrdering_METISND()
72 PetscCall(MatRestoreRowIJ(B, 0, PETSC_TRUE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_METISND()
75 PetscCall(MatRestoreRowIJ(mat, 0, PETSC_TRUE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_METISND()
/petsc/src/vec/is/ao/impls/basic/
H A Daobasic.c47 static PetscErrorCode AOPetscToApplication_Basic(AO ao, PetscInt n, PetscInt *ia) in AOPetscToApplication_Basic() argument
54 if (ia[i] >= 0 && ia[i] < N) { in AOPetscToApplication_Basic()
55 ia[i] = aobasic->app[ia[i]]; in AOPetscToApplication_Basic()
57 ia[i] = -1; in AOPetscToApplication_Basic()
63 static PetscErrorCode AOApplicationToPetsc_Basic(AO ao, PetscInt n, PetscInt *ia) in AOApplicationToPetsc_Basic() argument
70 if (ia[i] >= 0 && ia[i] < N) { in AOApplicationToPetsc_Basic()
71 ia[i] = aobasic->petsc[ia[i]]; in AOApplicationToPetsc_Basic()
73 ia[i] = -1; in AOApplicationToPetsc_Basic()
158 PetscInt napp, *allpetsc, *allapp, ip, ia, N, i, *petsc = NULL, start; in AOCreate_Basic() local
226 ia = allapp[i]; in AOCreate_Basic()
[all …]
/petsc/src/mat/impls/aij/seq/aijmkl/
H A Daijmkl.h17 …#define mkl_cspblas_xcsrgemv(transa, m, a, ia, ja, x, y) mkl_cspblas_scsrgemv(transa, m, a, ia, ja… argument
19 …#define mkl_cspblas_xcsrgemv(transa, m, a, ia, ja, x, y) mkl_cspblas_dcsrgemv(transa, m, a, ia, ja… argument
23 …#define mkl_cspblas_xcsrgemv(transa, m, a, ia, ja, x, y) mkl_cspblas_ccsrgemv(transa, m, a, ia, ja… argument
25 …#define mkl_cspblas_xcsrgemv(transa, m, a, ia, ja, x, y) mkl_cspblas_zcsrgemv(transa, m, a, ia, ja… argument
/petsc/src/mat/impls/sell/seq/
H A Dfdsell.c10 …scBool symmetric, PetscBool inodecompressed, PetscInt *nn, const PetscInt *ia[], const PetscInt *j… in MatGetColumnIJ_SeqSELL_Color() argument
21 if (!ia) PetscFunctionReturn(PETSC_SUCCESS); in MatGetColumnIJ_SeqSELL_Color()
53 *ia = cia; in MatGetColumnIJ_SeqSELL_Color()
59 …tscBool symmetric, PetscBool inodecompressed, PetscInt *n, const PetscInt *ia[], const PetscInt *j… in MatRestoreColumnIJ_SeqSELL_Color() argument
62 if (!ia) PetscFunctionReturn(PETSC_SUCCESS); in MatRestoreColumnIJ_SeqSELL_Color()
63 PetscCall(PetscFree(*ia)); in MatRestoreColumnIJ_SeqSELL_Color()
/petsc/src/vec/is/ao/impls/memscalable/
H A Daomemscalable.c90 static PetscErrorCode AOMap_MemoryScalable_private(AO ao, PetscInt n, PetscInt *ia, const PetscInt … in AOMap_MemoryScalable_private() argument
118 …if (ia[i] < 0) owner[i] = -1; /* mark negative entries (which are not to be mapped) with a sp… in AOMap_MemoryScalable_private()
119 if (ia[i] >= ao->N) owner[i] = -2; /* mark out of range entries with special negative value */ in AOMap_MemoryScalable_private()
122 if (lastidx > (idx = ia[i])) j = 0; in AOMap_MemoryScalable_private()
168 else if (j == -2) { /* out of range entries get mapped to -1 */ ia[i] = -1; in AOMap_MemoryScalable_private()
171 sindices[start[j]++] = ia[i]; in AOMap_MemoryScalable_private()
173 ia[i] = maploc[ia[i] - owners[rank]]; in AOMap_MemoryScalable_private()
221 if (source == owner[i]) ia[i] = rbuf[count++]; in AOMap_MemoryScalable_private()
235 static PetscErrorCode AOPetscToApplication_MemoryScalable(AO ao, PetscInt n, PetscInt *ia) in AOPetscToApplication_MemoryScalable() argument
241 PetscCall(AOMap_MemoryScalable_private(ao, n, ia, app_loc)); in AOPetscToApplication_MemoryScalable()
[all …]

123