Mercurial > octave
view liboctave/numeric/Quad.h @ 32078:632f9b828de1
Avoid using file_stat in liboctave/util (bug #59711).
* cmd-edit.cc (looks_like_filename), cmd-hist.cc (gnu_history::do_append),
kpse.cc (kpse_element_dir), oct-glob.cc (glob, windows_glob),
url-transfer.cc (base_url_transfer::mget_directory): Use functions "dir_exists"
or "file_exists" instead of "file_stat".
* kpse.cc (dir_p), oct-glob.cc (single_match_exists): Remove unused static
functions.
author | Markus Mützel <markus.muetzel@gmx.de> |
---|---|
date | Sat, 06 May 2023 10:56:33 +0200 |
parents | 1daf8bfceac3 |
children | ae4e19c0a2b1 2e484f9f1f18 |
line wrap: on
line source
//////////////////////////////////////////////////////////////////////// // // Copyright (C) 1993-2023 The Octave Project Developers // // See the file COPYRIGHT.md in the top-level directory of this // distribution or <https://octave.org/copyright/>. // // 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 // <https://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); #include "Quad-opts.h" class OCTAVE_API Quad : public Quad_options { public: Quad () = delete; Quad (integrand_fcn fcn) : Quad_options (), m_f (fcn), m_ff () { } Quad (float_integrand_fcn fcn) : Quad_options (), m_f (), m_ff (fcn) { } OCTAVE_DEFAULT_COPY_MOVE (Quad) virtual ~Quad () = default; virtual double integrate () { octave_idx_type ier, neval; double abserr; return do_integrate (ier, neval, abserr); } virtual float float_integrate () { 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 m_f; float_integrand_fcn m_ff; }; class OCTAVE_API DefQuad : public Quad { public: DefQuad () = delete; DefQuad (integrand_fcn fcn) : Quad (fcn), m_lower_limit (0.0), m_upper_limit (1.0), m_singularities () { } DefQuad (integrand_fcn fcn, double ll, double ul) : Quad (fcn), m_lower_limit (ll), m_upper_limit (ul), m_singularities () { } DefQuad (integrand_fcn fcn, double ll, double ul, const ColumnVector& sing) : Quad (fcn), m_lower_limit (ll), m_upper_limit (ul), m_singularities (sing) { } DefQuad (integrand_fcn fcn, const ColumnVector& sing) : Quad (fcn), m_lower_limit (0.0), m_upper_limit (1.0), m_singularities (sing) { } OCTAVE_DEFAULT_COPY_MOVE_DELETE (DefQuad) 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 m_lower_limit; double m_upper_limit; ColumnVector m_singularities; }; class OCTAVE_API IndefQuad : public Quad { public: enum IntegralType { bound_to_inf, neg_inf_to_bound, doubly_infinite }; IndefQuad () = delete; IndefQuad (integrand_fcn fcn) : Quad (fcn), m_bound (0.0), m_type (bound_to_inf) { } IndefQuad (integrand_fcn fcn, double b, IntegralType t) : Quad (fcn), m_bound (b), m_type (t) { } OCTAVE_DEFAULT_COPY_MOVE_DELETE (IndefQuad) 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 m_bound; IntegralType m_type; }; class OCTAVE_API FloatDefQuad : public Quad { public: FloatDefQuad () = delete; FloatDefQuad (float_integrand_fcn fcn) : Quad (fcn), m_lower_limit (0.0), m_upper_limit (1.0), m_singularities () { } FloatDefQuad (float_integrand_fcn fcn, float ll, float ul) : Quad (fcn), m_lower_limit (ll), m_upper_limit (ul), m_singularities () { } FloatDefQuad (float_integrand_fcn fcn, float ll, float ul, const FloatColumnVector& sing) : Quad (fcn), m_lower_limit (ll), m_upper_limit (ul), m_singularities (sing) { } FloatDefQuad (float_integrand_fcn fcn, const FloatColumnVector& sing) : Quad (fcn), m_lower_limit (0.0), m_upper_limit (1.0), m_singularities (sing) { } OCTAVE_DEFAULT_COPY_MOVE_DELETE (FloatDefQuad) 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 m_lower_limit; float m_upper_limit; FloatColumnVector m_singularities; }; class OCTAVE_API FloatIndefQuad : public Quad { public: enum IntegralType { bound_to_inf, neg_inf_to_bound, doubly_infinite }; FloatIndefQuad () = delete; FloatIndefQuad (float_integrand_fcn fcn) : Quad (fcn), m_bound (0.0), m_type (bound_to_inf) { } FloatIndefQuad (float_integrand_fcn fcn, double b, IntegralType t) : Quad (fcn), m_bound (b), m_type (t) { } OCTAVE_DEFAULT_COPY_MOVE_DELETE (FloatIndefQuad) 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 m_bound; IntegralType m_type; }; #endif