Created attachment 265 [details]
Patch for the solution I came up with
SuperLUSupport.h does not allow a Eigen::Map class to act as 'x' or 'b' in the solve routine. I had a need to use already allocated data to act as 'x', so I simply added an additional Map structure (note the duplication of names here, I am referring to what was called MaptoSLU in 3.0.4) that can accept an Eigen::Map and convert it to the SLuMatrix datatype.
My solution has a lot of duplication with the Matrix conversion to SLuMatrix code, so there's probably a better way to do this, but as I am not versed in the Eigen source code I didn't want to mess about with pre-existing structures.
Indeed there exist a cleaner and more general solution. We should support any dense expressions and check the other solver for the same shortcoming.
date: 2012-06-06 09:37:59
summary: Fix bug 454: allow Block/Map objects for solving with SuperLU
-- GitLab Migration Automatic Message --
This bug has been migrated to gitlab.com's GitLab instance and has been closed from further activity.
You can subscribe and participate further through the new bug through this link to our GitLab instance: https://gitlab.com/libeigen/eigen/issues/454.