From a70b93fd016adb557d55db771fac488be20e8ff1 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Wed, 1 Aug 2018 17:03:52 +0200 Subject: Use a more efficient relaxation factor. --- ell_relax.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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++) { -- cgit v1.2.3