I replaced
Ch = C_n + dt/g_mid*(f- div(phi_mid))
with
expr = dolfinx.fem.Expression(C_n + dt/g_mid*(f- div(phi_mid)),V.element.interpolation_points())
Ch.interpolate(expr)
Is this on the right track?
I replaced
Ch = C_n + dt/g_mid*(f- div(phi_mid))
with
expr = dolfinx.fem.Expression(C_n + dt/g_mid*(f- div(phi_mid)),V.element.interpolation_points())
Ch.interpolate(expr)
Is this on the right track?