QuasiStaticModel.h

//  Copyright (c) 2019 Prashant K. Jha
//  Copyright (c) 2019 Patrick Diehl
//
//  Distributed under the Boost Software License, Version 1.0. (See accompanying
//  file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
 
#ifndef MODEL_QSMODEL_H
#define MODEL_QSMODEL_H
 
#include <hpx/config.hpp>
#include "model/model.h"
#include "inp/decks/modelDeck.h"
 
// forward declaration of class
namespace fe {
class MassMatrix;
class Mesh;
// class Quadrature;
}// namespace fe
 
namespace geometry {
class Fracture;
class InteriorFlags;
class Neighbor;
class VolumeCorrection;
} // namespace geometry
 
namespace inp {
struct ModelDeck;
struct RestartDeck;
struct OutputDeck;
struct SolverDeck;
class Input;
class Policy;
} // namespace inp
 
namespace loading {
class InitialCondition;
class ULoading;
class FLoading;
struct BCData;
} // namespace loading
 
namespace material {
namespace pd {
class BaseMaterial;
}
} // namespace material
 
namespace util {
class StateBasedHelperFunctions;
} // namespace util
 
namespace data {
class DataManager;
} // namespace data
 
namespace rw {
namespace writer {
class Writer;
}
}
 
namespace model {
 
template<class T>
class QuasiStaticModel: public Model {
 
public:
 
    QuasiStaticModel(inp::Input *deck);
 
    ~QuasiStaticModel();
 
private:
 
    void initHObjects();
 
    void computeForces(bool full=false);
 
    void computePertubatedForces(size_t thread); 
 
    void assembly_jacobian_matrix();
 
    void assembly_jacobian_matrix_part(size_t begin, size_t end, size_t thread);
 
    util::VectorXi newton_step(util::VectorXi &res);
 
    void solver();
 
    util::VectorXi computeResidual();
 
    void removeRow(util::Matrixij &matrix, size_t rowToRemove);
 
    void removeCol(util::Matrixij &matrix, size_t colToRemove);
 
    void removeRow(util::VectorXi &vector, size_t rowToRemove);
 
    size_t d_nnodes;
 
    double d_time;
 
    std::string d_name;
 
    size_t d_osThreads;
 
    util::Matrixij jacobian;
 
    std::vector<data::DataManager*> d_dataManagers;
 
    inp::ModelDeck *d_modelDeck_p;
 
    inp::OutputDeck *d_outputDeck_p;
 
    inp::Input *d_input_p;
 
    material::pd::BaseMaterial *d_material_p;
 
    data::DataManager *d_dataManager_p;
 
};
 
}
 
#endif