gsControlDisplacement.h

Go to the documentation of this file.

#include <typeinfo>
#include <gsStructuralAnalysis/src/gsStaticSolvers/gsContinuationBase.h>

#pragma once


namespace gismo
{

template <class T>
class gsControlDisplacement : public gsContinuationBase<T>
{
protected:

public:

    gsControlDisplacement(gsStaticBase<T> * solver) :
    m_solver(solver),
    first(true)
    {
        m_U = gsVector<T>::Zero(m_solver->numDofs());
        m_L = 0;
    }

    gsStatus step(T dL)
    {
        m_solver->setLoad(m_L + dL);
        if (!first)
            m_solver->setDisplacement(m_U);
        else first = false;
        m_status = m_solver->solve();

        if (m_status == gsStatus::Success)
        {
            m_L += dL;
            m_U += m_solver->update();
        }
        return m_status;
    }

    gsVector<T> & solutionU()
    {
        return m_U;
    }

    T solutionL()
    {
        return m_L;
    }

    void reset()
    {
        m_solver->reset();
    }

protected:
    mutable gsStaticBase<T> * m_solver;
    T m_L;
    bool first;

    gsVector<T> m_U;

    gsStatus m_status;
};


} // namespace gismo