mhd_stat.h

Go to the documentation of this file.

 
#ifndef MHD_STAT_H
#define MHD_STAT_H
 
#include "../ascot5.h"
#include "../error.h"
#include "../boozer.h"
#include "../spline/interp.h"
#include "../B_field.h"
 
typedef struct {
 
    int n_modes;                          
    int nrho;                             
    real rho_min;                         
    real rho_max;                         
    int nmode[MHD_MODES_MAX_NUM];         
    int mmode[MHD_MODES_MAX_NUM];         
    real amplitude_nm[MHD_MODES_MAX_NUM]; 
    real omega_nm[MHD_MODES_MAX_NUM];     
    real phase_nm[MHD_MODES_MAX_NUM];     
    int offload_array_length; 
} mhd_stat_offload_data;
 
typedef struct {
    int n_modes;                          
    real rho_min;                         
    real rho_max;                         
    int nmode[MHD_MODES_MAX_NUM];         
    int mmode[MHD_MODES_MAX_NUM];         
    real amplitude_nm[MHD_MODES_MAX_NUM]; 
    real omega_nm[MHD_MODES_MAX_NUM];     
    real phase_nm[MHD_MODES_MAX_NUM];     
    interp1D_data alpha_nm[MHD_MODES_MAX_NUM];
 
    interp1D_data phi_nm[MHD_MODES_MAX_NUM];
} mhd_stat_data;
 
int mhd_stat_init_offload(mhd_stat_offload_data* offload_data,
                          real** offload_array);
 
void mhd_stat_free_offload(mhd_stat_offload_data* offload_data,
                           real** offload_array);
 
void mhd_stat_init(mhd_stat_data* mhdata, mhd_stat_offload_data* offload_data,
                   real* offload_array);
DECLARE_TARGET_SIMD_UNIFORM(boozerdata, mhddata, includemode)
a5err mhd_stat_eval(real mhd_dmhd[10], real r, real phi, real z, real t,
                    int includemode, boozer_data* boozerdata,
                    mhd_stat_data* mhddata, B_field_data* Bdata);
DECLARE_TARGET_SIMD_UNIFORM(boozerdata, mhddata, Bdata, pertonly,\
                                 includemode)
a5err mhd_stat_perturbations(real pert_field[7], real r, real phi, real z,
                             real t, int pertonly, int includemode,
                             boozer_data* boozerdata,
                             mhd_stat_data* mhddata, B_field_data* Bdata);
 
 
#endif