Home
last modified time | relevance | path

Searched refs:pivots (Results 1 – 22 of 22) sorted by relevance

/petsc/src/ksp/ksp/tutorials/output/
H A Dex52_5.out3 Mumps num of null pivots detected = 6
4 Mumps row with null pivots is = 0
5 Mumps row with null pivots is = 833
6 Mumps row with null pivots is = 834
7 Mumps row with null pivots is = 1666
8 Mumps row with null pivots is = 1667
9 Mumps row with null pivots is = 2499
H A Dex50_tut_2.out26 Replace tiny pivots FALSE
/petsc/src/vec/is/utils/
H A Dpsort.c96 PetscInt *pivots, *finalpivots, i; in PetscParallelSampleSelect() local
105 PetscCall(PetscMalloc1(size - 1, &pivots)); in PetscParallelSampleSelect()
113 pivots[i] = PETSC_INT_MAX; in PetscParallelSampleSelect()
120 pivots[i] = keysin[index]; in PetscParallelSampleSelect()
124 PetscCall(PetscParallelSortInt_Bitonic(mapin->comm, size - 1, pivots)); in PetscParallelSampleSelect()
128 PetscCall(PetscParallelSortedInt(mapin->comm, size - 1, pivots, &sorted)); in PetscParallelSampleSelect()
152 for (i = 0; i < keys_per[rank]; i++) pivots[i] = pivots[my_first + i * non_empty]; in PetscParallelSampleSelect()
153 for (i = keys_per[rank]; i < max_keys_per; i++) pivots[i] = PETSC_INT_MAX; in PetscParallelSampleSelect()
154 …PetscCallMPI(MPI_Allgather(pivots, max_keys_per, MPIU_INT, finalpivots, max_keys_per, MPIU_INT, ma… in PetscParallelSampleSelect()
164 PetscCall(PetscFree(pivots)); in PetscParallelSampleSelect()
[all …]
/petsc/src/dm/dt/fe/impls/composite/
H A Dfecomposite.c23 PetscBLASInt *pivots; in PetscFESetUp_Composite() local
64 PetscCall(PetscMalloc2(spdim, &pivots, spdim, &work)); in PetscFESetUp_Composite()
89 …PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&n, &n, &invVscalar[s * spdim * spdim], &n, pivots, &inf… in PetscFESetUp_Composite()
90 …PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&n, &invVscalar[s * spdim * spdim], &n, pivots, work, &n… in PetscFESetUp_Composite()
96 PetscCall(PetscFree2(pivots, work)); in PetscFESetUp_Composite()
/petsc/src/ksp/ksp/tests/output/
H A Dex81_1.out31 INFOG(25) (after factorization: number of pivots modified by static pivoting): 0
32 INFOG(28) (after factorization: number of null pivots encountered): 0
80 INFOG(25) (after factorization: number of pivots modified by static pivoting): 0
81 INFOG(28) (after factorization: number of null pivots encountered): 0
H A Dex81_2.out43 INFOG(25) (after factorization: number of pivots modified by static pivoting): 0
44 INFOG(28) (after factorization: number of null pivots encountered): 0
111 INFOG(25) (after factorization: number of pivots modified by static pivoting): 0
112 INFOG(28) (after factorization: number of null pivots encountered): 0
H A Dex63_2.out43 Replace tiny pivots FALSE
H A Dex63_2_alt.out43 Replace tiny pivots FALSE
/petsc/include/petsc/private/kernels/
H A Dblockinvert.h44pivots, W, allowzeropivot, zeropivotdetected) ((PetscErrorCode)(PetscLINPACKgefa((A), (bs), (pivot… argument
/petsc/src/snes/tutorials/output/
H A Dex19_superlu_dist_3d.out48 Replace tiny pivots FALSE
H A Dex19_superlu_dist_3ds.out48 Replace tiny pivots FALSE
/petsc/src/mat/impls/dense/seq/
H A Ddense.c44 PetscCheck(mat->pivots, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Pivots not present"); in MatSeqDenseInvertFactors_Private()
50 …PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&n, mat->v, &mat->lda, mat->pivots, mat->fwork, &mat->lf… in MatSeqDenseInvertFactors_Private()
61 PetscCheck(mat->pivots, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Pivots not present"); in MatSeqDenseInvertFactors_Private()
64 …PetscCallBLAS("LAPACKhetri", LAPACKhetri_("L", &n, mat->v, &mat->lda, mat->pivots, mat->fwork, &in… in MatSeqDenseInvertFactors_Private()
69 PetscCheck(mat->pivots, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Pivots not present"); in MatSeqDenseInvertFactors_Private()
72 …PetscCallBLAS("LAPACKsytri", LAPACKsytri_("L", &n, mat->v, &mat->lda, mat->pivots, mat->fwork, &in… in MatSeqDenseInvertFactors_Private()
409 …APACKgetrs", LAPACKgetrs_(T ? "T" : "N", &m, &nrhs, mat->v, &mat->lda, mat->pivots, x, &m, &info)); in MatSolve_SeqDense_Internal_LU()
433 …PetscCallBLAS("LAPACKhetrs", LAPACKhetrs_("L", &m, &nrhs, mat->v, &mat->lda, mat->pivots, x, &m, &… in MatSolve_SeqDense_Internal_Cholesky()
440 …PetscCallBLAS("LAPACKsytrs", LAPACKsytrs_("L", &m, &nrhs, mat->v, &mat->lda, mat->pivots, x, &m, &… in MatSolve_SeqDense_Internal_Cholesky()
797 if (!mat->pivots) PetscCall(PetscMalloc1(A->rmap->n, &mat->pivots)); in MatLUFactor_SeqDense()
[all …]
H A Ddense.h18 PetscBLASInt *pivots; /* pivots in LU factorization */ member
/petsc/src/dm/dt/interface/
H A Ddt.c390 PetscBLASInt *pivots; in PetscDTJacobianInverse_Internal() local
393 PetscCall(PetscMalloc2(m, &pivots, m, &W)); in PetscDTJacobianInverse_Internal()
396 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&bm, &bm, Jinvs, &bm, pivots, &info)); in PetscDTJacobianInverse_Internal()
398 PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&bm, Jinvs, &bm, pivots, W, &bm, &info)); in PetscDTJacobianInverse_Internal()
400 PetscCall(PetscFree2(pivots, W)); in PetscDTJacobianInverse_Internal()
403 PetscBLASInt *pivots; in PetscDTJacobianInverse_Internal() local
407 PetscCall(PetscMalloc2(m, &pivots, m, &W)); in PetscDTJacobianInverse_Internal()
417 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&bm, &bm, JJT, &bm, pivots, &info)); in PetscDTJacobianInverse_Internal()
419 PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&bm, JJT, &bm, pivots, W, &bm, &info)); in PetscDTJacobianInverse_Internal()
429 PetscCall(PetscFree2(pivots, W)); in PetscDTJacobianInverse_Internal()
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaijfact.c569 PetscInt *pivots; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local
582 PetscCall(PetscMalloc1(bs, &pivots)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
689 …PetscCall(PetscKernel_A_gets_inverse_A(bs, diag, pivots, work, allowzeropivot, &zeropivotdetected)… in MatCholeskyFactorNumeric_SeqSBAIJ_N()
716 PetscCall(PetscFree(pivots)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
742 PetscInt *pivots; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local
752 PetscCall(PetscMalloc1(bs, &pivots)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
821 …PetscCall(PetscKernel_A_gets_inverse_A(bs, diag, pivots, work, allowzeropivot, &zeropivotdetected)… in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
848 PetscCall(PetscFree(pivots)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
/petsc/src/ksp/pc/impls/bddc/
H A Dbddcschurs.c396 PetscBLASInt B_N, B_ierr, B_lwork, *pivots; in PCBDDCSubSchursSetUp() local
639 PetscCall(PetscMalloc2(B_lwork, &Bwork, B_N, &pivots)); in PCBDDCSubSchursSetUp()
642 pivots = NULL; in PCBDDCSubSchursSetUp()
1724 …PetscCallBLAS("LAPACKsytrf", LAPACKsytrf_("L", &B_N, S_data, &B_N, pivots, Bwork, &B_lwork, &B_ier… in PCBDDCSubSchursSetUp()
1726 … PetscCallBLAS("LAPACKsytri", LAPACKsytri_("L", &B_N, S_data, &B_N, pivots, Bwork, &B_ierr)); in PCBDDCSubSchursSetUp()
1729 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&B_N, &B_N, S_data, &B_N, pivots, &B_ierr)); in PCBDDCSubSchursSetUp()
1731 … PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&B_N, S_data, &B_N, pivots, Bwork, &B_lwork, &B_ierr)); in PCBDDCSubSchursSetUp()
2008 …PetscCallBLAS("LAPACKsytrf", LAPACKsytrf_("L", &B_N, array + cum, &B_N, pivots, Bwork, &B_lwork, &… in PCBDDCSubSchursSetUp()
2010 … PetscCallBLAS("LAPACKsytri", LAPACKsytri_("L", &B_N, array + cum, &B_N, pivots, Bwork, &B_ierr)); in PCBDDCSubSchursSetUp()
2013 … PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&B_N, &B_N, array + cum, &B_N, pivots, &B_ierr)); in PCBDDCSubSchursSetUp()
[all …]
H A Dbddcfetidp.c498 PetscBLASInt *pivots, B_lwork, B_N, B_ierr; in PCBDDCSetupFETIDPMatContext() local
533 PetscCall(PetscMalloc3(mss * mss, &W, mss, &pivots, B_lwork, &Bwork)); in PCBDDCSetupFETIDPMatContext()
545 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&B_N, &B_N, W, &B_N, pivots, &B_ierr)); in PCBDDCSetupFETIDPMatContext()
547 PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&B_N, W, &B_N, pivots, Bwork, &B_lwork, &B_ierr)); in PCBDDCSetupFETIDPMatContext()
560 PetscCall(PetscFree3(W, pivots, Bwork)); in PCBDDCSetupFETIDPMatContext()
H A Dbddcscalingbasic.c37 …PetscCallBLAS("LAPACKgetrs", LAPACKgetrs_("T", &m, &nrhs, mat->v, &mat->lda, mat->pivots, x, &m, &… in PCBDDCMatTransposeMatSolve_SeqDense()
/petsc/include/petsc/private/
H A Dpetscscalapack.h23 PetscBLASInt *pivots; /* pivots in LU factorization */ member
/petsc/src/mat/impls/scalapack/
H A Dmatscalapack.c741 …trs", SCALAPACKgetrs_("N", &a->M, &nrhs, a->loc, &one, &one, a->desc, a->pivots, x2d, &one, &one, … in MatSolve_ScaLAPACK()
791 …trs", SCALAPACKgetrs_("N", &a->M, &x->N, a->loc, &one, &one, a->desc, a->pivots, x->loc, &one, &on… in MatMatSolve_ScaLAPACK()
815 if (!a->pivots) PetscCall(PetscMalloc1(a->locr + a->mb, &a->pivots)); in MatLUFactor_ScaLAPACK()
816 …AS("SCALAPACKgetrf", SCALAPACKgetrf_(&a->M, &a->N, a->loc, &one, &one, a->desc, a->pivots, &info)); in MatLUFactor_ScaLAPACK()
1256 PetscCall(PetscFree(a->pivots)); in MatDestroy_ScaLAPACK()
/petsc/src/dm/dt/fe/impls/basic/
H A Dfebasic.c49 PetscBLASInt *pivots; in PetscFESetUp_Basic() local
77 PetscCall(PetscMalloc2(pdim, &pivots, pdim, &work)); in PetscFESetUp_Basic()
79 PetscCallBLAS("LAPACKgetrf", LAPACKREALgetrf_(&n, &n, fem->invV, &n, pivots, &info)); in PetscFESetUp_Basic()
81 PetscCallBLAS("LAPACKgetri", LAPACKREALgetri_(&n, fem->invV, &n, pivots, work, &n, &info)); in PetscFESetUp_Basic()
83 PetscCall(PetscFree2(pivots, work)); in PetscFESetUp_Basic()
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A DPC.pyx1408 The size at which smaller pivots are treated as zero.