add wrapper for Mathematica's PolyLog[2,_] function

FlexibleSUSY · Sep 10, 2015 · a052ef5 · a052ef5
1 parent 69ef10a
commit a052ef5
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/src/wrappers.hpp b/src/wrappers.hpp
@@ -29,6 +29,8 @@
 #include <Eigen/Core>
 #include <boost/lexical_cast.hpp>
 
+#include "dilog.hpp"
+
 namespace flexiblesusy {
 
 static const double Pi = M_PI;
@@ -291,6 +293,13 @@ inline int Sign(int x)
    return (x >= 0 ? 1 : -1);
 }
 
+template <typename T>
+T PolyLog(int n, T z) {
+   if (n == 2)
+      return dilog(z);
+   assert(false && "PolyLog(n!=2) not implemented");
+}
+
 template <typename Base, typename Exponent>
 Base Power(Base base, Exponent exp)
 {