#ifndef _MTL_MATRIX_TRAITS_#define _MTL_MATRIX_TRAITS_#include "mtl/mtl_complex.h"#include "mtl/mtl_config.h"namespace mtl {  enum { RECT, BAND, TRI, SYMM, HERM, ROW_MAJOR, COL_MAJOR, DIAG, DENSE,	 PACKED, BAND_VIEW, SPARSE, ARRAY, ENVELOPE,	 COMPRESSED, SPARSE_PAIR, TREE, LINKED_LIST,	 upper, lower, unit_upper, unit_lower, dynamic_uplo,         internal = 0, external,         index_from_one = -1, index_from_zero = 0 };//:  The "traits" class for MTL matrices.//!component: type//!category: containers, tagstemplate <class Matrix>struct matrix_traits {  //: The shape of the matrix, either rectangle_tag, banded_tag, diagonal_tag, triangle_tag, or symmetric_tag  typedef typename Matrix::shape shape;  //: The orientation, either row_tag or column_tag  typedef typename Matrix::orientation orientation;  //: The sparsity, either dense_tag or sparse_tag  typedef typename Matrix::sparsity sparsity;  //: Used by the trans helper function  typedef typename Matrix::transpose_type transpose_type;  //: Used by the rows and columns helper functions  typedef typename Matrix::strided_type strided_type;  //: Whether the rows and columns functions can be used with this Matrix  typedef typename Matrix::strideability strideability;  //: The Matrix type resulting from wrapping a scaled adator around this Matrix  typedef typename Matrix::scaled_type scaled_type;  //: Whether the Matrix owns its data, either external_tag or internal_tag  typedef typename Matrix::storage_loc storage_loc;  //: A OneD part of a Matrix. This could be a Row, a Column or a Diagonal depending on the type of Matrix.  typedef typename Matrix::OneD OneD;  //: The element type of the matrix  typedef typename Matrix::value_type value_type;  typedef typename Matrix::reference reference;  typedef typename Matrix::const_reference const_reference;  typedef typename Matrix::pointer pointer;  //: A NonNegativeIntegral type  typedef typename Matrix::size_type size_type;  typedef typename Matrix::difference_type difference_type;};//: Row Matrix Traits//!component: type//!category: containers, tagstemplate <class Matrix>struct row_matrix_traits {  typedef typename Matrix::Row Row;};//: Column Matrix Traits//!component: type//!category: containers, tagstemplate <class Matrix>struct column_matrix_traits {  typedef typename Matrix::Column Column;};//: Diagonal Matrix Traits//!component: type//!category: containers, tagstemplate <class Matrix>struct diagonal_matrix_traits {  typedef typename Matrix::Diagonal Diagonal;};/* * Shape Tags * add enum for shape tags. -- Rich *///: Identifies rectangular matrices//!component: type//!category: containers, tagsclass rectangle_tag { public: enum{ id = RECT}; }; //mh//: Identifies banded matrices//!component: type//!category: containers, tagsclass banded_tag { public: enum{ id = BAND};}; //mh//: Identifies triangular matrices//!component: type//!category: containers, tagsclass triangle_tag : public banded_tag { public: enum{id = TRI};}; //mh//: Identifies symmetric matrices//!component: type//!category: containers, tagsclass symmetric_tag : public banded_tag { public: enum{id=SYMM};}; //mh//: Identifies hermitian matrices//!component: type//!category: containers, tagsclass hermitian_tag : public banded_tag{ public: enum{id=HERM};}; //mh//: Identifies diagonal matrices//!component: type//!category: containers, tagsclass diagonal_tag : public banded_tag { public: enum{id=DIAG};}; //mh/* * Storage Tags *///: Identifies the Matrix as not owning its data//!component: type//!category: containers, tagsclass external_tag { };//: Identifies the Matrix as owning its data//!component: type//!category: containers, tagsclass internal_tag { };/* * Orientation Tags *///: Identifies a row-major Matrix//!component: type//!category: containers, tagsstruct row_tag { enum { id = ROW_MAJOR }; };//: Identifies a column-major Matrix//!component: type//!category: containers, tagsstruct column_tag { enum { id = COL_MAJOR }; };/* * Sparsity Tags *///: Identifies a dense Matrix or Vector//!component: type//!category: containers, tagsstruct dense_tag { enum { id = DENSE }; };//: Identifies a sparse Matrix or Vector//!component: type//!category: containers, tagsstruct sparse_tag { enum { id = SPARSE }; };/* * Traits for linear algebra objects (matrices and vectors) *///: Identifies linear algebra objects as 1D (Vector)//!component: type//!category: containers, tagsstruct oned_tag { };//: Identifies linear algebra objects as 2D (Matrix)//!component: type//!category: containers, tagsstruct twod_tag { };template <class number_type>struct number_traits {  typedef number_type magnitude_type;};#if MTL_PARTIAL_SPECtemplate <class T>struct number_traits< std::complex<T> > {  typedef T magnitude_type;};#elsetemplate <>struct number_traits< std::complex<double> > {  typedef double magnitude_type;};template <>struct number_traits< std::complex<float> > {  typedef float magnitude_type;};#endif//: Linear Algebra Object (Matrix and Vector) Traits //!component: type//!category: containers, tagstemplate <class Linalg>struct linalg_traits {  /*  enum { dimension = Linalg::dimension }; 1 for vectors, 2 for matrices */  //: Whether the object is a 1D or 2D container  typedef typename Linalg::dimension dimension;  //: The element type within the container  typedef typename Linalg::value_type value_type;  //: Either sparse or dense  typedef typename Linalg::sparsity sparsity;  //: The return type for abs(value_type)  typedef typename number_traits<value_type>::magnitude_type magnitude_type;};//: Identifies matrices that can be used with the rows and columns functions//!component: type//!category: containers, tagsstruct strideable { };//: Identifies matrices that can not be used with the rows and columns functions//!component: type//!category: containers, tagsstruct not_strideable { };/* the following are used in constructors to avoid ambiguity and   compiler errors *///: blah//!noindex:struct do_transpose { };//: blah//!noindex:struct do_strided { };//: blah//!noindex:struct do_scaled { };//: blah//!noindex:struct do_stream { };} /* namespace mtl */#endif /* _MTL_MATRIX_TRAITS_ */