View Issue Details

IDProjectCategoryView StatusLast Update
0002363OpenFOAM[All Projects] Bugpublic2016-12-01 20:00
ReporterThalenberghenAssigned Tohenry 
PrioritynormalSeveritycrashReproducibilityalways
Status resolvedResolutionfixed 
PlatformGNU/LinuxOSCentOS OS Version7.0
Product Version 
Fixed in Version 
Summary0002363: mapFields and mapFieldsPar crash with parallelSource
DescriptionWhen using mapFields with the flag -parallelSource it crashes with a segmentation fault. This can be reproduced using i.e. two instances of the motorbike tutorial and mapping them with -consistent.Using mapFieldsPar instead gives me the following error>

Creating mesh-to-mesh addressing for region0 and region0 regions using cellVolumeWeight
[3]
[3]
[3] --> FOAM FATAL ERROR:
[3] Plane normal defined with zero length
Bad points:(2 -0.999998 3) (2 -1.49999 3) (2 -1.25 3)
[3]
[3] From function void plane::calcPntAndVec
(
    const point&,
    const point&,
    const point&
)
[3]
[3] in file meshes/primitiveShapes/plane/plane.C at line 116.
[3]
FOAM parallel run aborting
[3]
[3] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #1 Foam::error::abort() at ??:?
[3] #2 Foam::plane::calcPntAndVec(Foam::Vector<float> const&, Foam::Vector<float> const&, Foam::Vector<float> const&) at ??:?
[3] #3 Foam::tetOverlapVolume::tetTetOverlapVol(Foam::tetPoints const&, Foam::tetPoints const&) const at ??:?
[3] #4 Foam::tetOverlapVolume::cellCellOverlapVolumeMinDecomp(Foam::primitiveMesh const&, long, Foam::primitiveMesh const&, long, Foam::treeBoundBox const&) const[0]
[0]
Steps To Reproduce1. Create two meshed and decomposed instances of the motorbike tutorial
2. Solve the first for some timeSteps
3. Try to map the results to the second instance
TagsNo tags attached.

Activities

henry

2016-12-01 12:28

manager   ~0007365

Have you tested OpenFOAM-4.1 or better still OpenFOAM-dev?

Thalenberghen

2016-12-01 14:31

reporter   ~0007367

I just tested OpenFOAM-4.1 with the same result.

MattijsJ

2016-12-01 17:45

reporter   ~0007371

In src/sampling/meshToMesh0/calculateMeshToMesh0Addressing.C it was not resizing the mapping for the faces (when reusing the mapping structure). Attached a fix.

calculateMeshToMesh0Addressing.C (10,938 bytes)

henry

2016-12-01 20:00

manager   ~0007374

Resolved in OpenFOAM-4.x by commit 0b3c8e6fb43f9127de22ce41ca6a51ee5f77193d

Resolved in OpenFOAM-dev by commit 2cff1e084debdd317afb8c70feae48d6583d5e32

Issue History

Date Modified Username Field Change
2016-12-01 12:11 Thalenberghen New Issue
2016-12-01 12:28 henry Note Added: 0007365
2016-12-01 14:31 Thalenberghen Note Added: 0007367
2016-12-01 17:45 MattijsJ File Added: calculateMeshToMesh0Addressing.C
2016-12-01 17:45 MattijsJ Note Added: 0007371
2016-12-01 20:00 henry Assigned To => henry
2016-12-01 20:00 henry Status new => resolved
2016-12-01 20:00 henry Resolution open => fixed
2016-12-01 20:00 henry Fixed in Version => 4.x
2016-12-01 20:00 henry Note Added: 0007374