def eosBaryonDensity(self,pressure): return np.where(pressure<self.p_match,EOS_BPS.eosBaryonDensity(pressure),self.eosBaryonDensity_match_Quarkyonic(pressure))
def eosBaryonDensity(self, pressure): return np.where( pressure < self.p_match, EOS_BPS.eosBaryonDensity(pressure), np.where(pressure < self.p_max, self.eosBaryonDensity_matchPNM(pressure), self.eosCSS.eosBaryonDensity(pressure)))