I have to assign a scalar function that "blurs" a boundary between a lot of subdomains e.g. takes a value of 1.0 for the nodes shared by two subdomains and vanishes rapidly - kind of "smeared interface". Shall I use Expression() projected on CG space or a MeshFunction()? I'll use these values in my LHS. No idea how to do that...
