anandrdbz

Hi @boombiubla, I was able to run your test case for 2500 time steps with dt = 2.3e-7 (I used 10 times coarser grid so that it runs fast) which...

Hi @boombiubla, I have managed to make a 3D shock droplet case that works for this case. You can adjust the Mach number by changing the value of Min in...

@wilfonba check s_rectangle subroutine in m_patches.fpp, it only defines smooth_patch_id when ib is not true (lines 893-894), and you added a further check at line 928 where that variable is...

I'll redo the QBMM tests (using nb=1) along with some of the bug fixes for QBMM + non-polytropic that are yet to be integrated in a separate PR

Is there a reason why these are integer values and not float ?

@wilfonba, is there a example case that can be used to test this on your fork ?

@sbryngelson , but I presume this difference is only present in the body-force problems and not the regular cases like shock droplet or bubble screen ?

> @anandrdbz I just made a push with new two example cases, 2D_rayleigh_taylor and 2D_rising_bubble Okay, I guess this is the STRevert branch right ? Also if other cases had...

Okay sounds good, I'll take a look at the rising bubble problem tonight

The immersed boundary problem at least could potentially be related to the bug I observed with Cray compilers (round-off being significant)