PropPowerlaw.H

Go to the documentation of this file.

 
#ifndef _PROPPOWERLAW_H_
#define _PROPPOWERLAW_H_
 
#include "Propagation.H"
#include "../IO/ParmParser.H"
 
namespace criptic {
  namespace propagation {
 
    class PropPowerlaw : public Propagation {
 
    public:
 
      PropPowerlaw(const ParmParser& pp);
 
      virtual ~PropPowerlaw() { };
      
      virtual FieldQtyNeedType fieldQtyNeed() const {
#ifdef TRACK_PITCH_ANGLE
    return noFieldQty;
#else
    if (varStreamDir) return needFieldQtyGrad;
    else return noFieldQty;
#endif
      }
 
      virtual PropagationData
      operator()(const RealVec& x,
         const Real t,
         const gas::GasData& gd,
         const CRPacket& packet,
         const FieldQty& qty,
         const FieldQtyGrad& qtyGrad) const;
      
    private:
      
      Real kPar0;     
      Real kParIdx;   
      Real kPerp0;    
      Real kPerpIdx;  
      Real kPP0;      
      Real kPPIdx;    
#ifndef TRACK_PITCH_ANGLE
      Real vStr0;     
      Real vStrIdx;   
      bool varStreamDir; 
      bool vAStream;  
#else
      Real kMu0;      
      Real kMuIdx;    
#endif
    };
  }
}
 
#endif
// _PROPPOWERLAW_H_