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
Cell indices in subdomain case when running parallelly General	3	374	October 11, 2022
MPI demo in dolfinx	2	758	March 9, 2023
What is the more succint way to get the "cell" and "facet" tags? mesh	5	300	October 27, 2023
How to define the cell index in Parallel General	1	256	November 1, 2022
Save MeshFunction in parallel with parameters['ghost_mode'] = 'shared_facet' I/O mesh , io , mpi	5	15	August 16, 2024

Indices of cells and facets in parallel

Related Topics