| /petsc/src/mat/tests/ |
| H A D | ex51.c | 7 Mat A, B, E, Bt, *submatA, *submatB; in main() local
69 PetscCall(MatTranspose(B, MAT_INITIAL_MATRIX, &Bt)); in main()
70 PetscCall(MatAXPY(Bt, 1.0, B, DIFFERENT_NONZERO_PATTERN)); in main()
71 PetscCall(MatSetOption(Bt, MAT_SYMMETRIC, PETSC_TRUE)); in main()
72 PetscCall(MatConvert(Bt, MATSBAIJ, MAT_INITIAL_MATRIX, &E)); in main()
74 PetscCall(MatDestroy(&Bt)); in main()
|
| H A D | ex54.c | 7 Mat E, A, B, Bt, *submatA, *submatB; in main() local
90 PetscCall(MatTranspose(A, MAT_INITIAL_MATRIX, &Bt)); in main()
91 PetscCall(MatAXPY(Bt, 1.0, B, DIFFERENT_NONZERO_PATTERN)); in main()
92 PetscCall(MatSetOption(Bt, MAT_SYMMETRIC, PETSC_TRUE)); in main()
93 PetscCall(MatConvert(Bt, MATSBAIJ, MAT_INITIAL_MATRIX, &E)); in main()
95 PetscCall(MatDestroy(&Bt)); in main()
|
| H A D | ex70.c | 204 Mat X, B, A, Bt, T, T2, PtAP = NULL, RARt = NULL, R = NULL; in main() local
304 PetscCall(MatCreateDense(PETSC_COMM_WORLD, k, n, K, N, dataBt, &Bt)); in main()
305 if (local) PetscCall(MatDenseSetLDA(Bt, k + ldr)); in main()
310 PetscCall(MatDenseGetArrayRead(Bt, (const PetscScalar **)&dataBt)); in main()
314 PetscCall(MatDenseRestoreArrayRead(Bt, (const PetscScalar **)&dataBt)); in main()
325 PetscCall(MatSetRandom(Bt, NULL)); in main()
327 PetscCall(CheckLocal(Bt, B, aBt, aB)); in main()
332 PetscCall(MatConvert(Bt, MATDENSECUDA, MAT_INPLACE_MATRIX, &Bt)); in main()
431 PetscCall(MatMatTransposeMult(A, Bt, MAT_REUSE_MATRIX, PETSC_DETERMINE, &X)); in main()
432 PetscCall(CheckLocal(Bt, X, aBt, aX)); in main()
[all …]
|
| H A D | ex154.c | 45 Mat Bt; in main() local
47 PetscCall(MatCreateTranspose(B, &Bt)); in main()
49 B = Bt; in main()
|
| H A D | ex53.c | 7 Mat A, B, C, At, Bt; in main() local
163 PetscCall(MatTranspose(B, MAT_INITIAL_MATRIX, &Bt)); in main()
167 PetscCall(MatMult(Bt, xx, s2)); in main()
174 PetscCall(MatDestroy(&Bt)); in main()
|
| /petsc/src/mat/utils/ |
| H A D | multequal.c | 116 …atMatMultEqual_Private(Mat A, Mat B, Mat C, PetscInt n, PetscBool *flg, PetscBool At, PetscBool Bt) in MatMatMultEqual_Private() argument
137 PetscValidLogicalCollectiveBool(B, Bt, 7); in MatMatMultEqual_Private()
146 if (Bt) { in MatMatMultEqual_Private()
154 sop = sops[(At ? 1 : 0) + 2 * (Bt ? 1 : 0)]; in MatMatMultEqual_Private()
158 if (Bt) { in MatMatMultEqual_Private()
173 if (Bt) { in MatMatMultEqual_Private()
|
| /petsc/src/mat/impls/nest/ |
| H A D | matnest.c | 1818 Mat Bt; in MatConvert_Nest_SeqAIJ_fast() local
1820 PetscCall(MatTransposeGetMat(B, &Bt)); in MatConvert_Nest_SeqAIJ_fast()
1821 PetscCall(MatTranspose(Bt, MAT_INITIAL_MATRIX, &trans[i * nest->nc + j])); in MatConvert_Nest_SeqAIJ_fast()
1826 Mat Bt; in MatConvert_Nest_SeqAIJ_fast() local
1828 PetscCall(MatHermitianTransposeGetMat(B, &Bt)); in MatConvert_Nest_SeqAIJ_fast()
1829 PetscCall(MatHermitianTranspose(Bt, MAT_INITIAL_MATRIX, &trans[i * nest->nc + j])); in MatConvert_Nest_SeqAIJ_fast()
2042 Mat Bt; in MatConvert_Nest_AIJ() local
2044 PetscCall(MatTransposeGetMat(B, &Bt)); in MatConvert_Nest_AIJ()
2045 PetscCall(PetscObjectTypeCompare((PetscObject)Bt, MATSEQAIJ, &fast)); in MatConvert_Nest_AIJ()
2049 Mat Bt; in MatConvert_Nest_AIJ() local
[all …]
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | fsarkimex.c | 82 const PetscReal *Bt = tab->binterpt, *B = tab->binterp; in TSExtrapolate_ARKIMEX_FastSlowSplit() local
86 …PetscCheck(Bt && B, PetscObjectComm((PetscObject)ts), PETSC_ERR_SUP, "TSARKIMEX %s does not have a… in TSExtrapolate_ARKIMEX_FastSlowSplit()
93 bt[i] += h * Bt[i * pinterp + j] * tt; in TSExtrapolate_ARKIMEX_FastSlowSplit()
|
| H A D | arkimex.c | 1677 const PetscReal *Bt = tab->binterpt, *B = tab->binterp; in TSInterpolate_ARKIMEX() local
1680 …PetscCheck(Bt && B, PetscObjectComm((PetscObject)ts), PETSC_ERR_SUP, "%s %s does not have an inter… in TSInterpolate_ARKIMEX()
1697 bt[i] += h * Bt[i * pinterp + j] * tt; in TSInterpolate_ARKIMEX()
1718 const PetscReal *Bt = tab->binterpt, *B = tab->binterp; in TSExtrapolate_ARKIMEX() local
1721 …PetscCheck(Bt && B, PetscObjectComm((PetscObject)ts), PETSC_ERR_SUP, "TSARKIMEX %s does not have a… in TSExtrapolate_ARKIMEX()
1728 bt[i] += h * Bt[i * pinterp + j] * tt; in TSExtrapolate_ARKIMEX()
|
| /petsc/src/mat/impls/shell/ |
| H A D | shell.c | 489 Mat Bt; member
500 PetscCall(MatDestroy(&mmdata->Bt)); in MatProductCtxDestroy_MatMatShell()
601 MatReuse reuse = mdata->Bt ? MAT_REUSE_MATRIX : MAT_INITIAL_MATRIX; in MatProductNumeric_Shell_X()
602 MatStructure str = mdata->Bt ? SUBSET_NONZERO_PATTERN : DIFFERENT_NONZERO_PATTERN; in MatProductNumeric_Shell_X()
604 PetscCall(MatTranspose(mdata->B, reuse, &mdata->Bt)); in MatProductNumeric_Shell_X()
605 PetscCall(MatDiagonalScale(mdata->Bt, shell->left_work, NULL)); in MatProductNumeric_Shell_X()
606 PetscCall(MatAXPY(D, 1.0, mdata->Bt, str)); in MatProductNumeric_Shell_X()
|
| /petsc/src/mat/impls/aij/mpi/mumps/ |
| H A D | mumps.c | 2048 Mat Bt = NULL; in MatMatSolve_MUMPS() local
2086 PetscCall(MatTransposeGetMat(B, &Bt)); in MatMatSolve_MUMPS()
2107 PetscCall(MatSeqAIJGetArray(Bt, &aa)); in MatMatSolve_MUMPS()
2108 …PetscCall(MatGetRowIJ(Bt, 1, PETSC_FALSE, PETSC_FALSE, &spnr, (const PetscInt **)&ia, (const Petsc… in MatMatSolve_MUMPS()
2137 PetscCall(MatSeqAIJRestoreArray(Bt, &aa)); in MatMatSolve_MUMPS()
2138 …PetscCall(MatRestoreRowIJ(Bt, 1, PETSC_FALSE, PETSC_FALSE, &spnr, (const PetscInt **)&ia, (const P… in MatMatSolve_MUMPS()
2223 b = (Mat_MPIAIJ *)Bt->data; in MatMatSolve_MUMPS()
2303 b = (Mat_MPIAIJ *)Bt->data; in MatMatSolve_MUMPS()
2331 static PetscErrorCode MatMatTransposeSolve_MUMPS(Mat A, Mat Bt, Mat X) in MatMatTransposeSolve_MUMPS() argument
2337 PetscCall(PetscObjectTypeCompareAny((PetscObject)Bt, &flg, MATSEQAIJ, MATMPIAIJ, NULL)); in MatMatTransposeSolve_MUMPS()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | matmatmult.c | 1215 Mat Bt; in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ() local
1225 PetscCall(MatTransposeSymbolic(B, &Bt)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1226 PetscCall(MatSetBlockSizes(Bt, A->cmap->bs, B->cmap->bs)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1227 PetscCall(MatSetType(Bt, ((PetscObject)A)->type_name)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1232 PetscCall(MatMatMultSymbolic_SeqAIJ_SeqAIJ(A, Bt, fill, C)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
1294 PetscCall(MatDestroy(&Bt)); in MatMatTransposeMultSymbolic_SeqAIJ_SeqAIJ()
|
| /petsc/src/ts/impls/rosw/ |
| H A D | rosw.c | 1247 const PetscReal *Bt = ros->tableau->binterpt; in TSInterpolate_RosW() local
1253 …PetscCheck(Bt, PetscObjectComm((PetscObject)ts), PETSC_ERR_SUP, "TSRosW %s does not have an interp… in TSInterpolate_RosW()
1271 for (i = 0; i < s; i++) bt[i] += Bt[i * pinterp + j] * tt; in TSInterpolate_RosW()
|
| /petsc/src/mat/interface/ |
| H A D | matrix.c | 3835 PetscErrorCode MatMatTransposeSolve(Mat A, Mat Bt, Mat X) in MatMatTransposeSolve() argument
3840 PetscValidHeaderSpecific(Bt, MAT_CLASSID, 2); in MatMatTransposeSolve()
3842 PetscCheckSameComm(A, 1, Bt, 2); in MatMatTransposeSolve()
3845 …PetscCheck(X != Bt, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_IDN, "X and B must be different… in MatMatTransposeSolve()
3847 …>N == Bt->cmap->N, PetscObjectComm((PetscObject)A), PETSC_ERR_ARG_SIZ, "Mat A,Mat Bt: global dim %… in MatMatTransposeSolve()
3848 …PetscCheck(X->cmap->N >= Bt->rmap->N, PetscObjectComm((PetscObject)X), PETSC_ERR_ARG_SIZ, "Solutio… in MatMatTransposeSolve()
3853 PetscCall(PetscLogEventBegin(MAT_MatTrSolve, A, Bt, X, 0)); in MatMatTransposeSolve()
3854 PetscUseTypeMethod(A, mattransposesolve, Bt, X); in MatMatTransposeSolve()
3855 PetscCall(PetscLogEventEnd(MAT_MatTrSolve, A, Bt, X, 0)); in MatMatTransposeSolve()
|
| /petsc/src/mat/impls/aij/seq/seqcusparse/ |
| H A D | aijcusparse.cu | 2598 PetscScalar *Bt; member
2626 PetscCallCUDA(cudaFree(mmdata->Bt)); in MatProductCtxDestroy_MatMatCusparse()
2794 …>n, &PETSC_CUSPARSE_ONE, barray, blda, &PETSC_CUSPARSE_ZERO, barray, blda, mmdata->Bt, B->cmap->n); in MatProductNumeric_SeqAIJCUSPARSE_SeqDENSECUDA()
2803 …offsets->data().get(), csrmat->column_indices->data().get(), mmdata->Bt ? mmdata->Bt : barray, bld… in MatProductNumeric_SeqAIJCUSPARSE_SeqDENSECUDA()
2884 …oduct->type == MATPRODUCT_RARt) PetscCallCUDA(cudaMalloc((void **)&mmdata->Bt, (size_t)B->rmap->n … in MatProductSymbolic_SeqAIJCUSPARSE_SeqDENSECUDA()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiaij.c | 6003 Mat At, Bt, Ct; in MatMatMultNumeric_MPIDense_MPIAIJ() local
6007 PetscCall(MatTranspose(B, MAT_INITIAL_MATRIX, &Bt)); in MatMatMultNumeric_MPIDense_MPIAIJ()
6008 PetscCall(MatMatMult(Bt, At, MAT_INITIAL_MATRIX, PETSC_CURRENT, &Ct)); in MatMatMultNumeric_MPIDense_MPIAIJ()
6010 PetscCall(MatDestroy(&Bt)); in MatMatMultNumeric_MPIDense_MPIAIJ()
|
| /petsc/src/mat/impls/aij/seq/seqhipsparse/ |
| H A D | aijhipsparse.hip.cxx | 2174 PetscScalar *Bt; member
2197 PetscCallHIP(hipFree(mmdata->Bt)); in MatProductCtxDestroy_MatMatHipsparse()
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 2353 Mat Bt; in MatSeqAIJCompress() local
2385 PetscCall(MatCreateSeqAIJWithArrays(PetscObjectComm((PetscObject)A), n, m, bii, bij, bdata, &Bt)); in MatSeqAIJCompress()
2388 Mat_SeqAIJ *b = (Mat_SeqAIJ *)Bt->data; in MatSeqAIJCompress()
2393 *B = Bt; in MatSeqAIJCompress()
|