Integrator.hpp

Go to the documentation of this file.

//==============================================================================
//
//                       Seismo Virtual Laboratory
//             Module for Serial and Parallel Analysis of seismic 
//         wave propagation and soil-structure interaction simulation
//         Copyright (C) 2018-2021, The California Institute of Technology
//                         All Rights Reserved.
//
// Commercial use of this program without express permission of the California
// Institute of Technology, is strictly  prohibited. See  file "COPYRIGHT"  in
// main  directory  for  information on  usage  and  redistribution, and for a
// DISCLAIMER OF ALL WARRANTIES.
//
//==============================================================================
//
// Written by:
//   Danilo S. Kusanovic (dkusanov@caltech.edu)
//   Elnaz E. Seylabi    (elnaze@unr.edu)
//
// Supervised by:
//   Domniki M. Asimaki  (domniki@caltech.edu)
//
// References : 
//   [1]
//
// Description:
//This file contains the "Integrator object" declarations, and defines how the 
//dynamic solution between two steps is going to be performed.
//------------------------------------------------------------------------------

#ifndef _INTEGRATOR_HPP_
#define _INTEGRATOR_HPP_

#include <memory>
#include <Eigen/Dense>
#include <Eigen/SparseCore>

#include "Mesh.hpp"
#include "LoadCombo.hpp"

class Algorithm;


class Integrator{

    public:
        Integrator(const std::shared_ptr<Mesh> &mesh);

        virtual ~Integrator() = 0;

        virtual void Initialize(std::shared_ptr<Mesh> &mesh) = 0;

        virtual void SetLoadCombination(std::shared_ptr<LoadCombo> &combo) = 0;

        virtual void SetAlgorithm(std::shared_ptr<Algorithm> &algorithm) = 0;

        virtual const Eigen::VectorXd& GetDisplacements() = 0;    

        virtual const Eigen::VectorXd& GetVelocities() = 0;

        virtual const Eigen::VectorXd& GetAccelerations() = 0;

        virtual const Eigen::VectorXd& GetPMLHistoryVector() = 0;

        virtual bool ComputeNewStep(std::shared_ptr<Mesh> &mesh, unsigned int k=0) = 0;

        virtual Eigen::VectorXd ComputeReactionForce(std::shared_ptr<Mesh> &mesh, unsigned int k=0) = 0;

        virtual Eigen::VectorXd ComputeProgressiveForce(std::shared_ptr<Mesh> &mesh, unsigned int k=0) = 0;

        virtual void ComputeSupportMotionVector(std::shared_ptr<Mesh> &mesh, Eigen::VectorXd &Feff, double factor=1.00, unsigned int k=0) = 0;

        virtual void ComputeEffectiveForce(std::shared_ptr<Mesh> &mesh, Eigen::VectorXd &Feff, double factor=1.00, unsigned int k=0) = 0;

        virtual void ComputeEffectiveStiffness(std::shared_ptr<Mesh> &mesh, Eigen::SparseMatrix<double> &Keff) = 0;

    protected:
        Eigen::VectorXd SupportMotion;

        Eigen::SparseMatrix<double> Total2FreeMatrix;

    private:
        //No private variables.
};

#endif