// RefMatrix.h copyright (C) 2007  adolfo@di-mare.com

/** \file  RefMatrix.h
    \brief Muestra cómo implementar una matriz usando referencias \c lkptr<X>.

    \author Adolfo Di Mare <adolfo@di-mare.com>
    \date   2007

    - Why English names??? ==> http://www.di-mare.com/adolfo/binder/c01.htm#sc04
*/

#ifndef   RefMatrix_h
#define   RefMatrix_h

#include <cassert> // assert()

#include "lkptr.h" // javify
#include "Matrix.h"

/// Definido por la biblioteca C++ estándar
namespace std {} // Está acá para que Doxygen lo documente

/// Matriz chirrisquitica de adolfo@di-mare.com
namespace Mx {

/// Empaque alrededor de \c Matrix que usa \c lkptr<X> para mejorar la implementación de los operadores aritméticos.
template <class E>
class RefMatrix : public Mx::Matrix<E> {
public:
    RefMatrix(unsigned m = 1, unsigned n = 1) : Matrix<E>(m,n) { } ///< Constructor base
    RefMatrix(const RefMatrix& o)  : Matrix<E>(o) { } ///< Constructor de copia
    RefMatrix(const E& V) : Matrix<E>(V) { } ///< Construye una matriz unitaria
    ~RefMatrix() {}///< Destructor
public:
    /// Retorna un puntero a una copia de la instancia.
    RefMatrix<E>* clone() const { return new RefMatrix<E>(*this); }
    friend lkptr<RefMatrix> operator + (const lkptr<RefMatrix>& A, const lkptr<RefMatrix>& B)
    { lkptr<RefMatrix> Res ( A->clone() ); Res->add(*B); return Res; } ///< Retorna \c A+B
    friend lkptr<RefMatrix> operator - (const lkptr<RefMatrix>& A, const lkptr<RefMatrix>& B)
    { lkptr<RefMatrix> Res ( A->clone() ); Res->substract(*B); return Res; } ///< Retorna \c A-B
    friend lkptr<RefMatrix> operator * (const lkptr<RefMatrix>& A, const lkptr<RefMatrix>& B)
    { RefMatrix *Res = new RefMatrix( A->rows(), B->cols() );
      Res->multiply(*A, *B); return lkptr<RefMatrix>( Res ); } ///< Retorna \c A*B
#if 0
    {   // Hay que usar: [ V->at(i,j) ] o [ (*V)(i,i) ] PERO...
        lkptr< RefMatrix<unsigned> > rV ( new RefMatrix<unsigned>(M,N) );
        RefMatrix<unsigned>& VV = *rV; // VV(i,j) == (*V)(i,j)
        // ...
        for (i=0; i < rV->rows(); ++i) {
            for (j=0; j < rV->cols(); ++j,++k) {
                assert(  rV->at(i,j) ==  VV(i,j) );
                assert( &rV->at(i,j) == &VV(i,j) );
            }
        }
    }
#else
//  friend E& operator()( lkptr<RefMatrix> L, unsigned i, unsigned j );
//  { return (*L)(i,j); } // Error: operator()() must be a nonstatic member function
#endif
}; // RefMatrix

} // namespace Mx

#endif // RefMatrix_h
// EOF: RefMatrix.h