| /petsc/src/mat/tests/ |
| H A D | ex49.c | 7 Mat mat, tmat = 0; in main() local
61 tmat = mat; in main()
64 PetscCall(MatTranspose(mat, MAT_INITIAL_MATRIX, &tmat)); in main()
68 PetscCall(MatGetInfo(tmat, MAT_GLOBAL_SUM, &info)); in main()
70 PetscCall(MatNorm(tmat, NORM_FROBENIUS, &normf)); in main()
71 PetscCall(MatNorm(tmat, NORM_1, &norm1)); in main()
72 PetscCall(MatNorm(tmat, NORM_INFINITY, &normi)); in main()
74 PetscCall(MatView(tmat, PETSC_VIEWER_STDOUT_WORLD)); in main()
81 PetscCall(MatAXPY(tmat, alpha, mat, DIFFERENT_NONZERO_PATTERN)); in main()
82 PetscCall(MatView(tmat, PETSC_VIEWER_STDOUT_WORLD)); in main()
[all …]
|
| H A D | ex2.c | 80 Mat mat, tmat = 0; in main() local
131 tmat = mat; in main()
134 PetscCall(MatTranspose(mat, MAT_INITIAL_MATRIX, &tmat)); in main()
139 PetscCall(MatNorm(tmat, NORM_FROBENIUS, &normf)); in main()
140 PetscCall(MatNorm(tmat, NORM_1, &norm1)); in main()
141 PetscCall(MatNorm(tmat, NORM_INFINITY, &normi)); in main()
143 PetscCall(MatView(tmat, PETSC_VIEWER_STDOUT_WORLD)); in main()
150 PetscCall(MatAXPY(tmat, alpha, mat, DIFFERENT_NONZERO_PATTERN)); in main()
151 PetscCall(MatView(tmat, PETSC_VIEWER_STDOUT_WORLD)); in main()
154 PetscCall(MatAYPX(tmat, alpha, mat, DIFFERENT_NONZERO_PATTERN)); in main()
[all …]
|
| /petsc/src/mat/impls/composite/ |
| H A D | mcomposite.c | 656 Mat tmat, newmat; in MatCompositeMerge_Composite() local
670 PetscCall(MatDuplicate(next->mat, MAT_COPY_VALUES, &tmat)); in MatCompositeMerge_Composite()
671 if (shell->scalings) PetscCall(MatScale(tmat, shell->scalings[i++])); in MatCompositeMerge_Composite()
672 …while ((next = next->next)) PetscCall(MatAXPY(tmat, shell->scalings ? shell->scalings[i++] : 1.0, … in MatCompositeMerge_Composite()
675 PetscCall(MatDuplicate(prev->mat, MAT_COPY_VALUES, &tmat)); in MatCompositeMerge_Composite()
676 if (shell->scalings) PetscCall(MatScale(tmat, shell->scalings[i--])); in MatCompositeMerge_Composite()
677 …while ((prev = prev->prev)) PetscCall(MatAXPY(tmat, shell->scalings ? shell->scalings[i--] : 1.0, … in MatCompositeMerge_Composite()
681 PetscCall(MatDuplicate(next->mat, MAT_COPY_VALUES, &tmat)); in MatCompositeMerge_Composite()
683 PetscCall(MatMatMult(next->mat, tmat, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &newmat)); in MatCompositeMerge_Composite()
684 PetscCall(MatDestroy(&tmat)); in MatCompositeMerge_Composite()
[all …]
|
| /petsc/src/ksp/pc/impls/asm/ftn-custom/ |
| H A D | zasmf.c | 35 Mat *tmat; in pcasmgetlocalsubmatrices_() local
39 *ierr = PCASMGetLocalSubmatrices(*pc, &nloc, &tmat); in pcasmgetlocalsubmatrices_()
41 …if (mat) *ierr = F90Array1dCreate(tmat, MPIU_FORTRANADDR, 1, nloc, mat PETSC_F90_2PTR_PARAM(ptrd)); in pcasmgetlocalsubmatrices_()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex65.c | 112 Mat tmat; in CreateRestriction() local
117 PetscCall(DMCreateInterpolation(da1, da2, &tmat, NULL)); in CreateRestriction()
118 PetscCall(MatTranspose(tmat, MAT_INITIAL_MATRIX, mat)); in CreateRestriction()
119 PetscCall(MatDestroy(&tmat)); in CreateRestriction()
|
| /petsc/src/mat/graphops/order/ |
| H A D | sorder.c | 227 Mat tmat; in MatGetOrdering() local
228 PetscCall(MatPermute(mat, *rperm, *cperm, &tmat)); in MatGetOrdering()
229 PetscCall(MatViewFromOptions(tmat, (PetscObject)mat, "-mat_view_ordering")); in MatGetOrdering()
230 PetscCall(MatDestroy(&tmat)); in MatGetOrdering()
|
| /petsc/src/ksp/pc/impls/redistribute/ |
| H A D | redistribute.c | 87 Mat tmat; in PCSetUp_Redistribute() local
95 PetscCall(KSPGetOperators(red->ksp, NULL, &tmat)); in PCSetUp_Redistribute()
96 PetscCall(MatCreateSubMatrix(pc->pmat, red->is, red->is, MAT_REUSE_MATRIX, &tmat)); in PCSetUp_Redistribute()
97 PetscCall(KSPSetOperators(red->ksp, tmat, tmat)); in PCSetUp_Redistribute()
307 PetscCall(MatCreateSubMatrix(pc->pmat, red->is, red->is, MAT_INITIAL_MATRIX, &tmat)); in PCSetUp_Redistribute()
308 PetscCall(KSPSetOperators(red->ksp, tmat, tmat)); in PCSetUp_Redistribute()
309 PetscCall(MatDestroy(&tmat)); in PCSetUp_Redistribute()
|
| /petsc/src/mat/impls/dense/seq/ |
| H A D | dense.c | 1826 Mat tmat; in MatTranspose_SeqDense() local
1832 PetscCall(MatCreate(PetscObjectComm((PetscObject)A), &tmat)); in MatTranspose_SeqDense()
1833 PetscCall(MatSetSizes(tmat, A->cmap->n, A->rmap->n, A->cmap->n, A->rmap->n)); in MatTranspose_SeqDense()
1834 PetscCall(MatSetType(tmat, ((PetscObject)A)->type_name)); in MatTranspose_SeqDense()
1835 PetscCall(MatSeqDenseSetPreallocation(tmat, NULL)); in MatTranspose_SeqDense()
1836 } else tmat = *matout; in MatTranspose_SeqDense()
1839 PetscCall(MatDenseGetArray(tmat, &v2)); in MatTranspose_SeqDense()
1840 tmatd = (Mat_SeqDense *)tmat->data; in MatTranspose_SeqDense()
1845 PetscCall(MatDenseRestoreArray(tmat, &v2)); in MatTranspose_SeqDense()
1847 PetscCall(MatAssemblyBegin(tmat, MAT_FINAL_ASSEMBLY)); in MatTranspose_SeqDense()
[all …]
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 5352 Mat tmat; in PCBDDCSetUpCorrection() local
5354 …teSeqDense(PETSC_COMM_SELF, pcis->n_B, pcbddc->local_primal_size - pcbddc->benign_n, data, &tmat)); in PCBDDCSetUpCorrection()
5356 PetscCall(MatMatMult(C_B, tmat, MAT_INITIAL_MATRIX, 1.0, &CPHI)); in PCBDDCSetUpCorrection()
5357 PetscCall(MatDestroy(&tmat)); in PCBDDCSetUpCorrection()
7342 Mat tmat; in PCBDDCConstraintsSetUp() local
7346 PetscCall(MatDuplicate(pc->pmat, MAT_DO_NOT_COPY_VALUES, &tmat)); in PCBDDCConstraintsSetUp()
7347 PetscCall(MatISSetLocalMat(tmat, localChangeOfBasisMatrix)); in PCBDDCConstraintsSetUp()
7348 PetscCall(MatAssemblyBegin(tmat, MAT_FINAL_ASSEMBLY)); in PCBDDCConstraintsSetUp()
7349 PetscCall(MatAssemblyEnd(tmat, MAT_FINAL_ASSEMBLY)); in PCBDDCConstraintsSetUp()
7350 PetscCall(MatConvert(tmat, MATAIJ, MAT_INITIAL_MATRIX, &pcbddc->ChangeOfBasisMatrix)); in PCBDDCConstraintsSetUp()
[all …]
|