Point evaluation in parallel

[cwangaw]

Hi all,

I am running my code in parallel using mpiexec. However, I keep encountering the following error: ValueError: Points to evaluate are inconsistent among processes.

I have tried the following approaches to resolve the issue:

1. Using If COMM_WORLD.rank == 0:

However, when I include this condition, the program always freezes upon execution, regardless of the subsequent code.

2. Adjusting the tolerance as suggested in the documentation

Here is the corresponding part of my code:

def polylog(s,u):
   '''
      usage: given a function u(x), it returns a function f = polylog(s,u(x)) 
             defined in the same finite element space
   '''
   V = u.function_space()
   m = V.mesh()
   W = VectorFunctionSpace(m, V.ufl_element())
   X = assemble(interpolate(m.coordinates, W))
   a = X.dat.data_ro
   m.tolerance = 0
   eval = [float(np.real(mpmath.polylog(s,float(val)))) for val in u.at(a)]
   m.tolerance = 0.1
   f = Function(V)
   f.dat.data[:] = eval

   return f

However, it does not work.

I would appreciate any insights on how to fix this error or better handle the parallel execution in this context.

Thanks in advance!