mmc_raytrace.h

Go to the documentation of this file.

/***************************************************************************/
/***************************************************************************/
#ifndef _MMC_RAY_TRACING_H
#define _MMC_RAY_TRACING_H

#include "mmc_mesh.h"
#include "mmc_utils.h"

#define MAX_TRIAL          3          
#define FIX_PHOTON         1e-3f      
/***************************************************************************/
typedef struct MMC_ray {
    float3 p0;                    
    float3 vec;                   
    float3 pout;                  
    float4 bary0;                 
    int eid;                      
    int faceid;                   
    int isend;                    
    int nexteid;                  
    float weight;                 
    float photontimer;            
    float slen0;                  
    float slen;                   
    float Lmove;                  
    double Eabsorb;               
    unsigned int photonid;        
    float* partialpath;           
    void*  photonseed;            
    float focus;                  
    unsigned int posidx;          
    unsigned int oldidx;
    double oldweight;
    float* roisize;               
    int roitype;              
    int inroi;            
    int roiidx;           
    int refeid;           
} ray;

/***************************************************************************/
typedef struct MMC_visitor {
    float raytet;                 
    float raytet0;                
    float rtstep;                 
    int   detcount;               
    int   bufpos;                 
    int   reclen;                 
    float* partialpath;           
    void*  photonseed;            
    double* launchweight;         
    double* absorbweight;         
    double* kahanc0;              
    double* kahanc1;              
} visitor;

#ifdef  __cplusplus
extern "C" {
#endif
void interppos(float3* w, float3* p1, float3* p2, float3* p3, float3* pout);
void getinterp(float w1, float w2, float w3, float3* p1, float3* p2, float3* p3, float3* pout);
void fixphoton(float3* p, float3* nodes, int* ee);
void onephoton(size_t id, raytracer* tracer, tetmesh* mesh, mcconfig* cfg, RandType* ran, RandType* ran0, visitor* visit);
void launchphoton(mcconfig* cfg, ray* r, tetmesh* mesh, RandType* ran, RandType* ran0);
float reflectray(mcconfig* cfg, float3* c0, raytracer* tracer, int* oldeid, int* eid, int faceid, RandType* ran, int inroi);
float reflectrayroi(mcconfig* cfg, float3* c0, float3* ph, raytracer* tracer, int* eid, int* inroi, RandType* ran, int roitype, int roiidx, int refeid);
void save_scatter_events(ray* r, tetmesh* mesh, mcconfig* cfg, visitor* visit);
void albedoweight(ray* r, tetmesh* mesh, mcconfig* cfg, visitor* visit);
void visitor_init(mcconfig* cfg, visitor* visit);
void visitor_clear(visitor* visit);
void updateroi(int immctype, ray* r, tetmesh* mesh);
void traceroi(ray* r, raytracer* tracer, int roitype, int doinit);
void compute_distances_to_edge(ray* r, raytracer* tracer, int* ee, int edgeid, float d2d[2], float3 p2d[2], int* hitstatus);
void compute_distances_to_node(ray* r, raytracer* tracer, int* ee, int index, float nr, float3** center, int* hitstatus);
float ray_cylinder_intersect(ray* r, int index, float d2d[2], float3 p2d[2], int hitstatus);
float ray_sphere_intersect(ray* r, int index, float3* center, float nr, int hitstatus);
float ray_face_intersect(ray* r, raytracer* tracer, int* ee, int index, int baseid, int eid, int* hitstatus);

#ifdef MCX_CONTAINER
#ifdef __cplusplus
extern "C"
#endif
int mcx_throw_exception(const int id, const char* msg, const char* filename, const int linenum);
#endif

#ifdef  __cplusplus
}
#endif

#endif