Power with a real base and even exponent reports that it's non-negative

Signed-off-by: Vladimir V. Kisil <kisilv@maths.leeds.ac.uk>

diff --git a/ginac/power.cpp b/ginac/power.cpp
index b815d17a9e6e07f716748ff909c81dc080835994..87792b922099215e6b00b9744724eef7ca481ab6 100644 (file)
--- a/ginac/power.cpp
+++ b/ginac/power.cpp
@@ -242,6 +242,8 @@ bool power::info(unsigned inf) const
                        return (flags & status_flags::expanded);
                case info_flags::positive:
                        return basis.info(info_flags::positive) && exponent.info(info_flags::real);
                        return (flags & status_flags::expanded);
                case info_flags::positive:
                        return basis.info(info_flags::positive) && exponent.info(info_flags::real);

+               case info_flags::nonnegative:
+                       return basis.info(info_flags::real) && exponent.info(info_flags::integer) && exponent.info(info_flags::even);

                case info_flags::has_indices: {
                        if (flags & status_flags::has_indices)
                                return true;
                case info_flags::has_indices: {
                        if (flags & status_flags::has_indices)
                                return true;