Daniel Bauer
--- a/src/pystencils/sympyextensions/integer_functions.py
+++ b/src/pystencils/sympyextensions/integer_functions.py

-    def __new__(cls, integer, divisor):
-        if is_integer_sequence((integer, divisor)):
-            return integer // divisor if integer % divisor == 0 else (integer // divisor) + 1
-        else:
-            return super().__new__(cls, integer, divisor)
+    def _eval_op(self, arg1, arg2):
+        from ..utils import c_intdiv

-    #   TODO: Implement this in FreezeExpressions
-    # def to_c(self, print_func):
-    #     dtype = collate_types((get_type_of_expression(self.args[0]), get_type_of_expression(self.args[1])))
-    #     assert dtype.is_int()
-    #     code = "( ({0}) % ({1}) == 0 ? ({dtype})({0}) / ({dtype})({1}) : ( ({dtype})({0}) / ({dtype})({1}) ) +1 )"
-    #     return code.format(print_func(self.args[0]), print_func(self.args[1]), dtype=dtype)
+        return sp.sympify(c_intdiv(arg1 + arg2 - 1, arg2))