masa_puntual.C

Ir a la documentación de este archivo.

//===========================================================================
//= masa_puntual.cc                                           Julio de 1998 =
//=-------------------------------------------------------------------------=
//= Definiciones de la clase COSA                                           =
//=-------------------------------------------------------------------------=
//= 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 "jed_gl.h"
#include "core/cosas/cosa.h"
#include "lista.cc"
#include "toolkits/media/jed_img.h"
#include "toolkits/geom/mesh.h"
#include "toolkits/util/3dsload.h"
#include "matriz4.h"

#include <ctype.h>
#include <string.h>
#include <stdlib.h>

//===========================================================================
//= VARIABLES GLOBALES Y MACROS                                             =
//===========================================================================

#define ESPERO(tipo, msg) \
    if ( tipo_token != (tipo) ) { \
        fprintf(stderr, "<COSA/MASA_PUNTUAL> ERROR: Esperaba %s y recibi [%s].\n", \
          (msg), cad); fflush(stderr);  return FALSE; \
    }

//===========================================================================
//= Clase MASA_PUNTUAL                                                      =
//===========================================================================

MASA_PUNTUAL::MASA_PUNTUAL() : COSA()
{
    _tipo_de_cosa = CC_MASA_PUNTUAL;
    en_colision = FALSE;
}

MASA_PUNTUAL::~MASA_PUNTUAL()
{
    if ( _nombre_padre ) delete _nombre_padre;
}

BOOLEAN
MASA_PUNTUAL::consultar_variable(const char *nombre_variable, int &tipo, void **ref)
{
    // Deberia responder a solicitudes de orientacion con un quaternion por
    // defecto de "orientacion nula" o algo asi.  Eso permitiria que las
    // camaras pudiesen asociarse a masas puntuales...

    if ( strcmp(nombre_variable, "centro_de_masa") == 0 ) {
        tipo = T_VECTOR;    (*ref) = &_posicion;    return TRUE;
      }
      else {
        return COSA::consultar_variable(nombre_variable, tipo, ref);
      }
    ;
}

COSA *
MASA_PUNTUAL::crear_copia(void)
{
    MASA_PUNTUAL *m;

    m = new MASA_PUNTUAL;
    if ( !m ) return NULL;
    (*m) = (*this);
    return (COSA *)m;
}

GEOMETRIA *
MASA_PUNTUAL::geometria(void)
{
    return NULL;
}

void
MASA_PUNTUAL::minmax(VECTOR *min, VECTOR *max)
{
    min->x = _posicion.x - EPSILON;
    min->y = _posicion.y - EPSILON;
    min->z = _posicion.z - EPSILON;
    max->x = _posicion.x + EPSILON;
    max->y = _posicion.y + EPSILON;
    max->z = _posicion.z + EPSILON;
}

#ifdef GL_ENABLED
void
MASA_PUNTUAL::pintar_gl(CALIDAD_VISUAL * /*Calidad*/, CAMARA * /*Camara*/)
{
    glDisable(GL_LIGHTING);
    glColor3f(_color.r, _color.g, _color.b);
    glPointSize(4);
    glBegin(GL_POINTS);
        glVertex3d(_posicion.x, _posicion.y, _posicion.z);
    glEnd();

    #define L 0.05

    if ( en_colision ) {
        VECTOR mmin(-L, -L, -L), mmax(L, L, L);

        glPushMatrix();
        glTranslatef((float)_posicion.x,
                     (float)_posicion.y, (float)_posicion.z);
        pintar_paralelepipedo(mmin, mmax);
        en_colision = FALSE;
        glPopMatrix();
    }
}
#endif // GL_ENABLED

void
MASA_PUNTUAL::pintar_povray(FILE *fd)
{
    fprintf(fd, 
        "sphere {\n"
        "    <%f, %f, %f>, 0.05\n"
        "    texture { pigment { color rgb <%f, %f, %f> } }"
        "}\n",
        _posicion.x, _posicion.y, _posicion.z,
        _color.r, _color.g, _color.b
    );
}

