450       fextra[0] = pv2e * (p_current[0] - p_target[0])*3.0*scale*scale*vol0;  in min_energy()
453       fextra[0] = pv2e * (p_current[0] - p_target[0])*2.0*scale*vol0;  in min_energy()
465       fextra[0] = pv2e * (p_current[0] - p_hydro)*scaley*scalez*vol0;  in min_energy()
467       fextra[1] = pv2e * (p_current[1] - p_hydro)*scalex*scalez*vol0;  in min_energy()
469       fextra[2] = pv2e * (p_current[2] - p_hydro)*scalex*scaley*vol0;  in min_energy()
474         fextra[3] = pv2e*p_current[3]*scaley*yprdinit*scalex*xprdinit*yprdinit;  in min_energy()
476         fextra[4] = pv2e*p_current[4]*scalex*xprdinit*scaley*yprdinit*xprdinit;  in min_energy()
478         fextra[5] = pv2e*p_current[5]*scalex*xprdinit*scalez*zprdinit*xprdinit;  in min_energy()
709     p_current[0] = p_current[1] = p_current[2] = pressure->scalar;  in couple()
712     p_current[0] = p_current[1] = p_current[2] = ave;  in couple()
715     p_current[0] = p_current[1] = ave;  in couple()
716     p_current[2] = tensor[2];  in couple()
719     p_current[1] = p_current[2] = ave;  in couple()
720     p_current[0] = tensor[0];  in couple()
723     p_current[0] = p_current[2] = ave;  in couple()
724     p_current[1] = tensor[1];  in couple()
726     p_current[0] = tensor[0];  in couple()
727     p_current[1] = tensor[1];  in couple()
728     p_current[2] = tensor[2];  in couple()
731   if (!std::isfinite(p_current[0]) || !std::isfinite(p_current[1]) || !std::isfinite(p_current[2]))  in couple()
737     p_current[3] = tensor[5];  in couple()
738     p_current[4] = tensor[4];  in couple()
739     p_current[5] = tensor[3];  in couple()
741 …  if (!std::isfinite(p_current[3]) || !std::isfinite(p_current[4]) || !std::isfinite(p_current[5]))  in couple()