ResuelveSistemaLineal.hpp

Go to the documentation of this file.

// Clase para resolver un sistema lineal                                                    //
// Revisión: 1 de Junio del 2006                                                            //
//                                                                                          //
//                                                                                          //
// 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 __ResuelveSistemaLineal__
#define __ResuelveSistemaLineal__


#include <stdlib.h>
#include <stdio.h>
#include "Matriz.hpp"

// Método numérico
#define CGM            0
#define FACT_LU        1
#define GAUSSSEIDEL    2
#define INVERSA        3
#define JACOBI         4
#define TRIDIAGONAL    5
#define FACT_CHOLESKI  6


// Nombre del método
static const char Nmb[][50]= {"CGM","Factorización LU","Gauss-Seidel","Inversa","Jacobi","Tridiagonal","Choleski"};

// Requiere matriz del tipo Bandada, Dispersa, Densa
#define REQUIERE_MAT_BAND   0
#define REQUIERE_MAT_DISP   1
#define REQUIERE_MAT_DENS   2



class ResuelveSistemaLineal
{
protected:

   bool MetodoModificaMatriz;
   int NumIteraciones;
   int MetodoNumerico;
   int RequiereMatriz;

   Vector *X;
   Vector *B;
   Matriz *M;


   void error(const char *txt)
   {
      printf("\nError resolver el sistema Lineal: %s\n\n",txt);
      exit(1);
   }


public:


   void inicializa(Matriz *m, Vector *b, Vector *x)
   {
      M = m;
      X = x;
      B = b;
   }

   ResuelveSistemaLineal(void)
   {
      NumIteraciones = 0;
   }

   virtual void resuelve(void) = 0;


   virtual void resuelve(Vector *b, Vector *x) = 0;
   
   

   int retornaNumeroIteraciones(void)
   {
      return NumIteraciones;
   }


   bool elMetodoModificaMatriz(void)
   {
      return MetodoModificaMatriz;
   }


   int retornaMetodoNumerico(void)
   {
      return MetodoNumerico;
   }


   int retornaRequiereMatriz(void)
   {
      return RequiereMatriz;
   }

   void informacionMetodo(void)
   {
      if (MetodoModificaMatriz) printf("\nMétodo %s para resolver el sistema lineal\n",Nmb[MetodoNumerico]);
      else printf("\nMétodo %s, iteraciones para resolver el sistema lineal %d\n",Nmb[MetodoNumerico],NumIteraciones);
   }


};

#endif