Hi,

I would like to construct a series of matrix which similar to each other.
I tried to construct like below. However it cannot work. Would you pls help to give suggestion?
Thank you.

Code: Select all

SpMat C[n](Sum_nonzeros_C, Sum_nonzeros_C);
      vector<Tri> CT;
      CT.reserve(Sum_nonzeros_C);

for(int i=1; i<=n;i++)
{
for(t=0;t<Sum_nonzeros_C;t++)
         CT.push_back(Tri(IC[i][t], JC[i][t], PC[i][t]));
         C[i].setFromTriplets(CT.begin(), CT.end());
}

The ligne SpMat "C[n](Sum_nonzeros_C, Sum_nonzeros_C);" is not legal C++, better use a std::vector:

vector<SpMat> C(n, SpMat(Sum_nonzeros_C, Sum_nonzeros_C) );

Then, I guess the main loop should be:

for(int i=0; i<n;i++)

and you probably also need to clear CT at each iteration. It is also strange that the number of rows/columns are the same as the number of non-zeros.