NewtonRaphson.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] Finite Element Procedures, Bathe, K.J., Chapter 8: pages 755-759. 
//       Prentice-Hall, 1996. 
//
// Description: 
//------------------------------------------------------------------------------

#ifndef _NEWTONRAPHSON_HPP_
#define _NEWTONRAPHSON_HPP_

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

#include "Mesh.hpp"
#include "LoadCombo.hpp"
#include "Algorithm.hpp"
#include "Integrator.hpp"
#include "LinearSystem.hpp"

class NewtonRaphson : public Algorithm {

    public:
        NewtonRaphson(std::unique_ptr<LinearSystem> &solver, std::shared_ptr<Mesh> &mesh, double tol=1E-6, unsigned int nIters=50, unsigned int flag=1,  double ldFactor=0.02);

        ~NewtonRaphson();

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

        const Eigen::VectorXd& GetDisplacementIncrement();    

        void SetLoadFactor(double factor);

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

    protected:
        double Tolerance;

        double LoadFactor;

        unsigned int nMaxIterations;

        Eigen::VectorXd dU;
        
        std::unique_ptr<LinearSystem> theSolver;

        std::shared_ptr<Integrator> theIntegrator;
};

#endif