25   #undef PetscSqrtScalar
26   #define PetscSqrtScalar(x) sqrt(x)  macro
434   cl = PetscSqrtScalar(*gaml * *pl / *rl);  in riem1mdt()
435   cr = PetscSqrtScalar(*gamr * *pr / *rr);  in riem1mdt()
443   skpr1   = cr * (pst - 1.) * PetscSqrtScalar(2. / (*gamr * (*gamr - 1. + (*gamr + 1.) * pst)));  in riem1mdt()
447   skpr2   = cl * (pst - 1.) * PetscSqrtScalar(2. / (*gaml * (*gaml - 1. + (*gaml + 1.) * pst)));  in riem1mdt()
476       cstarl  = PetscSqrtScalar(*gaml * *pstar / *rstarl);  in riem1mdt()
482       cstarr  = PetscSqrtScalar(*gamr * *pstar / *rstarr);  in riem1mdt()
500 …ustarl  = *uxl - (pst - 1.) * cl * PetscSqrtScalar(2. / (*gaml * (*gaml - 1. + (*gaml + 1.) * pst)…  in riem1mdt()
501 …zl      = *pl / cl * PetscSqrtScalar(*gaml * 2. * (*gaml - 1. + (*gaml + 1.) * pst)) * (*gaml - 1.…  in riem1mdt()
505       cstarr  = PetscSqrtScalar(*gamr * *pstar / *rstarr);  in riem1mdt()
523 …ustarl = *uxl - (pst - 1.) * cl * PetscSqrtScalar(2. / (*gaml * (*gaml - 1. + (*gaml + 1.) * pst))…  in riem1mdt()
524 …zl     = *pl / cl * PetscSqrtScalar(*gaml * 2. * (*gaml - 1. + (*gaml + 1.) * pst)) * (*gaml - 1. …  in riem1mdt()
526 …ustarr = *uxr + (pst - 1.) * cr * PetscSqrtScalar(2. / (*gamr * (*gamr - 1. + (*gamr + 1.) * pst))…  in riem1mdt()
527 …zr     = *pr / cr * PetscSqrtScalar(*gamr * 2. * (*gamr - 1. + (*gamr + 1.) * pst)) * (*gamr - 1. …  in riem1mdt()
545       cstarl  = PetscSqrtScalar(*gaml * *pstar / *rstarl);  in riem1mdt()
549 …ustarr  = *uxr + (pst - 1.) * cr * PetscSqrtScalar(2. / (*gamr * (*gamr - 1. + (*gamr + 1.) * pst)…  in riem1mdt()
550 …zr      = *pr / cr * PetscSqrtScalar(*gamr * 2. * (*gamr - 1. + (*gamr + 1.) * pst)) * (*gamr - 1.…  in riem1mdt()
625   c0     = PetscSqrtScalar(*gam * p0 / r0);  in riemannsolver()
633   w0     = PetscSqrtScalar(w0);  in riemannsolver()
634   cstar  = PetscSqrtScalar(*gam * pstar / rstar);  in riemannsolver()
712     tmp  = PetscSqrtScalar(d__1 * d__1 + d__2 * d__2 + d__3 * d__3);  in godunovflux()