diff options
| author | Arun Isaac | 2021-02-23 16:24:48 +0530 |
|---|---|---|
| committer | Arun Isaac | 2021-02-23 16:24:48 +0530 |
| commit | 25343d3723ae70690f056b0a501f5941892fbf93 (patch) | |
| tree | bc7d465ba091587438d71ad5c258f195448aeda2 | |
| parent | 59c0f5d621d665a2ab48cc760b561a7b610ab5fd (diff) | |
| download | nsmc-25343d3723ae70690f056b0a501f5941892fbf93.tar.gz nsmc-25343d3723ae70690f056b0a501f5941892fbf93.tar.lz nsmc-25343d3723ae70690f056b0a501f5941892fbf93.zip | |
Provide test integrals.
* include/integrands.h, src/integrands.sc: New files.
| -rw-r--r-- | include/integrands.h | 15 | ||||
| -rw-r--r-- | src/integrands.sc | 19 |
2 files changed, 34 insertions, 0 deletions
diff --git a/include/integrands.h b/include/integrands.h new file mode 100644 index 0000000..b4735a8 --- /dev/null +++ b/include/integrands.h @@ -0,0 +1,15 @@ +#ifndef INTEGRANDS_H +#define INTEGRANDS_H + +#include <gsl/gsl_vector.h> + +typedef struct { + double *coefficients; + int degree; +} polynomial_integrand_params; + +double polynomial_integrand (double r, const gsl_vector* x, void *params); +double gaussian_integrand (double r, const gsl_vector* x, void *params); +double x_coordinate_integrand (double r, const gsl_vector* x, void *params); + +#endif diff --git a/src/integrands.sc b/src/integrands.sc new file mode 100644 index 0000000..750530c --- /dev/null +++ b/src/integrands.sc @@ -0,0 +1,19 @@ +(pre-include "gsl/gsl_math.h") +(pre-include "gsl/gsl_poly.h") +(pre-include "integrands.h") + +(define (polynomial-integrand r x -params) (double double (const gsl-vector*) void*) + (let* ((params polynomial-integrand-params* + (convert-type -params polynomial-integrand-params*))) + (return (gsl-poly-eval (: params coefficients) + (+ (: params degree) 1) + r)))) + +(define (gaussian-integrand r x -params) (double double (const gsl-vector*) void*) + ;; This is a workaround to fix negation bug in sph-sc. See + ;; https://github.com/sph-mn/sph-sc/issues/5 . TODO: Remove + ;; workaround once bug is fixed upstream. + (return (exp (- 0 (/ (gsl-pow-2 r) 2))))) + +(define (x-coordinate-integrand r x -params) (double double (const gsl-vector*) void*) + (return (fabs (* r (gsl-vector-get x 0))))) |