1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 | inline mxd::t m_mulV(mxd::t _m , vxd::t _v){ vxd::t v_(_m.rows()); v_.setZero(); i::t idx = 0; while( idx < _m.rows()){ co(idx, "idx :"); vxd::b y_top = v_.segment( idx, 1); mxd::b m_top = _m.block ( idx, 0, 1, _m.cols() ); y_top = m_top * _v; idx++; } return v_; } | cs |
matrix and vector multiplication with y = Ax + y notation which abstract the linear combination
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | //matrix and vector multiplication with Ax puls Y notation inline mxd::t m_mulV(mxd::t _m , vxd::t _v){ vxd::t v_(_m.rows()); v_.setZero(); i::t idx = 0; while(idx < _m.rows()){ v_ = v_axpy(_v(idx), _m.col(idx), v_); idx++; } return v_; } | cs |