Mercurial > hg > octave-lyh
view libinterp/octave-value/ov-int-traits.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) 2004-2012 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_value_int_traits_h) #define octave_value_int_traits_h 1 #include "ov-int8.h" #include "ov-int16.h" #include "ov-int32.h" #include "ov-int64.h" #include "ov-uint8.h" #include "ov-uint16.h" #include "ov-uint32.h" #include "ov-uint64.h" template <class T> class octave_value_int_traits { public: typedef T scalar_type; }; #define OCTAVE_VALUE_INT_TRAITS(MT, ST) \ template<> \ class \ octave_value_int_traits<MT> \ { \ public: \ typedef ST scalar_type; \ } OCTAVE_VALUE_INT_TRAITS(int8NDArray, octave_int8_scalar); OCTAVE_VALUE_INT_TRAITS(int16NDArray, octave_int16_scalar); OCTAVE_VALUE_INT_TRAITS(int32NDArray, octave_int32_scalar); OCTAVE_VALUE_INT_TRAITS(int64NDArray, octave_int64_scalar); OCTAVE_VALUE_INT_TRAITS(uint8NDArray, octave_uint8_scalar); OCTAVE_VALUE_INT_TRAITS(uint16NDArray, octave_uint16_scalar); OCTAVE_VALUE_INT_TRAITS(uint32NDArray, octave_uint32_scalar); OCTAVE_VALUE_INT_TRAITS(uint64NDArray, octave_uint64_scalar); #endif