Lines Matching refs:hmnorm
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()
974 …ment norm %10.2e, estimated error norms %10.2e %10.2e %10.2e\n", (double)hmnorm[0], (double)enorm[… in TSSolve_GLLE()
976 PetscCall(TSGLLEChooseNextScheme(ts, h, hmnorm, &next_scheme, &next_h, &final_step)); in TSSolve_GLLE()