B_2DS.h

Go to the documentation of this file.

 
#ifndef B_2DS_H
#define B_2DS_H
#include "../offload_acc_omp.h"
#include "../ascot5.h"
#include "../error.h"
#include "../spline/interp.h"
 
typedef struct {
    int n_r;                  
    int n_z;                  
    real r_min;               
    real r_max;               
    real z_min;               
    real z_max;               
    real psi0;                
    real psi1;                
    real axis_r;              
    real axis_z;              
    int offload_array_length; 
} B_2DS_offload_data;
 
typedef struct {
    real psi0;           
    real psi1;           
    real axis_r;         
    real axis_z;         
    interp2D_data psi;   
    interp2D_data B_r;   
    interp2D_data B_phi; 
    interp2D_data B_z;   
} B_2DS_data;
 
int B_2DS_init_offload(B_2DS_offload_data* offload_data, real** offload_array);
void B_2DS_free_offload(B_2DS_offload_data* offload_data, real** offload_array);
 
void B_2DS_init(B_2DS_data* Bdata, B_2DS_offload_data* offload_data,
                real* offload_array);
GPU_DECLARE_TARGET_SIMD_UNIFORM(Bdata)
a5err B_2DS_eval_psi(real* psi, real r, real phi, real z, B_2DS_data* Bdata);
DECLARE_TARGET_END
GPU_DECLARE_TARGET_SIMD_UNIFORM(Bdata)
a5err B_2DS_eval_psi_dpsi(real psi_dpsi[4], real r, real phi, real z,
                          B_2DS_data* Bdata);
DECLARE_TARGET_END
GPU_DECLARE_TARGET_SIMD_UNIFORM(Bdata)
a5err B_2DS_eval_rho_drho(real rho_drho[4], real r, real phi, real z,
                          B_2DS_data* Bdata);
DECLARE_TARGET_END
GPU_DECLARE_TARGET_SIMD_UNIFORM(Bdata)
a5err B_2DS_eval_B(real B[3], real r, real phi, real z, B_2DS_data* Bdata);
DECLARE_TARGET_END
GPU_DECLARE_TARGET_SIMD_UNIFORM(Bdata)
a5err B_2DS_eval_B_dB(real B_dB[12], real r, real phi, real z, B_2DS_data* Bdata);
DECLARE_TARGET_END
GPU_DECLARE_TARGET_SIMD_UNIFORM(Bdata)
a5err B_2DS_get_axis_rz(real rz[2], B_2DS_data* Bdata);
DECLARE_TARGET_END
#endif