/TrouNoir/trou_noir_MyFloat.c - Diff - Bench4GPU - Forge du Centre Blaise Pascal

Révision 212 TrouNoir/trou_noir_MyFloat.c

     	     MYFLOAT h,MYFLOAT r,MYFLOAT m,MYFLOAT bss)
+    {
       MYFLOAT flx;
       int fi;
       fi=(int)(rf*nbr/bss);
       flx=exp(q*log(r/m))*pow(rf,4)*b*db*h;
       return(flx);
+    }
-...
       MYFLOAT nu_rec,nu_em,qu,temp_em,flux_int;
       int fi,posfreq;
     #define  planck 6.62e-34
     #define planck 6.62e-34
     #define k 1.38e-23
     #define c2 9.e16
     #define temp 3.e7
     #define m_point 1.
     #define lplanck (log(6.62)-34.*log(10.))
     #define lk (log(1.38)-23.*log(10.))
     #define lc2 (log(9.)+16.*log(10.))
       MYFLOAT v=1.-3./r;
       qu=1./sqrt((1.-3./r)*r)*(sqrt(r)-sqrt(6.)+sqrt(3.)/2.*log((sqrt(r)+sqrt(3.))/(sqrt(r)-sqrt(3.))* 0.17157287525380988 ));
       temp_em=temp*sqrt(m)*exp(0.25*log(m_point)-0.75*log(r)-0.125*log(v)+0.25*log(fabs(qu)));
       flux_int=0;
       flx=0;
       flux_int=0.;
       flx=0.;
       for (fi=0;fi<nbr;fi++)
+        {
           nu_em=bss*fi/(MYFLOAT)nbr;
           nu_rec=nu_em*rf;
           nu_em=bss*(MYFLOAT)fi/(MYFLOAT)nbr;
           nu_rec=nu_em*rf;
           posfreq=(int)(nu_rec*(MYFLOAT)nbr/bss);
           if ((posfreq>0)&&(posfreq<nbr))
+    	{
     	  flux_int=2.*planck/c2*pow(nu_em,3)/(exp(planck*nu_em/(k*temp_em))-1.);
     	  flux_int*=pow(rf,3)*b*db*h;
     	  flx+=flux_int;
+    	}
+      	{
     	  flux_int=2.*planck/c2*pow(nu_em,3)/(exp(planck*nu_em/(k*temp_em))-1.)*exp(3.*log(rf))*b*db*h;
       	  flx+=flux_int;
+      	}
+        }
       return((MYFLOAT)flx);
-...
       if (raie==0)
+        {
           bss=1.e19;
           flx=spectre_cn(rf,b,db,h,r,m,bss);
+        }
       else
+        {
           bss=2.;
           flx=spectre(rf,q,b,db,h,r,m,bss);
+        }
-...
+        {
           fp[xi][yi]=flx;
+        }
+    }
     void sauvegarde_pgm(char nom[24],unsigned int **image,int dim)
-...
     int main(int argc,char *argv[])
+    {
       MYFLOAT m,rs,ri,re,tho,ro;
       MYFLOAT m,rs,ri,re,tho;
       int q;
       MYFLOAT bss,stp;
-...
+    	{
     	  printf("\nSimulation d'un disque d'accretion autour d'un trou noir\n");
     	  printf("\nParametres a definir :\n\n");
     	  printf("   1) Dimension de l'Image\n");
     	  printf("   2) Masse relative du trou noir\n");
     	  printf("   3) Dimension du disque exterieur\n");
     	  printf("   4) Distance de l'observateur\n");
     	  printf("   5) Inclinaison par rapport au disque (en degres)\n");
     	  printf("   6) Observation a distance FINIE ou INFINIE\n");
     	  printf("   7) Rayonnement de disque MONOCHROMATIQUE ou CORPS_NOIR\n");
     	  printf("   8) Normalisation des flux INTERNE ou EXTERNE\n");
     	  printf("   9) Normalisation de z INTERNE ou EXTERNE\n");
     	  printf("  10) Impression des images NEGATIVE ou POSITIVE\n");
     	  printf("  11) Nom de l'image des Flux\n");
     	  printf("  12) Nom de l'image des decalages spectraux\n");
     	  printf("  13) Valeur de normalisation des flux\n");
     	  printf("  14) Valeur de normalisation des decalages spectraux\n");
     	  printf("  1) Dimension de l'Image\n");
     	  printf("  2) Masse relative du trou noir\n");
     	  printf("  3) Dimension du disque exterieur\n");
     	  printf("  4) Inclinaison par rapport au disque (en degres)\n");
     	  printf("  5) Rayonnement de disque MONOCHROMATIQUE ou CORPS_NOIR\n");
     	  printf("  6) Impression des images NEGATIVE ou POSITIVE\n");
     	  printf("  7) Nom de l'image des Flux\n");
     	  printf("  8) Nom de l'image des decalages spectraux\n");
     	  printf("\nSi aucun parametre defini, parametres par defaut :\n\n");
     	  printf("   1) Dimension de l'image : 1024 pixels de cote\n");
     	  printf("   2) Masse relative du trou noir : 1\n");
     	  printf("   3) Dimension du disque exterieur : 12 \n");
     	  printf("   4) Distance de l'observateur : 100 \n");
     	  printf("   5) Inclinaison par rapport au disque (en degres) : 10\n");
     	  printf("   6) Observation a distance FINIE\n");
     	  printf("   7) Rayonnement de disque CORPS_NOIR\n");
     	  printf("   8) Normalisation des flux INTERNE\n");
     	  printf("   9) Normalisation des z INTERNE\n");
     	  printf("  10) Impression des images NEGATIVE ou POSITIVE\n");
            	  printf("  11) Nom de l'image des flux : flux.pgm\n");
     	  printf("  12) Nom de l'image des z : z.pgm\n");
     	  printf("  13) <non definie>\n");
     	  printf("  14) <non definie>\n");
     	  printf("  1) Dimension de l'image : 1024 pixels de cote\n");
     	  printf("  2) Masse relative du trou noir : 1\n");
     	  printf("  3) Dimension du disque exterieur : 12 \n");
     	  printf("  4) Inclinaison par rapport au disque (en degres) : 10\n");
     	  printf("  5) Rayonnement de disque CORPS_NOIR\n");
     	  printf("  6) Impression des images NEGATIVE ou POSITIVE\n");
            	  printf("  7) Nom de l'image des flux : flux.pgm\n");
     	  printf("  8) Nom de l'image des z : z.pgm\n");
+    	}
+        }
       if ((argc==13)||(argc==15))
       if (argc==9)
