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" DoF. I don't know how that would perform. We didn't design the constraints data structures with "more constrained than unconstrained DoFs" in mind. But I suppose if we just use it for a solution to compare to, it doesn't matter. > I may have some time to investigate a patch to the DofMap code at some > point, and if the above approach works (to demonstrate a base case > solution) will start looking into it in the next couple of weeks... Thanks. I don't think fixing the distribute_dofs_... code to write per-subdomain variables will be hard, and I could do that if you'd like; I'm just still waffling on the user APIs. The hard (or at least tedious) part may be fixing our I/O classes to write out and read in solutions with per-subdomain variables. I'm not familiar with the nitty-gritty details of our output formats, but I wouldn't be surprised if they didn't all support such a thing. --- Roy ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
