Mercurial > hg > octave-lyh
view libinterp/octave-value/ov-oncleanup.h @ 15696:2fac72a256ce
Add complex erf,erfc,erfcx,erfi,dawson routines from Faddeeva package.
* libinterp/corefcn/mappers.cc: Add erfi and dawson mapper functions,
and add complex-argument test cases for erf, erfc, erfcx, erfi, and
dawson.
* libinterp/octave-value/ov-base.cc, libinterp/octave-value/ov-base.h:
Add erfi and dawson mapper functions.
* libinterp/octave-value/ov-complex.cc, libinterp/octave-value/ov-cx-mat.cc,
libinterp/octave-value/ov-cx-sparse.cc, libinterp/octave-value/ov-float.cc,
libinterp/octave-value/ov-flt-complex.cc,
libinterp/octave-value/ov-flt-cx-mat.cc,
libinterp/octave-value/ov-flt-re-mat.cc,
libinterp/octave-value/ov-re-mat.c,
libinterp/octave-value/ov-re-sparse.cc,
libinterp/octave-value/ov-scalar.cc,
libinterp/octave-value/ov.h:
Support erf, erfc, erfcx, erfi, and dawson mapper functions for
real and complex matrices and scalars.
* liboctave/cruft/Faddeeva/Faddeeva.cc, liboctave/cruft/Faddeeva/Faddeeva.hh:
liboctave/cruft/Faddeeva/module.mk, liboctave/cruft/Makefile.am:
Add Faddeeva package (from http://ab-initio.mit.edu/Faddeeva) to
libcruft, to provide the various complex-argument error functions.
* liboctave/numeric/lo-specfun.cc, liboctave/numeric/lo-specfun.h:
Add complex-argument erf, erfc, erfcx, erfi, and dawson functions
to liboctave API. Delete previous real-argument erfcx implementation
in favor of Faddeeva::erfcx (which seems to be slightly faster
in gcc/x86-64 benchmarks, with similar accuracy).
* doc/interpreter/arith.txi: Include erfi and dawson documentation.
author | Steven G. Johnson <stevenj@alum.mit.edu> |
---|---|
date | Tue, 27 Nov 2012 23:39:54 -0500 |
parents | 2fc554ffbc28 |
children |
line wrap: on
line source
/* Copyright (C) 2010-2012 VZLU Prague 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/>. */ #ifdef HAVE_CONFIG_H #include <config.h> #endif #include <iosfwd> #include "ov-base.h" #include "ov-struct.h" #include "ov.h" static void gripe_internal (void) { error ("onCleanup: internal error: cloning nonempty object"); } class octave_oncleanup : public octave_base_value { public: octave_oncleanup (void) : fcn () { } octave_oncleanup (const octave_value& fcn); octave_base_value *clone (void) const { if (fcn.is_defined ()) gripe_internal (); return empty_clone (); } octave_base_value *empty_clone (void) const { return new octave_oncleanup (); } ~octave_oncleanup (void); bool is_defined (void) const { return true; } bool is_object (void) const { return true; } // do we want this? octave_map map_value (void) const { return scalar_map_value (); } octave_scalar_map scalar_map_value (void) const; dim_vector dims (void) const { static dim_vector dv (1, 1); return dv; } bool save_ascii (std::ostream& os); bool load_ascii (std::istream& is); bool save_binary (std::ostream& os, bool& save_as_floats); bool load_binary (std::istream& is, bool swap, oct_mach_info::float_format fmt); #if defined (HAVE_HDF5) bool save_hdf5 (hid_t loc_id, const char *name, bool save_as_floats); bool load_hdf5 (hid_t loc_id, const char *name); #endif void print (std::ostream& os, bool pr_as_read_syntax = false) const; void print_raw (std::ostream& os, bool pr_as_read_syntax = false) const; private: DECLARE_OV_TYPEID_FUNCTIONS_AND_DATA protected: octave_value fcn; };