jblas.hpp

Go to the documentation of this file.

/* $Id$ */

#ifndef JMATH_JBLAS_HPP
#define JMATH_JBLAS_HPP

#ifdef USE_JMATH_SERIALIZATION
#include "jmath/serialize_vector.hpp"
#include "jmath/serialize_matrix.hpp"
#include "jmath/serialize_symmetric.hpp"
#include "jmath/serialize_banded.hpp"
#else
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/symmetric.hpp>
#include <boost/numeric/ublas/banded.hpp>
#endif
#include <boost/numeric/ublas/io.hpp>
#include <boost/numeric/ublas/triangular.hpp>
#include <boost/numeric/ublas/matrix_proxy.hpp>
#include <boost/numeric/ublas/vector_proxy.hpp>
#include <boost/numeric/ublas/matrix_sparse.hpp>

#include "boost/math/quaternion.hpp"

#include "jmath/ublasCompatibility.hpp"
#include "jmath/boundedSymmetricMatrix.hpp"


//namespace ublas = boost::numeric::ublas; // don't use, swig doesn't handle it correctly

namespace jblas {

  /*
   * scalar type double
   */

  typedef boost::numeric::ublas::vector<double> vec;
  typedef boost::numeric::ublas::vector_range<vec> vec_range; 
  typedef boost::numeric::ublas::vector_slice<vec> vec_slice;
  typedef boost::numeric::ublas::vector_indirect<vec> vec_indirect;
  typedef boost::numeric::ublas::vector_range<const vec> cvec_range; 
  typedef boost::numeric::ublas::zero_vector<double> zero_vec;
  typedef boost::numeric::ublas::unit_vector<double> unit_vec;
  typedef boost::numeric::ublas::scalar_vector<double> scalar_vec;

  typedef boost::numeric::ublas::bounded_vector<double,1> vec1;
  typedef boost::numeric::ublas::bounded_vector<double,2> vec2;
  typedef boost::numeric::ublas::bounded_vector<double,3> vec3;
  typedef boost::numeric::ublas::bounded_vector<double,4> vec4;
  typedef boost::numeric::ublas::vector_range<vec4> vec4_range; 
  typedef boost::numeric::ublas::bounded_vector<double,6> vec6;
  typedef boost::numeric::ublas::bounded_vector<double,7> vec7;
  typedef boost::numeric::ublas::bounded_vector<double,8> vec8;
  typedef boost::numeric::ublas::bounded_vector<double,9> vec9;
  typedef boost::numeric::ublas::bounded_vector<double,10> vec10;
  typedef boost::numeric::ublas::bounded_vector<double,11> vec11;


  typedef boost::numeric::ublas::matrix<double> mat;
  
  typedef boost::numeric::ublas::matrix_range<mat> mat_range;
  typedef boost::numeric::ublas::matrix_slice<mat> mat_slice;
  typedef boost::numeric::ublas::matrix_indirect<mat> mat_indirect;
  typedef boost::numeric::ublas::matrix_range<const mat> cmat_range;
  typedef boost::numeric::ublas::zero_matrix<double> zero_mat;
  typedef boost::numeric::ublas::identity_matrix<double> identity_mat;
  typedef boost::numeric::ublas::matrix_vector_range<mat> mat_vec_range;
  typedef boost::numeric::ublas::scalar_matrix<double> scalar_mat;
  typedef boost::numeric::ublas::matrix<double, boost::numeric::ublas::column_major> mat_column_major;
  typedef boost::numeric::ublas::matrix_range<const mat_column_major> cmat_column_major_range;
  typedef boost::numeric::ublas::matrix_range<mat_column_major> mat_column_major_range;

  typedef boost::numeric::ublas::diagonal_matrix<double> diag_mat;

  typedef boost::numeric::ublas::bounded_matrix<double,2,2> mat22; // 2x2
  typedef boost::numeric::ublas::bounded_matrix<double,2,3> mat23; // 2x3
  typedef boost::numeric::ublas::bounded_matrix<double,3,2> mat32; // 3x2
  typedef boost::numeric::ublas::bounded_matrix<double,3,3> mat33; // 3x3
  typedef boost::numeric::ublas::matrix_range<mat33> mat33_range; //  3x3
  typedef boost::numeric::ublas::bounded_matrix<double,3,4> mat34; // 3x4
  typedef boost::numeric::ublas::matrix_range<mat34> mat34_range; //  3x4
  typedef boost::numeric::ublas::bounded_matrix<double,4,3> mat43; // 4x3
  typedef boost::numeric::ublas::matrix_range<mat43> mat43_range; //  4x3
  typedef boost::numeric::ublas::bounded_matrix<double,4,4> mat44; // 4x4
  typedef boost::numeric::ublas::matrix_range<mat44> mat44_range; //  4x4
  typedef boost::numeric::ublas::bounded_matrix<double,6,6> mat66; // 6x6
  typedef boost::numeric::ublas::matrix_range<mat66> mat66_range; //  6x6
  typedef boost::numeric::ublas::bounded_matrix<double,8,8> mat88; // 8x8
  typedef boost::numeric::ublas::matrix_range<mat88> mat88_range; //  8x8
  typedef boost::numeric::ublas::bounded_matrix<double,9,9> mat99; // 9x9
  typedef boost::numeric::ublas::matrix_range<mat99> mat99_range; //  9x9
  
