checkGeometry.C

                            globalFaces,

                            mergedFaces,
                            mergedPoints
                        );
                        // Collect field
                        scalarField mergedWeights;
                        globalFaces().gather
                        (
                            UPstream::worldComm,
                            ListOps::create<label>
                            (
                                UPstream::procID(UPstream::worldComm),
                                labelOp<int>()  // int -> label
                            ),
                            ami.srcWeightsSum(),
                            mergedWeights
                        );

                        if (Pstream::master())
                        {
                            wr.write
                            (
                                functionObject::outputPrefix,
                                "src_" + tmName,
                                meshedSurfRef
                                (
                                    mergedPoints,
                                    mergedFaces
                                ),
                                "weightsSum",
                                mergedWeights,
                                false
                            );
                        }
                    }
                    {
                        // Collect geometry
                        labelList pointToGlobal;
                        labelList uniqueMeshPointLabels;
                        autoPtr<globalIndex> globalPoints;
                        autoPtr<globalIndex> globalFaces;
                        faceList mergedFaces;
                        pointField mergedPoints;
                        Foam::PatchTools::gatherAndMerge
                        (
                            mesh,
                            cpp.neighbPatch().localFaces(),
                            cpp.neighbPatch().meshPoints(),
                            cpp.neighbPatch().meshPointMap(),

                            pointToGlobal,
                            uniqueMeshPointLabels,
                            globalPoints,
                            globalFaces,

                            mergedFaces,
                            mergedPoints
                        );
                        // Collect field
                        scalarField mergedWeights;
                        globalFaces().gather
                        (
                            UPstream::worldComm,
                            ListOps::create<label>
                            (
                                UPstream::procID(UPstream::worldComm),
                                labelOp<int>()  // int -> label
                            ),
                            ami.tgtWeightsSum(),
                            mergedWeights
                        );

                        if (Pstream::master())
                        {
                            wr.write
                            (
                                functionObject::outputPrefix,
                                "tgt_" + tmName,
                                meshedSurfRef
                                (
                                    mergedPoints,
                                    mergedFaces
                                ),
                                "weightsSum",
                                mergedWeights,
                                false
                            );
                        }
                    }
                }
            }
        }
    }

    return noFailedChecks;
}