double distanceSolidSolid3D( const Solid& gA, const Solid& gB ) { //SFCGAL_DEBUG( boost::format("dispatch distancePolygonGeometry3D(%s,%s)") % gA.asText() % gB.asText() ); if ( gA.isEmpty() || gB.isEmpty() ) { return std::numeric_limits< double >::infinity() ; } if ( intersects( gA, gB, NoValidityCheck() ) ) { return 0.0 ; } double dMin = std::numeric_limits< double >::infinity() ; for ( size_t i = 0; i < gA.numShells(); i++ ) { for ( size_t j = 0; j < gB.numShells(); j++ ) { dMin = std::min( dMin, gA.shellN( i ).distance3D( gB.shellN( j ) ) ); } } return dMin ; }
double distanceTriangleSolid3D( const Triangle& gA, const Solid& gB ) { if ( gA.isEmpty() || gB.isEmpty() ) { return std::numeric_limits< double >::infinity() ; } if ( intersects3D( gA, gB, NoValidityCheck() ) ) { return 0.0 ; } double dMin = std::numeric_limits< double >::infinity() ; for ( size_t i = 0; i < gB.numShells(); i++ ) { dMin = std::min( dMin, gB.shellN( i ).distance3D( gA ) ); } return dMin ; }
double distancePointSolid3D( const Point& gA, const Solid& gB ) { if ( gA.isEmpty() || gB.isEmpty() ) { return std::numeric_limits< double >::infinity() ; } if ( intersects3D( gA, gB, NoValidityCheck() ) ) { return 0.0 ; } double dMin = std::numeric_limits< double >::infinity() ; for ( size_t i = 0; i < gB.numShells(); i++ ) { dMin = std::min( dMin, distanceGeometryCollectionToGeometry3D( gB.shellN( i ), gA ) ); } return dMin ; }