Commit eeb44fd1 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Make the ghost iterations abort based on h_old not on h_init to avoid cases...

Make the ghost iterations abort based on h_old not on h_init to avoid cases where the code can't converge due to very poor and much too small initial guess.
parent 9a99bdeb
......@@ -1402,7 +1402,7 @@ void runner_do_ghost(struct runner *r, struct cell *c, int timer) {
}
/* Check whether the particle has an inappropriate smoothing length */
if (fabsf(h_new - h_old) > eps * h_init) {
if (fabsf(h_new - h_old) > eps * h_old) {
/* Ok, correct then */
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment