Extract Diagonal sparse and diagonal matrix product in Eigen

Question

I have a sparse matrix and a diagonal matrix from which I would like to calculate a weighted matrix outer product:

typedef Eigen::SparseMatrix S