Merge branch 'master' of /home/dm4/OpenFOAM/OpenFOAM-dev

This commit is contained in:
andy
2013-08-29 16:38:34 +01:00

View File

@ -140,7 +140,7 @@ List<Vb::Point> rayShooting::initialPoints() const
geometryToConformTo().findSurfaceNearestIntersection geometryToConformTo().findSurfaceNearestIntersection
( (
fC - normStart[0]*SMALL, fC - normStart[0]*pert,
fC - normStart[0]*maxRayLength, fC - normStart[0]*maxRayLength,
surfHitEnd, surfHitEnd,
hitSurfaceEnd hitSurfaceEnd
@ -166,8 +166,8 @@ List<Vb::Point> rayShooting::initialPoints() const
pointIndexHit procIntersection = pointIndexHit procIntersection =
decomposition().findLine decomposition().findLine
( (
l.start() + l.vec()*SMALL, l.start(),
l.end() - l.vec()*maxRayLength l.end()
); );
if (procIntersection.hit()) if (procIntersection.hit())
@ -187,12 +187,6 @@ List<Vb::Point> rayShooting::initialPoints() const
minimumSurfaceDistanceCoeffSqr_ minimumSurfaceDistanceCoeffSqr_
*sqr(cellShapeControls().cellSize(midPoint)); *sqr(cellShapeControls().cellSize(midPoint));
if
(
magSqr(midPoint - l.start()) > minDistFromSurfaceSqr
&& magSqr(midPoint - l.end()) > minDistFromSurfaceSqr
)
{
if (randomiseInitialGrid_) if (randomiseInitialGrid_)
{ {
midPoint.x() += pert*(rndGen().scalar01() - 0.5); midPoint.x() += pert*(rndGen().scalar01() - 0.5);
@ -200,6 +194,12 @@ List<Vb::Point> rayShooting::initialPoints() const
midPoint.z() += pert*(rndGen().scalar01() - 0.5); midPoint.z() += pert*(rndGen().scalar01() - 0.5);
} }
if
(
magSqr(midPoint - l.start()) > minDistFromSurfaceSqr
&& magSqr(midPoint - l.end()) > minDistFromSurfaceSqr
)
{
initialPoints.append(toPoint(midPoint)); initialPoints.append(toPoint(midPoint));
} }
} }