[op] corrected bug in several definition domains

src/core/operators/codac2_abs.h

@@ -55,7 +55,7 @@ namespace codac2
   {
     return {
       fwd(x1.a),
-      x1.def_domain && x1.a != 0. // def domain of the derivative of abs
+      x1.def_domain && !x1.a.contains(0.) // def domain of the derivative of abs
     };
   }
 
@@ -72,7 +72,7 @@ namespace codac2
       fwd(x1.m),
       fwd(x1.a),
       d,
-      x1.def_domain && x1.a != 0. // def domain of the derivative of abs
+      x1.def_domain && !x1.a.contains(0.) // def domain of the derivative of abs
     };
   }
 

src/core/operators/codac2_acos.h

@@ -56,7 +56,7 @@ namespace codac2
     return {
       fwd(x1.a),
       x1.a.is_subset({-1,1}) // def domain of acos
-      && x1.a != 1. // def domain of the derivative of acos
+      && !x1.a.contains(1.) // def domain of the derivative of acos
       && x1.def_domain
     };
   }
@@ -75,7 +75,7 @@ namespace codac2
       fwd(x1.a),
       d,
       x1.a.is_subset({-1,1}) // def domain of acos
-      && x1.a != 1. // def domain of the derivative of acos
+      && !x1.a.contains(1.) // def domain of the derivative of acos
       && x1.def_domain
     };
   }

src/core/operators/codac2_arith_div.h

@@ -82,7 +82,7 @@ namespace codac2
   {
     return {
       fwd(x1.a, x2.a),
-      x1.def_domain && x2.def_domain && x2.a != 0. // def domain of the derivative of div
+      x1.def_domain && x2.def_domain && !x2.a.contains(0.) // def domain of the derivative of div
     };
   }
 
@@ -101,7 +101,7 @@ namespace codac2
       fwd(x1.m, x2.m),
       fwd(x1.a, x2.a),
       d,
-      x1.def_domain && x2.def_domain && x2.a != 0. // def domain of the derivative of div
+      x1.def_domain && x2.def_domain && !x2.a.contains(0.) // def domain of the derivative of div
     };
   }
 
@@ -128,7 +128,7 @@ namespace codac2
   {
     return {
       fwd(x1.a, x2.a),
-      x1.def_domain && x2.def_domain && x2.a != 0. // def domain of the derivative of div
+      x1.def_domain && x2.def_domain && !x2.a.contains(0.) // def domain of the derivative of div
     };
   }
 
@@ -149,7 +149,7 @@ namespace codac2
       fwd(x1.m, x2.m),
       fwd(x1.a, x2.a),
       d,
-      x1.def_domain && x2.def_domain && x2.a != 0. // def domain of the derivative of div
+      x1.def_domain && x2.def_domain && !x2.a.contains(0.) // def domain of the derivative of div
     };
   }
 
@@ -169,7 +169,7 @@ namespace codac2
   {
     return {
       fwd(x1.a, x2.a),
-      x1.def_domain && x2.def_domain && x2.a != 0. // def domain of the derivative of div
+      x1.def_domain && x2.def_domain && !x2.a.contains(0.) // def domain of the derivative of div
     };
   }
 
@@ -191,7 +191,7 @@ namespace codac2
       fwd(x1.m, x2.m),
       fwd(x1.a, x2.a),
       d, 
-      x1.def_domain && x2.def_domain && x2.a != 0. // def domain of the derivative of div
+      x1.def_domain && x2.def_domain && !x2.a.contains(0.) // def domain of the derivative of div
     };
   }
 

src/core/operators/codac2_asin.h

@@ -56,7 +56,7 @@ namespace codac2
     return {
       fwd(x1.a),
       x1.a.is_subset({-1,1}) // def domain of asin
-      && x1.a != 1. // def domain of the derivative of asin
+      && !x1.a.contains(1.) // def domain of the derivative of asin
       && x1.def_domain
     };
   }

