Agena

3.11.5 Converse, March 24, 2024 Tuned, extended and fixed calculus functions that compute integrals: - `calc.intdei` and `calc.intdeo` now remember the last epsilon setting given by the user and run 40 % faster with succeeding calls done with the same epsilon value. This also benefits the wrapper function `calc.integ` (see below). - `calc.gauleg`, `calc.intde`, calc.intdei`, `calc.intdeo` and `calc.intcc` now accept multivariate functions. See the Primer and Reference on how to pass the second, third, etc. arguments in the respective calls. With univariate functions, as before, you do not have to change the code. - `calc.integ` has been rewritten to support multivariate functions. Optional epsilon, omega and sample values must now be given as explicit options, so you may have to adapt your code. Example: > calc.integ(<< x, a -> x + a >>, 1, 2, 0, eps=hEps, omega=1, samples=100): 1.5 2.999999997239e-015 - If the third argument to `calc.intdeo` was zero, the interpreter could crash on some platforms. This has been fixed. - If the number of sample points is zero, `calc.gauleg` now automatically sets it to 20 points internally. This prevents segmentation faults on some platforms and also false results of exactly zero. - If you pass zero as an epsilon value to `calc.intde`, calc.intdei`, `calc.intdeo` and `calc.intcc`, then it is automatically reset to the default setting 1e-15 to avert any infinite loops. - This release has been Valgrind-checked on x86 and AMD64 Linux to ensure there are no internal errors or memory leaks.