From 4b714867ece34b96fa3332c687b91d8bf66e48a8 Mon Sep 17 00:00:00 2001
From: Martin Bauer <martin.bauer@fau.de>
Date: Wed, 5 Dec 2018 13:42:23 +0100
Subject: [PATCH] Conversion from order-parameters to concentration for n-phase

-> last concentration is automatically computed as 1-others
-> makes code more general
---
 fd/spatial.py | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/fd/spatial.py b/fd/spatial.py
index 6b147fe96..43e88eaea 100644
--- a/fd/spatial.py
+++ b/fd/spatial.py
@@ -86,6 +86,8 @@ def discretize_spatial(expr, dx, stencil=fd_stencils_standard):
             stencil = fd_stencils_standard
         elif stencil == 'isotropic':
             stencil = fd_stencils_isotropic
+        elif stencil == 'isotropic_hd':
+            stencil = fd_stencils_isotropic_high_density_code
         else:
             raise ValueError("Unknown stencil. Supported 'standard' and 'isotropic'")
 
-- 
GitLab