exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Adjust round-off tolerances for optimising compilers

Browse source
This commit is contained in:
Hrvoje Jasak 2011-10-12 18:16:42 +01:00
parent c906570c2d
commit 388a75ed4f
1 changed files with 8 additions and 14 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

22
src/OpenFOAM/algorithms/octree/octree/treeBoundBox.C
View file

@ -333,8 +333,8 @@ bool Foam::treeBoundBox::intersects
//HJ, experimental: problems with round-off on Intel 12 compiler.
// HJ, 29/Sep/2011
const scalar smallNumber = SMALL;
const scalar kSmall = 1000*SMALL;
const scalar smallNumber = VSMALL;
const scalar kSmall = 1000*smallNumber;
while (true)
{
@ -362,8 +362,7 @@ bool Foam::treeBoundBox::intersects
if (Foam::mag(overallVec.x()) > kSmall)
{
s = (min().x() - overallStart.x())/
overallVec.x();
// stabilise(overallVec.x(), smallNumber);
stabilise(overallVec.x(), smallNumber);
pt.x() = min().x();
pt.y() = overallStart.y() + overallVec.y()*s;
@ -382,8 +381,7 @@ bool Foam::treeBoundBox::intersects
if (Foam::mag(overallVec.x()) > kSmall)
{
s = (max().x() - overallStart.x())/
overallVec.x();
// stabilise(overallVec.x(), smallNumber);
stabilise(overallVec.x(), smallNumber);
pt.x() = max().x();
pt.y() = overallStart.y() + overallVec.y()*s;
@ -400,8 +398,7 @@ bool Foam::treeBoundBox::intersects
if (Foam::mag(overallVec.y()) > kSmall)
{
s = (min().y() - overallStart.y())/
overallVec.y();
// stabilise(overallVec.y(), smallNumber);
stabilise(overallVec.y(), smallNumber);
pt.x() = overallStart.x() + overallVec.x()*s;
pt.y() = min().y();
@ -418,8 +415,7 @@ bool Foam::treeBoundBox::intersects
if (Foam::mag(overallVec.y()) > kSmall)
{
s = (max().y() - overallStart.y())/
overallVec.y();
// stabilise(overallVec.y(), smallNumber);
stabilise(overallVec.y(), smallNumber);
pt.x() = overallStart.x() + overallVec.x()*s;
pt.y() = max().y();
@ -436,8 +432,7 @@ bool Foam::treeBoundBox::intersects
if (Foam::mag(overallVec.z()) > kSmall)
{
s = (min().z() - overallStart.z())/
overallVec.z();
// stabilise(overallVec.z(), smallNumber);
stabilise(overallVec.z(), smallNumber);
pt.x() = overallStart.x() + overallVec.x()*s;
pt.y() = overallStart.y() + overallVec.y()*s;
@ -454,8 +449,7 @@ bool Foam::treeBoundBox::intersects
if (Foam::mag(overallVec.z()) > kSmall)
{
s = (max().z() - overallStart.z())/
overallVec.z();
// stabilise(overallVec.z(), smallNumber);
stabilise(overallVec.z(), smallNumber);
pt.x() = overallStart.x() + overallVec.x()*s;
pt.y() = overallStart.y() + overallVec.y()*s;