Home
last modified time | relevance | path

Searched refs:bn (Results 1 – 25 of 31) sorted by relevance

12

/petsc/src/ksp/ksp/impls/gmres/
H A Dgmreig.c8 PetscBLASInt bn, bN, lwork, idummy, lierr; in KSPComputeExtremeSingularValues_GMRES() local
13 PetscCall(PetscBLASIntCast(n, &bn)); in KSPComputeExtremeSingularValues_GMRES()
30 …PetscCallBLAS("LAPACKgesvd", LAPACKgesvd_("N", "N", &bn, &bn, R, &bN, realpart, &sdummy, &idummy, … in KSPComputeExtremeSingularValues_GMRES()
32 …PetscCallBLAS("LAPACKgesvd", LAPACKgesvd_("N", "N", &bn, &bn, R, &bN, realpart, &sdummy, &idummy, … in KSPComputeExtremeSingularValues_GMRES()
47 PetscBLASInt bn, bN, lwork, idummy, lierr = -1; in KSPComputeEigenvalues_GMRES() local
52 PetscCall(PetscBLASIntCast(n, &bn)); in KSPComputeEigenvalues_GMRES()
66 …PetscCallBLAS("LAPACKgeev", LAPACKgeev_("N", "N", &bn, R, &bN, realpart, imagpart, &sdummy, &idumm… in KSPComputeEigenvalues_GMRES()
81 PetscBLASInt bn, bN, lwork, idummy, lierr = -1; in KSPComputeEigenvalues_GMRES()
84 PetscCall(PetscBLASIntCast(n, &bn)); in KSPComputeEigenvalues_GMRES()
98 …PetscCallBLAS("LAPACKgeev", LAPACKgeev_("N", "N", &bn, R, &bN, eigs, &sdummy, &idummy, &sdummy, &i… in KSPComputeEigenvalues_GMRES()
[all …]
/petsc/src/vec/vec/impls/seq/
H A Dbvec1.c18 PetscBLASInt bn; in VecXDot_Seq_Private() local
21 PetscCall(PetscBLASIntCast(n, &bn)); in VecXDot_Seq_Private()
27 PetscCallBLAS("BLASdot", *z = BLASfn(&bn, ya, &one, xa, &one)); in VecXDot_Seq_Private()
58 PetscBLASInt bn; in VecScale_Seq() local
61 PetscCall(PetscBLASIntCast(xin->map->n, &bn)); in VecScale_Seq()
62 PetscCall(PetscLogFlops(bn)); in VecScale_Seq()
64 PetscCallBLAS("BLASscal", BLASscal_(&bn, &alpha, xarray, &one)); in VecScale_Seq()
78 PetscBLASInt bn; in VecAXPY_Seq() local
80 PetscCall(PetscBLASIntCast(yin->map->n, &bn)); in VecAXPY_Seq()
81 PetscCall(PetscLogFlops(2.0 * bn)); in VecAXPY_Seq()
[all …]
H A Dbvec2.c173 PetscBLASInt bn; in VecSwap_Seq() local
175 PetscCall(PetscBLASIntCast(xin->map->n, &bn)); in VecSwap_Seq()
178 PetscCallBLAS("BLASswap", BLASswap_(&bn, xa, &one, ya, &one)); in VecSwap_Seq()
195 PetscBLASInt bn = 0; in VecNorm_Seq() local
197 PetscCall(PetscBLASIntCast(n, &bn)); in VecNorm_Seq()
202 PetscCallBLAS("BLASnrm2", ztmp[type == NORM_1_AND_2] = BLASnrm2_(&bn, xx, &one)); in VecNorm_Seq()
204 …ASdot", ztmp[type == NORM_1_AND_2] = PetscSqrtReal(PetscRealPart(BLASdot_(&bn, xx, &one, xx, &one)… in VecNorm_Seq()
223 PetscCallBLAS("BLASasum", ztmp[0] = BLASasum_(&bn, xx, &one)); in VecNorm_Seq()
/petsc/src/ksp/ksp/impls/cg/
H A Dcgeig.c15 PetscBLASInt bn, lierr = 0, ldz = 1; in KSPComputeEigenvalues_CG() local
33 PetscCall(PetscBLASIntCast(n, &bn)); in KSPComputeEigenvalues_CG()
35 PetscCallBLAS("LAPACKREALstev", LAPACKREALstev_("N", &bn, r, ee, NULL, &ldz, NULL, &lierr)); in KSPComputeEigenvalues_CG()
48 PetscBLASInt bn, lierr = 0, ldz = 1; in KSPComputeExtremeSingularValues_CG() local
66 PetscCall(PetscBLASIntCast(n, &bn)); in KSPComputeExtremeSingularValues_CG()
68 PetscCallBLAS("LAPACKREALstev", LAPACKREALstev_("N", &bn, dd, ee, NULL, &ldz, NULL, &lierr)); in KSPComputeExtremeSingularValues_CG()
/petsc/src/ksp/ksp/utils/lmvm/blas_cyclic/
H A Dblas_cyclic.c260 PetscBLASInt bn = i_next - i_oldest; in TRSVCyclic_Private() local
263 if (x != y) PetscCall(PetscArraycpy(&y[i_oldest], &x[i_oldest], bn)); in TRSVCyclic_Private()
264 …PetscCallBLAS("BLAStrsv", BLAStrsv_("U", trans, "N", &bn, &A[i_oldest * (lda + 1)], &blda, &y[i_ol… in TRSVCyclic_Private()
266 PetscBLASInt bn = bm - i_oldest; in TRSVCyclic_Private() local
269 PetscCall(PetscArraycpy(&y[i_oldest], &x[i_oldest], bn)); in TRSVCyclic_Private()
273 …if (i_next > 0 && bn > 0) PetscCallBLAS("BLASgemv", BLASgemv_("N", &bn, &i_next, &minus_one, &A[i_… in TRSVCyclic_Private()
274 …if (bn > 0) PetscCallBLAS("BLAStrsv", BLAStrsv_("U", "N", "N", &bn, &A[i_oldest * (lda + 1)], &bld… in TRSVCyclic_Private()
276 …if (bn > 0) PetscCallBLAS("BLAStrsv", BLAStrsv_("U", "C", "N", &bn, &A[i_oldest * (lda + 1)], &bld… in TRSVCyclic_Private()
277 …if (i_next > 0 && bn > 0) PetscCallBLAS("BLASgemv", BLASgemv_("C", &bn, &i_next, &minus_one, &A[i_… in TRSVCyclic_Private()
355 PetscBLASInt bn = i_next - i_oldest; in HEMVCyclic_Private() local
[all …]
/petsc/src/mat/utils/
H A Dmultequal.c14 PetscInt am, an, bm, bn, k; in MatMultEqual_Private() local
30 PetscCall(MatGetLocalSize(B, &bm, &bn)); in MatMultEqual_Private()
31 …= bn, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Mat A,Mat B: local dim %" PetscInt_FMT " %" PetscInt_FM… in MatMultEqual_Private()
121 PetscInt am, an, bm, bn, cm, cn, k; in MatMatMultEqual_Private() local
139 PetscCall(MatGetLocalSize(B, &bm, &bn)); in MatMatMultEqual_Private()
147 PetscInt tt = bn; in MatMatMultEqual_Private()
148 bn = bm; in MatMatMultEqual_Private()
151 …bn == cn, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Mat A, B, C local dim %" PetscInt_FMT " %" PetscInt… in MatMatMultEqual_Private()
442 PetscInt i, am, an, bm, bn, cm, cn; in MatProjMultEqual_Private() local
447 PetscCall(MatGetLocalSize(B, &bm, &bn)); in MatProjMultEqual_Private()
[all …]
/petsc/src/mat/tests/
H A Dex120.c20 PetscBLASInt nn, nevs, il, iu, *iwork, *ifail, lwork, lierr, bn, one = 1; in main() local
134 PetscCall(PetscBLASIntCast(n, &bn)); in main()
142 LAPACKsyev_("V", "U", &bn, arrayA, &bn, evals, work, &lwork, rwork, &lierr); in main()
163 …LAPACKsyevx_("V", "I", "U", &bn, arrayA, &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_arr… in main()
172 LAPACKsygv_(&one, "V", "U", &bn, arrayA, &bn, arrayB, &bn, evals, work, &lwork, rwork, &lierr); in main()
192 …LAPACKsygvx_(&one, "V", "I", "U", &bn, arrayA, &bn, arrayB, &bn, &vl, &vu, &il, &iu, &abstol, &nev… in main()
H A Dex116.c28 PetscBLASInt *iwork, *ifail, lwork, lierr, bn; in main() local
66 PetscCall(PetscBLASIntCast(n, &bn)); in main()
73 LAPACKsyev_("V", "U", &bn, arrayA, &bn, evals, work, &lwork, &lierr); in main()
90 …LAPACKsyevx_("V", "I", "U", &bn, arrayA, &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_arr… in main()
/petsc/src/sys/utils/
H A Dsorti.c1038 PetscErrorCode PetscMergeIntArray(PetscInt an, const PetscInt aI[], PetscInt bn, const PetscInt bI[… in PetscMergeIntArray() argument
1044 PetscCall(PetscMalloc1(an + bn, L)); in PetscMergeIntArray()
1048 while (ak < an && bk < bn) { in PetscMergeIntArray()
1068 if (bk < bn) { in PetscMergeIntArray()
1069 PetscCall(PetscArraycpy(L_ + k, bI + bk, bn - bk)); in PetscMergeIntArray()
1070 k += (bn - bk); in PetscMergeIntArray()
1103 …tArrayPair(PetscInt an, const PetscInt aI[], const PetscInt aJ[], PetscInt bn, const PetscInt bI[]… in PetscMergeIntArrayPair() argument
1110 n_ = an + bn; in PetscMergeIntArrayPair()
1117 while (ak < an && bk < bn) { in PetscMergeIntArrayPair()
1135 if (bk < bn) { in PetscMergeIntArrayPair()
[all …]
/petsc/src/ksp/ksp/impls/gmres/dgmres/
H A Ddgmres.c650 PetscBLASInt bn; in KSPDGMRESComputeSchurForm_DGMRES() local
671 PetscCall(PetscBLASIntCast(n, &bn)); in KSPDGMRESComputeSchurForm_DGMRES()
673 ihi = ldQ = bn; in KSPDGMRESComputeSchurForm_DGMRES()
697 PetscCall(PetscMalloc1(bn * bn, &Ht)); in KSPDGMRESComputeSchurForm_DGMRES()
698 for (i = 0; i < bn; i++) { in KSPDGMRESComputeSchurForm_DGMRES()
699 for (j = 0; j < bn; j++) Ht[i * bn + j] = dgmres->hes_origin[j * ldA + i]; in KSPDGMRESComputeSchurForm_DGMRES()
703 PetscCall(PetscCalloc1(bn, &t)); in KSPDGMRESComputeSchurForm_DGMRES()
704 t[bn - 1] = dgmres->hes_origin[(bn - 1) * ldA + bn]; /* Pick the last element H(m+1,m) */ in KSPDGMRESComputeSchurForm_DGMRES()
705 PetscCall(PetscMalloc1(bn, &ipiv)); in KSPDGMRESComputeSchurForm_DGMRES()
710 PetscCallBLAS("LAPACKgesv", LAPACKgesv_(&bn, &nrhs, Ht, &bn, ipiv, t, &bn, &info)); in KSPDGMRESComputeSchurForm_DGMRES()
[all …]
/petsc/src/snes/tutorials/
H A Dex18.c297 PetscScalar dn, ds, de, dw, an, as, ae, aw, bn, bs, be, bw, gn, gs, ge, gw; in FormJacobian() local
355 bn = PetscPowScalar(an, bm1); in FormJacobian()
358 gn = coef * bn * (tn - t0); in FormJacobian()
397 bn = PetscPowScalar(an, bm1); in FormJacobian()
400 gn = coef * bn * (tn - t0); in FormJacobian()
424 bn = PetscPowScalar(an, bm1); in FormJacobian()
427 gn = coef * bn * (tn - t0); in FormJacobian()
483 bn = PetscPowScalar(an, bm1); in FormJacobian()
486 gn = coef * bn * (tn - t0); in FormJacobian()
510 bn = PetscPowScalar(an, bm1); in FormJacobian()
[all …]
/petsc/src/ts/tutorials/
H A Dex11.h683 …static PetscScalar bn[3], fn, ft, tg[3], pl, rl, pm, pr, rr, xp, ubl, ubm, ubr, dtt, unm, tmp, utl… in godunovflux() local
691 bn[k - 1] = 0.; in godunovflux()
703 bn[0] = -nn[2] * tg[1]; in godunovflux()
704 bn[1] = nn[2] * tg[0]; in godunovflux()
705 bn[2] = nn[0] * tg[1] - nn[1] * tg[0]; in godunovflux()
707 d__1 = bn[0]; in godunovflux()
709 d__2 = bn[1]; in godunovflux()
711 d__3 = bn[2]; in godunovflux()
714 for (k = 1; k <= i__1; ++k) bn[k - 1] /= tmp; in godunovflux()
720 bn[0] = 0.; in godunovflux()
[all …]
/petsc/src/mat/impls/aij/seq/
H A Dmatmatmult.c148 PetscInt am = A->rmap->N, bn = B->cmap->N, bm = B->rmap->N; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() local
163 PetscCall(PetscHMapICreateWithSize(bn, &ta)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
168 PetscCall(PetscLLCondensedCreate(Crmax, bn, &lnk, &lnkbt)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
198 PetscCall(PetscLLCondensedClean(bn, cnzi, current_space->array, lnk, lnkbt)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
215 …PetscCall(MatSetSeqAIJWithArrays_private(PetscObjectComm((PetscObject)A), am, bn, ci, cj, NULL, ((… in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
384 PetscInt am = A->rmap->N, bn = B->cmap->N, bm = B->rmap->N; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast() local
398 PetscCall(PetscHMapICreateWithSize(bn, &ta)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast()
452 …PetscCall(MatSetSeqAIJWithArrays_private(PetscObjectComm((PetscObject)A), am, bn, ci, cj, ca, ((Pe… in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast()
487 PetscInt am = A->rmap->N, bn = B->cmap->N, bm = B->rmap->N; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable() local
501 PetscCall(PetscHMapICreateWithSize(bn, &ta)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable()
[all …]
H A Dmattransposematmult.c109 PetscInt bn = B->cmap->n; in MatTMatTMultNumeric_SeqAIJ_SeqDense() local
113 for (j = 0; j < bn; j++) { in MatTMatTMultNumeric_SeqAIJ_SeqDense()
/petsc/src/mat/graphops/order/
H A Dspectral.c133 PetscBLASInt bn, bN, lwork = 0, lierr, idummy; in MatGetOrdering_Spectral() local
139 PetscCall(PetscBLASIntCast(n, &bn)); in MatGetOrdering_Spectral()
145 …PetscCallBLAS("LAPACKgeev", LAPACKgeev_("N", "V", &bn, a, &bN, realpart, imagpart, &sdummy, &idumm… in MatGetOrdering_Spectral()
/petsc/src/ksp/ksp/interface/
H A Deige.c168 PetscBLASInt bn; in KSPComputeEigenvaluesExplicitly() local
170 PetscCall(PetscBLASIntCast(n, &bn)); in KSPComputeEigenvaluesExplicitly()
172 …PetscCallBLAS("LAPACKgeev", LAPACKgeev_("N", "N", &bn, array, &bn, realpart, imagpart, &sdummy, &i… in KSPComputeEigenvaluesExplicitly()
/petsc/src/snes/tests/
H A Dex20.c452 PetscScalar ***x, bn, bs, be, bw, bu, bd, gn, gs, ge, gw, gu, gd; in FormJacobian() local
513 bn = PetscPowScalar(an, bm1); in FormJacobian()
516 gn = coef * bn * (tn - t0); in FormJacobian()
582 bn = PetscPowScalar(an, bm1); in FormJacobian()
585 gn = coef * bn * (tn - t0); in FormJacobian()
792 bn = PetscPowScalar(an, bm1); in FormJacobian()
795 gn = coef * bn * (tn - t0); in FormJacobian()
922 bn = PetscPowScalar(an, bm1); in FormJacobian()
925 gn = coef * bn * (tn - t0); in FormJacobian()
1132 bn = PetscPowScalar(an, bm1); in FormJacobian()
[all …]
/petsc/src/vec/is/utils/
H A Disltog.c1527 PetscInt **bprocs = NULL, *bn_procs = NULL, bs, i, j, k, bn; in ISLocalToGlobalMappingGetNodeInfo() local
1532 PetscCall(ISLocalToGlobalMappingGetBlockNodeInfo(mapping, &bn, &bn_procs, &bprocs)); in ISLocalToGlobalMappingGetNodeInfo()
1537 PetscCall(PetscMalloc1(bn * bs, &tn_procs)); in ISLocalToGlobalMappingGetNodeInfo()
1538 for (i = 0, c = 0; i < bn; i++) { in ISLocalToGlobalMappingGetNodeInfo()
1542 if (n) *n = bn * bs; in ISLocalToGlobalMappingGetNodeInfo()
1545 PetscInt tn = bn * bs; in ISLocalToGlobalMappingGetNodeInfo()
1550 for (i = 0; i < bn; i++) { in ISLocalToGlobalMappingGetNodeInfo()
1560 if (n) *n = bn; in ISLocalToGlobalMappingGetNodeInfo()
/petsc/src/ksp/ksp/impls/minres/
H A Dminres.c878 PetscBLASInt bn, lierr = 0, ldz = 1; in KSPComputeEigenvalues_MINRES() local
896 PetscCall(PetscBLASIntCast(n, &bn)); in KSPComputeEigenvalues_MINRES()
898 PetscCallBLAS("LAPACKREALstev", LAPACKREALstev_("N", &bn, r, ee, NULL, &ldz, NULL, &lierr)); in KSPComputeEigenvalues_MINRES()
911 PetscBLASInt bn, lierr = 0, ldz = 1; in KSPComputeExtremeSingularValues_MINRES() local
929 PetscCall(PetscBLASIntCast(n, &bn)); in KSPComputeExtremeSingularValues_MINRES()
931 PetscCallBLAS("LAPACKREALstev", LAPACKREALstev_("N", &bn, dd, ee, NULL, &ldz, NULL, &lierr)); in KSPComputeExtremeSingularValues_MINRES()
/petsc/src/ksp/pc/tutorials/
H A Dex4.c166 …PetscInt boundary_indices_size, am, an, bm, bn, condensed_am, astart, aend, Dstart, Dend,… in main() local
237 PetscCall(MatGetLocalSize(Bcondensed, &bm, &bn)); in main()
/petsc/src/dm/dt/interface/
H A Ddt.c376 PetscBLASInt bm, bn, info; in PetscDTJacobianInverse_Internal() local
381 PetscCall(PetscBLASIntCast(n, &bn)); in PetscDTJacobianInverse_Internal()
447 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&bn, &bn, JTJ, &bn, pivots, &info)); in PetscDTJacobianInverse_Internal()
449 PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&bn, JTJ, &bn, pivots, W, &bn, &info)); in PetscDTJacobianInverse_Internal()
1427 PetscBLASInt bn, bm, ldz; /* bm will equal bn on exit */ in PetscDTSymmetricTridiagonalEigensolve() local
1440 PetscCall(PetscBLASIntCast(n, &bn)); in PetscDTSymmetricTridiagonalEigensolve()
1446 …PetscCallBLAS("LAPACKstegr", LAPACKstegr_(&jobz, &range, &bn, diag, subdiag, &VL, &VU, &IL, &IU, &… in PetscDTSymmetricTridiagonalEigensolve()
1452 …PetscCallBLAS("LAPACKstegr", LAPACKstegr_(&jobz, &range, &bn, diag, subdiag, &VL, &VU, &IL, &IU, &… in PetscDTSymmetricTridiagonalEigensolve()
1462 PetscCallBLAS("LAPACKsteqr", LAPACKsteqr_("I", &bn, diag, subdiag, V, &ldz, work, &info)); in PetscDTSymmetricTridiagonalEigensolve()
/petsc/config/
H A Dtestparse.py397 bn=re.sub("new_","",basename)
399 testname="run"+getlangsplit(bn)
/petsc/src/mat/impls/dense/mpi/
H A Dmpidense.c2270 PetscInt i, cK = A->cmap->N, sendsiz, recvsiz, k, j, bn; in MatMatTransposeMultNumeric_MPIDense_MPIDense_Cyclic() local
2294 bn = B->rmap->n; in MatMatTransposeMultNumeric_MPIDense_MPIDense_Cyclic()
2295 if (blda == bn) { in MatMatTransposeMultNumeric_MPIDense_MPIDense_Cyclic()
2300 for (j = 0; j < bn; j++, k++) sendbuf[k] = bv[i * blda + j]; in MatMatTransposeMultNumeric_MPIDense_MPIDense_Cyclic()
2319 sendsiz = cK * bn; in MatMatTransposeMultNumeric_MPIDense_MPIDense_Cyclic()
2334 bn = nextbn; in MatMatTransposeMultNumeric_MPIDense_MPIDense_Cyclic()
2355 PetscInt blda, i, cK = A->cmap->N, k, j, bn; in MatMatTransposeMultNumeric_MPIDense_MPIDense_Allgatherv() local
2374 bn = B->rmap->n; in MatMatTransposeMultNumeric_MPIDense_MPIDense_Allgatherv()
2377 for (k = 0, j = 0; j < bn; j++) { in MatMatTransposeMultNumeric_MPIDense_MPIDense_Allgatherv()
2381 PetscCall(PetscMPIIntCast(bn * cK, &ibn)); in MatMatTransposeMultNumeric_MPIDense_MPIDense_Allgatherv()
/petsc/src/dm/impls/plex/
H A Dplexgmsh.c2032 for (int bn = 0; bn < elem->numNodes + s; ++bn) { in DMPlexCreateGmsh() local
2033 if (lexorder[bn] < 0) continue; in DMPlexCreateGmsh()
2035 const PetscInt bnode = elem->nodes[lexorder[bn]]; in DMPlexCreateGmsh()
2036 …for (d = 0; d < coordDim; ++d) cellCoords[n * coordDim + d] += weights[bn] * (PetscReal)coords[bno… in DMPlexCreateGmsh()
/petsc/src/mat/impls/aij/mpi/
H A Dmpiov.c284 PetscInt is_id, tnz, an, bn, rstart, cstart, row, start, end, col, totalrows, *sbdata; in MatIncreaseOverlap_MPIAIJ_Send_Scalable() local
299 PetscCall(MatGetRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
302 tnz = ai[an] + bi[bn]; in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
398 PetscCall(MatRestoreRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
405 PetscInt tnz, an, bn, i, j, row, start, end, rstart, cstart, col, k, *indices_temp; in MatIncreaseOverlap_MPIAIJ_Local_Scalable() local
419 PetscCall(MatGetRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Local_Scalable()
422 tnz = ai[an] + bi[bn]; in MatIncreaseOverlap_MPIAIJ_Local_Scalable()
465 PetscCall(MatRestoreRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Local_Scalable()

12