RadBlastWave test

[kenshen04]

I'm starting to implement rad hydro for a project and tried to run the RadBlastWave test in 2D. After a few timesteps, I got a bunch of Riemann solve errors (e.g., "WARNING: (rho e)_l < 0 or pl < small_pres in Riemann: 624383.9679 -1498276.932 1e-100") and then "amrex::Abort::14::State has NaNs in the density component::check_for_nan() !!!"

Before I dig too far, I just wanted to check: should the test run successfully out of the box, in which case maybe I installed hypre incorrectly or something like that? Or has the test not been validated in a while?

Thanks,
Ken

[zingale]

what inputs file are you using? I just ran in 2D with inputs.2d for 128 steps, but then it had hypre solver issues.

[kenshen04]

I also used inputs.2d but encountered the same error with inputs.kp.2d. If it's helpful, I compiled with gcc 11.4 and hypre 2.32 and both with and without OpenMP.

[kenshen04]

I was able to get it run by cranking up the temperature in the ambient medium to 3e4 K, and the results look reasonable. I guess that was easier on the Riemann solver? But in any case, I've convinced myself that hypre is working correctly!