Hessian called on single volumes

PrincetonUniversity / SPEC

The Stepped-Pressure Equilibrium Code, an advanced MRxMHD equilibrium solver.

https://princetonuniversity.github.io/SPEC/

GNU General Public License v3.0

25 stars 6 forks source link

Hessian called on single volumes #61

Closed Spiffmeister closed 5 years ago

Spiffmeister commented 5 years ago

Just flagging that trying to run the single volume TestCases (i.e. G2V01L0Fi.001.sp) is failing because the .sp files have LHMatrix = T, setting LHMatrix = F fixes the issue.

The error message is:

** On entry to DGETRF parameter number  4 had an illegal value
Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

I don't know enough about the workings of SPEC but I'm guessing it shouldn't be trying to determine a H matrix on a single volume.

jloizu commented 5 years ago

Yes, that is correct, we should make this small change in the input file.

@abaillod can you do that?

abaillod commented 5 years ago

Done!