search

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

illegal value in interpolation #220

Closed gokhalen closed 1 year ago

gokhalen commented 1 year ago

I'm developing some parametric DMD code based on parametric heat conduction. When I have 4 training points, everything seem to run ok. But with 25 training points, I'm getting the following error

 ** On entry to DGEBAL parameter number  3 had an illegal value
 ** On entry to DGEHRD parameter number  2 had an illegal value
 ** On entry to DORGHR parameter number  2 had an illegal value
 ** On entry to DHSEQR parameter number  4 had an illegal value

from the call 

CAROM::getParametricDMD(dmd_interp,param_vectors,dmd_model_out,desired_param,
                                "G","LS",closest_rbf_val);

Any advice would be appreciated. Is there a choice of interpolation method that avoids these LAPACK/BLAS routines?

Edit: There seem to be no NaNs in my data, so I'm not sure where these illegal values are coming from. Edit 2: I tried all combinations of RBF type and interpolation type and I get the illegal values in all of them.

-Nachiket

gokhalen commented 1 year ago

This is my fault; I was feeding it illegal values because of a subtle error.