I am trying to run the following function:

Code: Select all

void ScalarField::LinTrans(SparseMatrix<double>* Mat, ScalarField** out){
  VectorXd tempvector(numPoints);                                                     
  tempvector=(*Mat)*(*value);
...
}


where (*value) is a VectorXd. It compiles fine, but when I run it SOMETIMES I get this error:

Code: Select all

.../Eigen/Eigen/src/Core/DenseCoeffsBase.h:391: Eigen::DenseCoeffsBase<Derived, 1>::Scalar& Eigen::DenseCoeffsBase<Derived, 1>::operator[](Eigen::DenseCoeffsBase<Derived, 1>::Index) [with Derived = Eigen::Matrix<double, -0x00000000000000001, 1>, Eigen::DenseCoeffsBase<Derived, 1>::Scalar = double, Eigen::DenseCoeffsBase<Derived, 1>::Index = long int]: Assertion `index >= 0 && index < size()' failed.


This happens about half the times I run it, and the other half it runs fine. Is this a bug? Or am I doing something I shouldn't?

Thanks

I presume you've ensured the dimensions are consistent, i.e. cols of Mat = rows of *value.

please be more specific. running your programming under a debugger to get the trace and the values of index, size(), etc.