solver.C

Ir a la documentación de este archivo.

//===========================================================================
//= solver.cc                                                  Mayo de 1999 =
//=-------------------------------------------------------------------------=
//= Definiciones para los SOLVERs numericos!                                =
//=-------------------------------------------------------------------------=
//= ADVERTENCIA: ESTE SOFTWARE NO ESTA CONCEBIDO NI DISENNADO PARA EL USO   =
//= EN EQUIPO DE CONTROL EN LINEA EN ENTORNOS PELIGROSOS QUE REQUIERAN UN   =
//= DESEMPENNO LIBRE DE FALLAS, COMO LA OPERACION DE PLANTAS NUCLEARES,     = 
//= SISTEMAS DE NAVEGACION O COMUNICACION EN AVIONES, TRAFICO AEREO,        =
//= EQUIPO MEDICO DEL CUAL DEPENDAN VIDAS HUMANAS O SISTEMAS DE ARMAMENTO,  =
//= EN LOS CUALES UNA FALLA EN EL SOFTWARE PUEDA IMPLICAR DIRECTAMENTE LA   =
//= MUERTE, DANNOS PERSONALES O DANNOS FISICOS Y/O AMBIENTALES GRAVES       =
//= ("ACTIVIDADES DE ALGO RIESGO").                                         =
//=-------------------------------------------------------------------------=
//= Autor original: Oscar J. Chavarro G.  A.K.A. JEDILINK. Copyright (c),   =
//= 1997 - 2003, oscarchavarro@hotmail.com                                  =
//= AQUYNZA es software libre, y se rige bajo los terminos de la licencia   =
//= LGPL de GNU (http://www.gnu.org). Para mayor informacion respecto a la  =
//= licencia de uso, consulte el archivo ./doc/LICENCIA en la distribucion. =
//===========================================================================

#include "core/simul/solver.h"

#include <stdio.h>
#include <stdlib.h>

#define VERIFICAR(x)                                                          \
  if ( !(x) ) {                                                               \
      fprintf(stderr,"<SOLVER> - ERROR FATAL: No hay memoria suficiente!\n"); \
      fflush(stderr);                                                         \
      exit(1);                                                                \
  }

//===========================================================================
//= CLASE SOLVER                                                            =
//===========================================================================

SOLVER::SOLVER(ECUACION_DIFERENCIAL *Mi_ecuacion)
{
    Ecuacion = Mi_ecuacion;

    temp1 = NULL;
    temp2 = NULL;
    num_variables_de_estado = Ecuacion->ODE_numero_de_variables_de_estado();
    if ( !num_variables_de_estado ) return;

    temp1 = new double[num_variables_de_estado];
    VERIFICAR(temp1);
    temp2 = new double[num_variables_de_estado];
    VERIFICAR(temp2);
}

SOLVER::~SOLVER()
{
    if ( temp1 ) {
        delete temp1;
        temp1 = NULL;
        delete temp2;
        temp2 = NULL;
        num_variables_de_estado = 0;
    }
}

void
SOLVER::escalar_vector(double factor, double *vector)
{
    int i;

    for ( i = 0; i < num_variables_de_estado; i++ ) {
        vector[i] *= factor;
    }
}

void
SOLVER::sumar_vectores(void)
{
    int i;

    for ( i = 0; i < num_variables_de_estado; i++ ) {
        temp2[i] += temp1[i];
    }
}

//===========================================================================
//= CLASE SOLVER_EULER                                                      =
//===========================================================================

SOLVER_EULER::SOLVER_EULER(ECUACION_DIFERENCIAL *Mi_ecuacion) :
                                                          SOLVER(Mi_ecuacion)
{
    ;
}

SOLVER_EULER::~SOLVER_EULER()
{
    // OJO: Aqui se llama al destructor de SOLVER o toca dejar esto?
    if ( temp1 ) {
        delete temp1;
        temp1 = NULL;
        delete temp2;
        temp2 = NULL;
        num_variables_de_estado = 0;
    }
}

void
SOLVER_EULER::ejecutar_solucion(double delta_t)
{
    Ecuacion->ODE_calcular_la_primera_derivada(temp1);
    escalar_vector(delta_t, temp1);
    Ecuacion->ODE_reportar_variables_de_estado(temp2);
    sumar_vectores(); // temp2 = temp2 + temp1;
    Ecuacion->ODE_actualizar_variables_de_estado(temp2);
}

//===========================================================================
//= EOF                                                                     =
//===========================================================================

---

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.