Assembly over interior facets is not yet working in parallel

Bug #922322 reported by Simon Funke
This bug affects 3 people
Affects Status Importance Assigned to Milestone

Bug Description

I am trying to write a parallel DG solver. Unfortunately, when I run the Python script with "mpirun -n 2 ... " I get following message:

*** -------------------------------------------------------------------------
*** Error: Unable to perform operation in parallel.
*** Reason: Assembly over interior facets is not yet working in parallel. Consider filing a bug report at
*** Where: This error was encountered inside log.cpp.
*** -------------------------------------------------------------------------

I searched the mailing list, but the last post about this issue seems to be from 2009:

The following code is a reduced example of my code. It runs fine in serial and produces the above error when run with mpirun on 2 processors.

from dolfin import *

def p1dgp2(mesh):
    "Return a function space U*H on mesh from the rt0 space."
    V = VectorFunctionSpace(mesh, 'DG', 1, dim=2)
    H = FunctionSpace(mesh, 'CG', 2)
    return W

def construct_div(W):
    (v, q)=TestFunctions(W)
    (u, h)=TrialFunctions(W)
    return C

mesh = Rectangle(0, 0, 1, 1, 10, 10)

solve(G==rhs, tmpstate)

Revision history for this message
Anders Logg (logg) wrote :

This is a known bug, or rather we have not yet implemented parallel assembly for DG methods. It should come in 1.1.

Changed in dolfin:
status: New → Confirmed
milestone: none → 1.1.0
importance: Undecided → High
Changed in dolfin:
milestone: 1.1.0 → trunk
To post a comment you must log in.
This report contains Public information  
Everyone can see this information.

Other bug subscribers

Related questions

Remote bug watches

Bug watches keep track of this bug in other bug trackers.