Algorithm 932: {PANG}: Software for Non-Matching Grid Projections in 2d and 3d with Linear Complexity

International audienceWe design and analyze an algorithm with linear complexity to perform projections between 2d and 3d non-matching grids. This algorithm, named the PANG algorithm, is based on an advancing front technique and neighboring information. Its implementation is surprisingly short, and we give the entire Matlab code. For computing the intersections, we use a direct and numerically robust approach. We show numerical experiments both for 2d and 3d grids, which illustrate the optimal complexity and negligible overhead of the algorithm. An outline of this algorithm has been presented already in a short proceedings paper of the 18th international conference on domain decomposition methods, (see Gander and Japhet [2008])