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