MatrizDensa.hpp

Go to the documentation of this file.

// Clase para el trabajar con matrices densas de punto flotante                             //
//                                                                                          //
// Análisis y Diseño y Programación:                                                        //
//                                                                                          //
// Nombre:   Antonio Carrillo Ledesma                                                       //
// E-mail:   acl@www.mmc.igeofcu.unam.mx                                                    //
// Página:   http://www.mmc.igeofcu.unam.mx/acl                                             //
//                                                                                          //
//                                                                                          //
// Este programa es software libre. Puede redistribuirlo y/o modificarlo                    //
// bajo los términos de la Licencia Pública General de GNU según es                         //
// publicada por la Free Software Foundation, bien de la versión 2 de                       //
// dicha Licencia o bien (según su elección) de cualquier versión                           //
// posterior.                                                                               //
//                                                                                          //
// Este programa se distribuye con la esperanza de que sea útil, pero SIN                   //
// NINGUNA GARANTÍA, incluso sin la garantía MERCANTIL implícita o sin                      //
// garantizar la CONVENIENCIA PARA UN PROPÓSITO PARTICULAR. Véase la                        //
// Licencia Pública General de GNU para más detalles.                                       //
//                                                                                          //
// Debería haber recibido una copia de la Licencia Pública General junto                    //
// con este programa. Si no ha sido así, escriba a la Free Software                         //
// Foundation, Inc., en 675 Mass Ave, Cambridge, MA 02139, EEUU.                            //
//                                                                                          //
//                                                                                          //



#ifndef __MatrizDensa__
#define __MatrizDensa__

#include "Matriz.hpp"




class MatrizDensa: public Matriz
{
private:

   ldouble **M;
   int Estado;



   void solicitaMemoria(const int ren, const int col);



public:

   MatrizDensa(void) : Matriz()
   {
      Estado = 0;
      solicitaMemoria(1, 1);
      asignaNombre("SinNombre");
   }


   MatrizDensa(MatrizDensa &mat) : Matriz()
   {
      solicitaMemoria(mat.renglones(),mat.columnas());
      asignaNombre(mat.nombre());
      int i, j;
      for (i = 0; i < Col; i++)
      {
         for (j = 0; j < Ren; j++) M[i][j] = mat.retorna(i,j);
      }
   }


   MatrizDensa(const int ren, const int col) : Matriz()
   {
      Estado = 0;
      solicitaMemoria(ren, col);
      asignaNombre("SinNombre");
   }


   MatrizDensa(const int ren, const int col, const char *nmb) : Matriz()
   {
      Estado = 0;
      solicitaMemoria(ren, col);
      asignaNombre(nmb);
   }

   ~MatrizDensa()
   {
      liberaMemoria();
   }

   void liberaMemoria(void);


   void redimensiona(int nren, int ncol)
   {
      liberaMemoria();
      solicitaMemoria(nren, ncol);
      inicializa(0.0);
   }


   void inicializa(const ldouble val);


#ifdef DEPURAR

   void asigna(const int ren, const int col, const ldouble val);


   ldouble retorna(const int ren, const int col);
#else

   inline void asigna(const int ren, const int col, const ldouble val)
   {
      M[ren][col] = val;
   }


   inline ldouble retorna(const int ren, const int col)
   {
      return M[ren][col];
   }
#endif



   void multiplica(ldouble esc);


   void multiplica(MatrizDensa *a, MatrizDensa *b);


   void multiplica(Vector *b, Vector *r);


   int retornaNumeroColumnasBanda(int ren)
   {
      printf("\nFuncion no implementada en esta clase %d\n",ren);
      return 0;
   }


   int retornaNumeroColumna(int ren, int ind)
   {
      printf("\nFuncion no implementada en esta clase: %d %d\n",ren, ind);
      return 0;
   }


   inline double retornaValorColumna(int ren, int ind)
   {
      printf("\nFuncion no implementada en esta clase: %d %d\n",ren, ind);
      return 0;
   }

};

#endif