Assigning dirichlet boundary conditions directly to a function

Question

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!

Answer 1

Try

bc->apply(u.vector())

Comments

Very nice, thank you!