| /petsc/src/mat/tests/ |
| H A D | ex98.c | 16 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 D | ex80.c | 17 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 D | ex82.c | 17 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 D | ex83.c | 17 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 D | mmio.c | 23 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 D | mmloader.c | 8 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 D | ex213.c | 16 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 D | ex79f.F90 | 18 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 D | ex89.c | 15 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 D | ex232.c | 8 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 D | ex233.c | 8 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 D | ex111.c | 48 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 D | ex96.c | 52 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 D | spnd.c | 10 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 D | sprcm.c | 10 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 D | sp1wd.c | 10 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 D | spqmd.c | 10 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 D | wbm.c | 35 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 D | aijmkl.h | 17 …#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 D | ex83f.F90 | 19 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 D | metisnd.c | 11 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 D | ij.c | 33 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 D | amd.c | 22 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 D | essl.c | 16 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 D | fdsell.c | 10 …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()
|