| /petsc/src/mat/impls/sell/mpi/ |
| H A D | mpisell.c | 692 static PetscErrorCode MatGetInfo_MPISELL(Mat matin, MatInfoType flag, MatInfo *info) in MatGetInfo_MPISELL() argument
694 Mat_MPISELL *mat = (Mat_MPISELL *)matin->data; in MatGetInfo_MPISELL()
722 …PIU_Allreduce(isend, irecv, 5, MPIU_PETSCLOGDOUBLE, MPI_MAX, PetscObjectComm((PetscObject)matin))); in MatGetInfo_MPISELL()
730 …PIU_Allreduce(isend, irecv, 5, MPIU_PETSCLOGDOUBLE, MPI_SUM, PetscObjectComm((PetscObject)matin))); in MatGetInfo_MPISELL()
1000 static PetscErrorCode MatDuplicate_MPISELL(Mat matin, MatDuplicateOption cpvalues, Mat *newmat) in MatDuplicate_MPISELL() argument
1003 Mat_MPISELL *a, *oldmat = (Mat_MPISELL *)matin->data; in MatDuplicate_MPISELL()
1007 PetscCall(MatCreate(PetscObjectComm((PetscObject)matin), &mat)); in MatDuplicate_MPISELL()
1008 PetscCall(MatSetSizes(mat, matin->rmap->n, matin->cmap->n, matin->rmap->N, matin->cmap->N)); in MatDuplicate_MPISELL()
1009 PetscCall(MatSetBlockSizesFromMats(mat, matin, matin)); in MatDuplicate_MPISELL()
1010 PetscCall(MatSetType(mat, ((PetscObject)matin)->type_name)); in MatDuplicate_MPISELL()
[all …]
|
| /petsc/src/mat/impls/cdiagonal/ |
| H A D | cdiagonal.c | 246 static PetscErrorCode MatSolve_ConstantDiagonal(Mat matin, Vec b, Vec x) in MatSolve_ConstantDiagonal() argument
248 Mat_ConstantDiagonal *ctx = (Mat_ConstantDiagonal *)matin->data; in MatSolve_ConstantDiagonal()
251 if (ctx->diag == 0.0) matin->factorerrortype = MAT_FACTOR_NUMERIC_ZEROPIVOT; in MatSolve_ConstantDiagonal()
252 else matin->factorerrortype = MAT_FACTOR_NOERROR; in MatSolve_ConstantDiagonal()
257 static PetscErrorCode MatSOR_ConstantDiagonal(Mat matin, Vec x, PetscReal omega, MatSORType flag, P… in MatSOR_ConstantDiagonal() argument
260 PetscCall(MatSolve_ConstantDiagonal(matin, x, y)); in MatSOR_ConstantDiagonal()
|
| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | mpisbaij.c | 1171 static PetscErrorCode MatGetRow_MPISBAIJ(Mat matin, PetscInt row, PetscInt *nz, PetscInt **idx, Pet… in MatGetRow_MPISBAIJ() argument
1173 Mat_MPISBAIJ *mat = (Mat_MPISBAIJ *)matin->data; in MatGetRow_MPISBAIJ()
1175 PetscInt bs = matin->rmap->bs, bs2 = mat->bs2, i, *cworkA, *cworkB, **pcA, **pcB; in MatGetRow_MPISBAIJ()
1176 PetscInt nztot, nzA, nzB, lrow, brstart = matin->rmap->rstart, brend = matin->rmap->rend; in MatGetRow_MPISBAIJ()
1438 static PetscErrorCode MatGetInfo_MPISBAIJ(Mat matin, MatInfoType flag, MatInfo *info) in MatGetInfo_MPISBAIJ() argument
1440 Mat_MPISBAIJ *a = (Mat_MPISBAIJ *)matin->data; in MatGetInfo_MPISBAIJ()
1445 info->block_size = (PetscReal)matin->rmap->bs; in MatGetInfo_MPISBAIJ()
1469 …PIU_Allreduce(isend, irecv, 5, MPIU_PETSCLOGDOUBLE, MPI_MAX, PetscObjectComm((PetscObject)matin))); in MatGetInfo_MPISBAIJ()
1477 …PIU_Allreduce(isend, irecv, 5, MPIU_PETSCLOGDOUBLE, MPI_SUM, PetscObjectComm((PetscObject)matin))); in MatGetInfo_MPISBAIJ()
2397 static PetscErrorCode MatDuplicate_MPISBAIJ(Mat matin, MatDuplicateOption cpvalues, Mat *newmat) in MatDuplicate_MPISBAIJ() argument
[all …]
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 1299 static PetscErrorCode MatGetRow_MPIBAIJ(Mat matin, PetscInt row, PetscInt *nz, PetscInt **idx, Pets… in MatGetRow_MPIBAIJ() argument
1301 Mat_MPIBAIJ *mat = (Mat_MPIBAIJ *)matin->data; in MatGetRow_MPIBAIJ()
1303 PetscInt bs = matin->rmap->bs, bs2 = mat->bs2, i, *cworkA, *cworkB, **pcA, **pcB; in MatGetRow_MPIBAIJ()
1304 PetscInt nztot, nzA, nzB, lrow, brstart = matin->rmap->rstart, brend = matin->rmap->rend; in MatGetRow_MPIBAIJ()
1402 static PetscErrorCode MatGetInfo_MPIBAIJ(Mat matin, MatInfoType flag, MatInfo *info) in MatGetInfo_MPIBAIJ() argument
1404 Mat_MPIBAIJ *a = (Mat_MPIBAIJ *)matin->data; in MatGetInfo_MPIBAIJ()
1409 info->block_size = (PetscReal)matin->rmap->bs; in MatGetInfo_MPIBAIJ()
1434 …PIU_Allreduce(isend, irecv, 5, MPIU_PETSCLOGDOUBLE, MPI_MAX, PetscObjectComm((PetscObject)matin))); in MatGetInfo_MPIBAIJ()
1442 …PIU_Allreduce(isend, irecv, 5, MPIU_PETSCLOGDOUBLE, MPI_SUM, PetscObjectComm((PetscObject)matin))); in MatGetInfo_MPIBAIJ()
1449 …} else SETERRQ(PetscObjectComm((PetscObject)matin), PETSC_ERR_ARG_WRONG, "Unknown MatInfoType argu… in MatGetInfo_MPIBAIJ()
[all …]
|
| /petsc/src/mat/impls/diagonal/ |
| H A D | diagonal.c | 532 static PetscErrorCode MatSolve_Diagonal(Mat matin, Vec b, Vec x) in MatSolve_Diagonal() argument
534 Mat_Diagonal *ctx = (Mat_Diagonal *)matin->data; in MatSolve_Diagonal()
537 PetscCall(MatDiagonalSetUpInverseDiagonal(matin)); in MatSolve_Diagonal()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiaij.c | 1399 static PetscErrorCode MatSOR_MPIAIJ(Mat matin, Vec bb, PetscReal omega, MatSORType flag, PetscReal … in MatSOR_MPIAIJ() argument
1401 Mat_MPIAIJ *mat = (Mat_MPIAIJ *)matin->data; in MatSOR_MPIAIJ()
1471 PetscCall(MatCreateVecs(matin, &mat->diag, NULL)); in MatSOR_MPIAIJ()
1472 PetscCall(MatGetDiagonal(matin, mat->diag)); in MatSOR_MPIAIJ()
1474 PetscCall(MatHasOperation(matin, MATOP_MULT_DIAGONAL_BLOCK, &hasop)); in MatSOR_MPIAIJ()
1476 PetscCall(MatMultDiagonalBlock(matin, xx, bb1)); in MatSOR_MPIAIJ()
1488 } else SETERRQ(PetscObjectComm((PetscObject)matin), PETSC_ERR_SUP, "Parallel SOR not supported"); in MatSOR_MPIAIJ()
1492 matin->factorerrortype = mat->A->factorerrortype; in MatSOR_MPIAIJ()
1612 static PetscErrorCode MatGetInfo_MPIAIJ(Mat matin, MatInfoType flag, MatInfo *info) in MatGetInfo_MPIAIJ() argument
1614 Mat_MPIAIJ *mat = (Mat_MPIAIJ *)matin->data; in MatGetInfo_MPIAIJ()
[all …]
|