Comment 4 for bug 1077912
Revision history for this message
| Johan Hake (johan-hake) wrote Re: [Bug 1077912] Re: Incorrect DirichletBC behaviour with assemble_system | #4 |
© 2004
Canonical Ltd.
•
Terms of use
•
Data privacy
•
Contact Launchpad Support
•
Blog
•
Careers
•
System status
•
bbfa235
(Get the code!)
On 11/12/2012 02:27 PM, Martin Sandve Alnæs wrote:
> The zero rows is only a symptom of the underlying bug. Looking at the
> sparsity pattern of B you can see there are dramatic changes between
> setting enable_bug=1/0, so this is not an assemble_system specific
> issue.
Sure but then this would be a separate bug, compared to what you mention
below.
> If DirichletBC with method=topological only applies to dofs on the
> boundary, then the markers in the interior facets should not matter.
> However they do get used, causing the boundary condition to be applied
> to interior dofs as well as boundary dofs in the enable_bug=1 case in
> the provided code.
>
> While this may very well be useful for some situations, it is contrary
> to the description in the documentation as quoted, and contrary to the
> name "boundary conditions".
Agree. Not sure how to deal with it as setting dirichlet dofs inside the
domain can be very useful.
Johan