Fix RoundOff problems

When printing Rationals with pystencils at the moment the number of digits which will be printed is only 15. This is not enough for DP and thus Error is introduced in the printed code...

The reason for this is the following line: str(expr.evalf().num)

Here it is important to know what is happening. First of all the Rational is evaluated with Sympy with 15 digits, because SymPy defaults to 15 digits when using evalf: https://docs.sympy.org/latest/modules/evalf.html

but the story continues: .num is not a SymPy function !!!! This comes from mpmath so the Rational is converted to an sp.Float and then to a mpmath.ctx_mp_python.mpf. And yet again mpmath provides also a notion of precision. This could be set with mpmath.mp.dps = 30 for example. However, why should this be done in that way here anyway?

In this MR the .num function is removed and all evalf functions in pystencils are set to evalf(17). Note here that the SymPy printer on the other hand knows very well which precision it needs to provide. Thus when using a float inside the Assignment SymPy prints this correctly.

pystencils/backends/cbackend.py

@@ -441,7 +441,7 @@ class CustomSympyPrinter(CCodePrinter):
     def _print_Pow(self, expr):
         """Don't use std::pow function, for small integer exponents, write as multiplication"""
         if not expr.free_symbols:
-            return self._typed_number(expr.evalf(), get_type_of_expression(expr.base))
+            return self._typed_number(expr.evalf(17), get_type_of_expression(expr.base))
 
         if expr.exp.is_integer and expr.exp.is_number and 0 < expr.exp < 8:
             return f"({self._print(sp.Mul(*[expr.base] * expr.exp, evaluate=False))})"
@@ -452,7 +452,7 @@ class CustomSympyPrinter(CCodePrinter):
 
     def _print_Rational(self, expr):
         """Evaluate all rationals i.e. print 0.25 instead of 1.0/4.0"""
-        res = str(expr.evalf().num)
+        res = str(expr.evalf(17))
         return res
 
     def _print_Equality(self, expr):