dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

The problem is that we have changed the numbering of Nedelec elements
to start at q = 1, not q = 0.

Strange that this was not detected when running the serial tests. I
get the same error here, both in serial and in parallel. I'll commit a
fix soon.

--
Anders


On Wed, Sep 16, 2009 at 08:38:31AM +0200, Johannes Ring wrote:
> The buildslaves still times out, now on the demo/pde/curl-curl/python
> demo. Here is the output:
>
> fenics-slave@byggmester:python$ mpirun -n 2 python demo.py
> Process 0: UnitSphere is experimental. It may be of poor quality mesh
> Process 0: Debug: creating local mesh data [at
> dolfin/mesh/MeshPartitioning.cpp:62 in partition()]
> Process 1: Debug: creating local mesh data [at
> dolfin/mesh/MeshPartitioning.cpp:62 in partition()]
> Process 1: Debug: check [at dolfin/mesh/LocalMeshData.cpp:228 in
> receive_mesh_data()]
> Process 0: Number of global vertices: 729
> Process 0: Number of global cells: 3072
> Process 0: Debug: check [at dolfin/mesh/LocalMeshData.cpp:158 in
> broadcast_mesh_data()]
> Process 0: Debug: check [at dolfin/mesh/LocalMeshData.cpp:172 in
> broadcast_mesh_data()]
> Process 0: Sending 365 vertices to process 0, range is (0, 365)
> Process 0: Sending 364 vertices to process 1, range is (365, 729)
> Process 1: Debug: check [at dolfin/mesh/LocalMeshData.cpp:240 in
> receive_mesh_data()]
> Process 1: Received 364 vertex coordinates
> Process 1: Debug: check [at dolfin/mesh/LocalMeshData.cpp:248 in
> receive_mesh_data()]
> Process 0: Received 365 vertex coordinates
> Process 0: Debug: check [at dolfin/mesh/LocalMeshData.cpp:191 in
> broadcast_mesh_data()]
> Process 0: Received 365 vertex indices
> Process 0: Debug: check [at dolfin/mesh/LocalMeshData.cpp:205 in
> broadcast_mesh_data()]
> Process 0: Sending 1536 cells to process 0, range is (0, 1536)
> Process 0: Sending 1536 cells to process 1, range is (1536, 3072)
> Process 1: Received 364 vertex indices
> Process 1: Debug: check [at dolfin/mesh/LocalMeshData.cpp:256 in
> receive_mesh_data()]
> Process 1: Received 1536 cell vertices
> Process 1: Debug: created local mesh data [at
> dolfin/mesh/MeshPartitioning.cpp:64 in partition()]
> Process 0: Received 1536 cell vertices
> Process 0: Debug: created local mesh data [at
> dolfin/mesh/MeshPartitioning.cpp:64 in partition()]
> Process 1: Partitioned mesh, edge cut is 207.
> Process 0: Partitioned mesh, edge cut is 207.
> Process 0: Calling FFC just-in-time (JIT) compiler, this may take some time.
> Traceback (most recent call last):
>   File "demo.py", line 25, in <module>
>     PN = FunctionSpace(mesh, "Nedelec 1st kind H(curl)", 0)
>   File "/opt/fenics/lib/python2.6/site-packages/dolfin/functionspace.py",
> line 190, in __init__
>     FunctionSpaceBase.__init__(self, mesh, element)
>   File "/opt/fenics/lib/python2.6/site-packages/dolfin/functionspace.py",
> line 43, in __init__
>     ufc_element, ufc_dofmap = jit(self._element)
>   File "/opt/fenics/lib/python2.6/site-packages/dolfin/jit.py", line
> 48, in mpi_jit
>     output = local_jit(*args, **kwargs)
>   File "/opt/fenics/lib/python2.6/site-packages/dolfin/jit.py", line 105, in jit
>     return jit_compile(form, options)
>   File "/opt/fenics/lib/python2.6/site-packages/ffc/jit/jit.py", line 56, in jit
>     return jit_element(object, options)
>   File "/opt/fenics/lib/python2.6/site-packages/ffc/jit/jit.py", line
> 122, in jit_element
>     (compiled_form, module, form_data) = jit_form(form, options)
>   File "/opt/fenics/lib/python2.6/site-packages/ffc/jit/jit.py", line
> 86, in jit_form
>     compile(form, signature, options)
>   File "/opt/fenics/lib/python2.6/site-packages/ffc/compiler/compiler.py",
> line 81, in compile
>     form_data = analyze_form(form.form_data(), options)
>   File "/opt/fenics/lib/python2.6/site-packages/ffc/compiler/compiler.py",
> line 118, in analyze_form
>     form_data.ffc_elements = [create_element(element) for element in
> form_data.elements]
>   File "/opt/fenics/lib/python2.6/site-packages/ffc/fem/createelement.py",
> line 48, in create_element
>     ffc_element = FFCFiniteElement(ufl_element.family(),
> ufl_element.cell().domain(), ufl_element.degree(), domain)
>   File "/opt/fenics/lib/python2.6/site-packages/ffc/fem/finiteelement.py",
> line 79, in __init__
>     (self.__fiat_element, self.__mapping) =
> self.__choose_element(family, shape, degree)
>   File "/opt/fenics/lib/python2.6/site-packages/ffc/fem/finiteelement.py",
> line 317, in __choose_element
>     return (Nedelec(fiat_shape, degree-1),
>   File "/opt/fenics/lib/python2.6/site-packages/FIAT/Nedelec.py", line
> 268, in __init__
>     U = NedelecSpace( shape , k )
>   File "/opt/fenics/lib/python2.6/site-packages/FIAT/Nedelec.py", line
> 92, in NedelecSpace
>     return NedelecSpace3D( degree )
>   File "/opt/fenics/lib/python2.6/site-packages/FIAT/Nedelec.py", line
> 37, in NedelecSpace3D
>     PkCrossX = polynomial.VectorPolynomialSet( Pkp1.base , PkCrossXcoeffs )
>   File "/opt/fenics/lib/python2.6/site-packages/FIAT/polynomial.py",
> line 242, in __init__
>     "Illegal coeff matrix: VectorPolynomialSet"
> RuntimeError: Illegal coeff matrix: VectorPolynomialSet
>
> Why does the demo hang and not just exit to the command line?
>
> Johannes
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev

Attachment: signature.asc
Description: Digital signature