Searched refs:hmnorm (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/ts/impls/implicit/glle/ |
| H A D | glle.c | 790 static PetscErrorCode TSGLLEChooseNextScheme(TS ts, PetscReal h, const PetscReal hmnorm[], PetscInt… in TSGLLEChooseNextScheme() argument
803 if (sc->p == cur_p - 1) errors[n] = PetscAbsScalar(sc->alpha[0]) * hmnorm[0]; in TSGLLEChooseNextScheme()
804 else if (sc->p == cur_p) errors[n] = PetscAbsScalar(sc->alpha[0]) * hmnorm[1]; in TSGLLEChooseNextScheme()
805 else if (sc->p == cur_p + 1) errors[n] = PetscAbsScalar(sc->alpha[0]) * (hmnorm[2] + hmnorm[3]); in TSGLLEChooseNextScheme()
872 PetscReal h, hmnorm[4], enorm[3], next_h; in TSSolve_GLLE() local
947 for (i = 0; i < 3; i++) PetscCall(TSGLLEVecNormWRMS(ts, gl->himom[i], &hmnorm[i + 1])); in TSSolve_GLLE()
948 enorm[0] = PetscRealPart(scheme->alpha[0]) * hmnorm[1]; in TSSolve_GLLE()
949 enorm[1] = PetscRealPart(scheme->beta[0]) * hmnorm[2]; in TSSolve_GLLE()
950 enorm[2] = PetscRealPart(scheme->gamma[0]) * hmnorm[3]; in TSSolve_GLLE()
971 PetscCall(TSGLLEVecNormWRMS(ts, gl->X[scheme->r - 1], &hmnorm[0])); in TSSolve_GLLE()
[all …]
|