Incorrect linear system size for mixed space
I am consistently getting an error from PETSC when I am trying to solve a variational problem of the following type (minimal Python example):
from dolfin import *
mesh=UnitSquare
V1=FunctionSpac
V2=FunctionSpac
V=MixedFunction
f,g=TrialFuncti
h,i=TestFunctio
a=(f*h+g*h+i*f)*dx
L=h*dx
U=Function(V)
solve(a==L,U)
I am on OSX running the nightly binary snapshot. The error message is (note that even the size of the system is incorrect):
Calling FFC just-in-time (JIT) compiler, this may take some time.
Solving linear variational problem.
[0]PETSC ERROR: -------
[0]PETSC ERROR: Arguments are incompatible!
[0]PETSC ERROR: Local size 513 not compatible with block size 2!
[0]PETSC ERROR: -------
[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 3, Wed Aug 29 11:26:24 CDT 2012
[0]PETSC ERROR: See docs/changes/
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: -------
[0]PETSC ERROR: Unknown Name on a darwin11. named leo.local by toxa Mon Jul 29 19:07:16 2013
[0]PETSC ERROR: Libraries linked from /Users/
[0]PETSC ERROR: Configure run at Tue May 14 10:31:34 2013
[0]PETSC ERROR: Configure options --prefix=
[0]PETSC ERROR: -------
[0]PETSC ERROR: PetscLayoutSetB
[0]PETSC ERROR: MatSetBlockSize() line 6704 in /Users/
Is this a known bug/is there a workaround?
Sincerely,
--Anton
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- DOLFIN Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Anton for more information if necessary.