eval: fix evaluating undifferentiated ψ in cylindrical coordinates

1 files changed, 1 insertions, 1 deletions

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;