Mercurial > octave
view liboctave/Array-fC.cc @ 8721:e9cb742df9eb
imported patch sort3.diff
author | Jaroslav Hajek <highegg@gmail.com> |
---|---|
date | Wed, 11 Feb 2009 15:25:53 +0100 |
parents | 314be237cd5b |
children | d5af326a3ede |
line wrap: on
line source
/* Copyright (C) 1994, 1995, 1996, 1997, 1998, 2001, 2003, 2004, 2005, 2006, 2007 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/>. */ #ifdef HAVE_CONFIG_H #include <config.h> #endif // Instantiate Arrays of FloatComplex values. #include "oct-cmplx.h" #include "lo-mappers.h" #include "Array.h" #include "Array.cc" #include "oct-sort.cc" template <> inline bool _sort_isnan (FloatComplex x) { return xisnan (x); } template <> bool octave_sort<FloatComplex>::ascending_compare (FloatComplex a, FloatComplex b) { return ((std::abs (a) < std::abs (b)) || ((std::abs (a) == std::abs (b)) && (arg (a) < arg (b)))); } template <> bool octave_sort<FloatComplex>::descending_compare (FloatComplex a, FloatComplex b) { return ((std::abs (a) > std::abs (b)) || ((std::abs (a) == std::abs (b)) && (arg (a) > arg (b)))); } static bool nan_ascending_compare (FloatComplex x, FloatComplex y) { return xisnan (y) ? ! xisnan (x) : ((std::abs (x) < std::abs (x)) || ((std::abs (x) == std::abs (x)) && (arg (x) < arg (x)))); } static bool nan_descending_compare (FloatComplex x, FloatComplex y) { return xisnan (x) ? ! xisnan (y) : ((std::abs (x) > std::abs (x)) || ((std::abs (x) == std::abs (x)) && (arg (x) > arg (x)))); } bool (*_sortrows_comparator (sortmode mode, const Array<FloatComplex>& a , bool allow_chk)) (FloatComplex, FloatComplex) { bool (*result) (FloatComplex, FloatComplex) = 0; if (allow_chk) { octave_idx_type k = 0; for (; k < a.numel () && ! xisnan (a(k)); k++) ; if (k == a.numel ()) { if (mode == ASCENDING) result = octave_sort<FloatComplex>::ascending_compare; else if (mode == DESCENDING) result = octave_sort<FloatComplex>::descending_compare; } } if (! result) { if (mode == ASCENDING) result = nan_ascending_compare; else if (mode == DESCENDING) result = nan_descending_compare; } return result; } INSTANTIATE_ARRAY_SORT (FloatComplex); INSTANTIATE_ARRAY (FloatComplex, OCTAVE_API); #include "Array2.h" template class OCTAVE_API Array2<FloatComplex>; #include "ArrayN.h" #include "ArrayN.cc" template class OCTAVE_API ArrayN<FloatComplex>; template OCTAVE_API std::ostream& operator << (std::ostream&, const ArrayN<FloatComplex>&); #include "DiagArray2.h" #include "DiagArray2.cc" template class OCTAVE_API DiagArray2<FloatComplex>; /* ;;; Local Variables: *** ;;; mode: C++ *** ;;; End: *** */