gsALMRiks.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsStructuralAnalysis/src/gsALMSolvers/gsALMBase.h>
 
namespace gismo
{
 
template <class T>
class gsALMRiks : public gsALMBase<T>
{
 
    typedef gsALMBase<T> Base;
 
    typedef typename Base::ALResidual_t  ALResidual_t;
    typedef typename Base::Jacobian_t    Jacobian_t;
    typedef typename Base::dJacobian_t   dJacobian_t;
 
public:
 
    using Base::setLength;
 
protected:
 
    using Base::defaultOptions;
    using Base::getOptions;
    using Base::computeJacobian;
    using Base::computeResidual;
    using Base::computeResidualNorms;
    using Base::computeUt;
    using Base::computeUbar;
    using Base::computeStability;
    using Base::computeLength;
 
public:
 
    gsALMRiks(  const Jacobian_t  &Jacobian,
                const ALResidual_t&ALResidual,
                const gsVector<T> &Force )
    : Base(Jacobian,ALResidual,Force)
    {
        defaultOptions();
        getOptions();
 
        initMethods();
    }
 
    gsALMRiks(  const dJacobian_t &dJacobian,
                const ALResidual_t&ALResidual,
                const gsVector<T> &Force )
    : Base(dJacobian,ALResidual,Force)
    {
        defaultOptions();
        getOptions();
 
        initMethods();
    }
 
public:
    T distance(const gsVector<T>& DeltaU, const T DeltaL) const
    {
        return math::pow(m_phi * math::pow(DeltaU.norm(),2.0) + (1.0-m_phi) * math::pow(DeltaL,2.0),0.5);
    }
 
protected:
 
    void initMethods();
    void initiateStep();
    void iterationFinish();
 
    void quasiNewtonPredictor();
    void quasiNewtonIteration();
 
    void predictor();
    void predictorGuess();
    void iteration();
 
    void initOutput();
    void stepOutput();
 
protected:
 
    // Number of degrees of freedom
    using Base::m_numDof;
 
    using Base::m_jacobian;
    using Base::m_djacobian;
    using Base::m_residualFun;
    using Base::m_forcing;
 
    using Base::m_options;
 
    using Base::m_numIterations;
 
    using Base::m_maxIterations;
 
    using Base::m_arcLength;
    using Base::m_arcLength_prev;
 
    using Base::m_verbose;
 
    using Base::m_note;
 
    using Base::m_converged;
 
    using Base::m_residueF;
 
    using Base::m_residueU;
 
    using Base::m_residueL;
 
    using Base::m_indicator;
    using Base::m_negatives;
 
    using Base::m_relax;
 
    // Previous update
    using Base::m_DeltaUold;
    using Base::m_DeltaLold;
    using Base::m_U;
    using Base::m_Uprev;
    using Base::m_Uguess;
    using Base::m_DeltaU;
    using Base::m_deltaUbar;
    using Base::m_deltaUt;
    using Base::m_deltaU;
 
    using Base::m_L;
    using Base::m_Lprev;
    using Base::m_Lguess;
    using Base::m_DeltaL;
    using Base::m_deltaL;
    using Base::m_deltaLs;
 
    using Base::m_jacMat;
    using Base::m_detKT;
 
    T m_phi;
};
 
} // namespace gismo
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsALMRiks.hpp)
#endif