#!/usr/bin/env python3
from pySecDec.integral_interface import IntegralLibrary
import sympy as sp
if __name__ == "__main__":
# load c++ library
box1L = IntegralLibrary('box1L/box1L_pylink.so')
# choose integrator
box1L.use_Vegas()
# integrate
str_integral_without_prefactor, str_prefactor, str_integral_with_prefactor = box1L(
real_parameters=[4.0, -0.75, 1.25, 1.0], verbose=True)
# convert complex numbers from c++ to sympy notation
str_integral_with_prefactor = str_integral_with_prefactor.replace(
',', '+I*')
str_prefactor = str_prefactor.replace(',', '+I*')
str_integral_without_prefactor = str_integral_without_prefactor.replace(
',', '+I*')
# convert result to sympy expressions
integral_with_prefactor = sp.sympify(
str_integral_with_prefactor.replace('+/-', '*value+error*'))
integral_with_prefactor_err = sp.sympify(
str_integral_with_prefactor.replace('+/-', '*value+error*'))
prefactor = sp.sympify(str_prefactor)
integral_without_prefactor = sp.sympify(
str_integral_without_prefactor.replace('+/-', '*value+error*'))
#!/usr/bin/env python3
from pySecDec.integral_interface import IntegralLibrary
import sympy as sp
if __name__ == "__main__":
# load c++ library
triangle2L = IntegralLibrary('triangle2L/triangle2L_pylink.so')
# choose integrator
triangle2L.use_Vegas() # ``flags=2``: verbose --> see Cuba manual
number_of_real_parameters = int(triangle2L.info['number_of_real_parameters'])
number_of_complex_parameters = int(triangle2L.info['number_of_complex_parameters'])
with open('kinematics.input') as f:
with open('results_P126.txt', 'w') as resultsfile:
for line in f:
point = line.split()
assert len(point) == 1+number_of_real_parameters+number_of_complex_parameters, "Invalid point: " + str(point)
# convert to float and complex
name = point[0]
vals_real_parameters = [float(point[1+i]) for i in range(number_of_real_parameters)]
vals_complex_parameters = [complex(point[1+number_of_real_parameters+i]) for i in range(number_of_complex_parameters)]
# compute the integral
str_integral_without_prefactor, str_prefactor, str_integral_with_prefactor = triangle2L(vals_real_parameters,vals_complex_parameters, verbose=True)
