Searched refs:M11 (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbrdn.c | 36 LMProducts M11 = lsb->products[LMVMModeMap(SYMBROYDEN_PRODUCTS_M11, mode)]; in SymBroydenKernel_Recursive_Inner() local
51 PetscCall(LMProductsMult(M11, oldest, next, 1.0, YtX, 1.0, V, PETSC_FALSE)); in SymBroydenKernel_Recursive_Inner()
252 …[], PetscBLASInt lda00, PetscScalar M01[], PetscBLASInt lda01, PetscScalar M11[], PetscBLASInt lda… in SymBroydenCompactDenseUpdateArrays() argument
270 PetscCall(PetscArrayzero(M11, m * lda11)); in SymBroydenCompactDenseUpdateArrays()
288 M11[i + i * lda11] = delta; in SymBroydenCompactDenseUpdateArrays()
302 …PetscCallBLAS("BLASgemv", BLASgemv_("N", &m, &m, &sone, M11, &lda11, q1, &ione, &sone, p1, &ione)); in SymBroydenCompactDenseUpdateArrays()
331 PetscCallBLAS("LAPACKgerc", LAPACKgerc_(&m, &m, &alpha, p1, &ione, p1, &ione, M11, &lda11)); in SymBroydenCompactDenseUpdateArrays()
332 for (PetscInt j = 0; j < m; j++) M11[j + i * lda11] = beta * p1[j]; in SymBroydenCompactDenseUpdateArrays()
333 for (PetscInt j = 0; j < m; j++) M11[i + j * lda11] = PetscConj(beta * p1[j]); in SymBroydenCompactDenseUpdateArrays()
334 M11[i + i * lda11] = delta; in SymBroydenCompactDenseUpdateArrays()
[all …]
|