submaps error

Asked by Chaffra Affouda

I am still have problem with submaps error. The example below does not work and I cannot specify boundary conditions for my mixed problem. I have the latest dolfin from the trunk. Can anyone confirm this?

mesh = UnitInterval(10)
Q = FunctionSpace(mesh, 'CG',1)
QV = MixedFunctionSpace([Q,Q,Q])

project(Constant(1.0), Q) #works
project(Constant(1.0),QV.sub(0)) #breaks

*** -------------------------------------------------------------------------
*** Error: Unable to access ownership range of degree of freedom mapping.
*** Reason: Cannot determine ownership range for submaps.
*** Where: This error was encountered inside DofMap.cpp.
*** -------------------------------------------------------------------------

Question information

Language:
English Edit question
Status:
Answered
For:
DOLFIN Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Joachim Haga (jobh) said :
#1

Yes, I can confirm it's an earlier error (of mine) that's snuck back
in through a merge. I'll post a merge request for the fix in the
morning (no access to the code right now).

Joachim.

On 5 March 2012 18:35, Chaffra Affouda
<email address hidden> wrote:
> New question #189750 on DOLFIN:
> https://answers.launchpad.net/dolfin/+question/189750
>
> I am still have problem with submaps error. The example below does not work and I cannot specify boundary conditions for my mixed problem. I have the latest dolfin from the trunk. Can anyone confirm this?
>
> mesh = UnitInterval(10)
> Q = FunctionSpace(mesh, 'CG',1)
> QV = MixedFunctionSpace([Q,Q,Q])
>
> project(Constant(1.0), Q) #works
> project(Constant(1.0),QV.sub(0)) #breaks
>
> *** -------------------------------------------------------------------------
> *** Error: Unable to access ownership range of degree of freedom mapping.
> *** Reason: Cannot determine ownership range for submaps.
> *** Where: This error was encountered inside DofMap.cpp.
> *** -------------------------------------------------------------------------
>
> --
> You received this question notification because you are a member of
> DOLFIN Team, which is an answer contact for DOLFIN.

Revision history for this message
Joachim Haga (jobh) said :
#2

Oops, I have to (partially) retract that. I'll fix the DirichletBC problem,
but the project() failure that you describe is elsewhere:

#3 0xaf56bf86 in dolfin::DofMap::ownership_range (this=0x901e518) at
../../dolfin/fem/DofMap.cpp:283
#4 0xaf577710 in dolfin::AssemblerTools::init_global_tensor (...)
    at ../../dolfin/fem/AssemblerTools.cpp:174
#5 0xaf590462 in dolfin::SystemAssembler::assemble (...)
    at ../../dolfin/fem/SystemAssembler.cpp:162

-j.

On 5 March 2012 22:49, Joachim Berdal Haga <email address hidden> wrote:

> Yes, I can confirm it's an earlier error (of mine) that's snuck back
> in through a merge. I'll post a merge request for the fix in the
> morning (no access to the code right now).
>
> Joachim.
>
> On 5 March 2012 18:35, Chaffra Affouda
> <email address hidden> wrote:
> > New question #189750 on DOLFIN:
> > https://answers.launchpad.net/dolfin/+question/189750
> >
> > I am still have problem with submaps error. The example below does not
> work and I cannot specify boundary conditions for my mixed problem. I have
> the latest dolfin from the trunk. Can anyone confirm this?
> >
> > mesh = UnitInterval(10)
> > Q = FunctionSpace(mesh, 'CG',1)
> > QV = MixedFunctionSpace([Q,Q,Q])
> >
> > project(Constant(1.0), Q) #works
> > project(Constant(1.0),QV.sub(0)) #breaks
> >
> > ***
> -------------------------------------------------------------------------
> > *** Error: Unable to access ownership range of degree of freedom mapping.
> > *** Reason: Cannot determine ownership range for submaps.
> > *** Where: This error was encountered inside DofMap.cpp.
> > ***
> -------------------------------------------------------------------------
> >
> > --
> > You received this question notification because you are a member of
> > DOLFIN Team, which is an answer contact for DOLFIN.
>

Revision history for this message
lvleph (erichlf) said :
#3

I have this exact same bug in the latest nightly
*** -------------------------------------------------------------------------
*** Warning: UnitInterval has been deprecated in DOLFIN version 1.1.0.
*** The class UnitInterval has been replaced by UnitIntervalMesh.
*** -------------------------------------------------------------------------

Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Solving linear system of size 11 x 11 (PETSc Krylov solver).
Traceback (most recent call last):
  File "./test.py", line 9, in <module>
    project(Constant(1.0),QV.sub(0)) #breaks
  File "/Applications/FEniCS.app/Contents/Resources/lib/python2.7/site-packages/dolfin/fem/projection.py", line 104, in project
    form_compiler_parameters=form_compiler_parameters)
  File "/Applications/FEniCS.app/Contents/Resources/lib/python2.7/site-packages/dolfin/fem/assembling.py", line 317, in assemble_system
    assembler.assemble(A_tensor, b_tensor)
RuntimeError:

*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at
***
*** <email address hidden>
***
*** Remember to include the error message listed below and, if possible,
*** include a *minimal* running example to reproduce the error.
***
*** -------------------------------------------------------------------------
*** Error: Unable to access ownership range of degree of freedom mapping.
*** Reason: Cannot determine ownership range for sub-dofmaps.
*** Where: This error was encountered inside DofMap.cpp.
*** Process: 0
*** DOLFIN version: 1.2.0+
*** Git changeset:
*** -------------------------------------------------------------------------

Revision history for this message
Chaffra Affouda (chaffra) said :
#4

fenics has moved. so try submitting on fenics bitbucket.

Revision history for this message
lvleph (erichlf) said :
#5

Thank you, I will.

Can you help with this problem?

Provide an answer of your own, or ask Chaffra Affouda for more information if necessary.

To post a message you must log in.