AbsSin.H

Go to the documentation of this file.

#ifndef MODEL_INTERFACE_GB_ABSSIN_H
#define MODEL_INTERFACE_GB_ABSSIN_H
 
#include <iostream>
#include <fstream>
 
#include "AMReX.H"
#include "GB.H"
#include "Set/Set.H"
#include "Util/Util.H"
 
namespace Model
{
namespace Interface
{
namespace GB
{
class AbsSin : public GB
{
public:
    static constexpr const char* name = "abssin";
 
AbsSin() {};
AbsSin(IO::ParmParse &pp) {pp_queryclass(*this);};
AbsSin(IO::ParmParse &pp, std::string name) {pp_queryclass(name,*this);};
 
AbsSin(amrex::Real _theta0, amrex::Real _sigma0, amrex::Real _sigma1)
{       
    Define(_theta0,_sigma0,_sigma1);
};
void Define(amrex::Real _theta0, amrex::Real _sigma0, amrex::Real _sigma1)
{
    theta0 = _theta0;
    sigma0 = _sigma0;
    sigma1 = _sigma1;
};
void Randomize()
{
    theta0 = Util::Random()*Set::Constant::Pi;
    sigma0 = Util::Random();
    sigma1 = Util::Random();
};
Set::Scalar W(const Set::Scalar theta) const
{
    //sigma(theta)=sigma0+sigma1*fabs(sin(n*(theta-theta0)))
    //n=2:  
    return sigma0+sigma1*fabs(sin(2*(theta-theta0)));  
  
};
Set::Scalar DW(const Set::Scalar theta) const
{
    //sigma'(theta)=(n*sigma1*sin(n*(theta-theta0))*cos(n*(theta-theta0)))/fabs(sin(n*(theta-theta0)))
    //n=2:
    return (2*sigma1*sin(2*(theta-theta0))*cos(2*(theta-theta0)))/fabs(sin(2*(theta-theta0)));
};
Set::Scalar DDW(const Set::Scalar theta) const 
{
    //sigma''(theta)=-n^2*sigma1*fabs(sin(n*(theta-theta0)))
    //n=2:
    return -4*sigma1*fabs(sin(2*(theta-theta0)));
};
 
private:
    Set::Scalar theta0 = NAN, sigma0 = NAN, sigma1 = NAN;
 
public:    
    static void Parse(AbsSin & value, amrex::ParmParse & pp)
    {
        pp_query("theta0",value.theta0); // Angle offset (degrees)
        value.theta0 *= 0.01745329251; // convert degrees into radians
        pp_query("sigma0",value.sigma0); // Minimum energy
        pp_query("sigma1",value.sigma1); // Energy multiplier
    }
 
};
}
}
}
 
#endif