Exterior facet assembly does not detect the sparsity pattern
Affects | Status | Importance | Assigned to | Milestone | |
---|---|---|---|---|---|
DOLFIN |
Fix Released
|
Medium
|
Garth Wells |
Bug Description
Hi,
I suspect that the exterior facet assembly does not detect the sparsity pattern.
The code to reproduce this is reported below.
I cannot find any reference to the case "exterior facets" in SparsityPattern
http://
May you please validate or invalidate this bug?
Thanks for your help.
Corrado
-------
from dolfin import *
set_log_
mesh = UnitSquare(4,4)
V = VectorFunctionS
u = TrialFunction(V)
v = TestFunction(V)
# This create a matrix with zeros
assemble(
# This create a matrix without zeros
assemble(
-------
Changed in dolfin: | |
assignee: | nobody → Garth Wells (garth-wells) |
milestone: | none → 1.1.0 |
Changed in dolfin: | |
status: | Confirmed → Fix Committed |
Changed in dolfin: | |
status: | Fix Committed → Fix Released |
A hackish work around is to assemble a a matrix using the first form and then reuse the matrix while assemble the second one.