404 static int riem1mdt(PetscScalar *gaml, PetscScalar *gamr, PetscScalar *rl, PetscScalar *pl, PetscSc…  in riem1mdt()  argument
434   cl = PetscSqrtScalar(*gaml * *pl / *rl);  in riem1mdt()
440   *pstar  = (wl * *pr + wr * *pl) / (wl + wr);  in riem1mdt()
442   pst     = *pl / *pr;  in riem1mdt()
446   pst     = *pr / *pl;  in riem1mdt()
451   if (PetscRealPart(*pr) < PetscRealPart(*pl)) {  in riem1mdt()
473       d__1    = *pstar / *pl;  in riem1mdt()
499       pst     = *pstar / *pl;  in riem1mdt()
501 …zl      = *pl / cl * PetscSqrtScalar(*gaml * 2. * (*gaml - 1. + (*gaml + 1.) * pst)) * (*gaml - 1.…  in riem1mdt()
522       pst    = *pstar / *pl;  in riem1mdt()
524 …zl     = *pl / cl * PetscSqrtScalar(*gaml * 2. * (*gaml - 1. + (*gaml + 1.) * pst)) * (*gaml - 1. …  in riem1mdt()
542       d__1    = *pstar / *pl;  in riem1mdt()
563   if (PetscRealPart(*pstar) > PetscRealPart(*pl)) {  in riem1mdt()
564     pst     = *pstar / *pl;  in riem1mdt()
582 …lar *xp, PetscScalar *dtt, PetscScalar *rl, PetscScalar *uxl, PetscScalar *pl, PetscScalar *utl, P…  in riemannsolver()  argument
592   if (*rl == *rr && *pr == *pl && *uxl == *uxr && *gaml == *gamr) {  in riemannsolver()
594     *px  = *pl;  in riemannsolver()
610   iwave = riem1mdt(gaml, gamr, rl, pl, uxl, rr, pr, uxr, &rstarl, &rstarr, &pstar, &ustar);  in riemannsolver()
618   p0     = cvmgm_(pl, pr, &sgn0);  in riemannsolver()
683 …static PetscScalar bn[3], fn, ft, tg[3], pl, rl, pm, pr, rr, xp, ubl, ubm, ubr, dtt, unm, tmp, utl…  in godunovflux()  local
756   pl   = (gamma - 1.) * (ul[ndim + 1] - rl * .5 * (d__1 * d__1 + d__2 * d__2 + d__3 * d__3));  in godunovflux()
767 …iwave = riemannsolver(&xcen, &xp, &dtt, &rl, &uxl, &pl, &utl, &ubl, &gaml, &rho1l, &rr, &uxr, &pr,…  in godunovflux()