From f2ff50e9ef476d2113290d07068cb527671bebd6 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Felix=20K=C3=B6hler?= <27728103+Ceyron@users.noreply.github.com> Date: Tue, 22 Oct 2024 13:38:55 +0200 Subject: [PATCH 1/2] Have general diffusion coefficient consisten with norm/diff --- exponax/stepper/generic/_polynomial.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/exponax/stepper/generic/_polynomial.py b/exponax/stepper/generic/_polynomial.py index 372611b..d9067e5 100644 --- a/exponax/stepper/generic/_polynomial.py +++ b/exponax/stepper/generic/_polynomial.py @@ -18,7 +18,7 @@ def __init__( num_points: int, dt: float, *, - coefficients: tuple[float, ...] = (10.0, 0.0, 0.01), + coefficients: tuple[float, ...] = (10.0, 0.0, 1.0), polynomial_scales: tuple[float, ...] = (0.0, 0.0, 10.0), order=2, dealiasing_fraction: float = 2 / 3, From d28c5c88506f8012002a641d3616c342da3fb5f8 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Felix=20K=C3=B6hler?= <27728103+Ceyron@users.noreply.github.com> Date: Tue, 22 Oct 2024 13:39:20 +0200 Subject: [PATCH 2/2] Use correct sign on quadratic nonlinearity --- exponax/stepper/generic/_polynomial.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/exponax/stepper/generic/_polynomial.py b/exponax/stepper/generic/_polynomial.py index d9067e5..2780569 100644 --- a/exponax/stepper/generic/_polynomial.py +++ b/exponax/stepper/generic/_polynomial.py @@ -19,7 +19,7 @@ def __init__( dt: float, *, coefficients: tuple[float, ...] = (10.0, 0.0, 1.0), - polynomial_scales: tuple[float, ...] = (0.0, 0.0, 10.0), + polynomial_scales: tuple[float, ...] = (0.0, 0.0, -10.0), order=2, dealiasing_fraction: float = 2 / 3, num_circle_points: int = 16,