Make curvilinear mesh work with domain parallelization
Description
Two routines were not working parallel in domains that make FD operators skew-adjoint or self-adjoint. They have been fixed by making them work only on the local part without the need to gather global arrays. This also makes it more efficient.
News snippet
Make curvilinear mesh work with domain parallelization
Checklist
-
I have checked that my code follows the Octopus coding standards -
I have added tests for all the new features added in this request.