Consider proper matrix interface for Eigen expressions::impl::AsDiagonalMatrixExpr<T> returned by as_diagonal_matrix(). At present it only returns diagonal values (m(i,j) fails for i != j). If it returned zeros for non-diagonal elements, it could be used as a regular matrix (e.g. in ops::product()).
@fnrizzi
Consider proper matrix interface for Eigen
expressions::impl::AsDiagonalMatrixExpr<T>returned byas_diagonal_matrix(). At present it only returns diagonal values (m(i,j)fails fori != j). If it returned zeros for non-diagonal elements, it could be used as a regular matrix (e.g. inops::product()).related to #304