src/core/operators/codac2_atan2.h

@@ -58,7 +58,7 @@ namespace codac2
     return {
       fwd(x1.a, x2.a),
       x1.def_domain && x2.def_domain
-      && !(x1.a == 0. && x2.a == 0.) // def domain of the derivative of atan2
+      && !(x1.a.contains(0.) && x2.a.contains(0.)) // def domain of the derivative of atan2
     };
   }
 
@@ -79,7 +79,7 @@ namespace codac2
       fwd(x1.a, x2.a),
       d,
       x1.def_domain && x2.def_domain
-      && !(x1.a == 0. && x2.a == 0.) // def domain of the derivative of atan2
+      && !(x1.a.contains(0.) && x2.a.contains(0.)) // def domain of the derivative of atan2
     };
   }
 

src/core/operators/codac2_log.h

@@ -59,7 +59,7 @@ namespace codac2
     return {
       fwd(x1.a),
       x1.a.is_subset({0,oo}) // def domain of log
-      && x1.a != 0. // def domain of the derivative of log
+      && !x1.a.contains(0.) // def domain of the derivative of log
       && x1.def_domain
     };
   }
@@ -78,7 +78,7 @@ namespace codac2
       fwd(x1.a),
       d,
       x1.a.is_subset({0,oo}) // def domain of log
-      && x1.a != 0. // def domain of the derivative of log
+      && !x1.a.contains(0.) // def domain of the derivative of log
       && x1.def_domain
     };
   }

src/core/operators/codac2_max.h

@@ -69,7 +69,7 @@ namespace codac2
       fwd(x1.a, x2.a),
       d,
       x1.def_domain && x2.def_domain
-        && (x1.a != x2.a) // def domain of the derivative of max
+        && !x1.a.intersects(x2.a) // def domain of the derivative of max
     };
   }
 

src/core/operators/codac2_sign.h

@@ -55,7 +55,7 @@ namespace codac2
   {
     return {
       fwd(x1.a),
-      x1.def_domain && x1.a != 0. // def domain of the derivative of sign
+      x1.def_domain && !x1.a.contains(0.) // def domain of the derivative of sign
     };
   }
 
@@ -68,7 +68,7 @@ namespace codac2
       fwd(x1.m),
       fwd(x1.a),
       IntervalMatrix::zero(1,x1.da.size()),
-      x1.def_domain && x1.a != 0. // def domain of the derivative of sign
+      x1.def_domain && !x1.a.contains(0.) // def domain of the derivative of sign
     };
   }
 

src/core/operators/codac2_sqrt.h

@@ -56,7 +56,7 @@ namespace codac2
     return {
       fwd(x1.a),
       x1.a.is_subset({0,oo}) // def domain of sqrt
-      && x1.a != 0. // def domain of the derivative of sqrt
+      && !x1.a.contains(0.) // def domain of the derivative of sqrt
       && x1.def_domain
     };
   }
@@ -75,7 +75,7 @@ namespace codac2
       fwd(x1.a),
       d,
       x1.a.is_subset({0,oo}) // def domain of sqrt
-      && x1.a != 0. // def domain of the derivative of sqrt
+      && !x1.a.contains(0.) // def domain of the derivative of sqrt
       && x1.def_domain
     };
   }

src/core/operators/codac2_tan.h

@@ -55,7 +55,7 @@ namespace codac2
   {
     return {
       fwd(x1.a),
-      x1.def_domain && cos(x1.a) != 0. // def domain of the derivative of tan
+      x1.def_domain && !cos(x1.a).contains(0.) // def domain of the derivative of tan
     };
   }
 
@@ -72,7 +72,7 @@ namespace codac2
       fwd(x1.m),
       fwd(x1.a),
       d,
-      x1.def_domain && cos(x1.a) != 0. // def domain of the derivative of tan
+      x1.def_domain && !cos(x1.a).contains(0.) // def domain of the derivative of tan
     };
   }