gsAdditiveOp.hpp

Go to the documentation of this file.

namespace gismo
{

template<typename T>
void gsAdditiveOp<T>::apply(const gsMatrix<T>& input, gsMatrix<T>& x) const
{
    GISMO_ASSERT( this->rows() == input.rows(), "The dimensions do not fit." );

    x.setZero( input.rows(), input.cols() );

    const index_t n = m_ops.size();
    gsMatrix<T> res_local, corr_local;

    for (index_t i=0; i<n; ++i)
    {
        res_local.noalias() = m_transfers[i]->transpose()*input;
        m_ops[i]->apply(res_local, corr_local);
        x.noalias() += *(m_transfers[i])*corr_local;
    }
}

} // namespace gismo