From f6d151bb8b919cd7e75ccb63eeb866a798f1ea0f Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Sebastian=20M=C3=BCller?= <mueller.seb@posteo.de>
Date: Thu, 13 Mar 2025 11:04:51 +0100
Subject: [PATCH] polynomial: correctly determine order of polynomial with only
zeros
---
src/mo_polynomial.f90 | 1 +
1 file changed, 1 insertion(+)
diff --git a/src/mo_polynomial.f90 b/src/mo_polynomial.f90
index e6aee44..3dd9894 100644
--- a/src/mo_polynomial.f90
+++ b/src/mo_polynomial.f90
@@ -82,6 +82,7 @@ contains
n = size(p)
if (n == 0_i4) call error_message("mo_polynomial : p is an empty array.")
+ poly_order = 0_i4
do i = 1_i4, n
if (abs(p(i)) > 0.0_dp) then
poly_order = (n - i)
--
GitLab