404 … PetscScalar *uxr, PetscScalar *rstarl, PetscScalar *rstarr, PetscScalar *pstar, PetscScalar *usta…  in riem1mdt()  argument
440   *pstar  = (wl * *pr + wr * *pl) / (wl + wr);  in riem1mdt()
441   *pstar  = PetscMax(PetscRealPart(*pstar), PetscRealPart(smallp));  in riem1mdt()
473       d__1    = *pstar / *pl;  in riem1mdt()
476       cstarl  = PetscSqrtScalar(*gaml * *pstar / *rstarl);  in riem1mdt()
479       d__1    = *pstar / *pr;  in riem1mdt()
482       cstarr  = PetscSqrtScalar(*gamr * *pstar / *rstarr);  in riem1mdt()
486       *pstar -= dpstar;  in riem1mdt()
487       *pstar = PetscMax(PetscRealPart(*pstar), PetscRealPart(smallp));  in riem1mdt()
488       if (PetscAbsScalar(dpstar) / PetscRealPart(*pstar) <= 1e-8) {  in riem1mdt()
499       pst     = *pstar / *pl;  in riem1mdt()
502       d__1    = *pstar / *pr;  in riem1mdt()
505       cstarr  = PetscSqrtScalar(*gamr * *pstar / *rstarr);  in riem1mdt()
509       *pstar -= dpstar;  in riem1mdt()
510       *pstar = PetscMax(PetscRealPart(*pstar), PetscRealPart(smallp));  in riem1mdt()
511       if (PetscAbsScalar(dpstar) / PetscRealPart(*pstar) <= 1e-8) {  in riem1mdt()
522       pst    = *pstar / *pl;  in riem1mdt()
525       pst    = *pstar / *pr;  in riem1mdt()
529       *pstar -= dpstar;  in riem1mdt()
530       *pstar = PetscMax(PetscRealPart(*pstar), PetscRealPart(smallp));  in riem1mdt()
531       if (PetscAbsScalar(dpstar) / PetscRealPart(*pstar) <= 1e-8) {  in riem1mdt()
542       d__1    = *pstar / *pl;  in riem1mdt()
545       cstarl  = PetscSqrtScalar(*gaml * *pstar / *rstarl);  in riem1mdt()
548       pst     = *pstar / *pr;  in riem1mdt()
552       *pstar -= dpstar;  in riem1mdt()
553       *pstar = PetscMax(PetscRealPart(*pstar), PetscRealPart(smallp));  in riem1mdt()
554       if (PetscAbsScalar(dpstar) / PetscRealPart(*pstar) <= 1e-8) {  in riem1mdt()
563   if (PetscRealPart(*pstar) > PetscRealPart(*pl)) {  in riem1mdt()
564     pst     = *pstar / *pl;  in riem1mdt()
567   if (PetscRealPart(*pstar) > PetscRealPart(*pr)) {  in riem1mdt()
568     pst     = *pstar / *pr;  in riem1mdt()
589 …static PetscScalar cstar, pstar, rstar, ustar, xstar, wspst, ushock, streng, rstarl, rstarr, rstar…  in riemannsolver()  local
610   iwave = riem1mdt(gaml, gamr, rl, pl, uxl, rr, pr, uxr, &rstarl, &rstarr, &pstar, &ustar);  in riemannsolver()
626   streng = pstar - p0;  in riemannsolver()
629   d__1   = p0 / pstar;  in riemannsolver()
634   cstar  = PetscSqrtScalar(*gam * pstar / rstar);  in riemannsolver()
660   *px  = cvmgm_(px, &pstar, &d__1);  in riemannsolver()