| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | borthog.c | 34 PetscScalar *hh, *hes; in KSPGMRESModifiedGramSchmidtOrthogonalization() local
39 hh = HH(0, it); in KSPGMRESModifiedGramSchmidtOrthogonalization()
43 PetscCall(VecDot(VEC_VV(it + 1), VEC_VV(j), hh)); in KSPGMRESModifiedGramSchmidtOrthogonalization()
44 KSPCheckDot(ksp, *hh); in KSPGMRESModifiedGramSchmidtOrthogonalization()
46 *hes++ = *hh; in KSPGMRESModifiedGramSchmidtOrthogonalization()
48 PetscCall(VecAXPY(VEC_VV(it + 1), -(*hh++), VEC_VV(j))); in KSPGMRESModifiedGramSchmidtOrthogonalization()
|
| H A D | borthog2.c | 39 PetscScalar *hh, *hes, *lhh; in KSPGMRESClassicalGramSchmidtOrthogonalization() local
49 hh = HH(0, it); in KSPGMRESClassicalGramSchmidtOrthogonalization()
54 hh[j] = 0.0; in KSPGMRESClassicalGramSchmidtOrthogonalization()
76 hh[j] -= lhh[j]; /* hh += <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
108 hh[j] -= lhh[j]; /* hh += <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
|
| H A D | gmres.c | 37 PetscInt hh, hes, rs, cc; in KSPSetUp_GMRES() local
43 hh = (max_k + 2) * (max_k + 1); in KSPSetUp_GMRES()
48 …PetscCall(PetscCalloc5(hh, &gmres->hh_origin, hes, &gmres->hes_origin, rs, &gmres->rs_origin, cc, … in KSPSetUp_GMRES()
349 PetscScalar *hh, *cc, *ss, tt; in KSPGMRESUpdateHessenberg() local
354 hh = HH(0, it); in KSPGMRESUpdateHessenberg()
361 tt = *hh; in KSPGMRESUpdateHessenberg()
362 *hh = PetscConj(*cc) * tt + *ss * *(hh + 1); in KSPGMRESUpdateHessenberg()
363 hh++; in KSPGMRESUpdateHessenberg()
364 *hh = *cc++ * *hh - (*ss++ * tt); in KSPGMRESUpdateHessenberg()
374 tt = PetscSqrtScalar(PetscConj(*hh) * *hh + PetscConj(*(hh + 1)) * *(hh + 1)); in KSPGMRESUpdateHessenberg()
[all …]
|
| /petsc/src/ksp/ksp/impls/gmres/pgmres/ |
| H A D | pgmres.c | 238 PetscScalar *hh, *cc, *ss, *rs; in KSPPGMRESUpdateHessenberg() local
244 hh = HH(0, it); /* pointer to beginning of column to update */ in KSPPGMRESUpdateHessenberg()
250 for (j = 0; j <= it + 1; j++) *HES(j, it) = hh[j]; in KSPPGMRESUpdateHessenberg()
253 hapbnd = PetscMin(PetscAbsScalar(hh[it + 1] / rs[it]), pgmres->haptol); in KSPPGMRESUpdateHessenberg()
254 if (PetscAbsScalar(hh[it + 1]) < hapbnd) { in KSPPGMRESUpdateHessenberg()
264 PetscScalar hhj = hh[j]; in KSPPGMRESUpdateHessenberg()
265 hh[j] = PetscConj(cc[j]) * hhj + ss[j] * hh[j + 1]; in KSPPGMRESUpdateHessenberg()
266 hh[j + 1] = -ss[j] * hhj + cc[j] * hh[j + 1]; in KSPPGMRESUpdateHessenberg()
282 …PetscReal delta = PetscSqrtReal(PetscSqr(PetscAbsScalar(hh[it])) + PetscSqr(PetscAbsScalar(hh[it +… in KSPPGMRESUpdateHessenberg()
288 cc[it] = hh[it] / delta; /* new cosine value */ in KSPPGMRESUpdateHessenberg()
[all …]
|
| /petsc/src/ksp/ksp/impls/gmres/pipefgmres/ |
| H A D | pipefgmres.c | 52 PetscScalar *hh, *hes, *lhh, shift = pipefgmres->shift; in KSPPIPEFGMRESCycle() local
170 hh = HH(0, loc_it); in KSPPIPEFGMRESCycle()
173 hh[j] = lhh[j]; in KSPPIPEFGMRESCycle()
176 hh[loc_it] = lhh[loc_it] + shift; in KSPPIPEFGMRESCycle()
201 hh[loc_it + 1] = tt; in KSPPIPEFGMRESCycle()
380 PetscScalar *hh, *cc, *ss, *rs; in KSPPIPEFGMRESUpdateHessenberg() local
386 hh = HH(0, it); /* pointer to beginning of column to update */ in KSPPIPEFGMRESUpdateHessenberg()
392 for (j = 0; j <= it + 1; j++) *HES(j, it) = hh[j]; in KSPPIPEFGMRESUpdateHessenberg()
395 hapbnd = PetscMin(PetscAbsScalar(hh[it + 1] / rs[it]), pipefgmres->haptol); in KSPPIPEFGMRESUpdateHessenberg()
396 if (PetscAbsScalar(hh[it + 1]) < hapbnd) { in KSPPIPEFGMRESUpdateHessenberg()
[all …]
|
| /petsc/src/ksp/ksp/impls/gmres/fgmres/ |
| H A D | fgmres.c | 304 PetscScalar *hh, *cc, *ss, tt; in KSPFGMRESUpdateHessenberg() local
309 hh = HH(0, it); /* pointer to beginning of column to update - so in KSPFGMRESUpdateHessenberg()
320 tt = *hh; in KSPFGMRESUpdateHessenberg()
321 *hh = PetscConj(*cc) * tt + *ss * *(hh + 1); in KSPFGMRESUpdateHessenberg()
322 hh++; in KSPFGMRESUpdateHessenberg()
323 *hh = *cc++ * *hh - (*ss++ * tt); in KSPFGMRESUpdateHessenberg()
340 tt = PetscSqrtScalar(PetscConj(*hh) * *hh + PetscConj(*(hh + 1)) * *(hh + 1)); in KSPFGMRESUpdateHessenberg()
346 *cc = *hh / tt; /* new cosine value */ in KSPFGMRESUpdateHessenberg()
347 *ss = *(hh + 1) / tt; /* new sine value */ in KSPFGMRESUpdateHessenberg()
352 *hh = PetscConj(*cc) * *hh + *ss * *(hh + 1); in KSPFGMRESUpdateHessenberg()
|
| /petsc/src/ksp/ksp/impls/gmres/lgmres/ |
| H A D | lgmres.c | 441 PetscScalar *hh, *cc, *ss, tt; in KSPLGMRESUpdateHessenberg() local
446 …hh = HH(0, it); /* pointer to beginning of column to update - so incrementing hh "steps down" the … in KSPLGMRESUpdateHessenberg()
455 tt = *hh; in KSPLGMRESUpdateHessenberg()
456 *hh = PetscConj(*cc) * tt + *ss * *(hh + 1); in KSPLGMRESUpdateHessenberg()
457 hh++; in KSPLGMRESUpdateHessenberg()
458 *hh = *cc++ * *hh - (*ss++ * tt); in KSPLGMRESUpdateHessenberg()
475 tt = PetscSqrtScalar(PetscConj(*hh) * *hh + PetscConj(*(hh + 1)) * *(hh + 1)); in KSPLGMRESUpdateHessenberg()
480 *cc = *hh / tt; /* new cosine value */ in KSPLGMRESUpdateHessenberg()
481 *ss = *(hh + 1) / tt; /* new sine value */ in KSPLGMRESUpdateHessenberg()
486 *hh = PetscConj(*cc) * *hh + *ss * *(hh + 1); in KSPLGMRESUpdateHessenberg()
|
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmres.c | 366 PetscScalar *hh, *cc, *ss, tt; in KSPDGMRESUpdateHessenberg() local
371 hh = HH(0, it); in KSPDGMRESUpdateHessenberg()
378 tt = *hh; in KSPDGMRESUpdateHessenberg()
379 *hh = PetscConj(*cc) * tt + *ss * *(hh + 1); in KSPDGMRESUpdateHessenberg()
380 hh++; in KSPDGMRESUpdateHessenberg()
381 *hh = *cc++ * *hh - (*ss++ * tt); in KSPDGMRESUpdateHessenberg()
391 tt = PetscSqrtScalar(PetscConj(*hh) * *hh + PetscConj(*(hh + 1)) * *(hh + 1)); in KSPDGMRESUpdateHessenberg()
396 *cc = *hh / tt; in KSPDGMRESUpdateHessenberg()
397 *ss = *(hh + 1) / tt; in KSPDGMRESUpdateHessenberg()
400 *hh = PetscConj(*cc) * *hh + *ss * *(hh + 1); in KSPDGMRESUpdateHessenberg()
|
| /petsc/src/sys/classes/draw/impls/x/ |
| H A D | xinit.c | 274 unsigned int ww = 0, hh = 0, dummy; in PetscDrawXiGetGeometry() local
278 XGetGeometry(XiWin->disp, XiWin->win, &parent, &xx, &yy, &ww, &hh, &dummy, &dummy); in PetscDrawXiGetGeometry()
286 XGetGeometry(XiWin->disp, XiWin->drw, &root, &xx, &yy, &ww, &hh, &dummy, &dummy); in PetscDrawXiGetGeometry()
291 if (h) *h = (int)hh; in PetscDrawXiGetGeometry()
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex2.c | 94 PetscReal *hh, low[3], high[3]; in PerturbVertices() local
105 PetscCall(PetscCalloc1(cdim, &hh)); in PerturbVertices()
106 for (d = 0; d < cdim; ++d) hh[d] = (user->L[d]) / PetscPowReal(cEnd, 1. / cdim); in PerturbVertices()
116 coord[d] = PetscMax(low[d], PetscMin(high[d], PetscRealPart(coord[d] + value * hh[d]))); in PerturbVertices()
121 PetscCall(PetscFree(hh)); in PerturbVertices()
|
| /petsc/ |
| H A D | .gitignore | 79 *.hh.html
|
| /petsc/doc/manual/ |
| H A D | other.md | 892 -regex '.*\.\(cc\|hh\|cpp\|cxx\|C\|hpp\|c\|h\|cu\)$' \
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 463 (ignition/press-exp-hh 0.63)
464 (ignition/activate-hh 618840)
465 (ignition/cetane-hh 25)
466 (ignition/pre-exp-hh 0.36)
|