Line data    Source code

       1              : #ifndef MODEL_INTERFACE_GB_ABSSIN_H
       2              : #define MODEL_INTERFACE_GB_ABSSIN_H
       3              : 
       4              : #include <iostream>
       5              : #include <fstream>
       6              : 
       7              : #include "AMReX.H"
       8              : #include "GB.H"
       9              : #include "Set/Set.H"
      10              : #include "Util/Util.H"
      11              : 
      12              : namespace Model
      13              : {
      14              : namespace Interface
      15              : {
      16              : namespace GB
      17              : {
      18              : class AbsSin : public GB
      19              : {
      20              : public:
      21              :     static constexpr const char* name = "abssin";
      22              : 
      23              : AbsSin() {};
      24              : AbsSin(IO::ParmParse &pp) {pp_queryclass(*this);};
      25            0 : AbsSin(IO::ParmParse &pp, std::string name) {pp_queryclass(name,*this);};
      26              : 
      27              : AbsSin(amrex::Real _theta0, amrex::Real _sigma0, amrex::Real _sigma1)
      28              : {       
      29              :     Define(_theta0,_sigma0,_sigma1);
      30              : };
      31              : void Define(amrex::Real _theta0, amrex::Real _sigma0, amrex::Real _sigma1)
      32              : {
      33              :     theta0 = _theta0;
      34              :     sigma0 = _sigma0;
      35              :     sigma1 = _sigma1;
      36              : };
      37              : void Randomize()
      38              : {
      39              :     theta0 = Util::Random()*Set::Constant::Pi;
      40              :     sigma0 = Util::Random();
      41              :     sigma1 = Util::Random();
      42              : };
      43            0 : Set::Scalar W(const Set::Scalar theta) const
      44              : {
      45              :     //sigma(theta)=sigma0+sigma1*fabs(sin(n*(theta-theta0)))
      46              :     //n=2:  
      47            0 :     return sigma0+sigma1*fabs(sin(2*(theta-theta0)));  
      48              :   
      49              : };
      50            0 : Set::Scalar DW(const Set::Scalar theta) const
      51              : {
      52              :     //sigma'(theta)=(n*sigma1*sin(n*(theta-theta0))*cos(n*(theta-theta0)))/fabs(sin(n*(theta-theta0)))
      53              :     //n=2:
      54            0 :     return (2*sigma1*sin(2*(theta-theta0))*cos(2*(theta-theta0)))/fabs(sin(2*(theta-theta0)));
      55              : };
      56            0 : Set::Scalar DDW(const Set::Scalar theta) const 
      57              : {
      58              :     //sigma''(theta)=-n^2*sigma1*fabs(sin(n*(theta-theta0)))
      59              :     //n=2:
      60            0 :     return -4*sigma1*fabs(sin(2*(theta-theta0)));
      61              : };
      62              :  
      63              : private:
      64              :     Set::Scalar theta0 = NAN, sigma0 = NAN, sigma1 = NAN;
      65              : 
      66              : public:    
      67            0 :     static void Parse(AbsSin & value, amrex::ParmParse & pp)
      68              :     {
      69            0 :         pp_query("theta0",value.theta0); // Angle offset (degrees)
      70            0 :         value.theta0 *= 0.01745329251; // convert degrees into radians
      71            0 :         pp_query("sigma0",value.sigma0); // Minimum energy
      72            0 :         pp_query("sigma1",value.sigma1); // Energy multiplier
      73            0 :     }
      74              : 
      75              : };
      76              : }
      77              : }
      78              : }
      79              : 
      80              : #endif