Re: [eigen] Dead code in Pardiso

[Rasmus Munk Larsen <rmlarsen@xxxxxxxxxx>]

It looks like p_null is uninitialized, so it looks like undefined behavior to me, but I'm not entirely sure. Did you try to remove .twistedby(p_null); ?

On Tue, Mar 22, 2022 at 2:03 PM Edward Lam <edward@xxxxxxxxxx> wrote:

Hi,

At https://gitlab.com/libeigen/eigen/-/blob/master/Eigen/src/PardisoSupport/PardisoSupport.h#L474 (and also at line 538), we have code like this:

void getMatrix(const MatrixType& matrix)
{
  // PARDISO supports only upper, row-major matrices
  PermutationMatrix<Dynamic,Dynamic,StorageIndex> p_null;
  m_matrix.resize(matrix.rows(), matrix.cols());
  m_matrix.template selfadjointView<Upper>() = matrix.template selfadjointView<UpLo>().twistedBy(p_null);
  m_matrix.makeCompressed();
}

As far as I can tell, the call to twistedBy(p_null) in the code is a no-op right? Or does it serve some purpose that I'm missing?

Thanks,

-Edward