| /petsc/src/ksp/ksp/utils/lmvm/rescale/ |
| H A D | symbrdnrescale.c | 8 static PetscErrorCode SymBroydenRescaleUpdateScalar(Mat B, SymBroydenRescale ldb) in SymBroydenRescaleUpdateScalar() argument
16 next = ldb->k; in SymBroydenRescaleUpdateScalar()
17 oldest = PetscMax(0, ldb->k - ldb->sigma_hist); in SymBroydenRescaleUpdateScalar()
20 if (ldb->sigma_hist == 0) { in SymBroydenRescaleUpdateScalar()
24 if (ldb->alpha == 1.0) { in SymBroydenRescaleUpdateScalar()
25 for (PetscInt i = 0; i < next - oldest; ++i) signew += ldb->yts[i] / ldb->yty[i]; in SymBroydenRescaleUpdateScalar()
26 } else if (ldb->alpha == 0.5) { in SymBroydenRescaleUpdateScalar()
27 for (PetscInt i = 0; i < next - oldest; ++i) signew += ldb->sts[i] / ldb->yty[i]; in SymBroydenRescaleUpdateScalar()
29 } else if (ldb->alpha == 0.0) { in SymBroydenRescaleUpdateScalar()
30 for (PetscInt i = 0; i < next - oldest; ++i) signew += ldb->sts[i] / ldb->yts[i]; in SymBroydenRescaleUpdateScalar()
[all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/diagbrdn/ |
| H A D | diagbrdn.c | 29 SymBroydenRescale ldb = (SymBroydenRescale)lmvm->ctx; in MatUpdate_DiagBrdn() local
50 PetscCall(SymBroydenRescaleUpdate(B, ldb)); in MatUpdate_DiagBrdn()
53 PetscCall(SymBroydenRescaleReset(B, ldb, MAT_LMVM_RESET_HISTORY)); in MatUpdate_DiagBrdn()
79 SymBroydenRescale ldb = (SymBroydenRescale)lmvm->ctx; in MatView_DiagBrdn() local
83 PetscCall(SymBroydenRescaleView(ldb, pv)); in MatView_DiagBrdn()
90 SymBroydenRescale ldb = (SymBroydenRescale)lmvm->ctx; in MatSetFromOptions_DiagBrdn() local
94 PetscCall(SymBroydenRescaleSetFromOptions(B, ldb, PetscOptionsObject)); in MatSetFromOptions_DiagBrdn()
101 SymBroydenRescale ldb = (SymBroydenRescale)lmvm->ctx; in MatReset_DiagBrdn() local
104 PetscCall(SymBroydenRescaleReset(B, ldb, mode)); in MatReset_DiagBrdn()
121 SymBroydenRescale ldb = (SymBroydenRescale)lmvm->ctx; in MatSetUp_DiagBrdn() local
[all …]
|
| /petsc/include/petsc/private/ |
| H A D | dmpleximpl.h | 571 static inline void DMPlex_MatMult2D_Internal(const PetscScalar A[], PetscInt n, PetscInt ldb, const… in DMPlex_MatMult2D_Internal() argument
576 for (int d = 0; d < PLEX_DIM__; ++d) z[d] = B[d * ldb + j]; in DMPlex_MatMult2D_Internal()
578 for (int d = 0; d < PLEX_DIM__; ++d) C[d * ldb + j] = z[d]; in DMPlex_MatMult2D_Internal()
583 static inline void DMPlex_MatMult3D_Internal(const PetscScalar A[], PetscInt n, PetscInt ldb, const… in DMPlex_MatMult3D_Internal() argument
588 for (int d = 0; d < PLEX_DIM__; ++d) z[d] = B[d * ldb + j]; in DMPlex_MatMult3D_Internal()
590 for (int d = 0; d < PLEX_DIM__; ++d) C[d * ldb + j] = z[d]; in DMPlex_MatMult3D_Internal()
595 …ex_MatMultTranspose2D_Internal(const PetscScalar A[], PetscInt n, PetscInt ldb, const PetscScalar … in DMPlex_MatMultTranspose2D_Internal() argument
600 for (int d = 0; d < PLEX_DIM__; ++d) z[d] = B[d * ldb + j]; in DMPlex_MatMultTranspose2D_Internal()
602 for (int d = 0; d < PLEX_DIM__; ++d) C[d * ldb + j] = z[d]; in DMPlex_MatMultTranspose2D_Internal()
607 …ex_MatMultTranspose3D_Internal(const PetscScalar A[], PetscInt n, PetscInt ldb, const PetscScalar … in DMPlex_MatMultTranspose3D_Internal() argument
[all …]
|
| H A D | cupmsolverinterface.hpp | 109 …, const cupmScalar_t *A, cupmBlasInt_t lda, cupmScalar_t *B, cupmBlasInt_t ldb, cupmScalar_t * /* … in PETSC_CUPMSOLVER_ALIAS_BLAS_FUNCTION()
111 return cupmSolverXpotrs_p(handle, uplo, n, nrhs, A, lda, B, ldb, dev_info); in PETSC_CUPMSOLVER_ALIAS_BLAS_FUNCTION()
152 …cupmBlasInt_t lda, cupmBlasInt_t *dev_ipiv, cupmScalar_t *B, cupmBlasInt_t ldb, cupmScalar_t * /* … in PETSC_CUPMSOLVER_ALIAS_BLAS_FUNCTION()
154 return cupmSolverXgetrs_p(handle, op, n, nrhs, A, lda, dev_ipiv, B, ldb, dev_info); in PETSC_CUPMSOLVER_ALIAS_BLAS_FUNCTION()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpimattransposematmult.c | 60 PetscInt i, j, m = A->rmap->n, n = A->cmap->n, ldb, BN = B->cmap->N, ldc; in MatTransposeMatMultNumeric_MPIAIJ_MPIDense() local
78 PetscCall(MatDenseGetLDA(B, &ldb)); in MatTransposeMatMultNumeric_MPIAIJ_MPIDense()
81 for (i = 0; i < m; i++) btarray[i * BN + j] = Barray[j * ldb + i]; in MatTransposeMatMultNumeric_MPIAIJ_MPIDense()
|
| /petsc/src/snes/impls/ngmres/ |
| H A D | snesngmres.h | 48 PetscBLASInt ldb; /* the padded vector dimension */ member
|
| H A D | ngmresfunc.c | 85 …&ngmres->n, &ngmres->nrhs, ngmres->h, &ngmres->lda, ngmres->beta, &ngmres->ldb, ngmres->s, &ngmres… in SNESNGMRESFormCombinedSolution_Private()
87 …&ngmres->n, &ngmres->nrhs, ngmres->h, &ngmres->lda, ngmres->beta, &ngmres->ldb, ngmres->s, &ngmres… in SNESNGMRESFormCombinedSolution_Private()
|
| H A D | snesngmres.c | 66 PetscCall(PetscBLASIntCast(msize, &ngmres->ldb)); in SNESSetUp_NGMRES()
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex13.c | 256 PetscBLASInt m, n, k, lda, ldb, ldc; in test() local
264 PetscCall(PetscBLASIntCast(Nf * Nbpt, &ldb)); in test()
266 …_("N", "T", &m, &n, &k, &alpha, p_scalar, &lda, &M_moment_real[f * Nbpt], &ldb, &beta, &p_trimmed_… in test()
|
| /petsc/src/dm/dt/fv/interface/ |
| H A D | fv.c | 2011 PetscBLASInt M, N, K, lda, ldb, ldwork, info; in PetscFVLeastSquaresPseudoInverse_Static() local
2039 ldb = lda; in PetscFVLeastSquaresPseudoInverse_Static()
2040 …sm_("Right", "Upper", "ConjugateTranspose", "NotUnitTriangular", &M, &N, &Alpha, R, &lda, Q, &ldb); in PetscFVLeastSquaresPseudoInverse_Static()
2048 …BLASgemm_("ConjugateTranspose", "Normal", &N, &K, &M, &Alpha, Ainv, &lda, Aback, &ldb, &Beta, work… in PetscFVLeastSquaresPseudoInverse_Static()
2066 PetscBLASInt M, N, lda, ldb, ldwork; in PetscFVLeastSquaresPseudoInverseSVD_Static() local
2081 PetscCall(PetscBLASIntCast(maxmn, &ldb)); in PetscFVLeastSquaresPseudoInverseSVD_Static()
2090 …PetscCallBLAS("LAPACKgelss", LAPACKgelss_(&M, &N, &nrhs, A, &lda, Brhs, &ldb, rtau, &rcond, &irank… in PetscFVLeastSquaresPseudoInverseSVD_Static()
2098 …PetscCallBLAS("LAPACKgelss", LAPACKgelss_(&M, &N, &nrhs, A, &lda, Brhs, &ldb, tau, &rcond, &irank,… in PetscFVLeastSquaresPseudoInverseSVD_Static()
|
| /petsc/src/mat/tests/ |
| H A D | ex70.c | 206 PetscInt m, n, k, M = 10, N = 10, K = 5, ldx = 3, ldb = 5, ldr = 4; in main() local
224 PetscCall(PetscOptionsGetInt(NULL, NULL, "-ldb", &ldb, NULL)); in main()
287 PetscCall(PetscMalloc1((n + ldb) * K, &dataB)); in main()
289 for (i = 0; i < (n + ldb) * K; i++) dataB[i] = MAGIC_NUMBER; in main()
295 PetscCall(MatDenseSetLDA(B, n + ldb)); in main()
|
| /petsc/src/mat/impls/diagonal/ |
| H A D | diagonal.c | 618 PetscInt ldb, ldc; in MatProductNumeric_Diagonal_Dense() local
626 PetscCall(MatDenseGetLDA(B, &ldb)); in MatProductNumeric_Diagonal_Dense()
631 for (PetscInt i = 0; i < B->rmap->n; i++) c[i + j * ldc] = alpha[i] * b[i + j * ldb]; in MatProductNumeric_Diagonal_Dense()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 802 PetscInt nz, neq, ldb, ldx; in MatMatSolve_SeqAIJ_inplace() local
820 PetscCall(MatDenseGetLDA(B, &ldb)); in MatMatSolve_SeqAIJ_inplace()
848 b += ldb; in MatMatSolve_SeqAIJ_inplace()
866 PetscInt nz, neq, ldb, ldx; in MatMatSolve_SeqAIJ() local
883 PetscCall(MatDenseGetLDA(B, &ldb)); in MatMatSolve_SeqAIJ()
912 b += ldb; in MatMatSolve_SeqAIJ()
930 PetscInt nz, neq, ldb, ldx; in MatMatSolveTranspose_SeqAIJ() local
947 PetscCall(MatDenseGetLDA(B, &ldb)); in MatMatSolveTranspose_SeqAIJ()
980 b += ldb; in MatMatSolveTranspose_SeqAIJ()
|
| /petsc/src/ts/impls/implicit/glle/ |
| H A D | glle.c | 131 PetscBLASInt m, n, one = 1, *ipiv, lwork, info, ldb; in TSGLLESchemeCreate() local
236 PetscCall(PetscBLASIntCast(ss, &ldb)); in TSGLLESchemeCreate()
240 …PetscCallBLAS("LAPACKgelss", LAPACKgelss_(&m, &n, &m, H, &m, bmat, &ldb, sing, &rcond, &rank, work… in TSGLLESchemeCreate()
243 …PetscCallBLAS("LAPACKgelss", LAPACKgelss_(&m, &n, &m, H, &m, bmat, &ldb, sing, &rcond, &rank, work… in TSGLLESchemeCreate()
|
| /petsc/src/snes/impls/composite/ |
| H A D | snescomposite.c | 33 PetscBLASInt ldb; /* the padded vector dimension */ member
359 PetscCall(PetscBLASIntCast(jac->nsnes, &jac->ldb)); in SNESSetUp_Composite()
|
| /petsc/src/mat/impls/dense/seq/ |
| H A D | dense.c | 643 PetscBLASInt nrhs = 0, m = 0, k = 0, ldb = 0, ldx = 0, ldy = 0; in MatMatSolve_SeqDense_SetUp() local
656 PetscCall(PetscBLASIntCast(_ldb, &ldb)); in MatMatSolve_SeqDense_SetUp()
662 if (ldb == m) { in MatMatSolve_SeqDense_SetUp()
663 PetscCall(PetscArraycpy(y, b, ldb * nrhs)); in MatMatSolve_SeqDense_SetUp()
665 for (PetscInt j = 0; j < nrhs; j++) PetscCall(PetscArraycpy(&y[j * m], &b[j * ldb], m)); in MatMatSolve_SeqDense_SetUp()
670 if (ldb == ldx) { in MatMatSolve_SeqDense_SetUp()
676 for (PetscInt j = 0; j < nrhs; j++) PetscCall(PetscArraycpy(&y[j * ldx], &b[j * ldb], m)); in MatMatSolve_SeqDense_SetUp()
2555 PetscInt i, j, nrows, ncols, ldb; in MatCreateSubMatrix_SeqDense() local
2585 PetscCall(MatDenseGetLDA(newmat, &ldb)); in MatCreateSubMatrix_SeqDense()
2589 bv += ldb; in MatCreateSubMatrix_SeqDense()
|
| /petsc/src/mat/impls/scalapack/ |
| H A D | matscalapack.c | 978 PetscInt i, M = A->rmap->N, N = A->cmap->N, m, n, rstart, rend, nz, ldb; in MatConvert_ScaLAPACK_Dense() local
1012 PetscCall(MatDenseGetLDA(Bmpi, &ldb)); in MatConvert_ScaLAPACK_Dense()
1013 PetscCall(PetscBLASIntCast(PetscMax(ldb, 1), &lld)); /* local leading dimension */ in MatConvert_ScaLAPACK_Dense()
1051 PetscCall(MatDenseGetLDA(Bmpi, &ldb)); in MatConvert_ScaLAPACK_Dense()
1052 PetscCall(PetscBLASIntCast(PetscMax(ldb, 1), &lld)); /* local leading dimension */ in MatConvert_ScaLAPACK_Dense()
|
| /petsc/src/ksp/pc/impls/bjacobi/ |
| H A D | bjacobi.c | 1227 PetscInt m, N, lda, ldb; in PCMatApply_BJacobi_Multiproc() local
1234 PetscCall(MatDenseGetLDA(Y, &ldb)); in PCMatApply_BJacobi_Multiproc()
1240 PetscCall(MatDenseSetLDA(sY, ldb)); in PCMatApply_BJacobi_Multiproc()
|
| /petsc/src/mat/impls/nest/ |
| H A D | matnest.c | 100 PetscInt i, j, M, N, nr, nc, ldb, ldc; in MatProductNumeric_Nest_Dense() local
118 PetscCall(MatDenseGetLDA(B, &ldb)); in MatProductNumeric_Nest_Dense()
131 PetscCall(MatDenseSetLDA(viewB, ldb)); in MatProductNumeric_Nest_Dense()
173 PetscInt i, j, M, N, m, n, nr, nc, maxm = 0, ldb; in MatProductSymbolic_Nest_Dense() local
220 PetscCall(MatDenseGetLDA(B, &ldb)); in MatProductSymbolic_Nest_Dense()
227 PetscCall(MatDenseSetLDA(viewB, ldb)); in MatProductSymbolic_Nest_Dense()
|
| /petsc/src/mat/impls/dense/seq/cupm/ |
| H A D | matseqdensecupm.hpp | 723 cupmBlasInt_t nrhs, ldb, ldx, ldy; in MatMatSolve_Factored_Dispatch_() local
741 PetscCall(PetscCUPMBlasIntCast(n, &ldb)); in MatMatSolve_Factored_Dispatch_()
774 PetscCall(PetscCUPMMemcpy2DAsync(y, ldy, b, ldb, m, nrhs, copy_mode, stream)); in MatMatSolve_Factored_Dispatch_()
|
| /petsc/src/dm/dt/fe/impls/basic/ |
| H A D | febasic.c | 105 PetscBLASInt n_, p_, k_, lda, ldb, ldc; in TensorContract_Private() local
114 ldb = n_; in TensorContract_Private()
116 …mm", BLASREALgemm_("N", "T", &p_, &n_, &k_, &one, A + i * k * p, &lda, B, &ldb, &zero, C + i * n *… in TensorContract_Private()
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 2790 PetscBLASInt M, N, K, lda, ldb, ldwork, info; in PetscDTPseudoInverseQR() local
2827 ldb = lda; in PetscDTPseudoInverseQR()
2828 …m_("Right", "Upper", "ConjugateTranspose", "NotUnitTriangular", &M, &N, &Alpha, R, &lda, Q, &ldb)); in PetscDTPseudoInverseQR()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact2.c | 2074 PetscInt nz, i, j, neq, ldb, ldx; in MatMatSolve_SeqSBAIJ_1_NaturalOrdering() local
2086 PetscCall(MatDenseGetLDA(B, &ldb)); in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2109 b += ldb; in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
|