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