Thanks! This makes the code much leaner. I thought about it too but I was unsure whether passing actual values would interfere somehow—thank you for suggesting this!
EDIT: I can’t find the default value for 1e-8 according to Computing One Dimensional Integralsrelative_tolerance
. Is there any recommendation on what value to set?
EDIT2: This is the part in the User Guide that discouraged me to pass actual values to x_r, x_i
This may require passing in zero-length arrays for data or a zero-length vector for parameters if the integral does not involve data or parameters.