diff --git a/Orb_AutFN.jl b/Orb_AutFN.jl
index 0120ba8..a474786 100644
--- a/Orb_AutFN.jl
+++ b/Orb_AutFN.jl
@@ -228,7 +228,13 @@ function orbit_check_propertyT(logger, sett::Settings)
    info(logger, "minimum(P) = $(minimum(P))")
 
    if λ > 0
-      sgap = PropertyT.check_distance_to_positive_cone(Δ, λ, P, 2*sett.radius, tol=sett.tol, rational=false)
+
+      isapprox(eigvals(P), abs(eigvals(P)), atol=sett.tol) ||
+          warn("The solution matrix doesn't seem to be positive definite!")
+     #  @assert P == Symmetric(P)
+      Q = real(sqrtm(Symmetric(P)))
+
+      sgap = PropertyT.check_distance_to_positive_cone(Δ, λ, Q, 2*sett.radius, tol=sett.tol, rational=false)
       if isa(sgap, Interval)
            sgap = sgap.lo
       end