Mercurial > hg > octave-jordi

comparison libinterp/octave-value/ov-flt-cx-mat.cc @ 15696:2fac72a256ce

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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 d63878346099
comparison
equal deleted inserted replaced
15683:806ea52af230 15696:2fac72a256ce
742 ARRAY_MAPPER (arg, float, std::arg); 742 ARRAY_MAPPER (arg, float, std::arg);
743 ARRAY_MAPPER (asin, FloatComplex, ::asin); 743 ARRAY_MAPPER (asin, FloatComplex, ::asin);
744 ARRAY_MAPPER (asinh, FloatComplex, ::asinh); 744 ARRAY_MAPPER (asinh, FloatComplex, ::asinh);
745 ARRAY_MAPPER (atan, FloatComplex, ::atan); 745 ARRAY_MAPPER (atan, FloatComplex, ::atan);
746 ARRAY_MAPPER (atanh, FloatComplex, ::atanh); 746 ARRAY_MAPPER (atanh, FloatComplex, ::atanh);
747 ARRAY_MAPPER (erf, FloatComplex, ::erf);
748 ARRAY_MAPPER (erfc, FloatComplex, ::erfc);
749 ARRAY_MAPPER (erfcx, FloatComplex, ::erfcx);
750 ARRAY_MAPPER (erfi, FloatComplex, ::erfi);
751 ARRAY_MAPPER (dawson, FloatComplex, ::dawson);
747 ARRAY_MAPPER (ceil, FloatComplex, ::ceil); 752 ARRAY_MAPPER (ceil, FloatComplex, ::ceil);
748 ARRAY_MAPPER (cos, FloatComplex, std::cos); 753 ARRAY_MAPPER (cos, FloatComplex, std::cos);
749 ARRAY_MAPPER (cosh, FloatComplex, std::cosh); 754 ARRAY_MAPPER (cosh, FloatComplex, std::cosh);
750 ARRAY_MAPPER (exp, FloatComplex, std::exp); 755 ARRAY_MAPPER (exp, FloatComplex, std::exp);
751 ARRAY_MAPPER (expm1, FloatComplex, ::expm1); 756 ARRAY_MAPPER (expm1, FloatComplex, ::expm1);