From acdf00b58ae4c1b49406873fc225bebcd5a00032 Mon Sep 17 00:00:00 2001
From: Pedro Gonnet <pedro.gonnet@durham.ac.uk>
Date: Sat, 6 Apr 2013 20:38:31 +0000
Subject: [PATCH] higher-order exp approximation.
Former-commit-id: 7e832b8c860870f5bbba6ab5cb13d1f73781aa39
---
src/engine.c | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/engine.c b/src/engine.c
index b9ae4506ed..67ab021e0d 100644
--- a/src/engine.c
+++ b/src/engine.c
@@ -251,7 +251,7 @@ void engine_map_kick_first ( struct cell *c , void *data ) {
if ( pdt > dt_step ) {
// rho = p->rho *= expf( -3.0f * h_dt / h * dt );
float w = -3.0f * h_dt / h * dt;
- rho = p->rho *= 1.0f + w*( -1.0f + w*( 0.5f - 1.0f/6.0f*w ) );
+ rho = p->rho *= 1.0f + w*( -1.0f + w*( 0.5f + w*(-1.0f/6.0f + 1.0f/24.0*w ) ) );
p->force.POrho2 = u * ( const_gamma - 1.0f ) / ( rho + h * p->rho_dh / 3.0f );
}
else {
@@ -531,7 +531,7 @@ void engine_step ( struct engine *e , int sort_queues ) {
}
e->step += 1;
}
- }
+ }
/* Set the system time. */
e->time = e->dt * e->step;
--
GitLab