  typedef boost::numeric::ublas::symmetric_matrix<double> sym_mat;
  typedef boost::numeric::ublas::matrix_range<sym_mat> sym_mat_range; 
  typedef boost::numeric::ublas::matrix_indirect<sym_mat> sym_mat_indirect;
  typedef boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::upper> up_sym_mat;
  typedef boost::numeric::ublas::matrix_range<up_sym_mat> up_sym_mat_range;
  typedef boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::lower> lo_sym_mat;

  typedef boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::column_major> sym_mat_column_major;
    
  typedef bounded_symmetric_matrix<double, 2 > sym_mat22;
  typedef bounded_symmetric_matrix<double, 3 > sym_mat33;
  typedef bounded_symmetric_matrix<double, 4 > sym_mat44;
  
  typedef boost::numeric::ublas::banded_matrix<double> banded_mat;

  typedef boost::numeric::ublas::matrix_vector_slice<mat> mat_vec_slice;
  typedef boost::numeric::ublas::matrix_vector_slice<sym_mat> sym_mat_vec_slice;

  typedef boost::numeric::ublas::matrix_column<mat> mat_column;

  typedef boost::numeric::ublas::compressed_matrix<double> sp_mat;
  
  typedef boost::numeric::ublas::indirect_array<> ind_array;

  typedef boost::math::quaternion<double> quat;

  /*
   * Scalar type float
   */

  typedef boost::numeric::ublas::vector<float> vecf;
  typedef boost::numeric::ublas::scalar_vector<float> scalar_vecf;

  typedef boost::numeric::ublas::bounded_vector<float,1> vecf1;
  typedef boost::numeric::ublas::bounded_vector<float,2> vecf2;
  typedef boost::numeric::ublas::bounded_vector<float,3> vecf3;
  typedef boost::numeric::ublas::bounded_vector<float,4> vecf4;
  typedef boost::numeric::ublas::bounded_vector<float,6> vecf6;
  typedef boost::numeric::ublas::bounded_vector<float,8> vecf8;

  /*
   * scalar type int
   */
  typedef boost::numeric::ublas::vector<int> veci;
  typedef boost::numeric::ublas::bounded_vector<int,2> veci2;
  typedef boost::numeric::ublas::bounded_vector<int,3> veci3;
  typedef boost::numeric::ublas::bounded_vector<int,4> veci4;
  
  typedef boost::numeric::ublas::zero_vector<int> zero_veci;
  typedef boost::numeric::ublas::matrix<int> mati;
  typedef boost::numeric::ublas::zero_matrix<int> zero_mati;
  
  /*
   * scalar type bool
   */
  
  typedef boost::numeric::ublas::matrix<bool> bool_mat;

  typedef boost::numeric::ublas::matrix<bool> matb;
  typedef boost::numeric::ublas::vector<bool> vecb;

  typedef boost::numeric::ublas::upper upper;
  typedef boost::numeric::ublas::lower lower;
  
  typedef boost::numeric::ublas::symmetric_adaptor<const jblas::mat       > sym_adapt;
  typedef boost::numeric::ublas::symmetric_adaptor<const jblas::mat, upper> up_sym_adapt;
  typedef boost::numeric::ublas::symmetric_adaptor<const jblas::mat, lower> lo_sym_adapt;
  typedef boost::numeric::ublas::symmetric_adaptor<jblas::mat_column_major        > sym_adapt_column_major; 
  typedef boost::numeric::ublas::symmetric_adaptor<jblas::mat_column_major, upper > up_sym_adapt_column_major; 
  typedef boost::numeric::ublas::symmetric_adaptor<jblas::mat_column_major, lower > lo_sym_adapt_column_major; 
  
  typedef boost::numeric::ublas::matrix_row< jblas::mat_column_major > mat_colum_major_row;
} // namespace jblas

/* End of Doxygen group */

#endif // JMATH_JBLAS_HPP