+        {
           printf("# Utilisation les valeurs definies par l'utilisateur\n");
           dim=atoi(argv[1]);
           m=atof(argv[2]);
           re=atof(argv[3]);
           ro=atof(argv[4]);
           tho=PI/180.*(90-atof(argv[5]));
           tho=PI/180.*(90-atof(argv[4]));
           rs=2.*m;
           ri=3.*rs;
           if (strcmp(argv[6],"FINIE")==0)
           if (strcmp(argv[5],"CORPS_NOIR")==0)
+    	{
     	  dist=0;
+    	}
           else
+    	{
     	  dist=1;
+    	}
           if (strcmp(argv[7],"CORPS_NOIR")==0)
+    	{
     	  raie=0;
+    	}
           else
-...
     	  raie=1;
+    	}
           if (strcmp(argv[8],"EXTERNE")==0)
+    	{
     	  fen=atof(argv[14]);
+    	}
           if (strcmp(argv[9],"EXTERNE")==0)
+    	{
     	  zen=atof(argv[15]);
+    	}
+        }
       else
+        {
-...
           rs=2.*m;
           ri=3.*rs;
           re=12.;
           ro=100.;
           tho=PI/180.*80;
           // Distance finie
           dist=0;
           // Corps noir
           raie=0;
+        }
-...
+        }
       else
+        {
           bss=1e19;
           bss=1.e19;
           q=-0.75;
+        }
-...
           printf("# Rayon singularite : %f\n",rs);
           printf("# Rayon interne : %f\n",ri);
           printf("# Rayon externe : %f\n",re);
           printf("# Distance de l'observateur : %f\n",ro);
           printf("# Inclinaison a la normale en radian : %f\n",tho);
       zp=(MYFLOAT**)calloc(dim,sizeof(MYFLOAT*));
-...
       stp=dim/(2.*nmx);
       bmx=1.25*re;
       b=0.;
       thx=asin(bmx/ro);
       pc=0;
       struct timeval tv1,tv2;
-...
           h=4.*PI/(MYFLOAT)TRACKPOINTS;
           d=stp*n;
           if (dist==1)
+    	{
     	  db=bmx/(MYFLOAT)nmx;
     	  b=db*(MYFLOAT)n;
     	  up=0.;
     	  vp=1.;
+    	}
           else
+    	{
     	  thi=thx/(MYFLOAT)nmx*(MYFLOAT)n;
     	  db=ro*sin(thi)-b;
     	  b=ro*sin(thi);
     	  up=sin(thi);
     	  vp=cos(thi);
+    	}
           db=bmx/(MYFLOAT)nmx;
           b=db*(MYFLOAT)n;
           up=0.;
           vp=1.;
           pp=0.;
           nh=1;
-...
       fmx=fp[0][0];
       zmx=zp[0][0];
       int zimx=0,zjmx=0,fimx=0,fjmx=0;
       for (i=0;i<dim;i++) for (j=0;j<dim;j++)
+        {
           if (fmx<fp[i][j])
+    	{
     	  fimx=i;
     	  fjmx=j;
     	  fmx=fp[i][j];
+    	}
           if (zmx<zp[i][j])
+    	{
     	  zimx=i;
     	  zjmx=j;
     	  zmx=zp[i][j];
+    	}
+        }
       printf("\nElapsed Time : %lf",(double)elapsed);
       printf("\nFlux max : %f",fmx);
       printf("\nZ max : %f\n\n",zmx);
       printf("\nZ max @(%i,%i) : %f",zimx,zjmx,zmx);
       printf("\nFlux max @(%i,%i) : %f\n\n",fimx,fjmx,fmx);
       if (strcmp(argv[8],"EXTERNE")==0)
+        {
           fmx=fen;
+        }
       if (strcmp(argv[9],"EXTERNE")==0)
+        {
           zmx=zen;
+        }
       for (i=0;i<dim;i++) for (j=0;j<dim;j++)
+        {
           zcl=(int)(255/zmx*zp[i][dim-1-j]);
           fcl=(int)(255/fmx*fp[i][dim-1-j]);
           if (strcmp(argv[8],"NEGATIVE")==0)
           if (strcmp(argv[6],"NEGATIVE")==0)
+    	{
     	  if (zcl>255)
+    	    {
-...
+        }
       if ((argc==14)||(argc==16))
       if (argc==9)
+       {
          sauvegarde_pgm(argv[11],ifp,dim);
          sauvegarde_pgm(argv[12],izp,dim);
          sauvegarde_pgm(argv[7],ifp,dim);
          sauvegarde_pgm(argv[8],izp,dim);
+       }
       else
+        {

Formats disponibles : Unified diff

Centre Blaise Pascal » Bench4GPU

Révision 212 TrouNoir/trou_noir_MyFloat.c