mhd.c

Go to the documentation of this file.

 
#include <stdlib.h>
#include "ascot5.h"
#include "error.h"
#include "print.h"
#include "mhd.h"
#include "B_field.h"
#include "boozer.h"
#include "mhd/mhd_stat.h"
#include "mhd/mhd_nonstat.h"
 
void mhd_free(mhd_data* data) {
    switch(data->type) {
        case mhd_type_stat:
            mhd_stat_free(&data->stat);
            break;
        case mhd_type_nonstat:
            mhd_nonstat_free(&data->nonstat);
            break;
    }
}
 
void mhd_offload(mhd_data* data) {
    switch(data->type) {
        case mhd_type_stat:
            mhd_stat_offload(&data->stat);
            break;
        case mhd_type_nonstat:
            mhd_nonstat_offload(&data->nonstat);
            break;
    }
}
 
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) {
    a5err err = 0;
 
    switch(mhddata->type) {
 
        case mhd_type_stat:
            err = mhd_stat_eval(mhd_dmhd, r, phi, z, t, includemode,
                                boozerdata, &(mhddata->stat), Bdata);
            break;
 
        case mhd_type_nonstat:
            err = mhd_nonstat_eval(mhd_dmhd, r, phi, z, t, includemode,
                                   boozerdata, &(mhddata->nonstat), Bdata);
            break;
 
        default:
            /* Unregonized input. Produce error. */
            err = error_raise( ERR_UNKNOWN_INPUT, __LINE__, EF_MHD );
            break;
    }
 
    return err;
}
 
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) {
    a5err err = 0;
 
    switch(mhddata->type) {
 
        case mhd_type_stat:
            err =  mhd_stat_perturbations(pert_field, r, phi, z, t, pertonly,
                                          includemode, boozerdata,
                                          &(mhddata->stat), Bdata);
            break;
 
        case mhd_type_nonstat:
            err =  mhd_nonstat_perturbations(pert_field, r, phi, z, t, pertonly,
                                             includemode, boozerdata,
                                             &(mhddata->nonstat), Bdata);
            break;
 
        default:
            /* Unregonized input. Produce error. */
            err = error_raise( ERR_UNKNOWN_INPUT, __LINE__, EF_MHD );
            break;
    }
 
    return err;
}
 
int mhd_get_n_modes(mhd_data* mhddata) {
    int val = 0;
    switch(mhddata->type) {
        case mhd_type_stat:
            val = mhddata->stat.n_modes;
            break;
        case mhd_type_nonstat:
            val = mhddata->nonstat.n_modes;
            break;
    }
    return val;
}
 
const int* mhd_get_nmode(mhd_data* mhddata) {
    const int* val = NULL;
    switch(mhddata->type) {
        case mhd_type_stat:
            val = mhddata->stat.nmode;
            break;
        case mhd_type_nonstat:
            val = mhddata->nonstat.nmode;
            break;
    }
    return val;
}
 
const int* mhd_get_mmode(mhd_data* mhddata) {
    const int* val = NULL;
    switch(mhddata->type) {
        case mhd_type_stat:
            val = mhddata->stat.mmode;
            break;
        case mhd_type_nonstat:
            val = mhddata->nonstat.mmode;
            break;
    }
    return val;
}
 
const real* mhd_get_amplitude(mhd_data* mhddata) {
    const real* val = NULL;
    switch(mhddata->type) {
        case mhd_type_stat:
            val = mhddata->stat.amplitude_nm;
            break;
        case mhd_type_nonstat:
            val = mhddata->nonstat.amplitude_nm;
            break;
    }
    return val;
}
 
const real* mhd_get_frequency(mhd_data* mhddata) {
    const real* val = NULL;
    switch(mhddata->type) {
        case mhd_type_stat:
            val = mhddata->stat.omega_nm;
            break;
        case mhd_type_nonstat:
            val = mhddata->nonstat.omega_nm;
            break;
    }
    return val;
}
 
const real* mhd_get_phase(mhd_data* mhddata) {
    const real* val = NULL;
    switch(mhddata->type) {
        case mhd_type_stat:
            val = mhddata->stat.phase_nm;
            break;
        case mhd_type_nonstat:
            val = mhddata->nonstat.phase_nm;
            break;
    }
    return val;
}