On Fri, 1 Aug 2008, Roy Stogner wrote:
> On Fri, 1 Aug 2008, Adam C Powell IV wrote:
>> This is a neat approach, which I could use myself in a separate problem.
>> One issue though: how do I get, say, T to be 0 on domain A, when it is
>> non-zero at nodes on the A-B boundary? Do you have any examples?
> No examples. I'd try to add a new constraint equation for each "fake"
Oh, wait, were you asking how to determine which DoFs are "fake" and
which aren't? Yeah, that's a bit trickier, since you can't just
assume that a variable B DoF is fake just because it's also on an
element from domain A. I'd create a vector<bool> indexed by dof
number, loop over all elements and all variables while setting entries
to "true" if you find a DoF from variable A on an element from domain
A (or B on B, etc), then loop over all DoFs adding a constraint for
every DoF whose entry is still set to "false".