Re: [eigen] SparseQR compilation error

[Gael Guennebaud <gael.guennebaud@xxxxxxxxx>]

not really because it fails on that line of SparseQR:

m_R = tempR * m_pivotperm;

where m_R and tempR are SparseMatrix<double> and m_pivotperm is a permutation matrix. The block expressions are triggered in the permutation product code. It would be interesting to see if the above _expression_ alone leads to the same error.

gael

On Mon, Aug 12, 2013 at 1:56 PM, Pavel Holoborodko <pavel@xxxxxxxxxxxxxxx> wrote:

Of course they are different :).

I meant the way they use components which compiler doesn't like.. 

I assume they both using Eigen::Block (?), in QR this leads to error but in LU it is fine.

Maybe this fact could help in adjusting QR to avoid problem with compiler (as in LU). 

That was just an idea.

On Mon, Aug 12, 2013 at 8:45 PM, Gael Guennebaud <gael.guennebaud@xxxxxxxxx> wrote:

On Mon, Aug 12, 2013 at 1:04 PM, Pavel Holoborodko <pavel@xxxxxxxxxxxxxxx> wrote:

I reported this issue on Intel C++ forum.  

From my experience Windows version of ICC seems to have many problems (as newer as more buggy).

Will probably try to roll back to older version.

ok, let's see.

 

How much SparseLU different from SparseQR architecture-wise? 

SparseLU has same interface and works just fine. Maybe we can do SparseQR in similar way ;)? 

beside the interface and the elimination tree, they have absolutely nothing in common! 

gael

 

On Mon, Aug 12, 2013 at 7:51 PM, Gael Guennebaud <gael.guennebaud@xxxxxxxxx> wrote:

Hi,

as already observed on the forum, this error seems to show up only with ICC on Windows. No problem with ICC 11, 12, 13 on Linux. Moreover, I don't understand this issue as the base class "Eigen::BlockImpl<const Eigen::SparseMatrix<double, 0, int>, -1, 1, 1, Eigen::Sparse>" does has an operator= line 194 in src/SparseCore/SparseBlock..h. So I'm clueless here.

gael

On Sat, Aug 10, 2013 at 11:23 AM, Pavel Holoborodko <pavel@xxxxxxxxxxxxxxx> wrote:

Hi, 

I am testing solvers for sparse systems.

Intel C++ compiler chokes on this code:

typedef Eigen::SparseMatrix< double > SparseDoubleMatrix;

SparseQR< SparseDoubleMatrix, COLAMDOrdering<int> > solver;

SparseDoubleMatrix A(10,10);

solver.compute(A);

Compiler points to solver.compute(A) with error message:

class "Eigen::BlockImpl<const Eigen::SparseMatrix<double, 0, int>, -1, 1, 1, Eigen::Sparse>" has no member "operator="

Error persist for any OrderingType.

Interestingly, error only shows for SparseQR. If I use SparseLU - everything compiles fine.

Full template instantiation chain lead to error is in attachment.

Would appreciate any help & pointers how to solve the problem. 

Thank you,

Pavel.