Fix all the normal double precision stuff

mcx @ ec2a7273

-Subproject commit 39471f3ba1cc6c3846c80aeec90a07787be195b0
+Subproject commit ec2a72739b884032e7bbd87a9b06aa0dadfbe372

include/teqp/models/cubics.hpp

@@ -141,7 +141,7 @@ public:
             throw std::invalid_argument("Sizes do not match");
         }
         auto b = get_b(T, molefrac);
-        auto Psiminus = -1.0*log(1.0 - b * rho);
+        auto Psiminus = -log(1.0 - b * rho);
         auto Psiplus = log((Delta1 * b * rho + 1.0) / (Delta2 * b * rho + 1.0)) / (b * (Delta1 - Delta2));
         auto val = Psiminus - get_a(T, molefrac) / (Ru * T) * Psiplus;
         return forceeval(val);