02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz/*
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * Matrix, MatrixView, ConstMatrixView classes wrap the gsl matrix routines;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * "views" mimic the semantic of C++ references: any operation performed
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * on a "view" is actually performed on the "viewed object"
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz *
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * Authors:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * Marco Cecchetti <mrcekets at gmail.com>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz *
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * Copyright 2008 authors
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz *
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * This library is free software; you can redistribute it and/or
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * modify it either under the terms of the GNU Lesser General Public
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * License version 2.1 as published by the Free Software Foundation
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * (the "LGPL") or, at your option, under the terms of the Mozilla
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * Public License Version 1.1 (the "MPL"). If you do not alter this
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * notice, a recipient may use your version of this file under either
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * the MPL or the LGPL.
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz *
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * You should have received a copy of the LGPL along with this library
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * in the file COPYING-LGPL-2.1; if not, write to the Free Software
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * You should have received a copy of the MPL along with this library
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * in the file COPYING-MPL-1.1
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz *
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * The contents of this file are subject to the Mozilla Public License
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * Version 1.1 (the "License"); you may not use this file except in
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * compliance with the License. You may obtain a copy of the License at
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * http://www.mozilla.org/MPL/
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz *
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * This software is distributed on an "AS IS" basis, WITHOUT WARRANTY
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * OF ANY KIND, either express or implied. See the LGPL or the MPL for
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz * the specific language governing rights and limitations.
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz */
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#ifndef _NL_MATRIX_H_
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#define _NL_MATRIX_H_
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen#include <2geom/exception.h>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#include <2geom/numeric/vector.h>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#include <cassert>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#include <utility> // for std::pair
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#include <algorithm> // for std::swap
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#include <sstream>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#include <string>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#include <gsl/gsl_matrix.h>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#include <gsl/gsl_linalg.h>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruznamespace Geom { namespace NL {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruznamespace detail
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzclass BaseMatrixImpl
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz public:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz virtual ~BaseMatrixImpl()
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz ConstVectorView row_const_view(size_t i) const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return ConstVectorView(gsl_matrix_const_row(m_matrix, i));
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz ConstVectorView column_const_view(size_t i) const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return ConstVectorView(gsl_matrix_const_column(m_matrix, i));
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz const double & operator() (size_t i, size_t j) const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return *gsl_matrix_const_ptr(m_matrix, i, j);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz const gsl_matrix* get_gsl_matrix() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return m_matrix;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz bool is_zero() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return gsl_matrix_isnull(m_matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz bool is_positive() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti for ( unsigned int i = 0; i < rows(); ++i )
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti {
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti for ( unsigned int j = 0; j < columns(); ++j )
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti {
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti if ( (*this)(i,j) <= 0 ) return false;
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti }
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti }
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti return true;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz bool is_negative() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti for ( unsigned int i = 0; i < rows(); ++i )
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti {
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti for ( unsigned int j = 0; j < columns(); ++j )
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti {
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti if ( (*this)(i,j) >= 0 ) return false;
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti }
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti }
40f90655923d3f6ab288055fdb6641b45066e8c9mcecchetti return true;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz bool is_non_negative() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz for ( unsigned int i = 0; i < rows(); ++i )
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz for ( unsigned int j = 0; j < columns(); ++j )
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz if ( (*this)(i,j) < 0 ) return false;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return true;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz double max() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return gsl_matrix_max(m_matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz double min() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return gsl_matrix_min(m_matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz std::pair<size_t, size_t>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz max_index() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz std::pair<size_t, size_t> indices;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_max_index(m_matrix, &(indices.first), &(indices.second));
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return indices;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz std::pair<size_t, size_t>
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz min_index() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz std::pair<size_t, size_t> indices;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_min_index(m_matrix, &(indices.first), &(indices.second));
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return indices;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz size_t rows() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return m_rows;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz size_t columns() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return m_columns;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz std::string str() const;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz protected:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz size_t m_rows, m_columns;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix* m_matrix;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}; // end class BaseMatrixImpl
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzinline
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzbool operator== (BaseMatrixImpl const& m1, BaseMatrixImpl const& m2)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz if (m1.rows() != m2.rows() || m1.columns() != m2.columns()) return false;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz for (size_t i = 0; i < m1.rows(); ++i)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz for (size_t j = 0; j < m1.columns(); ++j)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz if (m1(i,j) != m2(i,j)) return false;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return true;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruztemplate< class charT >
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzinline
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzstd::basic_ostream<charT> &
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzoperator<< (std::basic_ostream<charT> & os, const BaseMatrixImpl & _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz if (_matrix.rows() == 0 || _matrix.columns() == 0) return os;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz os << "[[" << _matrix(0,0);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz for (size_t j = 1; j < _matrix.columns(); ++j)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz os << ", " << _matrix(0,j);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz os << "]";
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz for (size_t i = 1; i < _matrix.rows(); ++i)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz os << ", [" << _matrix(i,0);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz for (size_t j = 1; j < _matrix.columns(); ++j)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz os << ", " << _matrix(i,j);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz os << "]";
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz os << "]";
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return os;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzinline
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzstd::string BaseMatrixImpl::str() const
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz std::ostringstream oss;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz oss << (*this);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return oss.str();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzclass MatrixImpl : public BaseMatrixImpl
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz public:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz typedef BaseMatrixImpl base_type;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz void set_all( double x )
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_set_all(m_matrix, x);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz void set_identity()
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_set_identity(m_matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen using base_type::operator(); // VSC legacy support
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen const double & operator() (size_t i, size_t j) const
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen {
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen return base_type::operator ()(i, j);
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen double & operator() (size_t i, size_t j)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
98a704c566ce5a750d76d5fc9675ccc804ac65f5johanengelen return *gsl_matrix_ptr(m_matrix, i, j);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz using base_type::get_gsl_matrix;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix* get_gsl_matrix()
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return m_matrix;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz VectorView row_view(size_t i)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return VectorView(gsl_matrix_row(m_matrix, i));
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz VectorView column_view(size_t i)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return VectorView(gsl_matrix_column(m_matrix, i));
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz void swap_rows(size_t i, size_t j)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_swap_rows(m_matrix, i, j);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz void swap_columns(size_t i, size_t j)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_swap_columns(m_matrix, i, j);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixImpl & transpose()
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz assert(columns() == rows());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_transpose(m_matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return (*this);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixImpl & scale(double x)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_scale(m_matrix, x);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return (*this);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixImpl & translate(double x)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_add_constant(m_matrix, x);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return (*this);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixImpl & operator+=(base_type const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_add(m_matrix, _matrix.get_gsl_matrix());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return (*this);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixImpl & operator-=(base_type const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_sub(m_matrix, _matrix.get_gsl_matrix());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return (*this);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}; // end class MatrixImpl
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz} // end namespace detail
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzusing detail::operator==;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzusing detail::operator<<;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelentemplate <size_t N>
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelenclass ConstBaseSymmetricMatrix;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzclass Matrix: public detail::MatrixImpl
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz public:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz typedef detail::MatrixImpl base_type;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz public:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz // the matrix is not inizialized
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix(size_t n1, size_t n2)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = n1;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = n2;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = gsl_matrix_alloc(n1, n2);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix(size_t n1, size_t n2, double x)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = n1;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = n2;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = gsl_matrix_alloc(n1, n2);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_set_all(m_matrix, x);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix(Matrix const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz : base_type()
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = _matrix.rows();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = _matrix.columns();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = gsl_matrix_alloc(rows(), columns());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_memcpy(m_matrix, _matrix.get_gsl_matrix());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz explicit
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix(base_type::base_type const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = _matrix.rows();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = _matrix.columns();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = gsl_matrix_alloc(rows(), columns());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_memcpy(m_matrix, _matrix.get_gsl_matrix());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen template <size_t N>
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen explicit
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen Matrix(ConstBaseSymmetricMatrix<N> const& _smatrix)
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen {
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen m_rows = N;
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen m_columns = N;
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen m_matrix = gsl_matrix_alloc(N, N);
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen for (size_t i = 0; i < N; ++i)
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen for (size_t j = 0; j < N ; ++j)
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen (*gsl_matrix_ptr(m_matrix, i, j)) = _smatrix(i,j);
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen }
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix & operator=(Matrix const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz assert( rows() == _matrix.rows() && columns() == _matrix.columns() );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_memcpy(m_matrix, _matrix.get_gsl_matrix());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return *this;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix & operator=(base_type::base_type const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz assert( rows() == _matrix.rows() && columns() == _matrix.columns() );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_memcpy(m_matrix, _matrix.get_gsl_matrix());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return *this;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen template <size_t N>
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen Matrix & operator=(ConstBaseSymmetricMatrix<N> const& _smatrix)
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen {
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen assert (rows() == N && columns() == N);
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen for (size_t i = 0; i < N; ++i)
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen for (size_t j = 0; j < N ; ++j)
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen (*this)(i,j) = _smatrix(i,j);
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen return *this;
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen }
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz virtual ~Matrix()
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_free(m_matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix & transpose()
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<Matrix &>( base_type::transpose() );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix & scale(double x)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<Matrix &>( base_type::scale(x) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix & translate(double x)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<Matrix &>( base_type::translate(x) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix & operator+=(base_type::base_type const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<Matrix &>( base_type::operator+=(_matrix) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Matrix & operator-=(base_type::base_type const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<Matrix &>( base_type::operator-=(_matrix) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz friend
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz void swap(Matrix & m1, Matrix & m2);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz friend
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz void swap_any(Matrix & m1, Matrix & m2);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}; // end class Matrix
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz// warning! this operation invalidates any view of the passed matrix objects
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzinline
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzvoid swap(Matrix & m1, Matrix & m2)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński assert(m1.rows() == m2.rows() && m1.columns() == m2.columns());
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński using std::swap;
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński swap(m1.m_matrix, m2.m_matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosińskiinline void swap_any(Matrix &m1, Matrix &m2)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński using std::swap;
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński swap(m1.m_matrix, m2.m_matrix);
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński swap(m1.m_rows, m2.m_rows);
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński swap(m1.m_columns, m2.m_columns);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzclass ConstMatrixView : public detail::BaseMatrixImpl
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz public:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz typedef detail::BaseMatrixImpl base_type;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz public:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz ConstMatrixView(const base_type & _matrix, size_t k1, size_t k2, size_t n1, size_t n2)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz : m_matrix_view( gsl_matrix_const_submatrix(_matrix.get_gsl_matrix(), k1, k2, n1, n2) )
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = n1;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = n2;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = const_cast<gsl_matrix*>( &(m_matrix_view.matrix) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz ConstMatrixView(const ConstMatrixView & _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz : base_type(),
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix_view(_matrix.m_matrix_view)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = _matrix.rows();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = _matrix.columns();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = const_cast<gsl_matrix*>( &(m_matrix_view.matrix) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz ConstMatrixView(const base_type & _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz : m_matrix_view(gsl_matrix_const_submatrix(_matrix.get_gsl_matrix(), 0, 0, _matrix.rows(), _matrix.columns()))
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = _matrix.rows();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = _matrix.columns();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = const_cast<gsl_matrix*>( &(m_matrix_view.matrix) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz private:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_const_view m_matrix_view;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}; // end class ConstMatrixView
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzclass MatrixView : public detail::MatrixImpl
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz public:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz typedef detail::MatrixImpl base_type;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz public:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView(base_type & _matrix, size_t k1, size_t k2, size_t n1, size_t n2)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = n1;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = n2;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix_view
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz = gsl_matrix_submatrix(_matrix.get_gsl_matrix(), k1, k2, n1, n2);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = &(m_matrix_view.matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView(const MatrixView & _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz : base_type()
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = _matrix.rows();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = _matrix.columns();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix_view = _matrix.m_matrix_view;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = &(m_matrix_view.matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView(Matrix & _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_rows = _matrix.rows();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_columns = _matrix.columns();
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix_view
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz = gsl_matrix_submatrix(_matrix.get_gsl_matrix(), 0, 0, rows(), columns());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz m_matrix = &(m_matrix_view.matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView & operator=(MatrixView const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz assert( rows() == _matrix.rows() && columns() == _matrix.columns() );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_memcpy(m_matrix, _matrix.m_matrix);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return *this;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView & operator=(base_type::base_type const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz assert( rows() == _matrix.rows() && columns() == _matrix.columns() );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_memcpy(m_matrix, _matrix.get_gsl_matrix());
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return *this;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView & transpose()
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<MatrixView &>( base_type::transpose() );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView & scale(double x)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<MatrixView &>( base_type::scale(x) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView & translate(double x)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<MatrixView &>( base_type::translate(x) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView & operator+=(base_type::base_type const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<MatrixView &>( base_type::operator+=(_matrix) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz MatrixView & operator-=(base_type::base_type const& _matrix)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz {
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz return static_cast<MatrixView &>( base_type::operator-=(_matrix) );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz }
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz friend
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz void swap_view(MatrixView & m1, MatrixView & m2);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz private:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz gsl_matrix_view m_matrix_view;
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}; // end class MatrixView
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzinline
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzvoid swap_view(MatrixView & m1, MatrixView & m2)
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz{
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński assert(m1.rows() == m2.rows() && m1.columns() == m2.columns());
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński using std::swap;
76addc201c409e81eaaa73fe27cc0f79c4db097cKrzysztof Kosiński swap(m1.m_matrix_view, m2.m_matrix_view);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz}
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzVector operator*( detail::BaseMatrixImpl const& A,
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz detail::BaseVectorImpl const& v );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzMatrix operator*( detail::BaseMatrixImpl const& A,
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz detail::BaseMatrixImpl const& B );
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruzMatrix pseudo_inverse(detail::BaseMatrixImpl const& A);
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelendouble trace (detail::BaseMatrixImpl const& A);
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelendouble det (detail::BaseMatrixImpl const& A);
d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz} } // end namespaces
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz#endif /*_NL_MATRIX_H_*/
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz/*
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz Local Variables:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz mode:c++
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz c-file-style:"stroustrup"
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz c-file-offsets:((innamespace . 0)(inline-open . 0)(case-label . +))
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz indent-tabs-mode:nil
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz fill-column:99
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz End:
02c55b36fc7069d6bc02a8f1bbf2e1c58428e21ajoncruz*/
a4030d5ca449e7e384bc699cd249ee704faaeab0Chris Morgan// vim: filetype=cpp:expandtab:shiftwidth=4:tabstop=8:softtabstop=4:fileencoding=utf-8:textwidth=99 :