| /petsc/src/mat/tests/ |
| H A D | ex248.c | 7 Mat A, B, C, K, Ad, Bd; in main() local
15 PetscCall(MatCreateDense(PETSC_COMM_SELF, p, q, p, q, NULL, &Bd)); in main()
17 PetscCall(MatSetRandom(Bd, NULL)); in main()
19 PetscCall(MatFilter(Bd, 0.2, PETSC_FALSE, PETSC_FALSE)); in main()
21 PetscCall(MatConvert(Bd, MATAIJ, MAT_INITIAL_MATRIX, &B)); in main()
26 PetscCall(MatDenseGetArrayRead(Bd, &Bv)); in main()
28 PetscCall(MatDenseRestoreArrayRead(Bd, &Bv)); in main()
33 PetscCall(MatScale(Bd, 0.3)); in main()
35 PetscCall(MatDenseGetArrayRead(Bd, &Bv)); in main()
37 PetscCall(MatDenseRestoreArrayRead(Bd, &Bv)); in main()
[all …]
|
| /petsc/src/ksp/ksp/utils/schurm/ |
| H A D | schurm.c | 481 Mat P, B, C, D, E = NULL, Bd, AinvBd, sub = NULL; in MatSchurComplementComputeExplicitOperator() local
506 PetscCall(MatDuplicate(AinvBd, MAT_DO_NOT_COPY_VALUES, &Bd)); in MatSchurComplementComputeExplicitOperator()
507 … PetscCall(MatDenseGetSubMatrix(Bd, PETSC_DECIDE, PETSC_DECIDE, A->cmap->N, AinvBd->cmap->N, s)); in MatSchurComplementComputeExplicitOperator()
511 PetscCall(MatDenseRestoreSubMatrix(Bd, s)); in MatSchurComplementComputeExplicitOperator()
512 PetscCall(MatDenseGetSubMatrix(Bd, PETSC_DECIDE, PETSC_DECIDE, 0, A->cmap->N, &sub)); in MatSchurComplementComputeExplicitOperator()
514 PetscCall(MatDenseRestoreSubMatrix(Bd, &sub)); in MatSchurComplementComputeExplicitOperator()
515 } else PetscCall(MatConvert(B, mtype, MAT_INITIAL_MATRIX, &Bd)); in MatSchurComplementComputeExplicitOperator()
520 PetscCall(MatConvert(B, mtype, MAT_INITIAL_MATRIX, &Bd)); in MatSchurComplementComputeExplicitOperator()
522 PetscCall(KSPMatSolve(ksp, Bd, AinvBd)); in MatSchurComplementComputeExplicitOperator()
554 PetscCall(MatDestroy(&Bd)); in MatSchurComplementComputeExplicitOperator()
|
| /petsc/src/mat/impls/aij/mpi/kokkos/ |
| H A D | mpiaijkok.kokkos.cxx | 1054 KokkosCsrMatrix Adt, Aot, Ad, Ao, Bd, Bo; in MatProductSymbolic_MPIAIJKokkos_AtB() local
1064 PetscCall(MatSeqAIJKokkosGetKokkosCsrMatrix(bmpi->A, &Bd)); in MatProductSymbolic_MPIAIJKokkos_AtB()
1083 PetscCallCXX(KokkosSparse::spgemm_symbolic(mm->kh3, Aot, false, Bd, false, mm->C3)); in MatProductSymbolic_MPIAIJKokkos_AtB()
1087 PetscCallCXX(KokkosSparse::spgemm_numeric(mm->kh3, Aot, false, Bd, false, mm->C3)); in MatProductSymbolic_MPIAIJKokkos_AtB()
1101 PetscCallCXX(KokkosSparse::spgemm_symbolic(mm->kh1, Adt, false, Bd, false, mm->C1)); in MatProductSymbolic_MPIAIJKokkos_AtB()
1103 PetscCallCXX(KokkosSparse::spgemm_numeric(mm->kh1, Adt, false, Bd, false, mm->C1)); in MatProductSymbolic_MPIAIJKokkos_AtB()
1132 KokkosCsrMatrix Adt, Aot, Bd, Bo; in MatProductNumeric_MPIAIJKokkos_AtB() local
1139 PetscCall(MatSeqAIJKokkosGetKokkosCsrMatrix(bmpi->A, &Bd)); in MatProductNumeric_MPIAIJKokkos_AtB()
1143 PetscCallCXX(KokkosSparse::spgemm_numeric(mm->kh3, Aot, false, Bd, false, mm->C3)); in MatProductNumeric_MPIAIJKokkos_AtB()
1150 PetscCallCXX(KokkosSparse::spgemm_numeric(mm->kh1, Adt, false, Bd, false, mm->C1)); in MatProductNumeric_MPIAIJKokkos_AtB()
[all …]
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcschurs.c | 268 Mat B, C, D, Bd, Cd, AinvBd; in PCBDDCComputeExplicitSchur() local
296 PetscCall(MatConvert(B, MATSEQDENSE, MAT_INITIAL_MATRIX, &Bd)); in PCBDDCComputeExplicitSchur()
298 Bd = B; in PCBDDCComputeExplicitSchur()
305 PetscCall(MatDuplicate(Bd, MAT_DO_NOT_COPY_VALUES, &AinvBd)); in PCBDDCComputeExplicitSchur()
306 PetscCall(MatMatSolve(fact, Bd, AinvBd)); in PCBDDCComputeExplicitSchur()
314 PetscCall(MatMatMult(Ainvd, Bd, MAT_INITIAL_MATRIX, fill, &AinvBd)); in PCBDDCComputeExplicitSchur()
321 PetscCall(MatDuplicate(Bd, MAT_DO_NOT_COPY_VALUES, &AinvBd)); in PCBDDCComputeExplicitSchur()
322 PetscCall(MatGetSize(Bd, &n, &nrhs)); in PCBDDCComputeExplicitSchur()
323 PetscCall(MatDenseGetArray(Bd, &arrayrhs)); in PCBDDCComputeExplicitSchur()
334 PetscCall(MatDenseRestoreArray(Bd, &arrayrhs)); in PCBDDCComputeExplicitSchur()
[all …]
|