source: trunk/source/geometry/solids/test/X11RayTracer/main.c @ 1350

#include "main.h"
#include "torus.h"
#include <signal.h>
#include <unistd.h>
#include <assert.h>
#include <math.h>

double cos(double x);
double sin(double x);

#define EPSILON 1e-8
#define NBCERCLES 36
#define NBPOINTS  36

#define RAYON 100

int NbAff = 0;
double *px,*py,*pz;

void alarme (int sig)
{
                printf("NbAff = %d\n",NbAff);
                NbAff = 0;
                alarm (10);
}


int main (int argc,char *argv[]) {

  if (argc > 1) {
    printf("%s\n",argv[0]);
    printf("%s\n",argv[1]);
    if (strcmp(argv[1],"BVM") == 0) { 
      printf("BVM only\n");
      BVM_ONLY = 1;
    }
  }     

  /* Init X11 for 24 bpp */

  init_x ();

  /* Set up the alarm */
  
  signal (SIGALRM,alarme);
  alarm (10);

  while (event_x()) {                                           
    refresh ();
    aff ();
  }

  for(;;);              
  /* Free */

  free (px);
  free (py);
  free (pz);

  /* Close X */
  
  close_x();
  
  /* End */

  exit (0);
}

void aff () {
                NbAff ++;
                XPutImage(dis,win,gc,xim,0,0,0,0,W,H);

}

void refresh () {
  int i,j;
  long int V,R ;
  double d;
  static int t=1;

  double zoom = 0.20 + 0.01*t;
  double theta = 0.15 + 0.05*t;

  Intersect Inter ;

  memset ( buffer , 0 , 4*W*H );
  
  for (j=0;j<H;j++)
    for (i=0;i<W;i++)
      {
        //fprintf(stderr,"%d %d\n",i,j);
        
        R = 0;

        Inter.dx = -1.0*cos(theta);
        Inter.dy = 0.0;
        Inter.dz = -1.0*sin(theta);
                
        Inter.x = (65.0)*cos(theta) + 0.20*(i - W/2)*sin(theta); //45.0
        Inter.y = 0.20*(j - H/2) ;
        Inter.z = (65.0)*sin(theta) + 0.20*(i - W/2)*(-cos(theta));
                                                  
        Inter.R0 = 20.0; //20.0; //20.0 25
        Inter.R1 = 24.0; //14.0;  //4.0 19.0
        Inter.phi = 0.0;
        Inter.deltaphi = M_PI;
        d = DistanceToTorus (&Inter);
                
        if (d < 0)
          {     
            V = 0;
            R = 255 ;
          } else {
            V = d*2;//4 16 32
          }
        if (V > 255) {
          V = 255 ;
        }
        
        if (V < 0) V = 0 ;
        buffer[4*(i*H+j)+0] = V;
        buffer[4*(i*H+j)+1] = V;
        buffer[4*(i*H+j)+2] = V+R;
      }
  t ++;
}