mass matrix between two field PtrList #552
Conversation
To alleviate RAM issue with large meshes
Sparse mass matrix
| for (int i=0; i<V.size(); i++) | ||
| { | ||
| M += V(i) * F.transpose().topRows(Msize).col(i) * F.leftCols(Msize).row(i); | ||
| } |
There was a problem hiding this comment.
Why this change? I think the vectorized version should be considerably faster.
There was a problem hiding this comment.
I think that I had some memory issue for large mesh (mesh size ~ 10^6).
I will double check.
And anyway, I should probably do something that does not slow down the algorithm in the "small mesh case".
| if (consider_volumes) | ||
| { | ||
| Eigen::VectorXd V = getMassMatrixFV(modes[0]); | ||
| for (int i=0; i<V.size(); i++) |
There was a problem hiding this comment.
I think also this one can be written as the previous one
Changing M computation to handle RAM issues for high resolution runs
|
Florian corrected this merge request. |
|
I close this merge request and I will create a new one, from a dedicated branch of my fork (up to date with ITHACA-FV/master) |
No description provided.