Algorithm.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:
//------------------------------------------------------------------------------

#ifndef _ALGORITHM_HPP_
#define _ALGORITHM_HPP_

#include <map>
#include <mpi.h>
#include <memory>
#include <Eigen/Dense>
#include <Eigen/SparseCore>

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

class Integrator;

class Algorithm{

    public:
        Algorithm(const std::shared_ptr<Mesh> &mesh, unsigned int flag=1, double NormFactor=1.0);

        virtual ~Algorithm() = 0;

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

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

        virtual void SetIntegrator(std::shared_ptr<Integrator> &integrator) = 0;

        virtual void SetLoadFactor(double factor) = 0;

        void ReverseStatesIncrements(std::shared_ptr<Mesh> &mesh);

        void UpdateStatesIncrements(std::shared_ptr<Mesh> &mesh, const Eigen::VectorXd &dU);

        void ReducedParallelResidual(const Eigen::VectorXd &Feff, double &Residual);

        double ComputeConvergence(const Eigen::VectorXd &Force, const Eigen::VectorXd &Delta, const Eigen::VectorXd &delta, unsigned int k);

    protected:
        Eigen::SparseMatrix<double> Total2FreeMatrix;

    private:
        unsigned int flag;

        double NormFactor;
};

#endif