Implementation Pressure entropy SPH

Here is a new flavour of SPH that works and gives better result than the default Gadget-2 scheme.

The only change affecting the rest of the code is the addition of the condition in engine_init_particles() to create a second density loop over the particles for the special conversions that are required in this scheme. Other than that, it just fits in the rest of the framework.

The make check target is failing with:

testInteractions.c: In function ‘dump_indv_particle_fields’:
testInteractions.c:113:65: error: ‘struct part’ has no member named ‘rho_dh’
           p->a_hydro[0], p->a_hydro[1], p->a_hydro[2], p->rho, p->rho_dh,

That's with HOPKINS_PE_SPH defined.

But if I configure --disable-vec that works.

Other than that I'm happy to accept. Do you want to keep the branch?

Added 1 commit:

• a3c53f13 - Make testInteractions compile with PE-SPH.

I have just fixed the test-case to make it compile. It's not yet part of the test-suite as James is still working on it but is still being compiled. Just forgot about it.

The branch can go.

Thanks !

Status changed to merged

mentioned in commit 5c46ced3