404 static int riem1mdt(PetscScalar *gaml, PetscScalar *gamr, PetscScalar *rl, PetscScalar *pl, PetscSc…  in riem1mdt()  argument
426   gascl4 = 1. / (*gaml - 1.);  in riem1mdt()
434   cl = PetscSqrtScalar(*gaml * *pl / *rl);  in riem1mdt()
447   skpr2   = cl * (pst - 1.) * PetscSqrtScalar(2. / (*gaml * (*gaml - 1. + (*gaml + 1.) * pst)));  in riem1mdt()
448   d__1    = (*gaml - 1.) / (*gaml * 2.);  in riem1mdt()
474       d__2    = 1. / *gaml;  in riem1mdt()
476       cstarl  = PetscSqrtScalar(*gaml * *pstar / *rstarl);  in riem1mdt()
500 …ustarl  = *uxl - (pst - 1.) * cl * PetscSqrtScalar(2. / (*gaml * (*gaml - 1. + (*gaml + 1.) * pst)…  in riem1mdt()
501 …l * PetscSqrtScalar(*gaml * 2. * (*gaml - 1. + (*gaml + 1.) * pst)) * (*gaml - 1. + (*gaml + 1.) *…  in riem1mdt()
523 …ustarl = *uxl - (pst - 1.) * cl * PetscSqrtScalar(2. / (*gaml * (*gaml - 1. + (*gaml + 1.) * pst))…  in riem1mdt()
524 …l * PetscSqrtScalar(*gaml * 2. * (*gaml - 1. + (*gaml + 1.) * pst)) * (*gaml - 1. + (*gaml + 1.) *…  in riem1mdt()
543       d__2    = 1. / *gaml;  in riem1mdt()
545       cstarl  = PetscSqrtScalar(*gaml * *pstar / *rstarl);  in riem1mdt()
565     *rstarl = ((*gaml + 1.) * pst + *gaml - 1.) / ((*gaml - 1.) * pst + *gaml + 1.) * *rl;  in riem1mdt()
582 …r *uxl, PetscScalar *pl, PetscScalar *utl, PetscScalar *ubl, PetscScalar *gaml, PetscScalar *rho1l…  in riemannsolver()  argument
592   if (*rl == *rr && *pr == *pl && *uxl == *uxr && *gaml == *gamr) {  in riemannsolver()
596     *gam = *gaml;  in riemannsolver()
610   iwave = riem1mdt(gaml, gamr, rl, pl, uxl, rr, pr, uxr, &rstarl, &rstarr, &pstar, &ustar);  in riemannsolver()
620   *gam   = cvmgm_(gaml, gamr, &sgn0);  in riemannsolver()
683 …l, pm, pr, rr, xp, ubl, ubm, ubr, dtt, unm, tmp, utl, utm, uxl, utr, uxr, gaml, gamm, gamr, xcen, …  in godunovflux()  local
748   gaml = gamma;  in godunovflux()
767 …iwave = riemannsolver(&xcen, &xp, &dtt, &rl, &uxl, &pl, &utl, &ubl, &gaml, &rho1l, &rr, &uxr, &pr,…  in godunovflux()