fixed easy exponents

Author: reinhold-willcox

src/HG.cpp

@@ -517,7 +517,7 @@ double HG::CalculateRadialExtentConvectiveEnvelope() const {
 	BaseStar clone = *this;                         // clone this star so can manipulate without changes persisiting
 	clone.ResolveEnvelopeLoss(true);                // update clone's attributes after envelope is lost
 
-    return PPOW(m_Tau, 1.0 / 2.0) * (m_Radius - clone.Radius());
+    return sqrt(m_Tau) * (m_Radius - clone.Radius());
 }
 
 

src/MainSequence.cpp

@@ -348,7 +348,7 @@ double MainSequence::CalculateDeltaR(const double p_Mass) const {
     double deltaR;
 
     if (utils::Compare(p_Mass, massCutoffs(MHook)) <= 0) deltaR = 0.0;   // this really is supposed to be 0
-    else if (utils::Compare(p_Mass, a[42])         <= 0) deltaR = a[43] * PPOW(((p_Mass - massCutoffs(MHook)) / (a[42] - massCutoffs(MHook))), 0.5);
+    else if (utils::Compare(p_Mass, a[42])         <= 0) deltaR = a[43] * sqrt((p_Mass - massCutoffs(MHook)) / (a[42] - massCutoffs(MHook)));
     else if (utils::Compare(p_Mass, 2.0)            < 0) deltaR = a[43] + ((m_RConstants[static_cast<int>(R_CONSTANTS::B_DELTA_R)] - a[43]) * PPOW(((p_Mass - a[42]) / (2.0 - a[42])), a[44]));
     else {
         // pow() is slow - use multiplication (sqrt() is faster than pow())
@@ -483,7 +483,7 @@ double MainSequence::CalculateRadiusOnPhase(const double p_Mass, const double p_
  * @return                                      Radial extent of the star's convective envelope in Rsol
  */
 double MainSequence::CalculateRadialExtentConvectiveEnvelope() const {
-    return utils::Compare(m_Mass, 0.35) <= 0 ? m_Radius * PPOW(1.0 - m_Tau, 1.0 / 4.0) : 0.0;
+    return utils::Compare(m_Mass, 0.35) <= 0 ? m_Radius * sqrt(sqrt(1.0 - m_Tau)) : 0.0;
 }
 
 

src/WhiteDwarfs.cpp

@@ -31,5 +31,7 @@ double WhiteDwarfs::CalculateLuminosityOnPhase_Static(const double p_Mass, const
  * @return                                      Radius of a White Dwarf in Rsol (since WD is ~ Earth sized, expect answer around 0.009)
  */
 double WhiteDwarfs::CalculateRadiusOnPhase_Static(const double p_Mass) {
-    return std::max(NEUTRON_STAR_RADIUS, 0.0115 * sqrt(PPOW(MCH / p_Mass, 2.0 / 3.0) - PPOW(p_Mass / MCH, 2.0 / 3.0)));
+    double MCH_Mass_one_third  = std::cbrt(MCH / p_Mass); 
+    double MCH_Mass_two_thirds = MCH_Mass_one_third* MCH_Mass_one_third
+    return std::max(NEUTRON_STAR_RADIUS, 0.0115 * sqrt((MCH_Mass_two_thirds - 1.0/MCH_Mass_two_thirds )));
 }

src/ppowtest

@@ -172,7 +172,7 @@ namespace utils {
         double a_cubed_SI_top    = G * ((p_Mass1 * MSOL_TO_KG) + (p_Mass2 * MSOL_TO_KG)) * p_Period * p_Period * SECONDS_IN_DAY * SECONDS_IN_DAY;
         double a_cubed_SI_bottom = 4.0 * M_PI * M_PI;
         double a_cubed_SI        = a_cubed_SI_top / a_cubed_SI_bottom;
-        double a_SI              = PPOW(a_cubed_SI, 1.0 / 3.0);
+        double a_SI              = std::cbrt(a_cubed_SI); 
 
         return a_SI / AU;
     }