BUG: turbulentDFSEMInlet; Option = planarInterpolation; decomposePar causes Floating Point Exception for Fields R/L/U

Summary

When turbulentDFSEMInlet with planarInterpolation option is ON, decomposePar causes FPE with the following error message:

Time = 0
Turbulent DFSEM patch inlet: interpolating field R from "/home/snoopy2/kuta/OpenFOAM/kuta-v1806/run/channel395DFSEM/constant/boundaryData/inlet/0/R"
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib64/libc.so.6
#3  Foam::pointToPointPlanarInterpolation::calcCoordinateSystem(Foam::Field<Foam::Vector<double> > const&) const at ??:?
#4  Foam::pointToPointPlanarInterpolation::pointToPointPlanarInterpolation(Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, double, bool) at ??:?
#5  Foam::turbulentDFSEMInletFvPatchVectorField::patchMapper() const at ??:?
#6  Foam::tmp<Foam::Field<Foam::SymmTensor<double> > > Foam::turbulentDFSEMInletFvPatchVectorField::interpolateBoundaryData<Foam::SymmTensor<double> >(Foam::word const&) const at ??:?
#7  Foam::turbulentDFSEMInletFvPatchVectorField::turbulentDFSEMInletFvPatchVectorField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#8  Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::turbulentDFSEMInletFvPatchVectorField>::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#9  Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#10  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#11  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) at ??:?
#12  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields() at ??:?
#13  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&, bool) at ??:?
#14  ? at ??:?
#15  ? at ??:?
#16  __libc_start_main in /lib64/libc.so.6
#17  ? at /home/abuild/rpmbuild/BUILD/glibc-2.18/csu/../sysdeps/x86_64/start.S:125
Floating point exception

Steps to reproduce

Please consider tutorials/incompressible/pimpleFoam/LES/channel395DFSEM case.

Change the following under 0/U:

    inlet
    {   
        type            turbulentDFSEMInlet;
        ...
        mapMethod       nearestCell;
        ...
    }   

to

    inlet
    {   
        type            turbulentDFSEMInlet;
        ...
        mapMethod       planarInterpolation;
        ...
    }   

Then execute the following respectively:

blockMesh

decomposePar

PS: scotch and hierarchial were tested with various numbers of numberOfSubdomains.

@andy