From 8ad50ae962fea404c4ea7d74b59a51e25400c249 Mon Sep 17 00:00:00 2001
From: Anton Khirnov <anton@khirnov.net>
Date: Fri, 23 Jul 2021 10:28:30 +0200
Subject: eval: fix evaluating undifferentiated ψ in cylindrical coordinates
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 eval.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/eval.c b/eval.c
index b5f153e..9b94141 100644
--- a/eval.c
+++ b/eval.c
@@ -40,7 +40,7 @@ static int eval_psi_cylindric(const BDContext *bd,
 
     double *basis_val_rho = NULL, *basis_val_z = NULL;
 
-    double add = 0.0;
+    double add = (diff_order[0] || diff_order[1]) ? 0.0 : 1.0;
 
     int ret = 0;
 
-- 
cgit v1.2.3