public void PointSearchTest()
{
ElementSet elmtSet = ReadMesh(@"..\..\Spatial\Data\oresund.mesh");
XYElementSearchTree <int> tree = XYElementSearchTree <int> .BuildSearchTree(elmtSet);
// Point in element 746 (has index 745)
XYPoint point = new XYPoint(355273.227764, 6201136.0892);
XYExtent pointExtent = XYExtentUtil.GetExtent(point);
ICollection <int> pointCandidates = tree.FindElements(pointExtent);
Assert.Contains(745, (ICollection)pointCandidates);
XYPolygon elmt745 = ElementMapper.CreateFromXYPolygon(elmtSet, 745);
Assert.IsTrue(XYGeometryTools.IsPointInPolygon(point.X, point.Y, elmt745));
// Point in element 3583 (has index 3582)
XYPoint point2 = new XYPoint(354683.522377, 6167696.720773);
XYExtent point2Extent = XYExtentUtil.GetExtent(point2);
ICollection <int> point2Candidates = tree.FindElements(point2Extent);
Assert.Contains(3582, (ICollection)point2Candidates);
XYPolygon elmt3582 = ElementMapper.CreateFromXYPolygon(elmtSet, 3582);
Assert.IsTrue(XYGeometryTools.IsPointInPolygon(point2.X, point2.Y, elmt3582));
// Node 2001 is part of 7 elements
XYPoint pointNode2001 = new XYPoint(355451.29058143805, 6167084.759883712);
XYExtent pointNode2001Extent = XYExtentUtil.GetExtent(pointNode2001);
ICollection <int> pointNode2001Candidates = tree.FindElements(pointNode2001Extent);
Assert.AreEqual(7, pointNode2001Candidates.Count);
}
