412   PetscScalar d__1, d__2;  in riem1mdt()  local
444   d__1    = (*gamr - 1.) / (*gamr * 2.);  in riem1mdt()
445   rarepr2 = gascr4 * 2. * cr * (1. - PetscPowScalar(pst, d__1));  in riem1mdt()
448   d__1    = (*gaml - 1.) / (*gaml * 2.);  in riem1mdt()
449   rarepr1 = gascl4 * 2. * cl * (1. - PetscPowScalar(pst, d__1));  in riem1mdt()
473       d__1    = *pstar / *pl;  in riem1mdt()
475       *rstarl = *rl * PetscPowScalar(d__1, d__2);  in riem1mdt()
479       d__1    = *pstar / *pr;  in riem1mdt()
481       *rstarr = *rr * PetscPowScalar(d__1, d__2);  in riem1mdt()
502       d__1    = *pstar / *pr;  in riem1mdt()
504       *rstarr = *rr * PetscPowScalar(d__1, d__2);  in riem1mdt()
542       d__1    = *pstar / *pl;  in riem1mdt()
544       *rstarl = *rl * PetscPowScalar(d__1, d__2);  in riem1mdt()
585   PetscScalar d__1, d__2;  in riemannsolver()  local
613   d__1 = *xp - x2;  in riemannsolver()
614   sgn0 = sign(d__1);  in riemannsolver()
629   d__1   = p0 / pstar;  in riemannsolver()
631   rstarr = r0 * PetscPowScalar(d__1, d__2);  in riemannsolver()
648   d__1 = cx * cx / (*gam * s);  in riemannsolver()
649   *rx  = PetscPowScalar(d__1, gasc4);  in riemannsolver()
651   d__1 = sgn0 * (x0 - *xp);  in riemannsolver()
652   *rx  = cvmgp_(rx, &r0, &d__1);  in riemannsolver()
653   d__1 = sgn0 * (x0 - *xp);  in riemannsolver()
654   *px  = cvmgp_(px, &p0, &d__1);  in riemannsolver()
655   d__1 = sgn0 * (x0 - *xp);  in riemannsolver()
656   *uxm = cvmgp_(uxm, &u0, &d__1);  in riemannsolver()
657   d__1 = sgn0 * (xstar - *xp);  in riemannsolver()
658   *rx  = cvmgm_(rx, &rstar, &d__1);  in riemannsolver()
659   d__1 = sgn0 * (xstar - *xp);  in riemannsolver()
660   *px  = cvmgm_(px, &pstar, &d__1);  in riemannsolver()
661   d__1 = sgn0 * (xstar - *xp);  in riemannsolver()
662   *uxm = cvmgm_(uxm, &ustar, &d__1);  in riemannsolver()
679   PetscScalar d__1, d__2, d__3;  in godunovflux()  local
707     d__1 = bn[0];  in godunovflux()
712     tmp  = PetscSqrtScalar(d__1 * d__1 + d__2 * d__2 + d__3 * d__3);  in godunovflux()
751   d__1 = uxl;  in godunovflux()
756   pl   = (gamma - 1.) * (ul[ndim + 1] - rl * .5 * (d__1 * d__1 + d__2 * d__2 + d__3 * d__3));  in godunovflux()
758   d__1 = uxr;  in godunovflux()
763   pr    = (gamma - 1.) * (ur[ndim + 1] - rr * .5 * (d__1 * d__1 + d__2 * d__2 + d__3 * d__3));  in godunovflux()