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