polynomial_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 polynomial
{
struct Parameters
{
    double R_0, 
           pp, 
           pi, 
           a,  
           elongation, 
           triangularity; 
    unsigned M, 
             N; 
    std::vector<double> c;  
    std::string description;
    Parameters() = default;
    Parameters( const dg::file::WrappedJsonValue& js) {
        pp  = js.get( "PP", 1).asDouble();
        pi  = js.get( "PI", 1).asDouble();
        M = js.get( "M", 1).asUInt();
        N = js.get( "N", 1).asUInt();
        c.resize(M*N);
        for (unsigned i=0;i<M*N;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();
        description = js.get( "description", "standardX").asString();
    }
    dg::file::JsonType dump( ) const
    {
        dg::file::JsonType js;
        js["M"] = M;
        js["N"] = N;
        js["PP"] = pp;
        js["PI"] = pi;
        for (unsigned i=0;i<M*N;i++) js["c"][i] = c[i];
        js["R_0"] = R_0;
        js["inverseaspectratio"] = a/R_0;
        js["elongation"] = elongation;
        js["triangularity"] = triangularity;
        js[ "equilibrium"] = "polynomial";
        js[ "description"] = description;
        return js;
    }
    bool isToroidal() const{
        if( pp == 0)
            return true;
        return false;
    }
    void display( std::ostream& os = std::cout ) const
    {
        dg::file::WrappedJsonValue js = dump();
        os << "Polynomial Geometrical parameters are: \n"
            <<js.toStyledString();
        os << std::flush;
 
    }
};
} //namespace polynomial
} //namespace geo
} //namespace dg