1 files changed, 53 insertions, 0 deletions
diff --git a/src/zz_integral.h b/src/zz_integral.h
new file mode 100644
index 0000000..b23c88a
--- /dev/null
+++ b/src/zz_integral.h
@@ -0,0 +1,53 @@
+/*
+Copyright 2016 Ian Jauslin
+
+Licensed under the Apache License, Version 2.0 (the "License");
+you may not use this file except in compliance with the License.
+You may obtain a copy of the License at
+
+    http://www.apache.org/licenses/LICENSE-2.0
+
+Unless required by applicable law or agreed to in writing, software
+distributed under the License is distributed on an "AS IS" BASIS,
+WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+See the License for the specific language governing permissions and
+limitations under the License.
+*/
+
+/*
+  Computation of z1-z2 and z1+z2
+*/
+
+#ifndef ZZ_INTEGRAL_H
+#define ZZ_INTEGRAL_H
+
+#include "ss_integral.h"
+#include <libinum.h>
+#include <mpfr.h>
+
+// I for z1-z2
+int zz_I(mpfr_t out, mpfr_t p1, mpfr_t p2, mpfr_t q1, mpfr_t q2, mpfr_t F1, mpfr_t F2, mpfr_t t1, mpfr_t t2, mpfr_t phi, mpfr_t W, array_mpfr* tmps1, struct ss_cache cache);
+
+// I for z1+z2
+int ZZ_I(mpfr_t out, mpfr_t p1, mpfr_t p2, mpfr_t q1, mpfr_t q2, mpfr_t F1, mpfr_t F2, mpfr_t t1, mpfr_t t2, mpfr_t phi, mpfr_t W, array_mpfr* tmps1, struct ss_cache cache);
+
+// zeta(k1,k2)
+int zz_zeta(mpfr_t zeta, mpfr_t alpha1, mpfr_t t2, mpfr_t phi);
+
+// xi(k1,k2)
+int zz_xi(mpfr_t xi, mpfr_t m, mpfr_t alpha1, mpfr_t t1, mpfr_t tmp);
+
+// m(k1,k2)
+int zz_m(mpfr_t m, mpfr_t k1, mpfr_t k2, mpfr_t t2, mpfr_t phi, mpfr_t W, struct ss_cache cache, mpfr_t tmp1);
+
+// m(k1,k2)/xi(k1,k2)
+int zz_mxi(mpfr_t mxi, mpfr_t m, mpfr_t xi);
+
+// alpha1(k1,k2)
+int zz_alpha1(mpfr_t alpha1, mpfr_t k1, mpfr_t k2, struct ss_cache cache, mpfr_t tmp1);
+
+// alpha2(k1,k2)
+int zz_alpha2(mpfr_t alpha2, mpfr_t k1, mpfr_t k2, struct ss_cache cache, mpfr_t tmp1);
+
+#endif
+