diff options
author | Ian Jauslin <ian.jauslin@roma1.infn.it> | 2016-05-24 13:39:23 +0000 |
---|---|---|
committer | Ian Jauslin <ian.jauslin@roma1.infn.it> | 2016-05-24 13:39:23 +0000 |
commit | fa9b6f2b9bcb80778e63ef2aa4e17c7573de0015 (patch) | |
tree | 92b740d0736c9ed6f5bda051c224c8bb7196bb03 /src/zz_integral.h |
Diffstat (limited to 'src/zz_integral.h')
-rw-r--r-- | src/zz_integral.h | 53 |
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 + |