gmpFemTransientSolver.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND,
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR
** A PARTICULAR PURPOSE.
**
************************************************************************/
 
#ifndef _GEMA_PLUGIN_FEM_TRANSIENTSOLVER_H_
#define _GEMA_PLUGIN_FEM_TRANSIENTSOLVER_H_
 
#include "gmpFemSolver.h"
 
class GmpFemTransientSolver : public GmpFemSolver
{
public:
  GmpFemTransientSolver(GmElementMesh* mesh, GmSimulationData* simulation, const QVector<GmpFemPhysics*>& physics,
                        GmNumSolver* solver, bool nonlinear, const GmpFemSolverOptions& options, const GmLogCategory& logger);
 
  virtual ~GmpFemTransientSolver();
 
  virtual bool init();
  bool runStep(double dt, bool updateDof = false);
  bool runLinearStep(double dt, int iter, bool updateDof = false);
 
  bool calcLinearResidual(double dt, double* rnorm, double* maxNodeDiff, double* avgNodeDiff);
 
  bool nonlinear() const { return _nonlinear; }
 
  virtual bool addStateItemsToGroup(GmStateDump* state, int groupId);
 
protected:
  virtual bool initElementSets(GmNumSolver* solver);
  virtual bool cleanup();
 
  bool step(double dt, int iter, bool updateDof, bool updateBcs);
  bool istep(double dt, int iter, bool updateDof, bool updateBcs);
  bool estep(double dt, bool updateDof,bool updateBcs);
  bool prepareLinearSystem(double dt);
  bool prepareExplicitLinearSystem(double dt);
  bool GuyanReductionSolve();
 
private:
  bool     _implicit;  
  bool     _nonlinear; 
  GmVector _oldx;      
};
 
#endif