2d/matrix.hh

Go to the documentation of this file.

/* -*- mia-c++  -*-
 *
 * This file is part of MIA - a toolbox for medical image analysis 
 * Copyright (c) Leipzig, Madrid 1999-2017 Gert Wollny
 *
 * MIA 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.
 *
 * This program 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 MIA; if not, see <http://www.gnu.org/licenses/>.
 *
 */

#ifndef __mia_2d_matrix_hh
#define __mia_2d_matrix_hh

#include <mia/2d/vector.hh>

NS_MIA_BEGIN

template <typename T> 
struct T2DMatrix: public T2DVector< T2DVector<T> > {
        
        T2DMatrix(); 

        T2DMatrix(const T2DMatrix<T>& o); 

        T2DMatrix(const T2DVector< T2DVector<T> >& o); 
        
        T2DMatrix(const T2DVector< T >& x, const T2DVector< T >& y ); 
        
        T2DMatrix& operator = (const T2DMatrix<T>& o);
        
        T2DMatrix& operator -= (const T2DMatrix<T>& o);

        T2DMatrix<T>  transposed()const; 


        static const T2DMatrix _1; 
        static const T2DMatrix _0; 

}; 

template <typename T> 
const T2DMatrix<T> T2DMatrix<T>::_1(T2DVector< T >(1,0), 
                                    T2DVector< T >(0,1));

template <typename T> 
const T2DMatrix<T> T2DMatrix<T>::_0;


template <typename T> 
T2DMatrix<T>::T2DMatrix():
        T2DVector<T2DVector<T> >(T2DVector<T>(), T2DVector<T>())
{
}

template <typename T> 
T2DMatrix<T>::T2DMatrix(const T2DMatrix<T>& o):
        T2DVector<T2DVector<T> >(o.x, o.y)
{
}

template <typename T> 
T2DMatrix<T>::T2DMatrix(const T2DVector< T2DVector<T> >& o):
        T2DVector<T2DVector<T> >(o.x, o.y)
{
}

template <typename T> 
T2DMatrix<T>::T2DMatrix(const T2DVector< T >& x, const T2DVector< T >& y ):
        T2DVector<T2DVector<T> >(x, y)
{
}


template <typename T> 
T2DMatrix<T>& T2DMatrix<T>::operator = (const T2DMatrix<T>& o)
{
        this->x = o.x; 
        this->y = o.y; 
        return *this; 
}

template <typename T> 
T2DMatrix<T>& T2DMatrix<T>::operator -= (const T2DMatrix<T>& o)
{
        this->x -= o.x; 
        this->y -= o.y; 
        return *this; 
}

template <typename T> 
T2DMatrix<T>  T2DMatrix<T>::transposed()const
{
        return T2DMatrix<T>(T2DVector< T >(this->x.x, this->y.x), 
                            T2DVector< T >(this->y.y, this->y.y)); 
}

template <typename T> 
T2DVector<T> operator * (const T2DVector<T>& x, const T2DMatrix<T>& m)
{
        return T2DVector<T>(dot(m.x, x), dot(m.y, x)); 
}

template <typename T> 
T2DVector<T> operator * (const T2DMatrix<T>& m, const T2DVector<T>& x  )
{
        return T2DVector<T>(m.x.x * x.x + m.y.x * x.y, 
                            m.x.y * x.x + m.y.y * x.y);
}

typedef T2DMatrix<float> C2DFMatrix; 

NS_MIA_END

#endif