/ - Diff - Bench4GPU - Forge du Centre Blaise Pascal

Révision 226

     # - Luc Blanchet for his disponibility about my questions in General Relativity
     # - Pierre Lena for his passion about science and vulgarisation
     # If crash on OpenCL Intel implementation, disable VECTORIZER with
     # export CL_CONFIG_USE_VECTORIZER=0
     import pyopencl as cl
     import numpy
     import time,string
-...
         PhysicsList={'Einstein':0,'Newton':1}
         return(PhysicsList)
+    #
     # Blank space below to simplify debugging on OpenCL code
+    #
     BlobOpenCL= """
     #define PI (float)3.14159265359
-...
        // Perform trajectory for each pixel, exit on hit
       private float m,rs,ri,re,tho;
       private int q,raie;
       float m,rs,ri,re,tho;
       int q,raie;
       m=Mass;
       rs=2.*m;
-...
       q=-2;
       raie=Line;
       private float d,bmx,db,b,h;
       private float rp0,rpp,rps;
       private float phi,thi,phd,php,nr,r;
       private int nh;
       private float zp,fp;
       float bmx,db,b,h;
       float rp0,rpp,rps;
       float phi,phd;
       int nh;
       float zp=0.,fp=0.;
       // Autosize for image
       bmx=1.25*re;
       b=0.;
       h=4.e0f*PI/(float)TRACKPOINTS;
-...
       nh=0;
       rungekutta(&ps,&us,&vs,pp,up,vp,h,m,b);
       rps=fabs(b/us);
       rp0=rps;
-...
          pp=ps;
          up=us;
          vp=vs;
          rungekutta(&ps,&us,&vs,pp,up,vp,h,m,b);
          rpp=rps;
          //    rungekutta(&ps,&us,&vs,pp,up,vp,h,m,b);
          rps=fabs(b/us);
          ExitOnImpact = ((fmod(pp,PI)<fmod(phd,PI))&&(fmod(ps,PI)>fmod(phd,PI)))&&(rps>ri)&&(rps<re)?1:0;
          ExitOnImpact = ((fmod(pp,PI)<fmod(phd,PI))&&(fmod(ps,PI)>fmod(phd,PI)))&&(rps>=ri)&&(rps<=re)?1:0;
       } while ((rps>=rs)&&(rps<=rp0)&&(ExitOnImpact==0));
-...
        // Perform trajectory for each pixel
       float m,rs,ri,re,tho;
       float m,ri,re,tho;
       int q,raie;
       m=Mass;
       rs=2.*m;
       ri=InternalRadius;
       re=ExternalRadius;
       tho=Angle;
       q=-2;
       raie=Line;
       float d,bmx,db,b,h;
       float phi,thi,phd,php,nr,r;
       int nh;
       float zp=0,fp=0;
       float bmx,db,b,h;
       float phi,phd,php,nr,r;
       float zp=0.,fp=0.;
       // Autosize for image, 25% greater than external radius
       bmx=1.25*re;
-...
        // Perform trajectory for each pixel
       float m,rs,ri,re,tho;
       float m,ri,re,tho;
       int q,raie;
       m=Mass;
       rs=2.*m;
       ri=InternalRadius;
       re=ExternalRadius;
       tho=Angle;
       raie=Line;
       float bmx,db,b,h;
       float phi,thi,phd;
       int nh;
       float phi,phd;
       float zp=0,fp=0;
       // Autosize for image

     #!/bin/bash
     DEVICE=1
     DEVICE=2
     GPU=OpenCL
     SEQ="6 1 14"
     HOST=$(hostname)
-...
         SIZE=$((2**$POWER))
         seq 1 1 10 | xargs -I TOTO /usr/bin/time python ./TrouNoir.py -d $DEVICE -g $GPU -n -s $SIZE -t $METHOD >>$LOGFILE 2>&1
     done
     LINE=BB
     METHOD=EachCircle
     echo -e "Experience : $LINE $METHOD" >>$LOGFILE
     echo -e "Experience : $LINE $METHOD $DEVICE $GPU"
     seq $SEQ | while read POWER ; do
         SIZE=$((2**$POWER))
         seq 1 1 10 | xargs -I TOTO /usr/bin/time python ./TrouNoir.py -d $DEVICE -g $GPU -n -s $SIZE -t $METHOD >>$LOGFILE 2>&1
     done
     LINE=MONO
     METHOD=EachCircle
     echo -e "Experience : $LINE $METHOD" >>$LOGFILE
     echo -e "Experience : $LINE $METHOD $DEVICE $GPU"
     seq $SEQ | while read POWER ; do
         SIZE=$((2**$POWER))
         seq 1 1 10 | xargs -I TOTO /usr/bin/time python ./TrouNoir.py -d $DEVICE -g $GPU -n -s $SIZE -t $METHOD >>$LOGFILE 2>&1
     done

Formats disponibles : Unified diff

Centre Blaise Pascal » Bench4GPU

Révision 226