Hello, When I access facets and cells in parallel, their indices are not global but the one in each process. How can I access the global index? The code I run is: from dolfin import * parameters['ghost_mode'] = 'shared_facets' N = 2 mesh = UnitSquareMesh(N,N) for c in cells(mesh): for f in…

How about: from dolfin import * parameters['ghost_mode'] = 'shared_facets' N = 2 mesh = UnitSquareMesh(N, N) tdim = mesh.topology().dim() mesh.init(tdim, tdim - 1) c_to_f = mesh.topology()(tdim, tdim-1) for c in cells(mesh): facets = c_to_f(c.index()) if MPI.comm_world.size != 1: …

Indices of cells and facets in parallel

mesh

Yingqi_Jia July 10, 2021, 6:08am 2

Maybe this post will be helpful, although it is given in dolfinx. Gather solutions in parallel in FEniCSX - #2 by dokken

Topic		Replies	Views
MPI demo in dolfinx	2	635	March 9, 2023
How to define the cell index in Parallel General	1	219	November 1, 2022
Questions about local dofs and ordering General	8	774	February 15, 2023
Problem regarding reading mesh in parallel mesh	6	1201	February 12, 2020
Parallel code for extraction of submatrices General	0	111	November 16, 2023

Indices of cells and facets in parallel

Related Topics