I am implementing a space time simulation of the Cahn-Hilliard equation
with a stable FE method called AVS-FE. I am using the z dimension as the
time variable of the problem. I am running into issues when implementing a
Dirichlet boundary condition. I need to implement a piecewise continuous
function over one of my boundaries with more than 2 values which is why I
cannot use the conditional statement. I have searched around and found this
has not been done much and am getting errors I do not understand.
Attached is the python script that I run which functions for a simple
boundary condition but not the more complex one I am trying to implement
with an Expression call.
If I could get assistance with this error as to why my expression is not
working properly it would be appreciated.
The error I am getting is:
Traceback (most recent call last):
File "/home/austin/firedrake/src/firedrake/firedrake/interpolation.py",
line 76, in interpolate
assembled_interpolator = self.frozen_assembled_interpolator
AttributeError: 'Interpolator' object has no attribute
'frozen_assembled_interpolator'
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "cahn_hilliard_prototype.py", line 224, in <module>
u_source.interpolate(MyExpression_t0())
File "/home/austin/firedrake/src/firedrake/firedrake/function.py", line
349, in interpolate
return interpolation.interpolate(expression, self, subset=subset)
File "/home/austin/firedrake/src/firedrake/firedrake/interpolation.py",
line 34, in interpolate
return Interpolator(expr, V, subset=subset, access=access).interpolate()
File "/home/austin/firedrake/src/firedrake/firedrake/interpolation.py",
line 79, in interpolate
assembled_interpolator = self.callable()
File "/home/austin/firedrake/src/firedrake/firedrake/interpolation.py",
line 170, in callable
l()
File "/home/austin/firedrake/src/PyOP2/pyop2/base.py", line 3563, in
compute
self._compute(iterset.core_part, fun, *arglist)
File "/home/austin/firedrake/src/PyOP2/pyop2/pyparloop.py", line 114, in
_compute
raise NotImplementedError
NotImplementedError
I don't understand the frozen assembled interpolator part which is the
problem I believe.
Regards,
--
Austin Kaul