void
MASA_PUNTUAL::actualizar(double delta_t)
{
    SOLVER_EULER solver(this);

    solver.ejecutar_solucion(delta_t);
}

void
MASA_PUNTUAL::grabar(FILE *fd)
{
    fprintf(fd, "// COSA (MASA_PUNTUAL sin metodo de grabar!)\n");
}

BOOLEAN
MASA_PUNTUAL::leer(TOKENIZADOR *Sabiondo)
{
    char cad[1000];
    int tipo_token = TK_DESCONOCIDO, pos;
    BOOLEAN con_mi_color = FALSE;

    pos = 1;
    while ( tipo_token != TK_CERRAR ) {
        tipo_token = Sabiondo->siguiente_token(cad);
        switch ( pos ) {
            case 1:
              ESPERO(TK_CADENA, "una cadena");
              if ( strlen(cad) > MAX_CAD-1 ) cad[MAX_CAD-1] = '\0'; 
              des_comille(cad);
              set_nombre(cad);
              pos++;
              break;
            case 2:  ESPERO(TK_ABRIR, "\"{\"");  pos++; break;
            default:
              if ( tipo_token == TK_CERRAR ) break;
              ESPERO(TK_IDENTIFICADOR, "un identificador (1)");
              if ( strcmp(cad, "espacio") == 0 ) {
                  tipo_token = Sabiondo->siguiente_token(cad);
                  ESPERO(TK_CADENA, "una cadena");
                  // OJO: Es responsabilidad del UNIVERSO asignarme luego
                  // mi ESPACIO
                  //printf("  - espacio: %s\n", cad); 
                }
                else if ( !leer_basico(Sabiondo, cad, &con_mi_color) ) {
                  return FALSE;
                }
              ;
              break;
        }
    }

    if ( !con_mi_color ) _color = Material->difusa();

    return TRUE;
}

//= SERVICIOS DE ECUACION DIFERENCIAL! ======================================

int
MASA_PUNTUAL::ODE_numero_de_variables_de_estado(void)
{
    if ( _estoy_fijo ) return 0;
    return 6;
}

void
MASA_PUNTUAL::ODE_reportar_variables_de_estado(double *Datos)
{
    if ( _estoy_fijo ) return;

    Datos[0] = _posicion.x;
    Datos[1] = _posicion.y;
    Datos[2] = _posicion.z;
    Datos[3] = _velocidad.x;
    Datos[4] = _velocidad.y;
    Datos[5] = _velocidad.z;
}

void
MASA_PUNTUAL::ODE_actualizar_variables_de_estado(double *Datos)
{
    if ( _estoy_fijo ) return;

    _posicion.x = Datos[0];
    _posicion.y = Datos[1];
    _posicion.z = Datos[2];
    _velocidad.x = Datos[3];
    _velocidad.y = Datos[4];
    _velocidad.z = Datos[5];
}

void
MASA_PUNTUAL::ODE_calcular_la_primera_derivada(double *Datos)
{
    //- 1. CALCULE LA CONTRIBUCION TOTAL DE LAS FUERZAS ---------------------
    VECTOR f(0, 0, 0);
    int i;

    for ( i = 0; i < lista_fuerzas_externas.tam(); i++ ) {
        f = f + lista_fuerzas_externas[i]->fuerza;
    }
    limpiar_fuerzas();

    if ( _estoy_fijo ) return;

    //- 2. ACTUALICE LA POSICION Y LA VELOCIDAD DE LA PARTICULA -------------
    Datos[0] = _velocidad.x;
    Datos[1] = _velocidad.y;
    Datos[2] = _velocidad.z;
    Datos[3] = f.x / _masa;
    Datos[4] = f.y / _masa;
    Datos[5] = f.z / _masa;
}

//===========================================================================
//= 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.