Home
last modified time | relevance | path

Searched refs:ja (Results 1 – 25 of 50) sorted by relevance

12

/petsc/src/mat/tests/
H A Dex98.c16 PetscInt *ia, *ja; in main() local
26 PetscCall(PetscMalloc1(16, &ja)); in main()
28 ja[0] = 1; in main()
29 ja[1] = 4; in main()
30 ja[2] = 0; in main()
31 ja[3] = 2; in main()
32 ja[4] = 5; in main()
33 ja[5] = 1; in main()
34 ja[6] = 3; in main()
35 ja[7] = 6; in main()
[all …]
H A Dex80.c17 PetscInt *ia, *ja, row; in main() local
29 PetscCall(PetscMalloc1(16, &ja)); in main()
31 ja[0] = 1; in main()
32 ja[1] = 4; in main()
33 ja[2] = 0; in main()
34 ja[3] = 2; in main()
35 ja[4] = 5; in main()
36 ja[5] = 1; in main()
37 ja[6] = 3; in main()
38 ja[7] = 6; in main()
[all …]
H A Dex82.c17 PetscInt *ia, *ja; in main() local
31 PetscCall(PetscMalloc1(16, &ja)); in main()
33 ja[0] = 1; in main()
34 ja[1] = 4; in main()
35 ja[2] = 0; in main()
36 ja[3] = 2; in main()
37 ja[4] = 5; in main()
38 ja[5] = 1; in main()
39 ja[6] = 3; in main()
40 ja[7] = 6; in main()
[all …]
H A Dex83.c17 PetscInt *ia, *ja, *indices_sc, isrows_localsize; in main() local
32 PetscCall(PetscMalloc1(16, &ja)); in main()
34 ja[0] = 1; in main()
35 ja[1] = 4; in main()
36 ja[2] = 0; in main()
37 ja[3] = 2; in main()
38 ja[4] = 5; in main()
39 ja[5] = 1; in main()
40 ja[6] = 3; in main()
41 ja[7] = 6; in main()
[all …]
H A Dmmio.c23 int *ia, *ja; in mm_read_unsymmetric_sparse() local
53 ja = (int *)malloc(nz * sizeof(int)); in mm_read_unsymmetric_sparse()
58 *J_ = ja; 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()
70 ja[i]--; 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()
38 ja[i]--; /* adjust from 1-based to 0-based */ in MatCreateFromMTX()
41 if (ja[i] != ia[i]) rownz[ja[i]]++; in MatCreateFromMTX()
43 if (symmetric) rownz[ja[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()
83 PetscCall(PetscFree4(ia, ja, val, rownz)); in MatCreateFromMTX()
H A Dex213.c16 PetscInt *ia, *ja, bs = 2; in main() local
36 PetscCall(PetscMalloc1(3 * n, &ja)); in main()
47 ja[ia[i]] = col; in main()
54 ja[ia[i]] = col; in main()
61 ja[ia[i]] = col; in main()
70 PetscCall(MatMPIAIJSetPreallocationCSR(A, ia, ja, NULL)); in main()
77 PetscCall(MatMPIBAIJSetPreallocationCSR(A, bs, ia, ja, NULL)); in main()
82 PetscCall(PetscFree(ja)); in main()
H A Dex79f.F9018 PetscInt, pointer :: ia(:), ja(:), icol(:)
38 PetscCallA(MatGetRowIJ(Ad, one, bb, bb, n, ia, ja, done, ierr))
43 write (7 + rank, *) ' ', j, ja(j) + rstart, aa(j)
46 PetscCallA(MatRestoreRowIJ(Ad, one, bb, bb, n, ia, ja, done, ierr))
51 PetscCallA(MatGetRowIJ(Ao, one, bb, bb, n, ia, ja, done, ierr))
56 write (7 + rank, *) ' ', j, icol(ja(j)) + 1, aa(j)
60 PetscCallA(MatRestoreRowIJ(Ao, one, bb, bb, n, ia, ja, done, ierr))
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()
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()
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()
74 PetscCall(MatRestoreRowIJ(AB, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
H A Dex232.c8 PetscInt ja[4] = {0, 1, 0, 1}; in main() local
33 PetscCall(MatSeqSBAIJSetPreallocationCSR(ssbaij, 2, ia, ja, c)); in main()
34 PetscCall(MatMPISBAIJSetPreallocationCSR(msbaij, 2, ia, ja, c)); in main()
H A Dex233.c8 PetscInt ja[4] = {0, 1, 0, 1}; in main() local
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()
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()
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()
122 PetscCall(MatRestoreRowIJ(AB, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
H A Dex96.c52 const PetscInt *ia, *ja; in main() local
109 PetscCall(MatGetRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
115 PetscCall(MatRestoreRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
119 PetscCall(MatGetRowIJ(AA, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
125 PetscCall(MatRestoreRowIJ(AA, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
126 PetscCall(MatGetRowIJ(AB, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
132 PetscCall(MatRestoreRowIJ(AB, 0, PETSC_FALSE, PETSC_FALSE, &nrows, &ia, &ja, &flg)); in main()
/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 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()
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()
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()
/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/ksp/ksp/tutorials/
H A Dex83f.F9019 PetscInt, allocatable :: ja(:)
41 allocate (ia(n + 1), ja(nz), a(nz))
56 ja(1) = 0
60 ja(2 + 3*(i - 2)) = i - 2
63 ja(2 + 3*(i - 2) + 1) = i - 1
66 ja(2 + 3*(i - 2) + 2) = i
70 ja(nz) = 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/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()
52 if (ja[i] != j) adjncy[iptr++] = ja[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/mat/impls/aij/seq/
H A Dij.c33 PetscInt *work, *ia, *ja, *j, i, nz, row, col; in MatToSymmetricIJ_SeqAIJ() local
67 PetscCall(PetscMalloc1(nz, &ja)); in MatToSymmetricIJ_SeqAIJ()
68 *jja = ja; in MatToSymmetricIJ_SeqAIJ()
81 if (col != row) ja[work[col]++] = row + shiftout; in MatToSymmetricIJ_SeqAIJ()
82 ja[work[row]++] = col + shiftout; in MatToSymmetricIJ_SeqAIJ()
/petsc/src/mat/graphops/order/amd/
H A Damd.c22 const PetscInt *ia, *ja; in MatGetOrdering_AMD() local
33 PetscCall(MatGetRowIJ(mat, 0, PETSC_FALSE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_AMD()
52 status = amd_AMD_order(nrow, ia, ja, perm, Control, Info); in MatGetOrdering_AMD()
69 PetscCall(MatRestoreRowIJ(mat, 0, PETSC_FALSE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_AMD()
/petsc/src/mat/impls/aij/seq/essl/
H A Dessl.c16 int *ja; member
32 if (essl->CleanUpESSL) PetscCall(PetscFree4(essl->a, essl->aux, essl->ia, essl->ja)); in MatDestroy_Essl()
47 dgss(&zero, &nessl, essl->a, essl->ia, essl->ja, &essl->lna, xx, essl->aux, &essl->naux); in MatSolve_Essl()
62 for (i = 0; i < aa->nz; i++) essl->ja[i] = aa->j[i] + 1; in MatLUFactorNumeric_Essl()
76 …dgsf(&one, &nessl, &essl->nz, essl->a, essl->ia, essl->ja, &essl->lna, essl->iparm, essl->rparm, e… in MatLUFactorNumeric_Essl()
100 …cMalloc4(essl->lna, &essl->a, essl->naux, &essl->aux, essl->lna, &essl->ia, essl->lna, &essl->ja)); in MatLUFactorSymbolic_Essl()
/petsc/src/mat/impls/sell/seq/
H A Dfdsell.c10 …cBool inodecompressed, PetscInt *nn, const PetscInt *ia[], const PetscInt *ja[], PetscInt *spidx[]… in MatGetColumnIJ_SeqSELL_Color() argument
54 *ja = cja; in MatGetColumnIJ_SeqSELL_Color()
59 …scBool inodecompressed, PetscInt *n, const PetscInt *ia[], const PetscInt *ja[], PetscInt *spidx[]… in MatRestoreColumnIJ_SeqSELL_Color() argument
64 PetscCall(PetscFree(*ja)); in MatRestoreColumnIJ_SeqSELL_Color()

12