scripts/general/triplequad.m

@(
)
@(
)
@(
)
@(
)
@(
)
@(
)
author	Rik <octave@nomad.inbox5.com>
date	Sun, 17 Apr 2011 19:57:07 -0700
parents	d0b799dafede
children	0713ad019e53 f96b9b9f141b
rev	line source
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 9245 diff changeset	1 ## Copyright (C) 2008-2011 David Bateman
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	2 ##
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	3 ## This file is part of Octave.
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	4 ##
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	5 ## Octave is free software; you can redistribute it and/or modify it
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	6 ## under the terms of the GNU General Public License as published by
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	7 ## the Free Software Foundation; either version 3 of the License, or (at
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	8 ## your option) any later version.
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	9 ##
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	10 ## Octave is distributed in the hope that it will be useful, but
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	13 ## General Public License for more details.
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	14 ##
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	15 ## You should have received a copy of the GNU General Public License
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	16 ## along with Octave; see the file COPYING. If not, see
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	17 ## <http://www.gnu.org/licenses/>.
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	18
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	19 ## -- texinfo --
12612 16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	20 ## @deftypefn {Function File} {} triplequad (@var{f}, @var{xa}, @var{xb}, @var{ya}, @var{yb}, @var{za}, @var{zb})
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	21 ## @deftypefnx {Function File} {} triplequad (@var{f}, @var{xa}, @var{xb}, @var{ya}, @var{yb}, @var{za}, @var{zb}, @var{tol})
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	22 ## @deftypefnx {Function File} {} triplequad (@var{f}, @var{xa}, @var{xb}, @var{ya}, @var{yb}, @var{za}, @var{zb}, @var{tol}, @var{quadf})
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	23 ## @deftypefnx {Function File} {} triplequad (@var{f}, @var{xa}, @var{xb}, @var{ya}, @var{yb}, @var{za}, @var{zb}, @var{tol}, @var{quadf}, @dots{})
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	24 ## Numerically evaluate the triple integral of @var{f}.
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	25 ## @var{f} is a function handle, inline function, or string
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	26 ## containing the name of the function to evaluate. The function @var{f} must
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	27 ## have the form @math{w = f(x,y,z)} where either @var{x} or @var{y} is a
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	28 ## vector and the remaining inputs are scalars. It should return a vector of
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	29 ## the same length and orientation as @var{x} or @var{y}.
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	30 ##
12612 16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	31 ## @var{xa}, @var{ya}, @var{za} and @var{xb}, @var{yb}, @var{zb} are the lower
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	32 ## and upper limits of integration for x, y, and z respectively. The
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	33 ## underlying integrator determines whether infinite bounds are accepted.
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	34 ##
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	35 ## The optional argument @var{tol} defines the absolute tolerance used to
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	36 ## integrate each sub-integral. The default value is @math{1e^{-6}}.
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	37 ##
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	38 ## The optional argument @var{quadf} specifies which underlying integrator
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	39 ## function to use. Any choice but @code{quad} is available and the default
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	40 ## is @code{quadgk}.
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	41 ##
9051 1bf0ce0930be Grammar check TexInfo in all .m files Rik <rdrider0-list@yahoo.com> parents: 7771 diff changeset	42 ## Additional arguments, are passed directly to @var{f}. To use the default
12612 16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	43 ## value for @var{tol} or @var{quadf} one may pass an empty matrix ([]).
12575 d0b799dafede Grammarcheck files for 3.4.1 release. Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	44 ## @seealso{dblquad, quad, quadv, quadl, quadgk, quadcc, trapz}
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	45 ## @end deftypefn
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	46
12612 16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	47 function q = triplequad(f, xa, xb, ya, yb, za, zb, tol, quadf, varargin)
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	48 if (nargin < 7)
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	49 print_usage ();
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	50 endif
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	51 if (nargin < 8 \|\| isempty (tol))
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11562 diff changeset	52 tol = 1e-6;
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	53 endif
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	54 if (nargin < 9 \|\| isempty (quadf))
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11562 diff changeset	55 quadf = @quadgk;
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	56 endif
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	57
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	58 inner = @__triplequad_inner__;
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	59 if (ischar (f))
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	60 f = @(x,y,z) feval (f, x, y, z, varargin{:});
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	61 varargin = {};
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	62 endif
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	63
12612 16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	64 q = dblquad(@(y, z) inner (y, z, f, xa, xb, tol, quadf, varargin{:}),ya, yb, za, zb, tol);
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	65 endfunction
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	66
12612 16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	67 function q = __triplequad_inner__ (y, z, f, xa, xb, tol, quadf, varargin)
16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	68 q = zeros (size(y));
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	69 for i = 1 : length (y)
12612 16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	70 q(i) = feval (quadf, @(x) f (x, y(i), z, varargin{:}), xa, xb, tol);
7771 680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	71 endfor
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	72 endfunction
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	73
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	74 %% These tests are too expensive to run normally. Disable them
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	75 % !#assert (triplequad (@(x,y,z) exp(-x.^2 - y.^2 - z.^2) , -1, 1, -1, 1, -1, 1, [], @quadgk), pi ^ (3/2) * erf(1).^3, 1e-6)
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	76 % !#assert (triplequad (@(x,y,z) exp(-x.^2 - y.^2 - z.^2) , -1, 1, -1, 1, -1, 1, [], @quadl), pi ^ (3/2) * erf(1).^3, 1e-6)
680631e787aa Add quadv, quadgk, dblquad and triplequad functions David Bateman <dbateman@free.fr> parents: diff changeset	77 % !#assert (triplequad (@(x,y,z) exp(-x.^2 - y.^2 - z.^2) , -1, 1, -1, 1, -1, 1, [], @quadv), pi ^ (3/2) * erf(1).^3, 1e-6)
12612 16cca721117b doc: Update all documentation for chapter on Numerical Integration Rik <octave@nomad.inbox5.com> parents: 12575 diff changeset	78