Mercurial > hg > octave-lyh
view liboctave/lo-traits.h @ 8966:1bba53c0a38d
Implement diag + sparse, diag - sparse, sparse + diag, sparse - diag.
Date: Mon, 9 Mar 2009 17:45:22 -0400
This does not use the typical sparse-mx-ops generator. I suspect the
sematics of elementwise multiplication and division to be rather
controversial, so they are not included. If comparison operations are
added, the implementation should be shifted over to use the typical
generator.
The template in Sparse-diag-op-defs.h likely could use const bools
rather than functional argument operations. I haven't measured which
is optimized more effectively.
Also, the Octave binding layer in op-dm-scm.cc likely could use all
sorts of macro or template trickery, but it's far easier to let Emacs
handle it for now. That would be worth revisiting if further
elementwise sparse and diagonal operations are added.
author | Jason Riedy <jason@acm.org> |
---|---|
date | Mon, 09 Mar 2009 17:49:14 -0400 |
parents | 0f6683a8150a |
children | e793865ede63 |
line wrap: on
line source
/* Copyright (C) 2009 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_liboctave_traits_h) #define octave_liboctave_traits_h 1 // Ideas for these classes taken from C++ Templates, The Complete // Guide by David Vandevoorde and Nicolai M. Josuttis, Addison-Wesley // (2003). // Select a type based on the value of a constant expression. template <bool cond, typename T1, typename T2> class if_then_else; template<typename T1, typename T2> class if_then_else<true, T1, T2> { public: typedef T1 result; }; template<typename T1, typename T2> class if_then_else<false, T1, T2> { public: typedef T2 result; }; // Determine whether a template paramter is a class type. template<typename T1> class is_class_type { private: typedef char one; typedef struct { char c[2]; } two; // Classes can have pointers to members. template<typename T2> static one is_class_type_test (int T2::*); // Catch everything else. template<typename T2> static two is_class_type_test (...); public: enum { yes = sizeof (is_class_type_test<T1> (0)) == 1 }; enum { no = ! yes }; }; // Define typename ref_param<T>::type as T const& if T is a class // type. Otherwise, define it to be T. template<typename T> class ref_param { public: typedef typename if_then_else<is_class_type<T>::no, T, T const&>::result type; }; #endif /* ;;; Local Variables: *** ;;; mode: C *** ;;; page-delimiter: "^/\\*" *** ;;; End: *** */