NonLinSIM.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _NON_LIN_SIM_H
#define _NON_LIN_SIM_H
 
#include "MultiStepSIM.h"
 
 
class NonLinSIM : public MultiStepSIM
{
public:
  enum CNORM { NONE_UPTAN=-1, NONE=0, L2=1, L2SOL=2, ENERGY=3 };
 
  explicit NonLinSIM(SIMbase& sim, CNORM n = ENERGY);
  virtual ~NonLinSIM();
 
  void setConvNorm(CNORM n) { if ((iteNorm = n) <= NONE) fromIni = true; }
 
  void init(size_t nSol, const RealArray& initVal);
 
  virtual void initPrm();
 
  virtual bool advanceStep(TimeStep& param, bool updateTime = true);
 
  SIM::ConvStatus solve(double zero_tol = 1.0e-8, std::streamsize outPrec = 0);
 
  virtual SIM::ConvStatus solveStep(TimeStep& param,
                                    SIM::SolutionMode mode = SIM::STATIC,
                                    double zero_tolerance = 1.0e-8,
                                    std::streamsize outPrec = 0);
 
  SIM::ConvStatus solveIteration(TimeStep& param);
 
  int getMaxit() const { return maxit; }
 
  virtual bool isLinear() const { return iteNorm <= NONE; }
 
protected:
  void printWorst(utl::LogStream& os, double eps);
 
  virtual SIM::ConvStatus checkConvergence(TimeStep& param);
  virtual bool updateConfiguration(TimeStep& param);
  virtual bool lineSearch(TimeStep& param);
 
  virtual bool assembleSystem(const TimeDomain& time, const Vectors& pSol,
                              bool newLHSmatrix = true, bool poorConvg = false);
 
public:
  virtual bool parse(char* keyWord, std::istream& is);
 
  virtual bool parse(const tinyxml2::XMLElement* elem);
 
protected:
  bool   fromIni; 
  CNORM  iteNorm; 
  double rTol;    
  double aTol;    
  double divgLim; 
  double eta;     
  double alpha;   
  double alphaO;  
  int    maxit;   
  int    maxIncr; 
  int    nupdat;  
  int    prnSlow; 
  bool   saveExL; 
  bool   updNewN; 
 
  std::map<int,int> slowNodes; 
 
public:
  static const char* inputContext; 
};
 
#endif