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();
    }
 
    void setZero()
    {
        m_solver->reset();
        m_U = gsVector<T>::Zero(m_solver->numDofs());
        m_L = 0;
        m_solver->setLoad(m_L);
    }
 
protected:
    mutable gsStaticBase<T> * m_solver;
    T m_L;
    bool first;
 
    gsVector<T> m_U;
 
    gsStatus m_status;
};
 
 
} // namespace gismo