-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 On 05/02/15 09:44, Homolya, Miklós wrote:
Hello Eike,
It fails exactly before the Mesh factory function would identify whether to create a simplex or quadrilateral mesh object. It tries to retrieve how many facets a cell has, which fails the assertion pStart <= cStart < pEnd. The only explanation I have is that some MPI processes own 0 (zero) cells of the mesh.
Good crystal-balling:
Maybe try using less MPI nodes or a higher refinement (more cells). ...
Running on 24 MPI processes
...
ref_count_coarse = 0
This coarse mesh has 20 cells, so four processes (at least) will not own any cells on the coarse mesh. While most of firedrake and pyop2 will work in this situation, various bits of mesh setup won't. Note that even if the coarsest mesh went through, when you refine 4 of the processes won't have any work. Lawrence -----BEGIN PGP SIGNATURE----- Version: GnuPG v1 iQEcBAEBAgAGBQJU0z3JAAoJECOc1kQ8PEYvW1oH/31a0AY/r9O2Axu0xgsYNi9C 28O4e/VrHl4jBhqa1Og+a7if6Puole7/5+aP97P5tszMvA4FRUA/0BudUnGmLpPV +vg+cf6RF98R9TepTDoBnbIk78q2iWcHeDNKpJFFK1lMWiy2+Bg3yFm29LqsIZ4V y+vN1rYusIiyALD0san7bY//OZGyfdORTxGqmJOl2uKWkguzcD9ubolowNvRLWTN Ipavo7ayfGitIcbHSe2g/yZNJFGCHZMAHG6zVm+fuzTNRbDcdF8QIHU/FAgaIZlp v+O1iPRONZicwyS0RYrOwJDqdtsrSW8qcVqAJPrREmq2TPErsaAuFkAZ5ME38js= =o7zR -----END PGP SIGNATURE-----