gsAffineFunction.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsCore/gsLinearAlgebra.h>
#include <gsCore/gsFunction.h>
 
namespace gismo {
 
template <typename T>
class gsAffineFunction : public gsFunction<T>
{
protected:
    gsMatrix<T> m_mat;
    gsMatrix<T> m_box1, m_box2;
    gsVector<T> m_trans;
public:
    typedef memory::shared_ptr< gsAffineFunction > Ptr;
 
    typedef memory::unique_ptr< gsAffineFunction > uPtr;
 
    gsAffineFunction(const gsAffineFunction &other)
        : m_mat(other.m_mat), m_trans(other.m_trans)
    {}
 
    GISMO_CLONE_FUNCTION(gsAffineFunction)
 
    
    gsAffineFunction(const gsMatrix<T> mat, const gsVector<T> trans);
 
    gsAffineFunction(const gsVector<index_t> &directions, const gsVector<bool> &orientation, const gsMatrix<T> &box1, const gsMatrix<T> &box2);
 
    static uPtr make(const gsMatrix<T> mat, const gsVector<T> trans)
    { return uPtr(new gsAffineFunction(mat, trans)); }
 
    static uPtr make(const gsVector<index_t> &directions, const gsVector<bool> &orientation, const gsMatrix<T> &box1, const gsMatrix<T> &box2)
    { return uPtr(new gsAffineFunction(directions, orientation, box1, box2)); }
 
    const gsMatrix<T> & matrix() const { return m_mat; }
    const gsVector<T> & translation() const { return m_trans; }
    
    virtual short_t domainDim() const;
    virtual short_t targetDim() const;
    virtual void eval_into(const gsMatrix<T>& u, gsMatrix<T>& result) const;
    virtual void eval_component_into(const gsMatrix<T>& u,
                                     const index_t comp,
                                     gsMatrix<T>& result) const;
    virtual void deriv_into(const gsMatrix<T>& u, gsMatrix<T>& result) const;
    virtual void deriv2_into( const gsMatrix<T>& u, gsMatrix<T>& result ) const;
};
 
 
} // namespace gismo
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsAffineFunction.hpp)
#endif