neutral.h

Go to the documentation of this file.

 
#ifndef NEUTRAL_H
#define NEUTRAL_H
 
#include "ascot5.h"
#include "error.h"
#include "neutral/N0_1D.h"
#include "neutral/N0_3D.h"
 
typedef enum neutral_type {
    neutral_type_1D, 
    neutral_type_3D, 
} neutral_type;
 
typedef struct {
    neutral_type type;        
    N0_1D_offload_data N01D;  
    N0_3D_offload_data N03D;  
    int offload_array_length; 
} neutral_offload_data;
 
typedef struct {
    neutral_type type; 
    N0_1D_data N01D;   
    N0_3D_data N03D;   
} neutral_data;
 
int neutral_init_offload(neutral_offload_data* offload_data,
                         real** offload_array);
void neutral_free_offload(neutral_offload_data* offload_data,
                          real** offload_array);
 
int neutral_init(neutral_data* ndata, neutral_offload_data* offload_data,
                 real* offload_array);
DECLARE_TARGET_SIMD_UNIFORM(ndata)
a5err neutral_eval_n0(real* n0, real rho, real r, real phi, real z, real t,
                      neutral_data* ndata);
DECLARE_TARGET_SIMD_UNIFORM(ndata)
a5err neutral_eval_t0(real* t0, real rho, real r, real phi, real z, real t,
                      neutral_data* ndata);
DECLARE_TARGET_SIMD_UNIFORM(ndata)
int neutral_get_n_species(neutral_data* ndata);
#endif