esfera.h

Ir a la documentación de este archivo.

//===========================================================================
//= esfera.h                                                  Marzo de 1999 =
//=-------------------------------------------------------------------------=
//= Definiciones de la clase ESFERA.                                        =
//= NOTA: No incluya este encabezado en los modulos.  Incluya a             =
//=       geometria.h.                                                      =
//=-------------------------------------------------------------------------=
//= 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_defs.h"

#ifndef __GEOMETRIA__
    #error "No incluya a esfera.h, incluya a geometria.h!"
#endif

#ifdef VEL_ROSITA
    #include "toolkits/geom/geometria.h"
#endif

class ESFERA : public PRIMITIVA_GEOMETRICA 
{
  private:
    double _radio;
    unsigned long int _lod;
    int _paralelos;
    int _meridianos;
#ifdef GL_ENABLED
    GLUquadricObj *Textura_glu; // OJO: No usar GLU, copy-pastear algo mejor...
#endif
    IMAGEN *imagen;

    void activar_lod(void);
  public:
    ESFERA(double r);
    virtual ~ESFERA();

    // Servicios para raytracers y algoritmos especiales de geometria
    double interseccion(RAYO *Rayo, VECTOR &punto, VECTOR &normal);

    // Otros
    double radio(void);
    void set_radio(double r);
    void set_lod(unsigned long int lod);
    void anexar_textura(IMAGEN *img);
    int clasificar_punto(VECTOR p);
    virtual GEOMETRIA *crear_copia(void);

    void minmax(VECTOR *min, VECTOR *max);

#ifdef GL_ENABLED
    void
    pintar_gl(CALIDAD_VISUAL *Calidad, MATERIAL* Material, CAMARA *Camara);
#endif
    void pintar_povray(FILE *fd);
    void pintar_aqz(FILE *fd);
};

//===========================================================================
//= Metodos inline de la clase ESFERA                                       =
//===========================================================================

inline double
ESFERA::interseccion(RAYO *Rayo, VECTOR &punto, VECTOR &normal)
{
    double B, C, discriminante, t0, t1;

    //- No se necesita calcular el termino A ya que el rayo.d es unitario ---
    // OJO: Asercion valida si se cumple la precondicion Q
    ; // A = 1; (dx^2 + dy^2 + dz^2)

    //- Calcula el termino B ------------------------------------------------
    B = 2 *
        ((Rayo->direccion.x * Rayo->origen.x) +
         (Rayo->direccion.y * Rayo->origen.y) +
         (Rayo->direccion.z * Rayo->origen.z));

    //- Calcula el termino C ------------------------------------------------
    C = SQUARE(Rayo->origen.x) +
        SQUARE(Rayo->origen.y) +
        SQUARE(Rayo->origen.z) -
        SQUARE(_radio);           // OJO: Precalcular radio^2 por eficiencia!

    //- Calcula el discriminante. Si el discriminante no es positivo el -----
    //- rayo no intersecta la esfera. retorna t = 0                         -
    discriminante = SQUARE(B) - 4*C;
    if ( discriminante <= EPSILON ) return 0.0;

    //- Resuelve la ecuacion cuadratica para las raices de la ecuacion. -----
    //- (-B +/- sqrt(B^2 - 4*A*C)) / 2A.                                    -
    discriminante = sqrt(discriminante);
    t0 = (-B-discriminante) * 0.5;

    // OJO: Sera que siempre se cumple que t0 < t1 sin calcular a t1?

    //- Si t0 es > 0 listo. Si no debemos calcular la otra raiz t1. ---------
    if ( t0 > EPSILON ) {
        // OJO: Aqui va el calculo del punto y la normal!
        punto = Rayo->origen + Rayo->direccion * t0;
        normal = punto;
        normal.normalizar();
        return t0;
    }

    t1 = (-B+discriminante) * 0.5;
    if ( t1 > EPSILON ) {
        // OJO: Aqui va el calculo del punto y la normal!
        punto = Rayo->origen + Rayo->direccion * t1;
        normal = punto;
    normal.normalizar();
        // OJO: Aqui se retorna una normal que va hacia AFUERA de la
        //      esfera, y se tiene el caso en el que el rayo se origina
        //      dentro de la esfera.
        return t1;
    }

    return 0.0;
}

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