Thank you very much! It seems that grid (antik) is automatically loaded with gsll by quicklisp. But that doesn't seem to be the problem: both
(grid:copy-to (vector 0.0d0 1.0d0 (sqrt 2.0d0) 3.0d0)) and (vector 0.0d0 1.0d0 (sqrt 2.0d0) 3.0d0) seem to give the same result - which should be the endpoints and singularities for computation of a numerical integral by use of the GSL QUADPACK routine QAGP, implemented in gsll as integration-qagp. But the command (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)) fails with a slew of errors. The debugger's first backtrace is ((:METHOD NO-APPLICABLE-METHOD (T)) #<STANDARD-GENERIC-FUNCTION GRID:FOREIGN-POINTER (1)> #(0.0 1.0 1.4142135623730951 3.0)) [fast-method] Given that the grid command above works without any errors, I don't understand this at all. But then, as I say, I seem to have programmed in many languages over the years - with the lamentable absence of lisp. So if my questions seem trivial and naive, my apologies! Thanks again, Alasdair On Mon, Nov 2, 2015 at 12:15 AM, David Catteeuw <davidcatte...@gmail.com> wrote: > Hi Alasdair, > > I guess grid refers to Antik's grids, a data type for matrices: > https://www.common-lisp.net/project/antik/ > <https://www.common-lisp.net/project/antik/#documentation> > > Regards, > david. > > <https://www.common-lisp.net/project/antik/#documentation> > > On Sunday, 1 November 2015, Alasdair McAndrew <amc...@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 >> -- >> [image: http://www.facebook.com/alasdair.mcandrew] >> <http://www.facebook.com/alasdair.mcandrew> [image: >> https://plus.google.com/+AlasdairMcAndrew/posts] >> <https://plus.google.com/+AlasdairMcAndrew/posts> [image: >> https://www.linkedin.com/pub/alasdair-mcandrew/a/178/108] >> <https://www.linkedin.com/pub/alasdair-mcandrew/a/178/108> [image: >> https://twitter.com/amca01] <https://twitter.com/amca01> [image: >> http://numbersandshapes.net] <http://numbersandshapes.net> >> > -- [image: http://www.facebook.com/alasdair.mcandrew] <http://www.facebook.com/alasdair.mcandrew> [image: https://plus.google.com/+AlasdairMcAndrew/posts] <https://plus.google.com/+AlasdairMcAndrew/posts> [image: https://www.linkedin.com/pub/alasdair-mcandrew/a/178/108] <https://www.linkedin.com/pub/alasdair-mcandrew/a/178/108> [image: https://twitter.com/amca01] <https://twitter.com/amca01> [image: http://numbersandshapes.net] <http://numbersandshapes.net>
