| /petsc/src/ksp/ksp/impls/gmres/pgmres/ |
| H A D | pgmres.c | 60 *HH(it - 1, it - 2) = newnorm; in KSPPGMRESCycle()
63 PetscCall(VecMDotEnd(Zcur, it, &(VEC_VV(0)), HH(0, it - 1))); in KSPPGMRESCycle()
68 PetscCall(VecScale(VEC_VV(it - 1), 1. / *HH(it - 1, it - 2))); in KSPPGMRESCycle()
92 PetscCall(VecScale(Zcur, 1. / *HH(it - 1, it - 2))); in KSPPGMRESCycle()
94 PetscCall(VecScale(Znext, 1. / *HH(it - 1, it - 2))); in KSPPGMRESCycle()
97 for (k = 0; k < it; k++) *HH(k, it - 1) /= *HH(it - 1, it - 2); in KSPPGMRESCycle()
100 *HH(it - 1, it - 1) /= *HH(it - 1, it - 2); in KSPPGMRESCycle()
128 for (j = PetscMax(0, k - 1); j < it - 1; j++) work[k] -= *HES(k, j) * *HH(j, it - 1); in KSPPGMRESCycle()
131 PetscCall(VecAXPY(Znext, -*HH(it - 1, it - 1), Zcur)); in KSPPGMRESCycle()
134 for (k = 0; k < it; k++) work[k] = -*HH(k, it - 1); in KSPPGMRESCycle()
[all …]
|
| H A D | pgmresimpl.h | 10 #define HH(a, b) (pgmres->hh_origin + (b) * (pgmres->max_k + 2) + (a)) macro
|
| /petsc/src/ksp/ksp/impls/gmres/fgmres/ |
| H A D | fgmres.c | 137 *HH(loc_it + 1, loc_it) = tt; in KSPFGMRESCycle()
276 if (*HH(it, it) != 0.0) { in KSPFGMRESBuildSoln()
277 nrs[it] = *RS(it) / *HH(it, it); in KSPFGMRESBuildSoln()
284 for (j = k + 1; j <= it; j++) tt = tt - *HH(k, j) * nrs[j]; in KSPFGMRESBuildSoln()
285 nrs[k] = tt / *HH(k, k); in KSPFGMRESBuildSoln()
309 hh = HH(0, it); /* pointer to beginning of column to update - so in KSPFGMRESUpdateHessenberg()
|
| H A D | fgmresimpl.h | 23 #define HH(a, b) (fgmres->hh_origin + (b) * (fgmres->max_k + 2) + (a)) macro
|
| /petsc/src/ksp/ksp/impls/gmres/lgmres/ |
| H A D | lgmres.c | 191 *HH(loc_it + 1, loc_it) = tt; in KSPLGMRESCycle()
392 …PetscCheck(*HH(it, it) != 0.0, PETSC_COMM_SELF, PETSC_ERR_CONV_FAILED, "HH(it,it) is identically z… in KSPLGMRESBuildSoln()
393 if (*HH(it, it) != 0.0) { in KSPLGMRESBuildSoln()
394 nrs[it] = *GRS(it) / *HH(it, it); in KSPLGMRESBuildSoln()
402 for (j = k + 1; j <= it; j++) tt = tt - *HH(k, j) * nrs[j]; in KSPLGMRESBuildSoln()
403 nrs[k] = tt / *HH(k, k); in KSPLGMRESBuildSoln()
446 …hh = HH(0, it); /* pointer to beginning of column to update - so incrementing hh "steps down" the … in KSPLGMRESUpdateHessenberg()
|
| H A D | lgmresimpl.h | 35 #define HH(a, b) (lgmres->hh_origin + (b) * (lgmres->max_k + 2) + (a)) macro
|
| /petsc/src/ksp/ksp/impls/gmres/pipefgmres/ |
| H A D | pipefgmres.c | 170 hh = HH(0, loc_it); in KSPPIPEFGMRESCycle()
357 if (*HH(it, it) != 0.0) nrs[it] = *RS(it) / *HH(it, it); in KSPPIPEFGMRESBuildSoln()
362 for (j = k + 1; j <= it; j++) tt -= *HH(k, j) * nrs[j]; in KSPPIPEFGMRESBuildSoln()
363 nrs[k] = tt / *HH(k, k); in KSPPIPEFGMRESBuildSoln()
386 hh = HH(0, it); /* pointer to beginning of column to update */ in KSPPIPEFGMRESUpdateHessenberg()
397 …etscInt_FMT ") = %14.12e\n", (double)hapbnd, it + 1, it, (double)PetscAbsScalar(*HH(it + 1, it)))); in KSPPIPEFGMRESUpdateHessenberg()
|
| H A D | pipefgmresimpl.h | 26 #define HH(a, b) (pipefgmres->hh_origin + (b) * (pipefgmres->max_k + 2) + (a)) macro
|
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | borthog.c | 39 hh = HH(0, it); in KSPGMRESModifiedGramSchmidtOrthogonalization()
|
| H A D | gmres.c | 144 *HH(it + 1, it) = tt; in KSPGMRESCycle()
313 if (*HH(it, it) != 0.0) { in KSPGMRESBuildSoln()
314 nrs[it] = *GRS(it) / *HH(it, it); in KSPGMRESBuildSoln()
325 for (j = k + 1; j <= it; j++) tt = tt - *HH(k, j) * nrs[j]; in KSPGMRESBuildSoln()
326 if (*HH(k, k) == 0.0) { in KSPGMRESBuildSoln()
332 nrs[k] = tt / *HH(k, k); in KSPGMRESBuildSoln()
354 hh = HH(0, it); in KSPGMRESUpdateHessenberg()
|
| H A D | borthog2.c | 49 hh = HH(0, it); in KSPGMRESClassicalGramSchmidtOrthogonalization()
|
| H A D | gmresimpl.h | 78 #define HH(a, b) (gmres->hh_origin + (b) * (gmres->max_k + 2) + (a)) macro
|
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmresimpl.h | 53 #define HH(a, b) (dgmres->hh_origin + (b) * (dgmres->max_k + 2) + (a)) macro
|
| H A D | dgmres.c | 165 *HH(it + 1, it) = tt; in KSPDGMRESCycle()
333 …PetscCheck(*HH(it, it) != 0.0, PetscObjectComm((PetscObject)ksp), PETSC_ERR_CONV_FAILED, "Likely y… in KSPDGMRESBuildSoln()
334 if (*HH(it, it) != 0.0) nrs[it] = *GRS(it) / *HH(it, it); in KSPDGMRESBuildSoln()
340 for (j = k + 1; j <= it; j++) tt = tt - *HH(k, j) * nrs[j]; in KSPDGMRESBuildSoln()
341 …PetscCheck(*HH(k, k) != 0.0, PetscObjectComm((PetscObject)ksp), PETSC_ERR_CONV_FAILED, "Likely you… in KSPDGMRESBuildSoln()
342 nrs[k] = tt / *HH(k, k); in KSPDGMRESBuildSoln()
371 hh = HH(0, it); in KSPDGMRESUpdateHessenberg()
|
| /petsc/doc/ |
| H A D | petsc.bib | 415 author = {Dupont, Todd and Kendall, Richard P and Rachford, Jr, HH},
|