A question about numerical-integration.lisp

[Liam Healy]

On Sun, Nov 1, 2015 at 5:10 AM, Alasdair McAndrew <amca01 at gmail.com> wrote:

> Hello,
>
> I am using the excellent gsll package (in the first instance), to provide
> an interface to quadpack for the mathematics system FriCAS.  I'm slowly
> going through calculus/numerical-integration.lisp one function at a time
> and writing each one into FriCAS as I go.  This means writing functions in
> FriCAS's own language SPAD which interface with gsll.
>
> This is made harder by me being a lisp newbie.
>
> However, I've come to integration-QAGP, and I've tried to run the test
> command (having first defined the function it calls).  However, all I get
> is errors.
>
> I'm using SBCL in emacs-slime, and I have installed gsll with quicklisp,
> so that I can call an integration routine with
>
> * (gsll:integration-qng (lambda (x) (exp (- (* x x)))) 0.0 1.0)
>
> But the commands
>
> * (defun integration-test-f454 (x)
>   (* (expt x 3) (* (log (abs (* (- (expt x 2) 1.0d0) (- (expt x 2)
> 2.0d0)))))))
>
> * (gsll:integration-QAGp
>   'integration-test-f454
>   (grid:copy-to (vector 0.0d0 1.0d0 (sqrt 2.0d0) 3.0d0))
>   0.0d0 1.0d-3 1000)
>
> just produces a long list of errors. Is "grid" a standard library, or does
> it need to be loaded first?  (Told you I was a newbie...)
>
> Thanks
> Alasdair
> --
>

Try this

(setf grid:*default-grid-type* 'grid:foreign-array)
(gsll:integration-QAGp 'integration-test-f454 (grid:grid 0.0d0 1.0d0 (sqrt
2.0d0) 3.0d0) 0.0d0 1.0d-3 1000)

and let us know what happens.

Thanks,
Liam
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.common-lisp.net/pipermail/gsll-devel/attachments/20151101/c31f0517/attachment.html>