// This function defines all inline functions!

matrix::Type matrix::get( const size_t& i_row, 
			  const size_t& j_col ) const
{
  return gsl_matrix_get( m_, i_row, j_col ) ;
}


gsl_matrix* matrix::get_gsl_matrix() const
{
  return m_;
}  


void matrix::set( const size_t& i_row, 
		  const size_t& j_col, const Type& val )
{
  gsl_matrix_set( m_, i_row, j_col, val );
}


void matrix::set_all( const Type& val )
{
  gsl_matrix_set_all( m_, val );
}