analytic_I.hpp

Go to the documentation of this file.

#ifndef ANALYTIC_I_HPP
#define ANALYTIC_I_HPP

#include "globals.hpp"
#include "magnetic_field.hpp"

// Provides setup and get I for an analytic example of I
// get_analytic_I() will be replaced with interp eventually

namespace{

// Simple expression for I and dI/dpsi
double get_analytic_I(double psi, double magnitude){
    return magnitude/(psi+0.1);
}

double get_analytic_dIdpsi(double psi, double magnitude){
    return -magnitude/((psi+0.1)*(psi+0.1));
}

// Create an I grid based off the analytic I
View<double*, HostType> analytic_I_grid(double magnitude, const View<double*, HostType>& psi){
    View<double*, HostType> I_grid("I_grid",psi.size());
    for (int i=0; i<psi.size(); i++){
        I_grid(i) = get_analytic_I(psi(i), magnitude);
    }
    return I_grid;
}

}
#endif