LLNL / libROM

Model reduction library with an emphasis on large scale parallelism and linear subspace methods
https://www.librom.net
Other
198 stars 36 forks source link

StaticSVD hot fix for transposed SVD operation. #235

Closed dreamer2368 closed 11 months ago

dreamer2368 commented 11 months ago

When the sample size is larger than the system dimension, transposed SVD operation is executed in StaticSVD. In this case, ncolumns must be less than or equal to the system dimension. Previously, this is not automatically enforced, unless d_max_basis_dimension is manually specified. This causes the scalapack Bad_Submatrix error in StaticSVD::computeSVD.

Now this is fixed by enforcing ncolumns to be limited by the system dimension, when the transposed SVD is executed.