wall_3d.h

Go to the documentation of this file.

#ifndef WALL_3D_H
#define WALL_3D_H
#include "../ascot5.h"
#include "../offload.h"

#define WALL_EPSILON 1e-9

typedef struct {
    int n;               
    real xmin;           
    real xmax;           
    real xgrid;          
    real ymin;           
    real ymax;           
    real ygrid;          
    real zmin;           
    real zmax;           
    real zgrid;          
    int depth;           
    int ngrid;           
    real* wall_tris;     
    int* flag;           

    int* tree_array;
    int tree_array_size; /* Number of elements in tree_array */
} wall_3d_data;
 
int wall_3d_init(wall_3d_data* data, int nelements, real* x1x2x3, real* y1y2y3,
                 real* z1z2z3, int* flag);
void wall_3d_free(wall_3d_data* data);
void wall_3d_offload(wall_3d_data* data);
GPU_DECLARE_TARGET_SIMD_UNIFORM(w)
int wall_3d_hit_wall(real r1, real phi1, real z1, real r2, real phi2,
                     real z2, wall_3d_data* w, real* w_coll);
DECLARE_TARGET_END
GPU_DECLARE_TARGET_SIMD_UNIFORM(w)
int wall_3d_hit_wall_full(real r1, real phi1, real z1, real r2, real phi2,
                          real z2, wall_3d_data* w, real* w_coll);
DECLARE_TARGET_END
GPU_DECLARE_TARGET_SIMD
double wall_3d_tri_collision(real q1[3], real q2[3], real t1[3], real t2[3],
                             real t3[3]);
DECLARE_TARGET_END
 
DECLARE_TARGET
void wall_3d_init_tree(wall_3d_data* w, real* offload_array);
DECLARE_TARGET_END
DECLARE_TARGET
int wall_3d_tri_in_cube(real t1[3], real t2[3], real t3[3], real bb1[3],
                        real bb2[3]);
DECLARE_TARGET_END
DECLARE_TARGET
int wall_3d_quad_collision(real q1[3], real q2[3], real t1[3], real t2[3],
                           real t3[3], real t4[3]);
DECLARE_TARGET_END
 
#endif