(Temporarily) replacing project with interpolate in solve_forward (#39)

mesh-adaptation · Mar 5, 2024 · aa719f4 · aa719f4
1 parent 670a4e5
commit aa719f4
Showing 1 changed file with 6 additions and 3 deletions.
diff --git a/goalie/mesh_seq.py b/goalie/mesh_seq.py
@@ -8,7 +8,6 @@
 from firedrake.petsc import PETSc
 from firedrake.pyplot import triplot
 from .function_data import ForwardSolutionData
-from .interpolation import project
 from .log import pyrint, debug, warning, info, logger, DEBUG
 from .options import AdaptParameters
 from animate.quality import QualityMeasure
@@ -324,7 +323,9 @@ def get_checkpoints(
                 checkpoints.append(
                     AttrDict(
                         {
-                            field: project(sols[field], fs[i + 1])
+                            field: firedrake.Function(fs[i + 1]).interpolate(
+                                sols[field]
+                            )
                             for field, fs in self._fs.items()
                         }
                     )
@@ -585,7 +586,9 @@ def solve_forward(self, solver_kwargs: dict = {}) -> AttrDict:
             if i < num_subintervals - 1:
                 checkpoint = AttrDict(
                     {
-                        field: project(checkpoint[field], fs[i + 1])
+                        field: firedrake.Function(fs[i + 1]).interpolate(
+                            checkpoint[field]
+                        )
                         for field, fs in self._fs.items()
                     }
                 )