Commit 9c1abb95 by Liam M. Healy

Revision dates

parent 9c043926
;; Numerical integration
;; Liam Healy, Wed Jul 5 2006 - 23:14
;; Time-stamp: <2012-01-13 12:01:38EST numerical-integration.lisp>
;; Time-stamp: <2015-01-22 23:32:29EST numerical-integration.lisp>
;;
;; Copyright 2006, 2007, 2008, 2009, 2010, 2011 Liam M. Healy
;; Copyright 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2015 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -370,6 +370,8 @@
(grid:copy-to (vector 0.0d0 1.0d0 (sqrt 2.0d0) 3.0d0))
0.0d0 1.0d-3 1000)
(integration-QAWc 'integration-test-f459 -1.0d0 5.0d0 0.0d0 0.0d0 1.0d-3 1000)
;; Check that integration-QAG is reentrant, but note this is not the
;; recommended way to do multivariate integration (see Monte Carlo).
(integration-QAG (lambda (x)
(integration-QAG (lambda (y) (* (sin x) y)) 0d0 1d0 :gauss41))
0d0 pi :gauss41))
(integration-QAG (lambda (y) (* (sin x) y)) 0d0 1d0 :gauss41))
0d0 pi :gauss41))
;; Foreign callback functions.
;; Liam Healy
;; Time-stamp: <2010-07-11 19:01:17EDT callback.lisp>
;; Time-stamp: <2015-01-22 23:27:55EST callback.lisp>
;;
;; Copyright 2009 Liam M. Healy
;; Copyright 2009, 2010, 2015 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......
;; Regression test NUMERICAL-INTEGRATION for GSLL, automatically generated
;;
;; Copyright 2009, 2010 Liam M. Healy
;; Copyright 2009, 2010, 2015 Liam M. Healy
;; Distributed under the terms of the GNU General Public License
;;
;; This program is free software: you can redistribute it and/or modify
......@@ -233,11 +233,12 @@
(MULTIPLE-VALUE-LIST
(INTEGRATION-QAWC 'INTEGRATION-TEST-F459 -1.0d0 5.0d0
0.0d0 0.0d0 0.001d0 1000))))
(LISP-UNIT:ASSERT-NUMERICAL-EQUAL (LIST 1.0d0 1.1102230246251565d-14)
(MULTIPLE-VALUE-LIST
(INTEGRATION-QAG
(LAMBDA (X)
(INTEGRATION-QAG
(LAMBDA (Y) (* (SIN X) Y)) 0.0d0 1.0d0
:GAUSS41))
0.0d0 PI :GAUSS41))))
(LISP-UNIT:ASSERT-NUMERICAL-EQUAL
(LIST 1.0d0 1.1102230246251565d-14)
(MULTIPLE-VALUE-LIST
(INTEGRATION-QAG
(LAMBDA (X)
(INTEGRATION-QAG
(LAMBDA (Y) (* (SIN X) Y)) 0.0d0 1.0d0
:GAUSS41))
0.0d0 PI :GAUSS41))))
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment