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.
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