I have a vector of DirichletBCs,
std::vector<const DirichletBC*>
which work just fine when used with a solver. I would like to be able to update the Function with the DirichletBCs directly, without going through a solver. I have been trying to find a way to interpolate the Expression within each of the DirichletBCs, within its respective SubDomain, onto my Function.
How might I go about this?
I was hoping to find something like
u->interpolate(Expression, SubDomain);
Thanks!