Loading basis-splines.lisp +6 −6 Original line number Diff line number Diff line ;; Basis splines. ;; Liam Healy 2008-02-18 14:43:20EST basis-splines.lisp ;; Time-stamp: <2011-05-26 12:37:29EDT basis-splines.lisp> ;; Time-stamp: <2012-01-13 12:01:11EST basis-splines.lisp> ;; ;; Copyright 2008, 2009, 2011 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading @@ -25,7 +25,7 @@ ;;; Should be subclass of interpolation? (defmobject basis-spline "gsl_bspline" ((order sizet) (number-of-breakpoints sizet)) ((order :sizet) (number-of-breakpoints :sizet)) "basis spline" :gsl-version (1 9) :documentation ; FDL Loading Loading @@ -79,26 +79,26 @@ (defmfun number-of-coefficients (bspline) "gsl_bspline_ncoeffs" (((mpointer bspline) :pointer)) :c-return sizet :c-return :sizet :gsl-version (1 9)) (defmfun order ((bspline basis-spline)) "gsl_bspline_order" (((mpointer bspline) :pointer)) :definition :method :c-return sizet :c-return :sizet :gsl-version (1 9)) (defmfun number-of-breakpoints (bspline) "gsl_bspline_nbreak" (((mpointer bspline) :pointer)) :c-return sizet :c-return :sizet :gsl-version (1 9) :documentation "The number of breakpoints of the basis spline bspline.") (defmfun breakpoint (i bspline) "gsl_bspline_breakpoint" ((i sizet) ((mpointer bspline) :pointer)) ((i :sizet) ((mpointer bspline) :pointer)) :c-return :double :gsl-version (1 9) :documentation "The ith breakpoint of the basis spline bspline.") Loading calculus/monte-carlo-structs.lisp +7 −7 Original line number Diff line number Diff line ;; CFFI-Grovel definitions for unix systems. ;; Liam Healy 2009-06-06 09:32:30EDT monte-carlo-structs.lisp ;; Time-stamp: <2010-05-23 11:36:42EDT monte-carlo-structs.lisp> ;; Time-stamp: <2012-01-13 12:01:40EST monte-carlo-structs.lisp> ;; ;; Copyright 2009, 2010 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading Loading @@ -31,18 +31,18 @@ (include "gsl/gsl_monte_plain.h") (cstruct plain-state "gsl_monte_plain_state" (dim "dim" :type sizet) (dim "dim" :type :sizet) (x "x" :type :pointer)) (include "gsl/gsl_monte_miser.h") (cstruct miser-state "gsl_monte_miser_state" (min-calls "min_calls" :type sizet) (min-calls-per-bisection "min_calls_per_bisection" :type sizet) (min-calls "min_calls" :type :sizet) (min-calls-per-bisection "min_calls_per_bisection" :type :sizet) (dither "dither" :type :double) (estimate-frac "estimate_frac" :type :double) (alpha "alpha" :type :double) (dim "dim" :type sizet) (dim "dim" :type :sizet) (estimate-style "estimate_style" :type :int) (depth "depth" :type :int) (verbose "verbose" :type :int) Loading @@ -65,8 +65,8 @@ (cstruct vegas-state "gsl_monte_vegas_state" ;; grid (dim "dim" :type sizet) (bins-max "bins_max" :type sizet) (dim "dim" :type :sizet) (bins-max "bins_max" :type :sizet) (bins "bins" :type :uint) ; uint (boxes "boxes" :type :uint) ; these are both counted along the axes (xi "xi" :type :pointer) Loading calculus/monte-carlo.lisp +7 −7 Original line number Diff line number Diff line ;; Monte Carlo Integration ;; Liam Healy Sat Feb 3 2007 - 17:42 ;; Time-stamp: <2011-10-30 10:01:57EDT monte-carlo.lisp> ;; Time-stamp: <2012-01-13 12:01:39EST monte-carlo.lisp> ;; ;; Copyright 2007, 2008, 2009, 2011 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading Loading @@ -31,7 +31,7 @@ (defmobject monte-carlo-plain "gsl_monte_plain" ((dim sizet)) ((dim :sizet)) "plain Monte Carlo integration" :documentation ; FDL "Make and initialize a workspace for Monte Carlo integration in dimension dim." Loading @@ -52,7 +52,7 @@ "gsl_monte_plain_integrate" ((callback :pointer) ((grid:foreign-pointer lower-limits) :pointer) ((grid:foreign-pointer upper-limits) :pointer) ((dim0 lower-limits) sizet) (number-of-samples sizet) ((dim0 lower-limits) :sizet) (number-of-samples :sizet) ((mpointer generator) :pointer) ((mpointer state) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) Loading Loading @@ -83,7 +83,7 @@ (defmobject monte-carlo-miser "gsl_monte_miser" ((dim sizet)) ((dim :sizet)) "miser Monte Carlo integration" :documentation ; FDL "Make and initialize a workspace for Monte Carlo integration in Loading Loading @@ -113,7 +113,7 @@ "gsl_monte_miser_integrate" ((callback :pointer) ((grid:foreign-pointer lower-limits) :pointer) ((grid:foreign-pointer upper-limits) :pointer) ((dim0 lower-limits) sizet) (number-of-samples sizet) ((dim0 lower-limits) :sizet) (number-of-samples :sizet) ((mpointer generator) :pointer) ((mpointer state) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) Loading Loading @@ -144,7 +144,7 @@ (defmobject monte-carlo-vegas "gsl_monte_vegas" ((dim sizet)) ((dim :sizet)) "vegas Monte Carlo integration" :documentation ; FDL "Make and initialize a workspace for Monte Carlo integration in Loading Loading @@ -174,7 +174,7 @@ "gsl_monte_vegas_integrate" ((callback :pointer) ((grid:foreign-pointer lower-limits) :pointer) ((grid:foreign-pointer upper-limits) :pointer) ((dim0 lower-limits) sizet) (number-of-samples sizet) ((dim0 lower-limits) :sizet) (number-of-samples :sizet) ((mpointer generator) :pointer) ((mpointer state) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) Loading calculus/numerical-integration-with-tables.lisp +5 −5 Original line number Diff line number Diff line ;; Numerical integration techniques that require tables ;; Liam Healy 2009-04-04 15:24:05EDT ;; Time-stamp: <2011-10-30 10:27:27EDT numerical-integration-with-tables.lisp> ;; Time-stamp: <2012-01-13 12:01:39EST numerical-integration-with-tables.lisp> ;; ;; Copyright 2009, 2011 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading Loading @@ -62,7 +62,7 @@ (a :double) (b :double) ((mpointer table) :pointer) (absolute-error :double) (relative-error :double) (limit sizet) ((mpointer workspace) :pointer) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks (callback (:struct fnstruct) nil (function :double (:input :double) :slug)) Loading @@ -86,7 +86,7 @@ (defmobject qawo-table "gsl_integration_qawo_table" ((omega :double) (L :double) (trig integrate-sine-cosine) (n sizet)) ((omega :double) (L :double) (trig integrate-sine-cosine) (n :sizet)) "table for QAWO numerical integration method" :documentation ; FDL "Make a table describing a sine or cosine weight function W(x) with Loading Loading @@ -119,7 +119,7 @@ ((callback :pointer) (a :double) (absolute-error :double) (relative-error :double) (limit sizet) ((mpointer workspace) :pointer) (limit :sizet) ((mpointer workspace) :pointer) ((mpointer table) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading Loading @@ -162,7 +162,7 @@ ((callback :pointer) (a :double) (absolute-error :double) (limit sizet) (limit :sizet) ((mpointer workspace) :pointer) ((mpointer cycle-workspace) :pointer) ((mpointer table) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) Loading calculus/numerical-integration.lisp +11 −11 Original line number Diff line number Diff line ;; Numerical integration ;; Liam Healy, Wed Jul 5 2006 - 23:14 ;; Time-stamp: <2011-10-30 10:01:11EDT numerical-integration.lisp> ;; Time-stamp: <2012-01-13 12:01:38EST numerical-integration.lisp> ;; ;; Copyright 2006, 2007, 2008, 2009, 2010, 2011 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading Loading @@ -50,7 +50,7 @@ (a :double) (b :double) (absolute-error :double) (relative-error :double) (result (:pointer :double)) (abserr (:pointer :double)) (neval (:pointer sizet))) (neval (:pointer :sizet))) :callbacks (callback (:struct fnstruct) nil (function :double (:input :double) :slug)) :callback-dynamic (nil (function)) Loading @@ -70,7 +70,7 @@ ;;;;**************************************************************************** (defmobject integration-workspace "gsl_integration_workspace" ((size sizet)) "gsl_integration_workspace" ((size :sizet)) "integration workspace" :documentation ; FDL "Make a workspace sufficient to hold n double Loading @@ -89,7 +89,7 @@ ((callback :pointer) (a :double) (b :double) (absolute-error :double) (relative-error :double) (limit sizet) (method integrate-method) ((mpointer workspace) :pointer) (limit :sizet) (method integrate-method) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks (callback (:struct fnstruct) nil (function :double (:input :double) :slug)) Loading Loading @@ -125,7 +125,7 @@ "gsl_integration_qags" ((callback :pointer) (a :double) (b :double) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading Loading @@ -157,8 +157,8 @@ (limit 1000) (workspace (make-integration-workspace limit))) "gsl_integration_qagp" ((callback :pointer) ((grid:foreign-pointer points) :pointer) ((dim0 points) sizet) (absolute-error :double) (relative-error :double) (limit sizet) ((grid:foreign-pointer points) :pointer) ((dim0 points) :sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :inputs (points) Loading Loading @@ -189,7 +189,7 @@ (limit 1000) (workspace (make-integration-workspace limit))) "gsl_integration_qagi" ((callback :pointer) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading @@ -214,7 +214,7 @@ (limit 1000) (workspace (make-integration-workspace limit))) "gsl_integration_qagiu" ((callback :pointer) (a :double) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading @@ -235,7 +235,7 @@ (limit 1000) (workspace (make-integration-workspace limit))) "gsl_integration_qagil" ((callback :pointer) (b :double) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading @@ -261,7 +261,7 @@ "gsl_integration_qawc" ((callback :pointer) (a :double) (b :double) (c :double) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading Loading
basis-splines.lisp +6 −6 Original line number Diff line number Diff line ;; Basis splines. ;; Liam Healy 2008-02-18 14:43:20EST basis-splines.lisp ;; Time-stamp: <2011-05-26 12:37:29EDT basis-splines.lisp> ;; Time-stamp: <2012-01-13 12:01:11EST basis-splines.lisp> ;; ;; Copyright 2008, 2009, 2011 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading @@ -25,7 +25,7 @@ ;;; Should be subclass of interpolation? (defmobject basis-spline "gsl_bspline" ((order sizet) (number-of-breakpoints sizet)) ((order :sizet) (number-of-breakpoints :sizet)) "basis spline" :gsl-version (1 9) :documentation ; FDL Loading Loading @@ -79,26 +79,26 @@ (defmfun number-of-coefficients (bspline) "gsl_bspline_ncoeffs" (((mpointer bspline) :pointer)) :c-return sizet :c-return :sizet :gsl-version (1 9)) (defmfun order ((bspline basis-spline)) "gsl_bspline_order" (((mpointer bspline) :pointer)) :definition :method :c-return sizet :c-return :sizet :gsl-version (1 9)) (defmfun number-of-breakpoints (bspline) "gsl_bspline_nbreak" (((mpointer bspline) :pointer)) :c-return sizet :c-return :sizet :gsl-version (1 9) :documentation "The number of breakpoints of the basis spline bspline.") (defmfun breakpoint (i bspline) "gsl_bspline_breakpoint" ((i sizet) ((mpointer bspline) :pointer)) ((i :sizet) ((mpointer bspline) :pointer)) :c-return :double :gsl-version (1 9) :documentation "The ith breakpoint of the basis spline bspline.") Loading
calculus/monte-carlo-structs.lisp +7 −7 Original line number Diff line number Diff line ;; CFFI-Grovel definitions for unix systems. ;; Liam Healy 2009-06-06 09:32:30EDT monte-carlo-structs.lisp ;; Time-stamp: <2010-05-23 11:36:42EDT monte-carlo-structs.lisp> ;; Time-stamp: <2012-01-13 12:01:40EST monte-carlo-structs.lisp> ;; ;; Copyright 2009, 2010 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading Loading @@ -31,18 +31,18 @@ (include "gsl/gsl_monte_plain.h") (cstruct plain-state "gsl_monte_plain_state" (dim "dim" :type sizet) (dim "dim" :type :sizet) (x "x" :type :pointer)) (include "gsl/gsl_monte_miser.h") (cstruct miser-state "gsl_monte_miser_state" (min-calls "min_calls" :type sizet) (min-calls-per-bisection "min_calls_per_bisection" :type sizet) (min-calls "min_calls" :type :sizet) (min-calls-per-bisection "min_calls_per_bisection" :type :sizet) (dither "dither" :type :double) (estimate-frac "estimate_frac" :type :double) (alpha "alpha" :type :double) (dim "dim" :type sizet) (dim "dim" :type :sizet) (estimate-style "estimate_style" :type :int) (depth "depth" :type :int) (verbose "verbose" :type :int) Loading @@ -65,8 +65,8 @@ (cstruct vegas-state "gsl_monte_vegas_state" ;; grid (dim "dim" :type sizet) (bins-max "bins_max" :type sizet) (dim "dim" :type :sizet) (bins-max "bins_max" :type :sizet) (bins "bins" :type :uint) ; uint (boxes "boxes" :type :uint) ; these are both counted along the axes (xi "xi" :type :pointer) Loading
calculus/monte-carlo.lisp +7 −7 Original line number Diff line number Diff line ;; Monte Carlo Integration ;; Liam Healy Sat Feb 3 2007 - 17:42 ;; Time-stamp: <2011-10-30 10:01:57EDT monte-carlo.lisp> ;; Time-stamp: <2012-01-13 12:01:39EST monte-carlo.lisp> ;; ;; Copyright 2007, 2008, 2009, 2011 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading Loading @@ -31,7 +31,7 @@ (defmobject monte-carlo-plain "gsl_monte_plain" ((dim sizet)) ((dim :sizet)) "plain Monte Carlo integration" :documentation ; FDL "Make and initialize a workspace for Monte Carlo integration in dimension dim." Loading @@ -52,7 +52,7 @@ "gsl_monte_plain_integrate" ((callback :pointer) ((grid:foreign-pointer lower-limits) :pointer) ((grid:foreign-pointer upper-limits) :pointer) ((dim0 lower-limits) sizet) (number-of-samples sizet) ((dim0 lower-limits) :sizet) (number-of-samples :sizet) ((mpointer generator) :pointer) ((mpointer state) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) Loading Loading @@ -83,7 +83,7 @@ (defmobject monte-carlo-miser "gsl_monte_miser" ((dim sizet)) ((dim :sizet)) "miser Monte Carlo integration" :documentation ; FDL "Make and initialize a workspace for Monte Carlo integration in Loading Loading @@ -113,7 +113,7 @@ "gsl_monte_miser_integrate" ((callback :pointer) ((grid:foreign-pointer lower-limits) :pointer) ((grid:foreign-pointer upper-limits) :pointer) ((dim0 lower-limits) sizet) (number-of-samples sizet) ((dim0 lower-limits) :sizet) (number-of-samples :sizet) ((mpointer generator) :pointer) ((mpointer state) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) Loading Loading @@ -144,7 +144,7 @@ (defmobject monte-carlo-vegas "gsl_monte_vegas" ((dim sizet)) ((dim :sizet)) "vegas Monte Carlo integration" :documentation ; FDL "Make and initialize a workspace for Monte Carlo integration in Loading Loading @@ -174,7 +174,7 @@ "gsl_monte_vegas_integrate" ((callback :pointer) ((grid:foreign-pointer lower-limits) :pointer) ((grid:foreign-pointer upper-limits) :pointer) ((dim0 lower-limits) sizet) (number-of-samples sizet) ((dim0 lower-limits) :sizet) (number-of-samples :sizet) ((mpointer generator) :pointer) ((mpointer state) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) Loading
calculus/numerical-integration-with-tables.lisp +5 −5 Original line number Diff line number Diff line ;; Numerical integration techniques that require tables ;; Liam Healy 2009-04-04 15:24:05EDT ;; Time-stamp: <2011-10-30 10:27:27EDT numerical-integration-with-tables.lisp> ;; Time-stamp: <2012-01-13 12:01:39EST numerical-integration-with-tables.lisp> ;; ;; Copyright 2009, 2011 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading Loading @@ -62,7 +62,7 @@ (a :double) (b :double) ((mpointer table) :pointer) (absolute-error :double) (relative-error :double) (limit sizet) ((mpointer workspace) :pointer) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks (callback (:struct fnstruct) nil (function :double (:input :double) :slug)) Loading @@ -86,7 +86,7 @@ (defmobject qawo-table "gsl_integration_qawo_table" ((omega :double) (L :double) (trig integrate-sine-cosine) (n sizet)) ((omega :double) (L :double) (trig integrate-sine-cosine) (n :sizet)) "table for QAWO numerical integration method" :documentation ; FDL "Make a table describing a sine or cosine weight function W(x) with Loading Loading @@ -119,7 +119,7 @@ ((callback :pointer) (a :double) (absolute-error :double) (relative-error :double) (limit sizet) ((mpointer workspace) :pointer) (limit :sizet) ((mpointer workspace) :pointer) ((mpointer table) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading Loading @@ -162,7 +162,7 @@ ((callback :pointer) (a :double) (absolute-error :double) (limit sizet) (limit :sizet) ((mpointer workspace) :pointer) ((mpointer cycle-workspace) :pointer) ((mpointer table) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) Loading
calculus/numerical-integration.lisp +11 −11 Original line number Diff line number Diff line ;; Numerical integration ;; Liam Healy, Wed Jul 5 2006 - 23:14 ;; Time-stamp: <2011-10-30 10:01:11EDT numerical-integration.lisp> ;; Time-stamp: <2012-01-13 12:01:38EST numerical-integration.lisp> ;; ;; Copyright 2006, 2007, 2008, 2009, 2010, 2011 Liam M. Healy ;; Distributed under the terms of the GNU General Public License Loading Loading @@ -50,7 +50,7 @@ (a :double) (b :double) (absolute-error :double) (relative-error :double) (result (:pointer :double)) (abserr (:pointer :double)) (neval (:pointer sizet))) (neval (:pointer :sizet))) :callbacks (callback (:struct fnstruct) nil (function :double (:input :double) :slug)) :callback-dynamic (nil (function)) Loading @@ -70,7 +70,7 @@ ;;;;**************************************************************************** (defmobject integration-workspace "gsl_integration_workspace" ((size sizet)) "gsl_integration_workspace" ((size :sizet)) "integration workspace" :documentation ; FDL "Make a workspace sufficient to hold n double Loading @@ -89,7 +89,7 @@ ((callback :pointer) (a :double) (b :double) (absolute-error :double) (relative-error :double) (limit sizet) (method integrate-method) ((mpointer workspace) :pointer) (limit :sizet) (method integrate-method) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks (callback (:struct fnstruct) nil (function :double (:input :double) :slug)) Loading Loading @@ -125,7 +125,7 @@ "gsl_integration_qags" ((callback :pointer) (a :double) (b :double) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading Loading @@ -157,8 +157,8 @@ (limit 1000) (workspace (make-integration-workspace limit))) "gsl_integration_qagp" ((callback :pointer) ((grid:foreign-pointer points) :pointer) ((dim0 points) sizet) (absolute-error :double) (relative-error :double) (limit sizet) ((grid:foreign-pointer points) :pointer) ((dim0 points) :sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :inputs (points) Loading Loading @@ -189,7 +189,7 @@ (limit 1000) (workspace (make-integration-workspace limit))) "gsl_integration_qagi" ((callback :pointer) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading @@ -214,7 +214,7 @@ (limit 1000) (workspace (make-integration-workspace limit))) "gsl_integration_qagiu" ((callback :pointer) (a :double) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading @@ -235,7 +235,7 @@ (limit 1000) (workspace (make-integration-workspace limit))) "gsl_integration_qagil" ((callback :pointer) (b :double) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading @@ -261,7 +261,7 @@ "gsl_integration_qawc" ((callback :pointer) (a :double) (b :double) (c :double) (absolute-error :double) (relative-error :double) (limit sizet) (absolute-error :double) (relative-error :double) (limit :sizet) ((mpointer workspace) :pointer) (result (:pointer :double)) (abserr (:pointer :double))) :callbacks Loading