In most cases, you would have to redefine the bu_u object, as it applies caching when applying the first boundary condition, see: Bitbucket
This can also be found in the documentation: https://fenicsproject.org/olddocs/dolfin/latest/cpp/d1/d95/classdolfin_1_1DirichletBC.html#details
1 Like