Mercurial > hg > octave-nkf
annotate liboctave/numeric/Quad.h @ 17822:ebb3ef964372
maint: Use common #define syntax "octave_filename_h" in h_files.
* libgui/src/dialog.h, libgui/src/find-files-dialog.h,
libgui/src/find-files-model.h, libgui/src/m-editor/file-editor-tab.h,
libgui/src/m-editor/octave-qscintilla.h, libgui/src/main-window.h,
libgui/src/octave-dock-widget.h, libgui/src/workspace-model.h,
libgui/src/workspace-view.h, libinterp/corefcn/Cell.h,
libinterp/corefcn/gl-render.h, libinterp/corefcn/gl2ps-renderer.h,
libinterp/corefcn/graphics.in.h, libinterp/corefcn/oct-errno.h,
libinterp/corefcn/oct-fstrm.h, libinterp/corefcn/oct-handle.h,
libinterp/corefcn/oct-hdf5.h, libinterp/corefcn/oct-hist.h,
libinterp/corefcn/oct-iostrm.h, libinterp/corefcn/oct-lvalue.h,
libinterp/corefcn/oct-prcstrm.h, libinterp/corefcn/oct-procbuf.h,
libinterp/corefcn/oct-stdstrm.h, libinterp/corefcn/oct-stream.h,
libinterp/corefcn/oct-strstrm.h, libinterp/corefcn/octave-link.h,
libinterp/corefcn/pt-jit.h, libinterp/octave-value/ov-base-diag.h,
libinterp/octave-value/ov-base-int.h, libinterp/octave-value/ov-base-mat.h,
libinterp/octave-value/ov-base-scalar.h,
libinterp/octave-value/ov-base-sparse.h, libinterp/octave-value/ov-base.h,
libinterp/octave-value/ov-bool-mat.h, libinterp/octave-value/ov-bool-sparse.h,
libinterp/octave-value/ov-bool.h, libinterp/octave-value/ov-builtin.h,
libinterp/octave-value/ov-cell.h, libinterp/octave-value/ov-ch-mat.h,
libinterp/octave-value/ov-class.h, libinterp/octave-value/ov-colon.h,
libinterp/octave-value/ov-complex.h, libinterp/octave-value/ov-cs-list.h,
libinterp/octave-value/ov-cx-diag.h, libinterp/octave-value/ov-cx-mat.h,
libinterp/octave-value/ov-cx-sparse.h, libinterp/octave-value/ov-dld-fcn.h,
libinterp/octave-value/ov-fcn-handle.h, libinterp/octave-value/ov-fcn-inline.h,
libinterp/octave-value/ov-fcn.h, libinterp/octave-value/ov-float.h,
libinterp/octave-value/ov-flt-complex.h,
libinterp/octave-value/ov-flt-cx-diag.h,
libinterp/octave-value/ov-flt-cx-mat.h,
libinterp/octave-value/ov-flt-re-diag.h,
libinterp/octave-value/ov-flt-re-mat.h, libinterp/octave-value/ov-int-traits.h,
libinterp/octave-value/ov-int16.h, libinterp/octave-value/ov-int32.h,
libinterp/octave-value/ov-int64.h, libinterp/octave-value/ov-int8.h,
libinterp/octave-value/ov-java.h, libinterp/octave-value/ov-lazy-idx.h,
libinterp/octave-value/ov-mex-fcn.h, libinterp/octave-value/ov-null-mat.h,
libinterp/octave-value/ov-perm.h, libinterp/octave-value/ov-range.h,
libinterp/octave-value/ov-re-diag.h, libinterp/octave-value/ov-re-mat.h,
libinterp/octave-value/ov-re-sparse.h, libinterp/octave-value/ov-scalar.h,
libinterp/octave-value/ov-str-mat.h, libinterp/octave-value/ov-struct.h,
libinterp/octave-value/ov-typeinfo.h, libinterp/octave-value/ov-uint16.h,
libinterp/octave-value/ov-uint32.h, libinterp/octave-value/ov-uint64.h,
libinterp/octave-value/ov-uint8.h, libinterp/octave-value/ov-usr-fcn.h,
libinterp/octave-value/ov.h, libinterp/parse-tree/pt-all.h,
libinterp/parse-tree/pt-arg-list.h, libinterp/parse-tree/pt-array-list.h,
libinterp/parse-tree/pt-assign.h, libinterp/parse-tree/pt-binop.h,
libinterp/parse-tree/pt-bp.h, libinterp/parse-tree/pt-cbinop.h,
libinterp/parse-tree/pt-cell.h, libinterp/parse-tree/pt-check.h,
libinterp/parse-tree/pt-cmd.h, libinterp/parse-tree/pt-colon.h,
libinterp/parse-tree/pt-const.h, libinterp/parse-tree/pt-decl.h,
libinterp/parse-tree/pt-eval.h, libinterp/parse-tree/pt-except.h,
libinterp/parse-tree/pt-exp.h, libinterp/parse-tree/pt-fcn-handle.h,
libinterp/parse-tree/pt-id.h, libinterp/parse-tree/pt-idx.h,
libinterp/parse-tree/pt-jump.h, libinterp/parse-tree/pt-loop.h,
libinterp/parse-tree/pt-mat.h, libinterp/parse-tree/pt-misc.h,
libinterp/parse-tree/pt-pr-code.h, libinterp/parse-tree/pt-select.h,
libinterp/parse-tree/pt-stmt.h, libinterp/parse-tree/pt-unop.h,
libinterp/parse-tree/pt-walk.h, libinterp/parse-tree/pt.h,
liboctave/array/CColVector.h, liboctave/array/CDiagMatrix.h,
liboctave/array/CMatrix.h, liboctave/array/CNDArray.h,
liboctave/array/CRowVector.h, liboctave/array/boolMatrix.h,
liboctave/array/chMatrix.h, liboctave/array/chNDArray.h,
liboctave/array/dColVector.h, liboctave/array/dDiagMatrix.h,
liboctave/array/dMatrix.h, liboctave/array/dNDArray.h,
liboctave/array/dRowVector.h, liboctave/array/fCColVector.h,
liboctave/array/fCDiagMatrix.h, liboctave/array/fCMatrix.h,
liboctave/array/fCNDArray.h, liboctave/array/fCRowVector.h,
liboctave/array/fColVector.h, liboctave/array/fDiagMatrix.h,
liboctave/array/fMatrix.h, liboctave/array/fNDArray.h,
liboctave/array/fRowVector.h, liboctave/cruft/misc/lo-error.h,
liboctave/numeric/CmplxAEPBAL.h, liboctave/numeric/CmplxCHOL.h,
liboctave/numeric/CmplxGEPBAL.h, liboctave/numeric/CmplxHESS.h,
liboctave/numeric/CmplxLU.h, liboctave/numeric/CmplxQR.h,
liboctave/numeric/CmplxQRP.h, liboctave/numeric/CmplxSCHUR.h,
liboctave/numeric/CmplxSVD.h, liboctave/numeric/Quad.h,
liboctave/numeric/SparseCmplxCHOL.h, liboctave/numeric/SparseCmplxLU.h,
liboctave/numeric/SparseCmplxQR.h, liboctave/numeric/SparseQR.h,
liboctave/numeric/SparsedbleCHOL.h, liboctave/numeric/SparsedbleLU.h,
liboctave/numeric/bsxfun.h, liboctave/numeric/dbleAEPBAL.h,
liboctave/numeric/dbleCHOL.h, liboctave/numeric/dbleGEPBAL.h,
liboctave/numeric/dbleHESS.h, liboctave/numeric/dbleLU.h,
liboctave/numeric/dbleQR.h, liboctave/numeric/dbleQRP.h,
liboctave/numeric/dbleSCHUR.h, liboctave/numeric/dbleSVD.h,
liboctave/numeric/fCmplxAEPBAL.h, liboctave/numeric/fCmplxCHOL.h,
liboctave/numeric/fCmplxGEPBAL.h, liboctave/numeric/fCmplxHESS.h,
liboctave/numeric/fCmplxLU.h, liboctave/numeric/fCmplxQR.h,
liboctave/numeric/fCmplxQRP.h, liboctave/numeric/fCmplxSCHUR.h,
liboctave/numeric/fCmplxSVD.h, liboctave/numeric/fEIG.h,
liboctave/numeric/floatAEPBAL.h, liboctave/numeric/floatCHOL.h,
liboctave/numeric/floatGEPBAL.h, liboctave/numeric/floatHESS.h,
liboctave/numeric/floatLU.h, liboctave/numeric/floatQR.h,
liboctave/numeric/floatQRP.h, liboctave/numeric/floatSCHUR.h,
liboctave/numeric/floatSVD.h, liboctave/numeric/lo-mappers.h,
liboctave/numeric/lo-specfun.h, liboctave/numeric/oct-convn.h,
liboctave/numeric/oct-norm.h, liboctave/numeric/oct-rand.h,
liboctave/operators/Sparse-diag-op-defs.h,
liboctave/operators/Sparse-op-defs.h,
liboctave/operators/Sparse-perm-op-defs.h, liboctave/system/lo-sysdep.h,
liboctave/system/oct-env.h, liboctave/system/oct-group.h,
liboctave/system/oct-openmp.h, liboctave/system/oct-passwd.h,
liboctave/system/oct-syscalls.h, liboctave/system/oct-time.h,
liboctave/system/oct-uname.h, liboctave/util/lo-array-gripes.h,
liboctave/util/lo-cutils.h, liboctave/util/lo-ieee.h,
liboctave/util/lo-macros.h, liboctave/util/lo-math.h,
liboctave/util/lo-regexp.h, liboctave/util/lo-traits.h,
liboctave/util/lo-utils.h, liboctave/util/oct-base64.h,
liboctave/util/oct-binmap.h, liboctave/util/oct-glob.h,
liboctave/util/oct-inttypes.h, liboctave/util/oct-locbuf.h,
liboctave/util/oct-mem.h, liboctave/util/oct-mutex.h,
liboctave/util/oct-refcount.h, liboctave/util/oct-rl-edit.h,
liboctave/util/oct-rl-hist.h, liboctave/util/oct-shlib.h,
liboctave/util/oct-sort.h, liboctave/util/oct-sparse.h:
Use common #define syntax "octave_filename_h" in h_files.
author | Rik <rik@octave.org> |
---|---|
date | Thu, 31 Oct 2013 22:19:07 -0700 |
parents | 49a5a4be04a1 |
children | 4197fc428c7d |
rev | line source |
---|---|
3 | 1 /* |
2 | |
17744
d63878346099
maint: Update copyright notices for release.
John W. Eaton <jwe@octave.org>
parents:
15271
diff
changeset
|
3 Copyright (C) 1993-2013 John W. Eaton |
3 | 4 |
5 This file is part of Octave. | |
6 | |
7 Octave is free software; you can redistribute it and/or modify it | |
8 under the terms of the GNU General Public License as published by the | |
7016 | 9 Free Software Foundation; either version 3 of the License, or (at your |
10 option) any later version. | |
3 | 11 |
12 Octave is distributed in the hope that it will be useful, but WITHOUT | |
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
15 for more details. | |
16 | |
17 You should have received a copy of the GNU General Public License | |
7016 | 18 along with Octave; see the file COPYING. If not, see |
19 <http://www.gnu.org/licenses/>. | |
3 | 20 |
21 */ | |
22 | |
382 | 23 #if !defined (octave_Quad_h) |
24 #define octave_Quad_h 1 | |
25 | |
1536 | 26 #include <cfloat> |
27 | |
1296 | 28 #include "dColVector.h" |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
29 #include "fColVector.h" |
7231 | 30 #include "lo-math.h" |
3 | 31 |
32 typedef double (*integrand_fcn) (double x); | |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
33 typedef float (*float_integrand_fcn) (float x); |
3 | 34 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
35 // FIXME: would be nice to not have to have this global variable. |
260 | 36 // Nonzero means an error occurred in the calculation of the integrand |
37 // function, and the user wants us to quit. | |
6108 | 38 extern OCTAVE_API int quad_integration_error; |
260 | 39 |
3998 | 40 #include "Quad-opts.h" |
289 | 41 |
1859 | 42 class |
6108 | 43 OCTAVE_API |
1859 | 44 Quad : public Quad_options |
3 | 45 { |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
46 public: |
3 | 47 |
1859 | 48 Quad (integrand_fcn fcn) |
11501
331fcc41ca23
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
49 : Quad_options (), f (fcn), ff () { } |
1859 | 50 |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
51 Quad (float_integrand_fcn fcn) |
11501
331fcc41ca23
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
52 : Quad_options (), f (), ff (fcn) { } |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
53 |
2425 | 54 virtual ~Quad (void) { } |
2424 | 55 |
1528 | 56 virtual double integrate (void) |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
57 { |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
58 octave_idx_type ier, neval; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
59 double abserr; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
60 return do_integrate (ier, neval, abserr); |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
61 } |
3 | 62 |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
63 virtual float float_integrate (void) |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
64 { |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
65 octave_idx_type ier, neval; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
66 float abserr; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
67 return do_integrate (ier, neval, abserr); |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
68 } |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
69 |
5275 | 70 virtual double integrate (octave_idx_type& ier) |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
71 { |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
72 octave_idx_type neval; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
73 double abserr; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
74 return do_integrate (ier, neval, abserr); |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
75 } |
1528 | 76 |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
77 virtual float float_integrate (octave_idx_type& ier) |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
78 { |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
79 octave_idx_type neval; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
80 float abserr; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
81 return do_integrate (ier, neval, abserr); |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
82 } |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
83 |
5275 | 84 virtual double integrate (octave_idx_type& ier, octave_idx_type& neval) |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
85 { |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
86 double abserr; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
87 return do_integrate (ier, neval, abserr); |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
88 } |
1528 | 89 |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
90 virtual float float_integrate (octave_idx_type& ier, octave_idx_type& neval) |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
91 { |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
92 float abserr; |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
93 return do_integrate (ier, neval, abserr); |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
94 } |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
95 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
96 virtual double integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
97 double& abserr) |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
98 { |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
99 return do_integrate (ier, neval, abserr); |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
100 } |
3511 | 101 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
102 virtual float float_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
103 float& abserr) |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
104 { |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
105 return do_integrate (ier, neval, abserr); |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
106 } |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
107 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
108 virtual double do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
109 double& abserr) = 0; |
3 | 110 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
111 virtual float do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
112 float& abserr) = 0; |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
113 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
114 protected: |
3 | 115 |
116 integrand_fcn f; | |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
117 float_integrand_fcn ff; |
3 | 118 }; |
119 | |
1859 | 120 class |
6108 | 121 OCTAVE_API |
1859 | 122 DefQuad : public Quad |
3 | 123 { |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
124 public: |
3 | 125 |
1859 | 126 DefQuad (integrand_fcn fcn) |
127 : Quad (fcn), lower_limit (0.0), upper_limit (1.0), singularities () { } | |
1528 | 128 |
1859 | 129 DefQuad (integrand_fcn fcn, double ll, double ul) |
130 : Quad (fcn), lower_limit (ll), upper_limit (ul), singularities () { } | |
1528 | 131 |
132 DefQuad (integrand_fcn fcn, double ll, double ul, | |
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
133 const ColumnVector& sing) |
1859 | 134 : Quad (fcn), lower_limit (ll), upper_limit (ul), |
135 singularities (sing) { } | |
1528 | 136 |
1859 | 137 DefQuad (integrand_fcn fcn, const ColumnVector& sing) |
138 : Quad (fcn), lower_limit (0.0), upper_limit (1.0), | |
139 singularities (sing) { } | |
1528 | 140 |
2424 | 141 ~DefQuad (void) { } |
142 | |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
143 double do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
144 double& abserr); |
3 | 145 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
146 float do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
147 float& abserr); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
148 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
149 private: |
3 | 150 |
151 double lower_limit; | |
152 double upper_limit; | |
153 | |
1528 | 154 ColumnVector singularities; |
3 | 155 }; |
156 | |
1859 | 157 class |
6108 | 158 OCTAVE_API |
1859 | 159 IndefQuad : public Quad |
3 | 160 { |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
161 public: |
3 | 162 |
163 enum IntegralType { bound_to_inf, neg_inf_to_bound, doubly_infinite }; | |
164 | |
1859 | 165 IndefQuad (integrand_fcn fcn) |
11501
331fcc41ca23
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
166 : Quad (fcn), bound (0.0), type (bound_to_inf), integration_error (0) { } |
1528 | 167 |
1859 | 168 IndefQuad (integrand_fcn fcn, double b, IntegralType t) |
11501
331fcc41ca23
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
169 : Quad (fcn), bound (b), type (t), integration_error (0) { } |
1528 | 170 |
2424 | 171 ~IndefQuad (void) { } |
172 | |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
173 double do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
174 double& abserr); |
3 | 175 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
176 float do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
177 float& abserr); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
178 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
179 private: |
3 | 180 |
181 double bound; | |
182 IntegralType type; | |
1859 | 183 int integration_error; |
3 | 184 }; |
185 | |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
186 class |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
187 OCTAVE_API |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
188 FloatDefQuad : public Quad |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
189 { |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
190 public: |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
191 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
192 FloatDefQuad (float_integrand_fcn fcn) |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
193 : Quad (fcn), lower_limit (0.0), upper_limit (1.0), singularities () { } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
194 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
195 FloatDefQuad (float_integrand_fcn fcn, float ll, float ul) |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
196 : Quad (fcn), lower_limit (ll), upper_limit (ul), singularities () { } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
197 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
198 FloatDefQuad (float_integrand_fcn fcn, float ll, float ul, |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
199 const FloatColumnVector& sing) |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
200 : Quad (fcn), lower_limit (ll), upper_limit (ul), |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
201 singularities (sing) { } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
202 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
203 FloatDefQuad (float_integrand_fcn fcn, const FloatColumnVector& sing) |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
204 : Quad (fcn), lower_limit (0.0), upper_limit (1.0), |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
205 singularities (sing) { } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
206 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
207 ~FloatDefQuad (void) { } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
208 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
209 double do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
210 double& abserr); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
211 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
212 float do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
213 float& abserr); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
214 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
215 private: |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
216 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
217 float lower_limit; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
218 float upper_limit; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
219 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
220 FloatColumnVector singularities; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
221 }; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
222 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
223 class |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
224 OCTAVE_API |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
225 FloatIndefQuad : public Quad |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
226 { |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
227 public: |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
228 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
229 enum IntegralType { bound_to_inf, neg_inf_to_bound, doubly_infinite }; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
230 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
231 FloatIndefQuad (float_integrand_fcn fcn) |
11501
331fcc41ca23
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
232 : Quad (fcn), bound (0.0), type (bound_to_inf), integration_error (0) { } |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
233 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
234 FloatIndefQuad (float_integrand_fcn fcn, double b, IntegralType t) |
11501
331fcc41ca23
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
235 : Quad (fcn), bound (b), type (t), integration_error (0) { } |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
236 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
237 ~FloatIndefQuad (void) { } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
238 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
239 double do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
240 double& abserr); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
241 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
242 float do_integrate (octave_idx_type& ier, octave_idx_type& neval, |
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
243 float& abserr); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
244 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
245 private: |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
246 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
247 float bound; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
248 IntegralType type; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
249 int integration_error; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
250 }; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7231
diff
changeset
|
251 |
3 | 252 #endif |