AdaptiveSetup.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _ADAPTIVE_SETUP_H
#define _ADAPTIVE_SETUP_H
 
#include "MatVec.h"
 
#include <utility>
 
class ScalarFunc;
class SIMoutput;
namespace tinyxml2 { class XMLElement; }
namespace LR { struct RefineData; }
 
 
class AdaptiveSetup
{
public:
  explicit AdaptiveSetup(SIMoutput& sim, bool sa = true);
  virtual ~AdaptiveSetup();
 
  void setAdaptationNorm(size_t g, size_t i = 0) { adaptor = g; adNorm = i; }
 
  bool initPrm(size_t normGroup);
 
  int calcRefinement(LR::RefineData& prm, int iStep,
                     const Vectors& gNorm, const Vector& refIn,
                     int currDofs = -1) const;
 
  bool parse(char* keyWord, std::istream& is);
  bool parse(const tinyxml2::XMLElement* elem);
 
  void printNorms(const Vectors& gNorm, const Vectors& dNorm,
                  const Matrix& eNorm, size_t w = 36,
                  bool printModelNorms = true) const;
 
  bool writeMesh(int iStep) const;
 
  size_t eIdx() const { return eRow; }
 
protected:
  SIMoutput& model; 
 
  size_t adaptor; 
  size_t adNorm;  
  size_t eRow;    
  double rCond;   
 
  using DblIdx = std::pair<double,int>;
 
  struct ErrorInfo
  {
    double limit = 0.0; 
    double sumErr = 0.0; 
    double curErr = 0.0; 
  };
 
  enum Threshold { NONE=0, MAXIMUM, AVERAGE, MINIMUM, TRUE_BETA,
                   DORFEL, SYMMETRIZED };
 
  static ErrorInfo errorLimits(const std::vector<DblIdx>& errors,
                               const double beta,
                               const Threshold threshold);
 
private:
  bool   alone;      
  bool   linIndep;   
  double beta;       
  ScalarFunc* betaFunc;
  double errTol;     
  double condLimit;  
  int    maxStep;    
  int    maxDOFs;    
  int    knot_mult;  
  int    maxTjoints; 
  double maxAspect;  
  bool   closeGaps;  
  double symmEps;    
 
  enum RefScheme { FULLSPAN=0, MINSPAN=1,
                   ISOTROPIC_FUNCTION=2,
                   ISOTROPIC_ELEMENT=3 };
 
  Threshold   threshold; 
  RefScheme   scheme;    
 
  int         storeMesh; 
  std::string mshPrefix; 
  std::string errPrefix; 
};
 
#endif