CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

plateHole in solidMechanics/elasticSolidFoam fails checkMesh in Opt but runs in Debug

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 10, 2015, 01:54
Default plateHole in solidMechanics/elasticSolidFoam fails checkMesh in Opt but runs in Debug
  #1
Member
 
Ripudaman Manchanda
Join Date: May 2013
Posts: 50
Rep Power: 5
ripudaman is on a distinguished road
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 is offline   Reply With Quote

Old   November 10, 2015, 12:46
Default more information
  #2
Member
 
Ripudaman Manchanda
Join Date: May 2013
Posts: 50
Rep Power: 5
ripudaman is on a distinguished road
I dug deeper and found that the error occurs in the creation of the mesh object in createNamePolyMesh.H (lines 21-30)
Code:
    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 is offline   Reply With Quote

Old   November 12, 2015, 13:26
Default Works in Linux and Windows
  #3
Member
 
Ripudaman Manchanda
Join Date: May 2013
Posts: 50
Rep Power: 5
ripudaman is on a distinguished road
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:
Code:
# WM_MPLIB = SYSTEMOPENMPI | OPENMPI | MPICH | MPICH-GM | HPMPI | GAMMA | MPI | QSMPI | MVAPICH2
: ${WM_MPLIB:=MACPORTOPENMPI}; export WM_MPLIB
# : ${WM_MPLIB:=OPENMPI}; export WM_MPLIB
Code:
# 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 is offline   Reply With Quote

Old   December 1, 2015, 19:08
Default Change the GCC version to fix this....
  #4
Member
 
Ripudaman Manchanda
Join Date: May 2013
Posts: 50
Rep Power: 5
ripudaman is on a distinguished road
I was able to solve the issue by changing the GCC version

GCC48 on pre ElCapitan and GCC49 on ElCapitan
ripudaman is offline   Reply With Quote

Reply

Tags
checkmesh, foam-extend-3.2, platehole

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem in Implementing a non-linear SGS model huangxianbei OpenFOAM Programming & Development 30 September 20, 2015 13:26
could not open file .vtk in paraview ali_atrian OpenFOAM Paraview & paraFoam 8 August 27, 2014 10:31
CheckMesh fails when cyclic boundaries are defined - OF 2.2.0 caduqued OpenFOAM Meshing & Mesh Conversion 2 May 21, 2013 18:31
GGI and Topological changes - checkMesh fails strakakl OpenFOAM Programming & Development 3 June 20, 2012 05:30
Compiling checkMesh fails kjetil OpenFOAM Installation 1 February 8, 2010 12:05


All times are GMT -4. The time now is 01:34.