E_field.c

Go to the documentation of this file.

#include <stdio.h>
#include "ascot5.h"
#include "error.h"
#include "print.h"
#include "E_field.h"
#include "B_field.h"
#include "Efield/E_TC.h"
#include "Efield/E_1DS.h"
#include "Efield/E_2DS.h"

void E_field_free(E_field_data* data) {
    switch(data->type) {
        case E_field_type_2DS:
            E_2DS_free(&data->E2DS);
            break;
        case E_field_type_1DS:
            E_1DS_free(&data->E1DS);
            break;
        case E_field_type_TC:
            E_TC_free(&data->ETC);
            break;
    }
}

void E_field_offload(E_field_data* data) {
    switch(data->type) {
        case E_field_type_2DS:
            E_2DS_offload(&data->E2DS);
            break;
        case E_field_type_1DS:
            E_1DS_offload(&data->E1DS);
            break;
        case E_field_type_TC:
            E_TC_offload(&data->ETC);
            break;
    }
}

a5err E_field_eval_E(real E[3], real r, real phi, real z, real t,
                     E_field_data* Edata, B_field_data* Bdata) {
    a5err err = 0;
 
    switch(Edata->type) {
 
        case E_field_type_2DS:
            err = E_2DS_eval_E(E, r, z, &(Edata->E2DS));
            break;
 
        case E_field_type_1DS:
            err = E_1DS_eval_E(E, r, phi, z, &(Edata->E1DS), Bdata);
            break;
 
        case E_field_type_TC:
            err = E_TC_eval_E(E, r, phi, z, &(Edata->ETC), Bdata);
            break;
 
        default:
            /* Unregonized input. Produce error. */
            err = error_raise( ERR_UNKNOWN_INPUT, __LINE__, EF_E_FIELD );
            break;
    }
 
    return err;
}