solovev_parameters.h

Go to the documentation of this file.

#pragma once
#include <string>
#include <vector>
#include <dg/file/json_utilities.h>
 
namespace dg
{
namespace geo
{
namespace solovev
{
struct Parameters
{
    double A, 
           R_0, 
           pp, 
           pi, 
           a,  
           elongation, 
           triangularity; 
    std::vector<double> c;  
    std::string description;
    Parameters() = default;
    Parameters( const dg::file::WrappedJsonValue& js) {
        A   = js.get("A", 0).asDouble();
        pp  = js.get("PP", 1).asDouble();
        pi  = js.get("PI", 1).asDouble();
        c.resize(12);
        for (unsigned i=0;i<12;i++)
            c[i] = js["c"].get(i,0.).asDouble();
 
        R_0  =          js.get( "R_0", 0.).asDouble();
        a        = R_0* js.get( "inverseaspectratio", 0.).asDouble();
        elongation=     js.get( "elongation", 1.).asDouble();
        triangularity=  js.get( "triangularity", 0.).asDouble();
        try{
            description = js.get( "description", "standardX").asString();
        } catch ( std::exception& err)
        {
            if( isToroidal())
                description = "none";
            else if( !hasXpoint())
                description = "standardO";
            else
                description = "standardX";
        }
    }
    dg::file::JsonType dump( ) const
    {
        dg::file::JsonType js;
        js["A"] = A;
        js["PP"] = pp;
        js["PI"] = pi;
        for (unsigned i=0;i<12;i++) js["c"][i] = c[i];
        js["R_0"] = R_0;
        js["inverseaspectratio"] = a/R_0;
        js["elongation"] = elongation;
        js["triangularity"] = triangularity;
        js[ "equilibrium"] = "solovev";
        js[ "description"] = description;
        return js;
    }
    bool hasXpoint( ) const{
        bool Xpoint = false;
        for( int i=7; i<12; i++)
            if( fabs(c[i]) >= 1e-10)
                Xpoint = true;
        return Xpoint;
    }
    bool isToroidal() const{
        if( pp == 0)
            return true;
        return false;
    }
    void display( std::ostream& os = std::cout ) const
    {
        dg::file::WrappedJsonValue js = dump();
        os << "Solovev Geometrical parameters are: \n"
            <<js.toStyledString();
        os << std::flush;
    }
};
} //namespace solovev
} //namespace geo
} //namespace dg