将Eigen :: DiagonalMatrix添加到Eigen :: SparseMatrix的最有效方法是什么

Question

我很惊讶地发现无法简单地向Eigen::DiagonalMatrix添加Eigen::SparseMatrix：

#include <Eigen/Core>
#include <Eigen/Sparse>

int main()
{
  Eigen::Matrix4d M;
  M<<1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16;
  Eigen::SparseMatrix<double> S = M.sparseView();
  Eigen::DiagonalMatrix<double,Eigen::Dynamic> D = 
    Eigen::VectorXd::LinSpaced(4,1,4).asDiagonal();
  Eigen::SparseMatrix<double> T = D+S;
}

不会编译：

 error: invalid operands to binary expression ('Eigen::DiagonalMatrix<double, Eigen::Dynamic>' and 'Eigen::SparseMatrix<double>')
  Eigen::SparseMatrix<double> T = D+S;

似乎在Eigen 3.3中我可以使用：

  T.diagonal() = D.diagonal() + T.diagonal();

但这不能编译为较老的（更稳定的）Eigen 3.2。*：

error: no viable overloaded '='
  T.diagonal() = D.diagonal() + T.diagonal();

我是否真的必须迭代条目并逐个插入？

 Eigen::SparseMatrix<double> T = S;
 for(int i = 0;i<T.rows();i++) T.coeffRef(i,i) += D.diagonal()(i);