mcx_core.h

Go to the documentation of this file.

/***************************************************************************/
/***************************************************************************/
#ifndef _MCEXTREME_GPU_LAUNCH_H
#define _MCEXTREME_GPU_LAUNCH_H

#include "mcx_utils.h"

#ifdef  __cplusplus
extern "C" {
#endif


#define ABS(a)  ((a)<0?-(a):(a))                 
#define MAX_ACCUM           1000.f

#define ROULETTE_SIZE       10.f                  
#ifdef  MCX_DEBUG
#define GPUDEBUG(x)        printf x             
#else
#define GPUDEBUG(x)                             
#endif

typedef float4 MCXpos; 
typedef struct __align__(16) SplitVoxel {
    unsigned char issplit; 
    unsigned char lower;   
    unsigned char upper;   
    unsigned char isupper; 
} SVox;

typedef struct __align__(16) StokesVector {
    float i; 
    float q; 
    float u; 
    float v; 
} Stokes;

typedef struct __align__(16) MCXSplit {
    SVox   sv; 
    float3 rp; 
    float3 nv; 
} MCXsp;

typedef struct  __align__(16) MCXDir {
    float x; 
    float y; 
    float z; 
    float nscat; 
} MCXdir;

typedef struct  __align__(16) MCXTimer {
    float pscat; 
    float t;     
    float pathlen; 
    float ndone; 
} MCXtime;

typedef union  __align__(16) GPosition {
    MCXpos d;
    float4 v;
    float  f[4];
} Gpos;

typedef union  __align__(16) GDirection {
    MCXdir d;
    float4 v;
    float  f[4];
} Gdir;

typedef union  __align__(16) GLength {
    MCXtime d;
    float4 v;
    float  f[4];
} Glen;

typedef union  __align__(16) GProperty {
    Medium d; /*defined in mcx_utils.h*/
    float4 v;
    float  f[4];
} Gprop;

typedef unsigned char uchar;

typedef struct  __align__(16) KernelParams {
    float3 vsize;                      
    float  minstep;                    
    float  twin0;                      
    float  twin1;                      
    float  tmax;                       
    float  oneoverc0;                  
    unsigned int save2pt;              
    unsigned int doreflect;            
    unsigned int dorefint;             
    unsigned int savedet;              
    float  Rtstep;                     
    float4 ps;                         
    float4 c0;                         
    float4 s0;                         
    float3 maxidx;                     
    uint4  dimlen;                     
    uint3  cp0;                        
    uint3  cp1;                        
    uint2  cachebox;                   
    float  minenergy;                  
    float  skipradius2;                
    float  minaccumtime;               
    int    srctype;                    
    float4 srcparam1;                  
    float4 srcparam2;                  
    int voidtime;                      
    unsigned int maxdetphoton;         
    unsigned int maxmedia;             
    unsigned int detnum;               
    unsigned int maxpolmedia;          
    unsigned int maxgate;              
    unsigned int idx1dorig;            
    unsigned int mediaidorig;          
    unsigned int isatomic;             
    unsigned int maxvoidstep;          
    unsigned int issaveseed;           
    unsigned int issaveref;            
    unsigned int isspecular;           
    unsigned int seedoffset;           
    int seed;                          
    unsigned int outputtype;           
    int threadphoton;                  
    int oddphotons;                    
    int faststep;                      
    unsigned int debuglevel;           
    unsigned int savedetflag;          
    unsigned int reclen;               
    unsigned int partialdata;          
    unsigned int w0offset;             
    unsigned int mediaformat;          
    unsigned int maxjumpdebug;         
    unsigned int gscatter;             
    unsigned int is2d;                 
    int replaydet;                     
    unsigned int srcnum;               
    unsigned int nphase;               
    float omega;                       
    unsigned char bc[12];              
} MCXParam;

void mcx_run_simulation(Config* cfg, GPUInfo* gpu);
int  mcx_list_gpu(Config* cfg, GPUInfo** info);

#ifdef  __cplusplus
}
#endif

#endif