Coverage for src/meshpy/geometric_search/cython.py: 82%

1# The MIT License (MIT) 

2# 

3# Copyright (c) 2018-2025 MeshPy Authors 

4# 

5# Permission is hereby granted, free of charge, to any person obtaining a copy 

6# of this software and associated documentation files (the "Software"), to deal 

7# in the Software without restriction, including without limitation the rights 

8# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell 

9# copies of the Software, and to permit persons to whom the Software is 

10# furnished to do so, subject to the following conditions: 

11# 

12# The above copyright notice and this permission notice shall be included in 

13# all copies or substantial portions of the Software. 

14# 

15# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 

16# IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 

17# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 

18# AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 

19# LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, 

20# OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN 

21# THE SOFTWARE. 

22"""This file defines the interface to the Cython geometric search 

23functionality.""" 

24 

25import warnings as _warnings 

26 

27from meshpy.geometric_search.utils import cython_is_available as _cython_is_available 

28 

29if _cython_is_available(): 

30 from meshpy.geometric_search.cython_lib import ( 

31 find_close_points as _find_close_points, 

32 ) 

33 

34 

35def find_close_points_brute_force_cython( 

36 point_coordinates, tol, *, n_points_performance_warning=5000 

37): 

38 """Call the Cython brute force implementation of find close_points.""" 

39 if _cython_is_available(): 

40 n_points = len(point_coordinates) 

41 if n_points > n_points_performance_warning: 

42 _warnings.warn( 

43 "The function find_close_points is called with the brute force algorithm " 

44 + f"with {n_points} points, for performance reasons other algorithms should be used!" 

45 ) 

46 return _find_close_points(point_coordinates, tol) 

47 else: 

48 raise ModuleNotFoundError( 

49 "Cython geometric search functionality is not available" 

50 )