[ripudaman]

Hi,
I am using foam-extend-3.2 on Mac OS 10.10

The plateHole case in solidMechanics/elasticSolidFoam gives a floating point error when I try to run the case. The mesh looks okay in paraview. However, when I run checkMesh on the mesh, it gives the same Floating Point Error.
These errors go away in debug mode and the case runs fine.
This also happens in the slidingFrictionBall case in the same elasticSolidFoam solver.
I pulled the latest version of the master branch and compiled it to test this out again but to no avail.
Please help.
Thank you.
Ripu

P.S. This has been posted as a bug on the foam-extend bug tracker on source forge

[ripudaman]

I dug deeper and found that the error occurs in the creation of the mesh object in createNamePolyMesh.H (lines 21-30)

    Foam::polyMesh mesh
    (
        Foam::IOobject
        (
            regionName,
            runTime.timeName(),
            runTime,
            Foam::IOobject::MUST_READ
        )
    );

This part of the code is the same as the code in fe31. Thus I believe the numerous changes made in polyMesh are the root of the problem. I compared the fe31 polyMesh and the fe32 polyMesh directories but I was very quickly overwhelmed. However, I noticed that the constructor declaration and definition both had no changes.

I need help with this.

Also, is there any other way of getting an error stack for this? Since the code runs fine in Debug mode... I do not know how to proceed and get the source of the error.

Thank you.
Ripu

[ripudaman]

So this actually seems like a problem inherent to the Mac OS. It occurs both in Yosemite (10.10) as well as El Capitan (10.11). The plateHole case does work in Win 7 and Ubuntu 14.04.

Here are my Mac etc/bashrc file changes:

# WM_MPLIB = SYSTEMOPENMPI | OPENMPI | MPICH | MPICH-GM | HPMPI | GAMMA | MPI | QSMPI | MVAPICH2
: ${WM_MPLIB:=MACPORTOPENMPI}; export WM_MPLIB
# : ${WM_MPLIB:=OPENMPI}; export WM_MPLIB

# Compiler: set to Gcc or Icc (for Intel's icc)
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
: ${WM_COMPILER:=Gcc46}; export WM_COMPILER
# : ${WM_COMPILER:=Gcc}; export WM_COMPILER
#: ${WM_COMPILER:=Icc}; export WM_COMPILER

Kindly point out any issues that you may feel are causing the strange behavior with the mesh.

[ripudaman]

I was able to solve the issue by changing the GCC version

GCC48 on pre ElCapitan and GCC49 on ElCapitan