Mercurial > octave
view liboctave/numeric/Quad.h @ 23443:3f1bf237908b
maint: Eliminate <cfloat.h> header from liboctave files.
* build-aux/mk-opts.pl: Remove "#include <cfloat>" from script which creates *-opts.h
* Array-util.cc, Array-util.h, Array.cc, Array.h, CColVector.h,
CDiagMatrix.h, CMatrix.cc, CMatrix.h, CNDArray.cc, CNDArray.h, CRowVector.h,
CSparse.cc, CSparse.h, MSparse.h, Matrix.h, Range.cc, Sparse.h, boolMatrix.h,
boolNDArray.h, boolSparse.h, chMatrix.h, chNDArray.h, dColVector.h,
dDiagMatrix.h, dMatrix.cc, dMatrix.h, dNDArray.cc, dNDArray.h, dRowVector.h,
dSparse.cc, dSparse.h, dim-vector.cc, dim-vector.h, fCColVector.h,
fCDiagMatrix.h, fCMatrix.cc, fCMatrix.h, fCNDArray.cc, fCNDArray.h,
fCRowVector.h, fColVector.h, fDiagMatrix.h, fMatrix.cc, fMatrix.h, fNDArray.cc,
fNDArray.h, fRowVector.h, int16NDArray.h, int32NDArray.h, int64NDArray.h,
int8NDArray.h, intNDArray.h, uint16NDArray.h, uint32NDArray.h, uint64NDArray.h,
uint8NDArray.h, CollocWt.cc, DASPK.cc, DASPK.h, DASRT.cc, DASRT.h, DASSL.cc,
DASSL.h, LSODE.cc, LSODE.h, Quad.h, eigs-base.cc, eigs-base.h, lo-mappers.cc,
lo-mappers.h, oct-norm.cc, lo-ieee.cc, lo-utils.cc, lo-utils.h, oct-refcount.h:
Remove "#include <cfloat>". Add "#include <limits>" if necessary. Remove
redundant #includes. Alphabetize #includes.
author | Rik <rik@octave.org> |
---|---|
date | Wed, 26 Apr 2017 09:57:28 -0700 |
parents | 092078913d54 |
children | 194eb4bd202b |
line wrap: on
line source
/* Copyright (C) 1993-2017 John W. Eaton This file is part of Octave. Octave is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version. Octave is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with Octave; see the file COPYING. If not, see <http://www.gnu.org/licenses/>. */ #if ! defined (octave_Quad_h) #define octave_Quad_h 1 #include "octave-config.h" #include "dColVector.h" #include "fColVector.h" typedef double (*integrand_fcn) (double x); typedef float (*float_integrand_fcn) (float x); // FIXME: would be nice to not have to have this global variable. // Nonzero means an error occurred in the calculation of the integrand // function, and the user wants us to quit. extern OCTAVE_API int quad_integration_error; #include "Quad-opts.h" class OCTAVE_API Quad : public Quad_options { public: Quad (integrand_fcn fcn) : Quad_options (), f (fcn), ff () { } Quad (float_integrand_fcn fcn) : Quad_options (), f (), ff (fcn) { } virtual ~Quad (void) = default; virtual double integrate (void) { octave_idx_type ier, neval; double abserr; return do_integrate (ier, neval, abserr); } virtual float float_integrate (void) { octave_idx_type ier, neval; float abserr; return do_integrate (ier, neval, abserr); } virtual double integrate (octave_idx_type& ier) { octave_idx_type neval; double abserr; return do_integrate (ier, neval, abserr); } virtual float float_integrate (octave_idx_type& ier) { octave_idx_type neval; float abserr; return do_integrate (ier, neval, abserr); } virtual double integrate (octave_idx_type& ier, octave_idx_type& neval) { double abserr; return do_integrate (ier, neval, abserr); } virtual float float_integrate (octave_idx_type& ier, octave_idx_type& neval) { float abserr; return do_integrate (ier, neval, abserr); } virtual double integrate (octave_idx_type& ier, octave_idx_type& neval, double& abserr) { return do_integrate (ier, neval, abserr); } virtual float float_integrate (octave_idx_type& ier, octave_idx_type& neval, float& abserr) { return do_integrate (ier, neval, abserr); } virtual double do_integrate (octave_idx_type& ier, octave_idx_type& neval, double& abserr) = 0; virtual float do_integrate (octave_idx_type& ier, octave_idx_type& neval, float& abserr) = 0; protected: integrand_fcn f; float_integrand_fcn ff; }; class OCTAVE_API DefQuad : public Quad { public: DefQuad (integrand_fcn fcn) : Quad (fcn), lower_limit (0.0), upper_limit (1.0), singularities () { } DefQuad (integrand_fcn fcn, double ll, double ul) : Quad (fcn), lower_limit (ll), upper_limit (ul), singularities () { } DefQuad (integrand_fcn fcn, double ll, double ul, const ColumnVector& sing) : Quad (fcn), lower_limit (ll), upper_limit (ul), singularities (sing) { } DefQuad (integrand_fcn fcn, const ColumnVector& sing) : Quad (fcn), lower_limit (0.0), upper_limit (1.0), singularities (sing) { } ~DefQuad (void) = default; double do_integrate (octave_idx_type& ier, octave_idx_type& neval, double& abserr); OCTAVE_NORETURN float do_integrate (octave_idx_type& ier, octave_idx_type& neval, float& abserr); private: double lower_limit; double upper_limit; ColumnVector singularities; }; class OCTAVE_API IndefQuad : public Quad { public: enum IntegralType { bound_to_inf, neg_inf_to_bound, doubly_infinite }; IndefQuad (integrand_fcn fcn) : Quad (fcn), bound (0.0), type (bound_to_inf) { } IndefQuad (integrand_fcn fcn, double b, IntegralType t) : Quad (fcn), bound (b), type (t) { } ~IndefQuad (void) = default; double do_integrate (octave_idx_type& ier, octave_idx_type& neval, double& abserr); OCTAVE_NORETURN float do_integrate (octave_idx_type& ier, octave_idx_type& neval, float& abserr); private: double bound; IntegralType type; }; class OCTAVE_API FloatDefQuad : public Quad { public: FloatDefQuad (float_integrand_fcn fcn) : Quad (fcn), lower_limit (0.0), upper_limit (1.0), singularities () { } FloatDefQuad (float_integrand_fcn fcn, float ll, float ul) : Quad (fcn), lower_limit (ll), upper_limit (ul), singularities () { } FloatDefQuad (float_integrand_fcn fcn, float ll, float ul, const FloatColumnVector& sing) : Quad (fcn), lower_limit (ll), upper_limit (ul), singularities (sing) { } FloatDefQuad (float_integrand_fcn fcn, const FloatColumnVector& sing) : Quad (fcn), lower_limit (0.0), upper_limit (1.0), singularities (sing) { } ~FloatDefQuad (void) = default; OCTAVE_NORETURN double do_integrate (octave_idx_type& ier, octave_idx_type& neval, double& abserr); float do_integrate (octave_idx_type& ier, octave_idx_type& neval, float& abserr); private: float lower_limit; float upper_limit; FloatColumnVector singularities; }; class OCTAVE_API FloatIndefQuad : public Quad { public: enum IntegralType { bound_to_inf, neg_inf_to_bound, doubly_infinite }; FloatIndefQuad (float_integrand_fcn fcn) : Quad (fcn), bound (0.0), type (bound_to_inf) { } FloatIndefQuad (float_integrand_fcn fcn, double b, IntegralType t) : Quad (fcn), bound (b), type (t) { } ~FloatIndefQuad (void) = default; OCTAVE_NORETURN double do_integrate (octave_idx_type& ier, octave_idx_type& neval, double& abserr); float do_integrate (octave_idx_type& ier, octave_idx_type& neval, float& abserr); private: float bound; IntegralType type; }; #endif