// Membrane - A linear solver for membrane problem using FEM
// Copyright (C) 2010-2026 Eric Bechet
//
// See the LICENSE file for license information and contributions.
// Please report all bugs and problems to <bechet@cadxfem.org>.

#ifndef _MEMBRANE_SOLVER_H_
#define _MEMBRANE_SOLVER_H_

#include <memory>
#include "dofManager.h"
#include "genTerm.h"
#include "genSolver.h"
#include "membraneDomain.h"
#include "genCache.h"
#include "savedGenTermManager.h"
#include "genFSpace.h"


class membraneSolver : public genSolver
{
protected :
  dofManager<double>* pAssembler;
  // function spaces
  genFSpace<genTensor1<double> >::Handle FSpaceDisp;
  genTerm<genTensor2<double>,1>::Handle EpsilonDisp;

  // supports for elastic
  std::vector<membraneDomain*> EDomains;

  // model name
  std::string modelname;

  // genTerm manager
  SavedGenTermManager savedGenTermManager;

public :
  membraneSolver() : genSolver(), pAssembler(0) { }
  virtual ~membraneSolver();
  //read data file
  virtual void readInputFile(const std::string &fileName);
  // builds local supports
  virtual void CheckProperties();
  // build the problem function spaces
  virtual void BuildFunctionSpaces();
  // assemble the problem
  virtual void AssembleRHS();
  virtual void AssembleLHS();
  virtual void BuildLinearSystem();

  // assemble the residu
//  virtual void AssembleResidual(); Has to be implemented souldn't be GBDTAG
  // export K.txt b.txt in file
  virtual void exportKb();
  virtual void solve();
  virtual PView* buildDisplacementView(const std::string &postFileName);
  virtual PView* buildElasticEnergyView(const std::string &postFileName);
};

#endif //_MEMBRANE_SOLVER_H_