Home
last modified time | relevance | path

Searched refs:aij (Results 1 – 25 of 53) sorted by relevance

123

/petsc/src/mat/impls/aij/mpi/
H A Dmmaij.c10 Mat_MPIAIJ *aij = (Mat_MPIAIJ *)mat->data; in MatSetUpMultiply_MPIAIJ() local
11 Mat_SeqAIJ *B = (Mat_SeqAIJ *)aij->B->data; in MatSetUpMultiply_MPIAIJ()
25 if (!aij->garray) { in MatSetUpMultiply_MPIAIJ()
26 PetscCheck(aij->B, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Missing B mat"); in MatSetUpMultiply_MPIAIJ()
29 PetscCall(PetscHMapICreateWithSize(aij->B->rmap->n, &gid1_lid1)); in MatSetUpMultiply_MPIAIJ()
30 for (i = 0; i < aij->B->rmap->n; i++) { in MatSetUpMultiply_MPIAIJ()
55 for (i = 0; i < aij->B->rmap->n; i++) { in MatSetUpMultiply_MPIAIJ()
63 PetscCall(PetscLayoutDestroy(&aij->B->cmap)); in MatSetUpMultiply_MPIAIJ()
64 …PetscCall(PetscLayoutCreateFromSizes(PetscObjectComm((PetscObject)aij->B), ec, ec, 1, &aij->B->cma… in MatSetUpMultiply_MPIAIJ()
70 for (i = 0; i < aij->B->rmap->n; i++) { in MatSetUpMultiply_MPIAIJ()
[all …]
H A Dmpb_aij.c5 Mat_MPIAIJ *aij = (Mat_MPIAIJ *)mat->data; in MatGetMultiProcBlock_MPIAIJ() local
6 Mat_SeqAIJ *aijB = (Mat_SeqAIJ *)aij->B->data; in MatGetMultiProcBlock_MPIAIJ()
35 PetscCall(PetscCalloc1(aij->B->cmap->n, &garrayCMap)); in MatGetMultiProcBlock_MPIAIJ()
36 for (i = 0; i < aij->B->cmap->n; i++) { in MatGetMultiProcBlock_MPIAIJ()
37 col = aij->garray[i]; in MatGetMultiProcBlock_MPIAIJ()
49 PetscCall(PetscCalloc1(aij->B->rmap->n, &nnz)); in MatGetMultiProcBlock_MPIAIJ()
50 for (i = 0; i < aij->B->rmap->n; i++) { in MatGetMultiProcBlock_MPIAIJ()
59 ((Mat_MPIAIJ *)((*subMat)->data))->A = aij->A; in MatGetMultiProcBlock_MPIAIJ()
60 PetscCall(PetscObjectReference((PetscObject)aij->A)); in MatGetMultiProcBlock_MPIAIJ()
61 } else if (((Mat_MPIAIJ *)(*subMat)->data)->A != aij->A) { in MatGetMultiProcBlock_MPIAIJ()
[all …]
H A Dmpiaij.c18 Mat_MPIAIJ *aij = (Mat_MPIAIJ *)mat->data; in MatReset_MPIAIJ() local
23 PetscCall(VecDestroy(&aij->diag)); in MatReset_MPIAIJ()
24 PetscCall(MatDestroy(&aij->A)); in MatReset_MPIAIJ()
25 PetscCall(MatDestroy(&aij->B)); in MatReset_MPIAIJ()
27 PetscCall(PetscHMapIDestroy(&aij->colmap)); in MatReset_MPIAIJ()
29 PetscCall(PetscFree(aij->colmap)); in MatReset_MPIAIJ()
31 PetscCall(PetscFree(aij->garray)); in MatReset_MPIAIJ()
32 PetscCall(VecDestroy(&aij->lvec)); in MatReset_MPIAIJ()
33 PetscCall(VecScatterDestroy(&aij->Mvctx)); in MatReset_MPIAIJ()
34 PetscCall(PetscFree2(aij->rowvalues, aij->rowindices)); in MatReset_MPIAIJ()
[all …]
H A Dfdmpiaij.c483 Mat_MPIAIJ *aij = (Mat_MPIAIJ *)mat->data; in MatFDColoringSetUp_MPIXAIJ() local
485 A = aij->A; in MatFDColoringSetUp_MPIXAIJ()
488 B = aij->B; in MatFDColoringSetUp_MPIXAIJ()
492 if (!aij->colmap) { in MatFDColoringSetUp_MPIXAIJ()
497 colmap = aij->colmap; in MatFDColoringSetUp_MPIXAIJ()
504 …scObjectComm((PetscObject)mat), mat->cmap->n, PETSC_DETERMINE, B->cmap->n, aij->garray, &c->vscale… in MatFDColoringSetUp_MPIXAIJ()
696 Mat_MPIAIJ *aij = (Mat_MPIAIJ *)mat->data; in MatFDColoringCreate_MPIXAIJ() local
704 A = aij->A; in MatFDColoringCreate_MPIXAIJ()
706 B = aij->B; in MatFDColoringCreate_MPIXAIJ()
H A Dmpiov.c2848 Mat_MPIAIJ *aij; in MatSetSeqMats_MPIAIJ() local
2890 aij = (Mat_MPIAIJ *)C->data; in MatSetSeqMats_MPIAIJ()
2891 if (!aij->A) { in MatSetSeqMats_MPIAIJ()
2893 PetscCall(MatCreate(PETSC_COMM_SELF, &aij->A)); in MatSetSeqMats_MPIAIJ()
2894 PetscCall(MatSetSizes(aij->A, C->rmap->n, C->cmap->n, C->rmap->n, C->cmap->n)); in MatSetSeqMats_MPIAIJ()
2895 PetscCall(MatSetBlockSizesFromMats(aij->A, C, C)); in MatSetSeqMats_MPIAIJ()
2896 PetscCall(MatSetType(aij->A, MATSEQAIJ)); in MatSetSeqMats_MPIAIJ()
2899 PetscCall(MatSetSeqMat_SeqAIJ(aij->A, rowemb, dcolemb, pattern, A)); in MatSetSeqMats_MPIAIJ()
2901 PetscCall(MatSetUp(aij->A)); in MatSetSeqMats_MPIAIJ()
2918 PetscCall(PetscHMapIDestroy(&aij->colmap)); in MatSetSeqMats_MPIAIJ()
[all …]
/petsc/src/mat/impls/aij/seq/matlab/
H A Daijmatlab.c12 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)B->data; in MatSeqAIJToMatlab() local
17 mat = mxCreateSparse(B->cmap->n, B->rmap->n, aij->nz, mxREAL); in MatSeqAIJToMatlab()
18 if (PetscArraycpy(mxGetPr(mat), aij->a, aij->nz)) return NULL; in MatSeqAIJToMatlab()
21 for (i = 0; i < aij->nz; i++) jj[i] = aij->j[i]; in MatSeqAIJToMatlab()
23 for (i = 0; i < B->rmap->n + 1; i++) ii[i] = aij->i[i]; in MatSeqAIJToMatlab()
43 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)mat->data; in MatSeqAIJFromMatlab() local
63 if (nz != aij->nz) { in MatSeqAIJFromMatlab()
65 aij->nz = nz; in MatSeqAIJFromMatlab()
66 PetscCall(PetscShmgetAllocateArray(aij->nz, sizeof(PetscScalar), (void **)&aij->a)); in MatSeqAIJFromMatlab()
67 PetscCall(PetscShmgetAllocateArray(aij->nz, sizeof(PetscInt), (void **)&aij->j)); in MatSeqAIJFromMatlab()
[all …]
/petsc/src/mat/impls/baij/mpi/
H A Dmpb_baij.c5 Mat_MPIBAIJ *aij = (Mat_MPIBAIJ *)mat->data; in MatGetMultiProcBlock_MPIBAIJ() local
6 Mat_SeqBAIJ *aijB = (Mat_SeqBAIJ *)aij->B->data; in MatGetMultiProcBlock_MPIBAIJ()
40 PetscCall(PetscCalloc1(aij->B->cmap->n / bs, &garrayCMap)); in MatGetMultiProcBlock_MPIBAIJ()
41 for (i = 0; i < aij->B->cmap->n / bs; i++) { in MatGetMultiProcBlock_MPIBAIJ()
42 col = aij->garray[i]; /* blocked column index */ in MatGetMultiProcBlock_MPIBAIJ()
54 PetscCall(PetscCalloc1(aij->B->rmap->n / bs, &nnz)); in MatGetMultiProcBlock_MPIBAIJ()
55 for (i = 0; i < aij->B->rmap->n / bs; i++) { in MatGetMultiProcBlock_MPIBAIJ()
65 ((Mat_MPIBAIJ *)((*subMat)->data))->A = aij->A; in MatGetMultiProcBlock_MPIBAIJ()
67 PetscCall(PetscObjectReference((PetscObject)aij->A)); in MatGetMultiProcBlock_MPIBAIJ()
68 } else if (((Mat_MPIBAIJ *)(*subMat)->data)->A != aij->A) { in MatGetMultiProcBlock_MPIBAIJ()
[all …]
H A Dmpibaij.c119 Mat_MPIBAIJ *aij = (Mat_MPIBAIJ *)mat->data; in MatStoreValues_MPIBAIJ() local
122 PetscCall(MatStoreValues(aij->A)); in MatStoreValues_MPIBAIJ()
123 PetscCall(MatStoreValues(aij->B)); in MatStoreValues_MPIBAIJ()
129 Mat_MPIBAIJ *aij = (Mat_MPIBAIJ *)mat->data; in MatRetrieveValues_MPIBAIJ() local
132 PetscCall(MatRetrieveValues(aij->A)); in MatRetrieveValues_MPIBAIJ()
133 PetscCall(MatRetrieveValues(aij->B)); in MatRetrieveValues_MPIBAIJ()
1127 Mat_MPIBAIJ *aij = (Mat_MPIBAIJ *)mat->data; in MatView_MPIBAIJ_Binary() local
1128 Mat_SeqBAIJ *A = (Mat_SeqBAIJ *)aij->A->data; in MatView_MPIBAIJ_Binary()
1129 Mat_SeqBAIJ *B = (Mat_SeqBAIJ *)aij->B->data; in MatView_MPIBAIJ_Binary()
1130 const PetscInt *garray = aij->garray; in MatView_MPIBAIJ_Binary()
[all …]
/petsc/src/ksp/pc/impls/cp/
H A Dcp.c20 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)pc->pmat->data; in PCSetUp_CP() local
37 if (!cp->a) PetscCall(PetscMalloc3(aij->nz, &cp->a, cp->n + 1, &cp->i, aij->nz, &cp->j)); in PCSetUp_CP()
40 for (i = 0; i < aij->nz; i++) colcnt[aij->j[i]]++; in PCSetUp_CP()
45 for (j = aij->i[i]; j < aij->i[i + 1]; j++) { /* over columns in row */ in PCSetUp_CP()
46 cp->j[cp->i[aij->j[j]] + colcnt[aij->j[j]]] = i; in PCSetUp_CP()
47 cp->a[cp->i[aij->j[j]] + colcnt[aij->j[j]]++] = aij->a[j]; in PCSetUp_CP()
/petsc/src/mat/tests/
H A Dex101.c10 PetscInt aij[3][3] = { in main() local
47 PetscCall(MatSetValues(A, 3, aij[0], 3, aij[0], pa, ADD_VALUES)); in main()
48 PetscCall(MatSetValues(A, 3, aij[1], 3, aij[1], pa, ADD_VALUES)); in main()
49 PetscCall(MatSetValues(A, 3, aij[2], 3, aij[2], pa, ADD_VALUES)); in main()
H A Dex195.c15 Mat aij; in main() local
108 PetscCall(MatConvert(nest, MATAIJ, MAT_INITIAL_MATRIX, &aij)); in main()
109 PetscCall(MatView(aij, PETSC_VIEWER_STDOUT_WORLD)); in main()
134 PetscCall(MatEqual(nest, aij, &equal)); in main()
150 PetscCall(MatDestroy(&aij)); in main()
H A Dex267.c122 PetscBool aij, inplace = PETSC_FALSE; in main() local
150 PetscCall(PetscObjectTypeCompareAny((PetscObject)A, &aij, MATSEQAIJ, MATSEQBAIJ, "")); in main()
158 if (!aij) { in main()
/petsc/src/ksp/ksp/tutorials/
H A Dex61f.F9081 PetscScalar :: aij
146 aij = 4.0
148 aij = -1.0
150 alist(nz) = aij
185 aij = alist(ip)
186 PetscCallA(MatSetValue(col_f_mat, ii, jj, aij, INSERT_VALUES, ierr))
206 aij = alist(ip)
207 b(i) = b(i) + aij*x(j)
/petsc/src/mat/graphops/color/impls/jp/
H A Djp.c34 Mat_MPIAIJ *aij; in MCJPGreatestWeight_Private() local
55 aij = (Mat_MPIAIJ *)G->data; in MCJPGreatestWeight_Private()
56 dG = aij->A; in MCJPGreatestWeight_Private()
57 oG = aij->B; in MCJPGreatestWeight_Private()
68 PetscCall(PetscSFSetGraphLayout(sf, layout, on, NULL, PETSC_COPY_VALUES, aij->garray)); in MCJPGreatestWeight_Private()
141 Mat_MPIAIJ *aij; in MCJPInitialLocalColor_Private() local
158 aij = (Mat_MPIAIJ *)G->data; in MCJPInitialLocalColor_Private()
159 dG = aij->A; in MCJPInitialLocalColor_Private()
160 oG = aij->B; in MCJPInitialLocalColor_Private()
272 Mat_MPIAIJ *aij; in MCJPMinColor_Private() local
[all …]
/petsc/src/mat/impls/aij/seq/
H A Daij.c37 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatGetColumnReductions_SeqAIJ() local
43 for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] += PetscAbsScalar(aij->a[i] * aij->a[i]); in MatGetColumnReductions_SeqAIJ()
45 for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] += PetscAbsScalar(aij->a[i]); in MatGetColumnReductions_SeqAIJ()
47 …for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] = PetscMax(PetscAbsScalar(aij->a[i]), reduct… in MatGetColumnReductions_SeqAIJ()
49 for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] += PetscRealPart(aij->a[i]); in MatGetColumnReductions_SeqAIJ()
51 for (i = 0; i < aij->i[m]; i++) reductions[aij->j[i]] += PetscImaginaryPart(aij->a[i]); in MatGetColumnReductions_SeqAIJ()
2248 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data, *bij = (Mat_SeqAIJ *)B->data; in MatIsTranspose_SeqAIJ() local
2262 aii = aij->i; in MatIsTranspose_SeqAIJ()
2264 adx = aij->j; in MatIsTranspose_SeqAIJ()
2300 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data, *bij = (Mat_SeqAIJ *)B->data; in MatIsHermitianTranspose_SeqAIJ() local
[all …]
H A Dmatmatmult.c24 Mat_SeqAIJ *aij; in MatSetSeqAIJWithArrays_private() local
38 aij = (Mat_SeqAIJ *)mat->data; in MatSetSeqAIJWithArrays_private()
39 ofree_a = aij->free_a; in MatSetSeqAIJWithArrays_private()
40 ofree_ij = aij->free_ij; in MatSetSeqAIJWithArrays_private()
44 PetscCall(PetscFree(aij->ilen)); in MatSetSeqAIJWithArrays_private()
45 PetscCall(PetscFree(aij->imax)); in MatSetSeqAIJWithArrays_private()
46 PetscCall(PetscMalloc1(m, &aij->imax)); in MatSetSeqAIJWithArrays_private()
47 PetscCall(PetscMalloc1(m, &aij->ilen)); in MatSetSeqAIJWithArrays_private()
48 for (ii = 0, aij->nonzerorowcnt = 0, aij->rmax = 0; ii < m; ii++) { in MatSetSeqAIJWithArrays_private()
50 aij->nonzerorowcnt += !!rnz; in MatSetSeqAIJWithArrays_private()
[all …]
/petsc/src/mat/impls/aij/seq/cholmod/
H A Daijcholmod.c6 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)A->data; in MatWrapCholmod_seqaij() local
9 const PetscInt *ai = aij->i, *aj = aij->j, *adiag; in MatWrapCholmod_seqaij()
/petsc/src/mat/graphops/color/utils/
H A Dweights.c47 Mat_SeqAIJ *aij; in MatColoringGetDegrees() local
61 aij = (Mat_SeqAIJ *)lG->data; in MatColoringGetDegrees()
62 Gi = aij->i; in MatColoringGetDegrees()
63 Gj = aij->j; in MatColoringGetDegrees()
148 Mat_SeqAIJ *aij; in MatColoringCreateSmallestLastWeights() local
165 aij = (Mat_SeqAIJ *)lG->data; in MatColoringCreateSmallestLastWeights()
166 Gi = aij->i; in MatColoringCreateSmallestLastWeights()
167 Gj = aij->j; in MatColoringCreateSmallestLastWeights()
/petsc/src/mat/graphops/color/impls/greedy/
H A Dgreedy.c21 Mat_MPIAIJ *aij = (Mat_MPIAIJ *)m->data; in GreedyColoringLocalDistanceOne_Private() local
46 md = aij->A; in GreedyColoringLocalDistanceOne_Private()
48 mo = aij->B; in GreedyColoringLocalDistanceOne_Private()
70 PetscCall(VecGetSize(aij->lvec, &no)); in GreedyColoringLocalDistanceOne_Private()
83 PetscCall(PetscSFSetGraphLayout(sf, layout, no, NULL, PETSC_COPY_VALUES, aij->garray)); in GreedyColoringLocalDistanceOne_Private()
175 Mat_MPIAIJ *aij = (Mat_MPIAIJ *)m->data; in GreedyColoringLocalDistanceTwo_Private() local
200 md = aij->A; in GreedyColoringLocalDistanceTwo_Private()
202 mo = aij->B; in GreedyColoringLocalDistanceTwo_Private()
240 PetscCall(VecGetLocalSize(aij->lvec, &no)); in GreedyColoringLocalDistanceTwo_Private()
243 PetscCall(PetscSFSetGraphLayout(sf, layout, no, NULL, PETSC_COPY_VALUES, aij->garray)); in GreedyColoringLocalDistanceTwo_Private()
/petsc/src/dm/impls/sliced/
H A Dsliced.c20 PetscBool aij = PETSC_FALSE; in DMCreateMatrix_Sliced() local
33 PetscCall(PetscObjectHasFunction((PetscObject)*J, "MatMPIAIJSetPreallocation_C", &aij)); in DMCreateMatrix_Sliced()
34 if (!aij) PetscCall(PetscObjectHasFunction((PetscObject)*J, "MatSeqAIJSetPreallocation_C", &aij)); in DMCreateMatrix_Sliced()
35 if (aij) { in DMCreateMatrix_Sliced()
/petsc/lib/petsc/bin/maint/
H A Dxclude95 petsc-dist/src/mat/impls/aij/seq/mi25bfac.f
96 petsc-dist/src/mat/impls/aij/seq/mi15blas.f
97 petsc-dist/src/mat/impls/aij/seq/LUSOL.f
98 petsc-dist/src/mat/impls/aij/seq/LUSOL_BLAS.f
/petsc/src/mat/tests/output/
H A Dex72_1.out4 Writing matrix to binary file petscmat.aij using PETSc AIJ format ...
H A Dex72_4.out4 Writing matrix to binary file petscmat.aij using PETSc AIJ format ...
H A Dex72_3.out4 Writing matrix to binary file petscmat2.aij using PETSc AIJ format ...
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaij.c296 Mat aij; in MatView_SeqSBAIJ_ASCII() local
302 PetscCall(MatConvert(A, MATSEQAIJ, MAT_INITIAL_MATRIX, &aij)); in MatView_SeqSBAIJ_ASCII()
304 if (((PetscObject)A)->name) PetscCall(PetscObjectSetName((PetscObject)aij, matname)); in MatView_SeqSBAIJ_ASCII()
305 PetscCall(MatView_SeqAIJ(aij, viewer)); in MatView_SeqSBAIJ_ASCII()
306 PetscCall(MatDestroy(&aij)); in MatView_SeqSBAIJ_ASCII()
1132 Mat_SeqSBAIJ *aij = (Mat_SeqSBAIJ *)Y->data; in MatShift_SeqSBAIJ() local
1135 if (!Y->preallocated || !aij->nz) PetscCall(MatSeqSBAIJSetPreallocation(Y, Y->rmap->bs, 1, NULL)); in MatShift_SeqSBAIJ()
1334 Mat_SeqSBAIJ *aij = (Mat_SeqSBAIJ *)mat->data; in MatStoreValues_SeqSBAIJ() local
1335 PetscInt nz = aij->i[mat->rmap->N] * mat->rmap->bs * aij->bs2; in MatStoreValues_SeqSBAIJ()
1338 …PetscCheck(aij->nonew == 1, PETSC_COMM_SELF, PETSC_ERR_ORDER, "Must call MatSetOption(A,MAT_NEW_NO… in MatStoreValues_SeqSBAIJ()
[all …]

123