Make precision work

demo/MetaData.py

@@ -32,3 +32,5 @@
     + inner(c, c) * inner(grad(u), grad(v)) * dx(1, degree=4)\
     + inner(c, c) * inner(grad(u), grad(v)) * dx(1, degree=2)\
     + inner(grad(u), grad(v)) * dx(1, degree=-1)
+
+L = v * dx(0, metadata={"precision": 1})

ffcx/analysis.py

@@ -202,7 +202,7 @@ def _analyze_form(form: ufl.form.Form, options: typing.Dict) -> ufl.algorithms.f
             p = precisions.pop()
         elif len(precisions) == 0:
             # Default precision
-            p = np.finfo("double").precision + 1  # == 16
+            p = None
         else:
             raise RuntimeError("Only one precision allowed within integrals grouped by subdomain.")
 

ffcx/codegeneration/C/c_implementation.py

@@ -150,11 +150,11 @@ def __init__(self, scalar, precision=None) -> None:
             self.precision = precision
 
     def _format_number(self, x):
-        prec = self.precision
+        p = self.precision
         if isinstance(x, complex):
-            return "({:.{prec}}+I*{:.{prec}})".format(x.real, x.imag, prec=prec)
+            return f"({x.real:.{p}}+I*{x.imag:.{p}})"
         elif isinstance(x, float):
-            return "{:.{prec}}".format(x, prec=prec)
+            return f"{x:.{p}}"
         return str(x)
 
     def _build_initializer_lists(self, values):

ffcx/codegeneration/C/integrals.py

@@ -36,7 +36,7 @@ def generator(ir, options):
     parts = ig.generate()
 
     # Format code as string
-    CF = CFormatter(options["scalar_type"])
+    CF = CFormatter(options["scalar_type"], ir.precision)
     body = CF.c_format(parts)
 
     # Generate generic FFCx code snippets and add specific parts