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
Local cell and global edges mesh	2	337	December 2, 2020
Cell indices in subdomain case when running parallelly General	3	377	October 11, 2022
MPI demo in dolfinx	2	797	March 9, 2023
What is the more succint way to get the "cell" and "facet" tags? mesh	5	335	October 27, 2023
How to define the cell index in Parallel General	1	256	November 1, 2022

Indices of cells and facets in parallel

Related topics