Use a more efficient relaxation factor.

1 files changed, 1 insertions, 1 deletions

diff --git a/ell_relax.c b/ell_relax.c
index f51c671..ffc4522 100644
--- a/ell_relax.c
+++ b/ell_relax.c
@@ -171,7 +171,7 @@ static void boundaries_apply(EllRelaxContext *ctx)
 int mg2di_ell_relax_step(EllRelaxContext *ctx)
 {
     EllRelaxInternal *priv = ctx->priv;
-    const double cfl_fac = (1.0 / 4.0) * ctx->step[0] * ctx->step[1];
+    const double cfl_fac = (1.0 / 5.0) * ctx->step[0] * ctx->step[1];
 
     for (int idx1 = 0; idx1 < ctx->domain_size[1]; idx1++)
         for (int idx0 = 0; idx0 < ctx->domain_size[0]; idx0++) {