| /petsc/src/ksp/ksp/impls/gmres/agmres/ |
| H A D | agmresdeflation.c | 99 PetscScalar *wi = agmres->wi; in KSPAGMRESSchurForm() local
127 …APACKhgeqz_("S", "I", "I", &KspSize, &ilo, &KspSize, A, &ldA, B, &ldB, wr, wi, beta, Q, &N, Z, &N,… in KSPAGMRESSchurForm()
130 …", LAPACKgges_("V", "V", "N", NULL, &KspSize, A, &ldA, B, &ldB, &sdim, wr, wi, beta, Q, &N, Z, &N,… in KSPAGMRESSchurForm()
138 wi[i] /= beta[i]; in KSPAGMRESSchurForm()
143 PetscCall(KSPAGMRESQuickSort(wr, wi, KspSize, perm)); in KSPAGMRESSchurForm()
148 if (wi[r] != 0) r += 2; in KSPAGMRESSchurForm()
158 …APACKtgsen_(&ijob, &wantQ, &wantZ, select, &KspSize, A, &ldA, B, &ldB, wr, wi, beta, Q, &N, Z, &N,… in KSPAGMRESSchurForm()
170 …if (agmres->DeflPrecond) PetscCall(KSPAGMRESLejaOrdering(wr, wi, agmres->Rshift, agmres->Ishift, m… in KSPAGMRESSchurForm()
|
| H A D | agmres.c | 133 PetscCall(KSPAGMRESLejaOrdering(agmres->wr, agmres->wi, agmres->Rshift, agmres->Ishift, max_k)); in KSPComputeShifts_DGMRES()
136 PetscScalar *wr, *wi, *Rshift, *Ishift; in KSPComputeShifts_DGMRES() local
137 PetscCall(PetscMalloc4(2 * max_k, &wr, 2 * max_k, &wi, 2 * max_k, &Rshift, 2 * max_k, &Ishift)); in KSPComputeShifts_DGMRES()
140 wi[i] = agmres->wi[i]; in KSPComputeShifts_DGMRES()
158 wi[max_k + i] = agmres->wi[i]; in KSPComputeShifts_DGMRES()
160 PetscCall(KSPAGMRESLejaOrdering(wr, wi, Rshift, Ishift, 2 * max_k)); in KSPComputeShifts_DGMRES()
168 PetscCall(PetscFree(wi)); in KSPComputeShifts_DGMRES()
502 PetscCall(PetscFree(agmres->wi)); in KSPDestroy_AGMRES()
|
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | gmreig.c | 120 PetscScalar *wr, *wi; /* Real and imaginary part of the Ritz values */ in KSPComputeRitz_GMRES() local
140 PetscCall(PetscMalloc4(bN * bN, &H, bn * bn, &Q, bn, &wr, bn, &wi)); in KSPComputeRitz_GMRES()
187 …PetscCallBLAS("LAPACKgeev", LAPACKgeev_("N", "V", &bn, H, &bN, wr, wi, &sdummy, &idummy, Q, &bn, w… in KSPComputeRitz_GMRES()
202 for (i = 0; i < bn; i++) modul[i] = PetscSqrtReal(wr[i] * wr[i] + wi[i] * wi[i]); in KSPComputeRitz_GMRES()
226 if (!wi[perm[nb]]) nb += 1; in KSPComputeRitz_GMRES()
235 tetai[i] = wi[perm[i]]; in KSPComputeRitz_GMRES()
240 if (wi[perm[bn - nb - 1]] == 0) nb += 1; in KSPComputeRitz_GMRES()
249 tetai[i] = wi[perm[bn - nb + i]]; in KSPComputeRitz_GMRES()
255 PetscCall(PetscFree4(H, Q, wr, wi)); in KSPComputeRitz_GMRES()
|
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmres.c | 64 dgmres->wi = NULL; in KSPSetUp_DGMRES()
286 PetscCall(PetscFree(dgmres->wi)); in KSPDestroy_DGMRES()
661 PetscReal *wr, *wi, *modul; /* Real and imaginary part and modulus of the eigenvalues of A */ in KSPDGMRESComputeSchurForm_DGMRES() local
685 PetscCall(PetscMalloc1(n, &dgmres->wi)); in KSPDGMRESComputeSchurForm_DGMRES()
688 wi = dgmres->wi; in KSPDGMRESComputeSchurForm_DGMRES()
722 …PetscCallBLAS("LAPACKhseqr", LAPACKhseqr_("S", "I", &bn, &ilo, &ihi, A, &ldA, wr, wi, Q, &ldQ, wor… in KSPDGMRESComputeSchurForm_DGMRES()
728 for (i = 0; i < n; i++) modul[i] = PetscSqrtReal(wr[i] * wr[i] + wi[i] * wi[i]); in KSPDGMRESComputeSchurForm_DGMRES()
737 if (wi[perm[NbrEig]] != 0) NbrEig += 2; in KSPDGMRESComputeSchurForm_DGMRES()
758 …PetscCallBLAS("LAPACKtrsen", LAPACKtrsen_("B", "V", select, &bn, A, &ldA, Q, &ldQ, wr, wi, &NbrEig… in KSPDGMRESComputeSchurForm_DGMRES()
847 …PetscReal *wr, *wi, *beta, *modul; /* Real and imaginary part and modulus of the eigenvalues of… in KSPDGMRESImproveEig_DGMRES() local
[all …]
|
| H A D | dgmresimpl.h | 21 PetscReal *wr, *wi, *modul; /* Real and complex part and modulus of eigenvalues */ \
|
| /petsc/src/ts/utils/dmplexlandau/kokkos/ |
| H A D | landau.kokkos.cxx | 612 const PetscReal wi = d_w[ipidx], x = d_x[ipidx], y = d_y[ipidx]; in LandauKokkosJacobian() local
638 temp1[0] *= wi; in LandauKokkosJacobian()
639 temp1[1] *= wi; in LandauKokkosJacobian()
641 temp1[2] *= wi; in LandauKokkosJacobian()
643 temp2 *= wi; in LandauKokkosJacobian()
|
| /petsc/src/ts/utils/dmplexlandau/ |
| H A D | plexland.c | 353 const PetscReal wi = ww[ipidx], x = xx[ipidx], y = yy[ipidx]; in LandauFormJacobian_Internal() local
376 temp1[0] *= wi; in LandauFormJacobian_Internal()
377 temp1[1] *= wi; in LandauFormJacobian_Internal()
379 temp1[2] *= wi; in LandauFormJacobian_Internal()
381 temp2 *= wi; in LandauFormJacobian_Internal()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex70.c | 132 static PetscInt map_wIwDI_uJuDJ(PetscInt wi, PetscInt wd, PetscInt w_NPE, PetscInt w_dof, PetscInt … in map_wIwDI_uJuDJ() argument
137 r = w_dof * wi + wd; in map_wIwDI_uJuDJ()
|
| H A D | ex43.c | 470 static PetscInt ASS_MAP_wIwDI_uJuDJ(PetscInt wi, PetscInt wd, PetscInt w_NPE, PetscInt w_dof, Petsc… in ASS_MAP_wIwDI_uJuDJ() argument
476 r = w_dof * wi + wd; in ASS_MAP_wIwDI_uJuDJ()
|
| H A D | ex42.c | 544 static PetscInt ASS_MAP_wIwDI_uJuDJ(PetscInt wi, PetscInt wd, PetscInt w_NPE, PetscInt w_dof, Petsc… in ASS_MAP_wIwDI_uJuDJ() argument
552 r = w_dof * wi + wd; in ASS_MAP_wIwDI_uJuDJ()
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/ |
| H A D | dspacelagrange.c | 132 PetscReal wi = nodesets[degree][degree - tup[i]]; in PetscNodeRecursive_Internal() local
136 for (j = 0; j < dim + 1; j++) node[j + (j >= i)] += wi * node[dim + 1 + j]; in PetscNodeRecursive_Internal()
137 w += wi; in PetscNodeRecursive_Internal()
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 1711 PetscReal wi = w[i]; in PetscDTGaussJacobiQuadrature_Internal() local
1716 w[i] = w[j] = (wi + wj) / 2.; in PetscDTGaussJacobiQuadrature_Internal()
|