updating fortran wrapper

scarrazza · Feb 16, 2021 · a3c5b45 · a3c5b45
1 parent bea5af4
commit a3c5b45
Show file tree

Hide file tree

Showing 2 changed files with 74 additions and 1 deletion.
diff --git a/src/qcdloop/wrapper.h b/src/qcdloop/wrapper.h
@@ -55,6 +55,11 @@ extern"C" {
   qcomplex qli2q_(qdouble const& p1, qdouble const& m1, qdouble const& m2, qdouble const& mu2, int const& ep);
   qcomplex qli2qc_(qdouble const& p1, qcomplex const& m1, qcomplex const& m2, qdouble const& mu2, int const& ep);
 
+  complex  qli2p_(double const& p1, double const& m1, double const& m2, double const& mu2, int const& ep);
+  complex  qli2pc_(double const& p1, complex const& m1, complex const& m2, double const& mu2, int const& ep);
+  qcomplex qli2pq_(qdouble const& p1, qdouble const& m1, qdouble const& m2, qdouble const& mu2, int const& ep);
+  qcomplex qli2pqc_(qdouble const& p1, qcomplex const& m1, qcomplex const& m2, qdouble const& mu2, int const& ep);
+
   complex  qli3_(double const& p1, double const& p2, double const& p3, double const& m1, double const& m2, double const& m3, double const& mu2, int const& ep);
   complex  qli3c_(double const& p1, double const& p2, double const& p3, complex const& m1, complex const& m2, complex const& m3, double const& mu2, int const& ep);
   qcomplex qli3q_(qdouble const& p1, qdouble const& p2, qdouble const& p3, qdouble const& m1, qdouble const& m2, qdouble const& m3, qdouble const& mu2, int const& ep);

diff --git a/src/wrapper.cc b/src/wrapper.cc
@@ -533,7 +533,7 @@ void qlboxcq_(qcomplex (&out)[3], qdouble const& mu2, qcomplex const& m1, qcompl
     {
       std::cout << e.what() << std::endl;
       exit(-1);
-    }    
+    }
     return r[abs(ep)];
   }
 
@@ -588,6 +588,74 @@ void qlboxcq_(qcomplex (&out)[3], qdouble const& mu2, qcomplex const& m1, qcompl
     return rq[abs(ep)];
   }
 
+  complex qli2p_(double const& p1, double const& m1, double const& m2, double const& mu2, int const& ep)
+  {
+    try
+    {
+      mI2[0] = m1;
+      mI2[1] = m2;
+      pI2[0] = p1;
+      bb.derivative(r, mu2, mI2, pI2);
+    }
+    catch (std::exception &e)
+    {
+      std::cout << e.what() << std::endl;
+      exit(-1);
+    }
+    return r[abs(ep)];
+  }
+
+  complex qli2pc_(double const& p1, complex const& m1, complex const& m2, double const& mu2, int const& ep)
+  {
+    try
+    {
+      mI2c[0] = m1;
+      mI2c[1] = m2;
+      pI2[0] = p1;
+      bbc.derivative(r, mu2, mI2c, pI2);
+    }
+    catch (std::exception &e)
+    {
+      std::cout << e.what() << std::endl;
+      exit(-1);
+    }
+    return r[abs(ep)];
+  }
+
+  qcomplex qli2pq_(qdouble const& p1, qdouble const& m1, qdouble const& m2, qdouble const& mu2, int const& ep)
+  {
+    try
+    {
+      mI2q[0] = m1;
+      mI2q[1] = m2;
+      pI2q[0] = p1;
+      bbq.derivative(rq, mu2, mI2q, pI2q);
+    }
+    catch (std::exception &e)
+    {
+      std::cout << e.what() << std::endl;
+      exit(-1);
+    }
+    return rq[abs(ep)];
+  }
+
+  qcomplex qli2pqc_(qdouble const& p1, qcomplex const& m1, qcomplex const& m2, qdouble const& mu2, int const& ep)
+  {
+    try
+    {
+      mI2cq[0] = m1;
+      mI2cq[1] = m2;
+      pI2q[0] = p1;
+      bbcq.derivative(rq, mu2, mI2cq, pI2q);
+    }
+    catch (std::exception &e)
+    {
+      std::cout << e.what() << std::endl;
+      exit(-1);
+    }
+    return rq[abs(ep)];
+  }
+
   complex qli3_(double const& p1, double const& p2, double const& p3, double const& m1, double const& m2, double const& m3, double const& mu2, int const& ep)
   {
     try