[field] Implement Mesh.__eq__()

phi/geom/_mesh.py

@@ -78,6 +78,13 @@ def __variable_attrs__(self):
     def __value_attrs__(self):
         return '_vertices',
 
+    def __eq__(self, other):
+        if self is other:
+            return True
+        if not isinstance(other, Mesh):
+            return False
+        return self._vertices == other._vertices and math.always_close(self._polygons, other._polygons) and self._boundaries == other._boundaries
+
     @property
     def shape(self) -> Shape:
         return shape(self._polygons).non_spatial & channel(self._vertices)