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 QList<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