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(m_DeltaU.norm(),2.0) + (1.0-m_phi) * math::pow(m_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