| 122 | } |
| 123 | |
| 124 | template <typename U> static Matrix diag( int size, const U& diag_value ) { |
| 125 | Matrix res(size,size) ; |
| 126 | res.fill_diag( diag_value ) ; |
| 127 | return res ; |
| 128 | } |
| 129 | |
| 130 | inline Proxy operator[]( R_xlen_t i ) { |
| 131 | return static_cast< Vector<RTYPE>* >( this )->operator[]( i ) ; |