NewmarkSIM.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _NEWMARK_SIM_H
#define _NEWMARK_SIM_H
 
#include "MultiStepSIM.h"
 
 
class NewmarkSIM : public MultiStepSIM
{
public:
  explicit NewmarkSIM(SIMbase& sim);
 
  using MultiStepSIM::parse;
  virtual bool parse(const tinyxml2::XMLElement* elem);
 
  virtual void printProblem(bool stopInputTimer = false) const;
 
  virtual void initPrm();
 
  bool initAcc(double zero_tolerance = 1.0e-8, std::streamsize outPrec = 0);
 
  virtual bool advanceStep(TimeStep& param, bool updateTime = true);
 
  virtual SIM::ConvStatus solveStep(TimeStep& param,
                                    SIM::SolutionMode = SIM::STATIC,
                                    double zero_tolerance = 1.0e-8,
                                    std::streamsize outPrec = 0);
 
  SIM::ConvStatus solveIteration(TimeStep& param);
 
  int getMaxit() const { return maxit; }
 
  virtual bool isDynamic() const { return true; }
 
protected:
  virtual bool solutionNorms(const TimeDomain&, double zero_tolerance = 1.0e-8,
                             std::streamsize outPrec = 0);
 
  virtual SIM::ConvStatus checkConvergence(TimeStep& param);
  virtual bool predictStep(TimeStep& param);
  virtual bool correctStep(TimeStep& param, bool = false);
  virtual void finalizeRHSvector(bool);
 
public:
  const Vector& getVelocity() const { return solution[1]; }
  const Vector& getAcceleration() const { return solution[2]; }
 
  using MultiStepSIM::saveStep;
  bool saveStep(int iStep, TimeStep& param);
 
  virtual void dumpResults(double time, utl::LogStream& os,
                           std::streamsize precision = 3,
                           bool formatted = true) const;
 
protected:
  // Time integration parameters
  double alpha1; 
  double alpha2; 
  double beta;   
  double gamma;  
 
  // Solution algorithm parameters
  bool   solveDisp; 
  char   predictor; 
  int    maxit;     
  int    maxIncr;   
  int    nupdat;    
  int    saveIts;   
  char   saveVelAc; 
  double rTol;      
  double aTol;      
  double divgLim;   
 
  unsigned short int cNorm; 
 
  std::string factor_file;   
  bool write_factor;         
  bool read_factor;          
  LinAlg::StorageFormat fmt; 
 
public:
  static const char* inputContext; 
};
 
#endif