From 4a33d00d0de5baa851e3ac6457190c5dd259aca5 Mon Sep 17 00:00:00 2001
From: Nacime Bouziani <48448063+nbouziani@users.noreply.github.com>
Date: Sat, 9 Dec 2023 20:04:02 +0000
Subject: [PATCH] Alignment with new interpolation behaviour

---
 ngsPETSc/utils/firedrake.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/ngsPETSc/utils/firedrake.py b/ngsPETSc/utils/firedrake.py
index c4805a2..5bc065c 100644
--- a/ngsPETSc/utils/firedrake.py
+++ b/ngsPETSc/utils/firedrake.py
@@ -85,7 +85,7 @@ def curveField(self, order):
     low_order_element = self.coordinates.function_space().ufl_element().sub_elements[0]
     element = low_order_element.reconstruct(degree=order)
     space = fd.VectorFunctionSpace(self, fd.BrokenElement(element))
-    newFunctionCoordinates = fd.interpolate(self.coordinates, space)
+    newFunctionCoordinates = fd.assemble(fd.interpolate(self.coordinates, space))
 
     #Computing reference points using fiat
     fiat_element = newFunctionCoordinates.function_space().finat_element.fiat_equivalent