ResuelveTridiagonal.hpp

Go to the documentation of this file.

// Clase para resolver un sistema lineal tridiagonal                                        //
//                                                                                          //
// Análisis y Diseño y Programación:                                                        //
//                                                                                          //
// Nombre:   Antonio Carrillo Ledesma                                                       //
// E-mail:   acl@www.mmc.geofisica.unam.mx                                                  //
// Página:   http://www.mmc.geofisica.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 __ResuelveTridiagonal__
#define __ResuelveTridiagonal__


#include "ResuelveSistemaLineal.hpp"



class ResuelveTridiagonal: public ResuelveSistemaLineal
{


public:

   ResuelveTridiagonal(void) : ResuelveSistemaLineal()
   {
      M = NULL;
      X = NULL;
      B = NULL;
      MetodoModificaMatriz = true;
      MetodoNumerico = TRIDIAGONAL;
   }


   ResuelveTridiagonal(Matriz *A) : ResuelveSistemaLineal()
   {
      if (A->maximoTamanoColumnaOcupado() > 3)
      {
         printf("\nError el tamaño de la banda es mayor que 3\n");
         exit(1);
      }
      M = A;
      X = NULL;
      B = NULL;
      MetodoModificaMatriz = true;
      MetodoNumerico = TRIDIAGONAL;
   }


   ResuelveTridiagonal(Matriz *A, Vector *x, Vector *b) : ResuelveSistemaLineal()
   {
      if (A->maximoTamanoColumnaOcupado() > 3)
      {
         printf("\nError el tamaño de la banda es mayor que 3\n");
         exit(1);
      }
      M = A;
      X = x;
      B = b;
      MetodoModificaMatriz = true;
      MetodoNumerico = TRIDIAGONAL;
   }


   void resuelve(void);



   void resuelve(Vector *x, Vector *b)
   {
      X = x;
      B = b;
      resuelve();
   }

};

#endif