implement more power functions

to avoid std::pow

meta/CConversion.m

@@ -153,7 +153,8 @@
   https://stackoverflow.com/questions/3136604/evaluate-beyond-one-level-within-hold-in-mathematica
 ";
 
-{ Sqr, Cube, Quad, Power2, Power3, Power4, Power5, Power6, Power7, Power8 };
+{ Sqr, Cube, Quad, Cbrt, Power2, Power3, Power4, Power5, Power6,
+  Power7, Power8, Power9, Power10, Power11, Power12 };
 
 Begin["`Private`"];
 
@@ -684,6 +685,12 @@
 Unprotect[Power];
 Format[Power[E,z_],CForm] :=
     Format["Exp(" <> ToString[CForm[z]] <> ")", OutputForm];
+Format[Power[b_,2],CForm] :=
+    Format["Sqr(" <> ToString[CForm[b]] <> ")", OutputForm];
+Format[Power[b_,0.5 | 1/2],CForm] :=
+    Format["Sqrt(" <> ToString[CForm[b]] <> ")", OutputForm];
+Format[Power[b_,1/3],CForm] :=
+    Format["Cbrt(" <> ToString[CForm[b]] <> ")", OutputForm];
 Protect[Power];
 
 Unprotect[If];
@@ -884,6 +891,14 @@
                     Power[a_,-7]             :> 1/Power7[a] /.
                     Power[a_,8]              :> Power8[a] /.
                     Power[a_,-8]             :> 1/Power8[a] /.
+                    Power[a_,9]              :> Power9[a] /.
+                    Power[a_,-9]             :> 1/Power9[a] /.
+                    Power[a_,10]             :> Power10[a] /.
+                    Power[a_,-10]            :> 1/Power10[a] /.
+                    Power[a_,11]             :> Power11[a] /.
+                    Power[a_,-11]            :> 1/Power11[a] /.
+                    Power[a_,12]             :> Power12[a] /.
+                    Power[a_,-12]            :> 1/Power12[a] /.
                     Sqrt[x_]/Sqrt[y_]        :> Sqrt[x/y];
            result = Apply[Function[code, Hold[CForm[code]], HoldAll],
                           Hold[#] &[result /. { (p:(Dot | SARAH`MatMul))[a__] :> times @@ p[a],

src/wrappers.hpp

@@ -522,19 +522,19 @@ Base Power(Base base, Exponent exp) noexcept
 template <typename Base>
 constexpr Base Power2(Base b) noexcept
 {
-   return Sqr(b);
+   return b * b;
 }
 
 template <typename Base>
 constexpr Base Power3(Base b) noexcept
 {
-   return Cube(b);
+   return b * b * b;
 }
 
 template <typename Base>
 constexpr Base Power4(Base b) noexcept
 {
-   return Quad(b);
+   return b * b * b * b;
 }
 
 template <typename Base>
@@ -552,13 +552,45 @@ constexpr Base Power6(Base b) noexcept
 template <typename Base>
 constexpr Base Power7(Base b) noexcept
 {
-   return b * b * b * b * b * b * b;
+   return b * b * b * b * b *
+          b * b;
 }
 
 template <typename Base>
 constexpr Base Power8(Base b) noexcept
 {
-   return b * b * b * b * b * b * b * b;
+   return b * b * b * b * b *
+          b * b * b;
+}
+
+template <typename Base>
+constexpr Base Power9(Base b) noexcept
+{
+   return b * b * b * b * b *
+          b * b * b * b;
+}
+
+template <typename Base>
+constexpr Base Power10(Base b) noexcept
+{
+   return b * b * b * b * b *
+          b * b * b * b * b;
+}
+
+template <typename Base>
+constexpr Base Power11(Base b) noexcept
+{
+   return b * b * b * b * b *
+          b * b * b * b * b *
+          b;
+}
+
+template <typename Base>
+constexpr Base Power12(Base b) noexcept
+{
+   return b * b * b * b * b *
+          b * b * b * b * b *
+          b * b;
 }
 
 inline constexpr double Re(double x) noexcept