mhd.h

Go to the documentation of this file.

 
#ifndef MHD_H
#define MHD_H
 
#include "ascot5.h"
#include "error.h"
#include "B_field.h"
#include "boozer.h"
#include "mhd/mhd_stat.h"
#include "mhd/mhd_nonstat.h"
 
#define MHD_INCLUDE_ALL -1
 
typedef enum mhd_type {
    mhd_type_stat,   
    mhd_type_nonstat 
} mhd_type;
 
typedef struct {
    mhd_type type;                   
    mhd_stat_offload_data stat;      
    mhd_nonstat_offload_data nonstat;
    int offload_array_length;        
} mhd_offload_data;
 
typedef struct {
    mhd_type type;            
    mhd_stat_data stat;       
    mhd_nonstat_data nonstat; 
} mhd_data;
 
int mhd_init_offload(mhd_offload_data* offload_data,
                     real** offload_array);
void mhd_free_offload(mhd_offload_data* offload_data,
                      real** offload_array);
 
int mhd_init(mhd_data* data, mhd_offload_data* offload_data,
             real* offload_array);
DECLARE_TARGET_SIMD_UNIFORM(boozerdata, mhddata, Bdata, includemode)
a5err mhd_eval(real mhd_dmhd[10], real r, real phi, real z, real t,
               int includemode, boozer_data* boozerdata, mhd_data* mhddata,
               B_field_data* Bdata);
DECLARE_TARGET_SIMD_UNIFORM(boozerdata, mhddata, Bdata, pertonly,\
                                 includemode)
a5err mhd_perturbations(real pert_field[7], real r, real phi, real z,
                        real t, int pertonly, int includemode,
                        boozer_data* boozerdata, mhd_data* mhddata,
                        B_field_data* Bdata);
DECLARE_TARGET_SIMD_UNIFORM(mhddata)
int mhd_get_n_modes(mhd_data* mhddata);
DECLARE_TARGET_SIMD_UNIFORM(mhddata)
const int* mhd_get_nmode(mhd_data* mhddata);
DECLARE_TARGET_SIMD_UNIFORM(mhddata)
const int* mhd_get_mmode(mhd_data* mhddata);
DECLARE_TARGET_SIMD_UNIFORM(mhddata)
const real* mhd_get_amplitude(mhd_data* mhddata);
DECLARE_TARGET_SIMD_UNIFORM(mhddata)
const real* mhd_get_frequency(mhd_data* mhddata);
DECLARE_TARGET_SIMD_UNIFORM(mhddata)
const real*  mhd_get_phase(mhd_data* mhddata);
#endif