matriznm.h

Ir a la documentación de este archivo.

/*
 * Universidad de los Andes
 * Departamento de sistemas y computación
 * http://ella.uniandes.edu.co
 *
 * ---------------------------------------------------------------
 * Ver documento README en la base de la distribución de libMath.
 * Historia de este archivo disponible en el log de CVS.
 * ---------------------------------------------------------------
 * 
 * Autores:
 * Sebastián González (se-gonza@uniandes.edu.co)
 */

#ifndef _MATRIXNXM_H_
#define _MATRIXNXM_H_

// OJO: Esta vaina que?
//#include <cassert>
//#include "MemoryException.h"

template <class T>
class MATRIZ_NxM
{
        unsigned numrows, numcols;
        T *entries;

public:

        // constructoras / destructoras
        MATRIZ_NxM(unsigned numrows, unsigned numcols);
        ~MATRIZ_NxM();

        // analizadoras
        unsigned numRows(void);
        unsigned numCols(void);

        // operadores
        virtual MATRIZ_NxM<T> operator +(const MATRIZ_NxM<T> &op);
        virtual MATRIZ_NxM<T> operator -(const MATRIZ_NxM<T> &op);
        virtual MATRIZ_NxM<T> operator *(const MATRIZ_NxM<T> &op);
        
        // depuradoras
        void imprimir();
};


/******************
  Implementaciones
 ******************/


template <class T>
MATRIZ_NxM<T>::MATRIZ_NxM(unsigned n, unsigned m)
: numrows(n), numcols(m)                // efficient!
{
        assert(n > 0 && m > 0);
        entries = new T[n * m];

        if(entries == NULL)
                throw new libRuntime::MemoryException();
}


template <class T>
MATRIZ_NxM<T>::~MATRIZ_NxM()
{
        if(entries != NULL)
                delete[] entries;
}


template <class T>
unsigned MATRIZ_NxM<T>::numRows(void)
{
        return numrows;
}


template <class T>
unsigned MATRIZ_NxM<T>::numCols(void)
{
        return numcols;
}


template <class T>
MATRIZ_NxM<T> MATRIZ_NxM<T>::operator +(const MATRIZ_NxM<T> &op)
{
        assert(numrows == op.numrows);
        assert(numcols == op.numcols);

    MATRIZ_NxM<T> result(numrows, numcols);

        for(unsigned i = 0; i < numcols * numrows; i++)
                result.entries[i] = entries[i] + op.entries[i];

    return result;
}


template <class T>
MATRIZ_NxM<T> MATRIZ_NxM<T>::operator -(const MATRIZ_NxM<T> &op)
{
        assert(numrows == op.numrows);
        assert(numcols == op.numcols);

    MATRIZ_NxM<T> result(numrows, numcols);

        for(unsigned i = 0; i < numcols * numrows; i++)
                result.entries[i] = entries[i] - op.entries[i];

    return result;
}


template <class T>
MATRIZ_NxM<T> MATRIZ_NxM<T>::operator *(const MATRIZ_NxM<T> &op)
{
        assert(numcols == op.numrows);

    MATRIZ_NxM<T> result(numrows, op.numcols);

        for(unsigned i = 0; i < numrows; i++)
        {
                for(unsigned j = 0; j < op.numcols; j++)
                {
                        // para cada entrada del resultado...
                        // R[i, j] = suma( A[i,k] * B[k, j] )

                        T accum = 0;

                        for(unsigned k = 0; k < numcols; k++)
                                accum += entries[i * numcols + k] * op.entries[k * op.numcols + j];

                        result.entries[i * result.numcols + j] = accum;
                }
        }

    return result;
}


void
MATRIZ_NxM::imprimir(void)
{
    int fila, columna;

    printf("=================================\n");
    for( fila = 0; fila < numrows; fila++ ) {
        for( columna = 0; columna < numcols; columna++ ) {
            printf("(%d,%d)=%.2f\t", fila, columna, entries[fila*numcols + columna]);
        }
        printf("\n");
    }
    printf("=================================\n");
    fflush(stdout);
}

#endif

---

Este archivo HTML ha sido generado automáticamente a partir del código fuente AQUYNZA. NO LO EDITE. Para mayor información contacte al autor.