Adjoint solution with PointSource as source term (f) for Poisson equation

Wrap this with mesh = Mesh(generate_mesh(domain_geometry, resolution)), as mshr is not overloaded by dolfin-adjoint.

1 Like