diag_orb.h

Go to the documentation of this file.

 
#ifndef DIAG_ORB_H
#define DIAG_ORB_H
 
#include <stdio.h>
#include "../particle.h"
 
#define DIAG_ORB_POINCARE 0      
#define DIAG_ORB_INTERVAL 1      
#define DIAG_ORB_MAXPOINCARES 30 
#define DIAG_ORB_FOFIELDS     16 
#define DIAG_ORB_GCFIELDS     16 
#define DIAG_ORB_MLFIELDS     11 
#define DIAG_ORB_HYBRIDFIELDS 19 
#define DIAG_ORB_FO 1 
#define DIAG_ORB_GC 2 
#define DIAG_ORB_ML 3 
DECLARE_TARGET_SIMD_UNIFORM(ang0)
real diag_orb_check_plane_crossing(real fang, real iang, real ang0);
DECLARE_TARGET_SIMD_UNIFORM(r0)
real diag_orb_check_radial_crossing(real fr, real ir, real r0);
 
 
typedef struct{
    int record_mode;    
    int mode;           
    int Npnt;           
    int Nmrk;           
    int Nfld;           
    real writeInterval; 
    int ntoroidalplots; 
    int npoloidalplots; 
    int nradialplots;   
    real toroidalangles[DIAG_ORB_MAXPOINCARES]; 
    real poloidalangles[DIAG_ORB_MAXPOINCARES]; 
    real radialdistances[DIAG_ORB_MAXPOINCARES];   
}diag_orb_offload_data;
 
typedef struct{
 
    real* id;     
    real* mileage;
    real* r;      
    real* phi;    
    real* z;      
    real* p_r;    
    real* p_phi;  
    real* p_z;    
    real* ppar;   
    real* mu;     
    real* zeta;   
    real* weight; 
    real* charge; 
    real* rho;    
    real* theta;  
    real* B_r;    
    real* B_phi;  
    real* B_z;    
    real* simmode;
    real* pncrid; 
    real* pncrdi; 
    integer* mrk_pnt;     
    real* mrk_recorded;   
    int mode;           
    int Npnt;           
    int Nmrk;           
    real writeInterval; 
    int ntoroidalplots; 
    int npoloidalplots; 
    int nradialplots;   
    real toroidalangles[DIAG_ORB_MAXPOINCARES]; 
    real poloidalangles[DIAG_ORB_MAXPOINCARES]; 
    real radialdistances[DIAG_ORB_MAXPOINCARES];   
}diag_orb_data;
 
void diag_orb_init(diag_orb_data* data, diag_orb_offload_data* offload_data,
                   real* offload_array);
 
void diag_orb_free(diag_orb_data* data);
 
void diag_orb_update_fo(diag_orb_data* data,
                        particle_simd_fo* p_f, particle_simd_fo* p_i);
 
void diag_orb_update_gc(diag_orb_data* data,
                        particle_simd_gc* p_f, particle_simd_gc* p_i);
 
void diag_orb_update_ml(diag_orb_data* data,
                        particle_simd_ml* p_f, particle_simd_ml* p_i);
 
#endif