rge.h

Go to the documentation of this file.

 
#ifndef RGE_H
#define RGE_H
 
#include "utils.h"
class DoubleVector; 
#include "linalg.h"
#include "numerics.h"
 
namespace softsusy {
  class Approx {
  private: 
    int loops; 
    int thresholds; 
  public: 
    Approx(): loops(2), thresholds(0) {};
    virtual ~Approx() {};
 
    inline const Approx & displayApprox() const { return *this; }
    const Approx & operator=(const Approx &);
    Approx(const Approx & a);
    Approx(int l, int t); 
 
    void setLoops(int l) { loops = l; };
    void setThresholds(int t) { thresholds = t; };
    
    int displayLoops() const { return loops; };
    int displayThresholds() const { return thresholds; };  
  };
 
class RGE {
private:
  double mu; 
  int numPars; 
public:
  RGE(): mu(0.0), numPars(0) {};
  virtual ~RGE() {};
  
  void setMu(double e) { mu = e; };
  virtual void setmu(double f) { setMu(f); };
  void setPars(int i) { numPars = i; };
 
  double displayMu() const { return mu; };
  int howMany() const { return numPars; };
  inline const RGE & displayRGE() const { return *this; }
  virtual const DoubleVector display() const = 0;
 
  virtual void set(const DoubleVector & s) = 0;
  virtual DoubleVector beta() const = 0;
 
  virtual int run(double from, double to, double eps= -1.0);
  virtual int runto(double to, double eps = -1.0);
  
  int callRK(double x1, double x2, DoubleVector & v, 
             DoubleVector (*derivs)(double, const DoubleVector &), 
             double eps = -1.0);
 
};  
 
DoubleVector allDerivs(double, const DoubleVector &);
 
} // namespace softsusy
 
#endif