LCEquanti.h

Go to the documentation of this file.

 
#ifndef LCEQUANTI_H
#define LCEQUANTI_H
 
#include "LCEbreed.h"
#include "ttquanti.h"
 
//LCE_QuantiInit
//
class LCE_QuantiInit : public virtual LifeCycleEvent
{
  
  TMatrix _init_values, _init_freq;
  bool _doByTraitValue, _doByAlleleFreq;
  unsigned int _nTraits, _nLoci;
  TProtoQuanti* _QuantiProto;
  
public:
  
  LCE_QuantiInit ( );
  
  virtual ~LCE_QuantiInit ( ) { }
  
  virtual void execute ();
  void init_trait_value(sex_t SEX, age_idx AGE, unsigned int size, unsigned int deme, double *values);
  void init_allele_freq(sex_t SEX, age_idx AGE, unsigned int size, unsigned int deme, double *values, const TMatrix& all_val);
  
  virtual LifeCycleEvent* clone ( ) {return new LCE_QuantiInit();}
  
  virtual bool setParameters ();
  
  //SimComponent implementation:
  virtual void loadFileServices ( FileServices* loader ) {}
  virtual void loadStatServices ( StatServices* loader ) {}
  virtual bool resetParameterFromSource (std::string param, SimComponent* cmpt) {return false;}
  virtual age_t removeAgeClass ( ) {return 0;}
  virtual age_t addAgeClass ( ) {return 0;}
  virtual age_t requiredAgeClass () {return 0;}
  
};
 
//LCE_QuantiModifier
//
class LCE_QuantiModifier : public virtual LifeCycleEvent
{
 
  TProtoQuanti*        _QuantiProto;
  /*A vector storing the target heritability of the trait(s).*/
  vector<double>       _h2;
  /*An array to hold the multiple occurrences of the modifier.*/
  TMatrix              _occurrences;
  /*The index of the next time the modifier must be called.*/
  unsigned int         _next_occurrence;
  /*A flag to indicate whether the modifier is called at each generation.*/
  bool                 _is_permanent;
  /*A flag to indicate on which age class the modifier is applied.*/
  age_t                _age_flag;
 
  vector<double> (TTQuantiSH::* _get_VA_func_ptr) (Patch*, const age_idx);
 
 
public:
 
  LCE_QuantiModifier ( );
 
  virtual ~LCE_QuantiModifier ( ) { }
 
  virtual void execute ();
 
  void   setVefromVa (Patch* patch, age_idx AGE);
 
  void   setIndPhenotype (Patch* patch, sex_t SEX, age_idx AGE);
 
 
  virtual LifeCycleEvent* clone ( ) {return new LCE_QuantiModifier();}
 
  virtual bool setParameters ();
 
  //SimComponent implementation:
  virtual void loadFileServices ( FileServices* loader ) {}
  virtual void loadStatServices ( StatServices* loader ) {}
  virtual bool resetParameterFromSource (std::string param, SimComponent* cmpt) {return false;}
  virtual age_t removeAgeClass ( ) {return 0;}
  virtual age_t addAgeClass ( ) {return 0;}
  virtual age_t requiredAgeClass () {return 0;}
 
};
 
// Class LCE_Breed_Quanti
//
class LCE_Breed_Quanti : public virtual LCE_Breed_base
{
 
  TProtoQuanti*        _QuantiProto;
  vector<double>       _h2;
 
  vector<double> (TTQuantiSH::* _get_VA_func_ptr) (Patch*, const age_idx);
 
  void (LCE_Breed_Quanti::*_myPopModelFuncPtr)(void);
 
public:
 
  LCE_Breed_Quanti ( ) : LifeCycleEvent("breed_quanti","quant"), _QuantiProto(0),
  _get_VA_func_ptr(0), _myPopModelFuncPtr(0) { }
 
  virtual ~LCE_Breed_Quanti   ( ) { }
 
  void   setVefromVa (Patch* patch, age_idx AGE);
 
  void   setIndPhenotype (Patch* patch, sex_t SEX, age_idx AGE);
 
  void   NonWrightFisherPopulation ();
 
  void   WrightFisherPopulation    ();
 
  virtual bool setParameters ();
  virtual void  execute ();
 
  virtual LifeCycleEvent* clone ( ) {return new LCE_Breed_Quanti();}
 
  virtual void loadFileServices ( FileServices* loader ) {}
  virtual void loadStatServices ( StatServices* loader ) {}
  virtual bool resetParameterFromSource (std::string param, SimComponent* cmpt) {return false;}
  virtual age_t removeAgeClass ( ) {return 0;}
  virtual age_t addAgeClass ( ) {return OFFSPRG;}
  virtual age_t requiredAgeClass () {return ADULTS;}
};
 
 
#endif