gsIterative.h

Go to the documentation of this file.

#pragma once

#include <gsIO/gsOptionList.h>
#include <gsElasticity/gsBaseUtils.h>
#include <functional>

namespace gismo
{

template <class T>
class gsBaseAssembler;
// TODO correct
template <class T>
class gsIterative
{
public:
    typedef memory::shared_ptr<gsIterative> Ptr;
    typedef memory::unique_ptr<gsIterative> uPtr;

    gsIterative(gsBaseAssembler<T> & assembler_);

    gsIterative(gsBaseAssembler<T> & assembler_,
             const gsMatrix<T> & initSolutionVector);

    gsIterative(gsBaseAssembler<T> & assembler_,
             const gsMatrix<T> & initSolutionVector,
             const std::vector<gsMatrix<T> > & initFixedDoFs);

    static gsOptionList defaultOptions();

    gsOptionList & options() { return m_options; }

    void solve();

    bool compute();

    const gsMatrix<T> & solution() const { return solVector; }

    const std::vector<gsMatrix<T> > & allFixedDofs() const { return fixedDoFs; }

    std::string status();

    void reset();

    virtual void setFixedDofs(const std::vector<gsMatrix<T> > & ddofs);

    index_t numberIterations() const {return numIterations;}

    void setSolutionVector(const gsMatrix<T> & solutionVector) { solVector = solutionVector; }

    void saveState();

    void recoverState();

protected:
    gsBaseAssembler<T> & assembler;
    gsMatrix<T> solVector;
    std::vector<gsMatrix<T> > fixedDoFs;
    index_t numIterations; 
    solver_status m_status;  
    T residualNorm; 
    T initResidualNorm; 
    T updateNorm; 
    T initUpdateNorm; 
    gsOptionList m_options;

    gsMatrix<T> solVecSaved;
    std::vector<gsMatrix<T> > ddofsSaved;
};

} // namespace ends

#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsIterative.hpp)
#endif