SolverCD.hpp

Go to the documentation of this file.

// This file is part of GranOO, a workbench for DEM simulation.
//
// Author(s)    : - Damien Andre       IRCER/UNILIM, Limoges France
//                  <damien.andre@unilim.fr>
//                - Jean-luc Charles   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jean-luc.charles@ensam.eu>
//                - Jeremie Girardot   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jeremie.girardot@ensam.eu>
//                - Cedric Hubert      LAMIH/UPHF, Valenciennes France
//                  <cedric.hubert@uphf.fr>
//                - Ivan Iordanoff     Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <ivan.iordanoff@ensam.eu>
//
// Copyright (C) 2008-2019 D. Andre, JL. Charles, J. Girardot, C. Hubert, I. Iordanoff
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program.  If not, see <http://www.gnu.org/licenses/>.
 
 
#ifndef _LibPlugIn_SolverCD_HPP
#define _LibPlugIn_SolverCD_HPP
 
#include <string>
#include <iostream>
#include <fstream>
 
#include "GranOO3/FEM/Element.hpp"
 
#include "GranOO3/Core/NeedSetOf.hpp"
#include "GranOO3/Core/PlugIn.hpp"
 
namespace GranOO3
{
  namespace PlugIn
  {
    class SolverCD : public Core::PlugInInterface<SolverCD>,
             public GranOO3::Core::NeedSetOf<GranOO3::Physic::Node>
    {
    public:
      DECLARE_STD_GRANOO_PLUGIN(SOLVER-CD)
 
      SolverCD();
      ~SolverCD() {};
 
      void parse_xml();
      void init();
      void run();
 
    private:
      void assemble_mass_matrix();
 
    private:
      Eigen::SparseMatrix<double> _MMat;
      Eigen::SparseVector<double> _lumpedMMat;
      Eigen::VectorXd _ooDispVector;   // global displacement vector (at t-dt)
      Eigen::VectorXd _oDispVector;    // global displacement vector (at t)
//        Eigen::SparseVector<double> _ooDispVector; // global displacement vector (at t-dt)
//        Eigen::SparseVector<double> _oDispVector;  // global displacement vector (at t)
      unsigned int _totalDofsCount;
    };
  }
}
 
#endif