blob: 750530c4c2b9f347468db0390fec2dc647d11726 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
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)))))
|
