Peakaboo kernel

[what-the-diff[bot]]

PR Summary

• Enhanced unit handling in blank_hit() and dimension_scale.py Improved unit management for smoother operation
• Implemented faster and flexible PeakFinder class in hitsearch() Replaced find_peaks() with a more efficient algorithm allowing user-controlled kernel size
• Introduced GPU acceleration for peak finding Created a new file utilizing CUDA via cupy library, increasing performance
• Resolved issues in DimensionScale index() method Fixed a bug related to handling of quantities
• Adjusted threshold in hitsearch test Changed the value from 15 to 30 to ensure proper testing
• Improved clean-up of output files in find_et tests Updated tests to remove output files, preventing interference with other tests

[telegraphic]

Speed report:

  _     ._   __/__   _ _  _  _ _/_   Recorded: 12:15:24  Samples:  29542
 /_//_/// /_\ / //_// / //_'/ //     Duration: 67.808    CPU time: 68.356
/   _/                      v4.1.1

Program: test_mars.py

67.807 <module>  <string>:1
   [7 frames hidden]  <string>, runpy, <built-in>
      67.806 _run_code  runpy.py:64
      └─ 67.806 <module>  test_mars.py:2
         ├─ 56.293 wrapper  hyperseti/utils.py:45
         │  └─ 56.291 find_et  hyperseti/pipeline.py:238
         │     ├─ 45.196 wrapper  hyperseti/utils.py:45
         │     │  └─ 45.196 run  hyperseti/pipeline.py:185
         │     │     ├─ 37.054 wrapper  hyperseti/utils.py:45
         │     │     │  ├─ 29.460 dedoppler  hyperseti/pipeline.py:146
         │     │     │  │  └─ 29.458 dedoppler  hyperseti/dedoppler.py:165
         │     │     │  │     ├─ 27.278 apply_boxcar_drift  hyperseti/dedoppler.py:55
         │     │     │  │     │  ├─ 20.372 asnumpy  cupy/__init__.py:756
         │     │     │  │     │  │     [6 frames hidden]  cupy, <built-in>, numpy
         │     │     │  │     │  │        20.363 ndarray.get  <built-in>:0
         │     │     │  │     │  ├─ 1.781 where  cupy/_sorting/search.py:177
         │     │     │  │     │  │     [22 frames hidden]  cupy, <built-in>, numpy
         │     │     │  │     │  ├─ 1.597 _array_repr_implementation  numpy/core/arrayprint.py:1426
         │     │     │  │     │  │     [110 frames hidden]  numpy, <built-in>, <__array_function_...
         │     │     │  │     │  ├─ 1.536 uniform_filter1d  cupyx/scipy/ndimage/filters.py:210
         │     │     │  │     │  │     [126 frames hidden]  cupyx, cupy, <built-in>, genericpath,...
         │     │     │  │     │  └─ 1.470 [self]
         │     │     │  │     └─ 1.145 _array_str_implementation  numpy/core/arrayprint.py:1534
         │     │     │  │           [60 frames hidden]  numpy, <built-in>, <__array_function_...
         │     │     │  ├─ 6.411 preprocess  hyperseti/pipeline.py:105
         │     │     │  │  ├─ 4.856 normalize  hyperseti/normalize.py:13
         │     │     │  │  │  └─ 3.981 polyfit  cupy/lib/_routines_poly.py:123
         │     │     │  │  │        [143 frames hidden]  cupy, <built-in>, genericpath, posixp...
         │     │     │  │  └─ 1.004 sk_flag  hyperseti/kurtosis.py:33
         │     │     │  └─ 1.181 hitsearch  hyperseti/pipeline.py:159
         │     │     │     └─ 1.121 hitsearch  hyperseti/hits.py:95
         │     │     ├─ 6.702 merge_hits  hyperseti/hits.py:55
         │     │     │  ├─ 5.066 query  pandas/core/frame.py:3910
         │     │     │  │     [3047 frames hidden]  pandas, abc, <built-in>, <__array_fun...
         │     │     │  └─ 0.704 wrapper  pandas/util/_decorators.py:302
         │     │     │        [603 frames hidden]  pandas, <built-in>, numpy, <__array_f...
         │     │     └─ 1.419 blank_hits  hyperseti/blanking.py:86
         │     │        └─ 1.387 blank_hit  hyperseti/blanking.py:48
         │     │           └─ 1.269 [self]
         │     ├─ 9.179 iterate_through_data  hyperseti/data_array.py:226
         │     │  └─ 9.178 isel  hyperseti/data_array.py:194
         │     │     ├─ 6.491 __getitem__  h5py/_hl/dataset.py:476
         │     │     │     [31 frames hidden]  h5py, <__array_function__ internals>,...
         │     │     └─ 2.351 asarray  cupy/_creation/from_data.py:44
         │     │           [4 frames hidden]  cupy, <built-in>
         │     └─ 1.773 __init__  hyperseti/pipeline.py:62
         └─ 11.451 <module>  hyperseti/__init__.py:1
            └─ 11.347 <module>  hyperseti/pipeline.py:1
               ├─ 7.202 <module>  setigen/__init__.py:1
               │     [7179 frames hidden]  setigen, <__array_function__ internal...
               ├─ 1.048 <module>  cupy/__init__.py:1
               │     [1663 frames hidden]  cupy, pytest, _pytest, attr, <built-i...
               ├─ 0.807 <module>  hyperseti/hits.py:1
               │  └─ 0.804 <module>  hyperseti/peak.py:1
               │     └─ 0.803 <module>  cusignal/__init__.py:14
               │           [1043 frames hidden]  cusignal, numba, abc, <built-in>, col...
               ├─ 0.771 <module>  hyperseti/dedoppler.py:1
               │  └─ 0.738 <module>  hyperseti/data_array.py:1
               └─ 0.722 <module>  pandas/__init__.py:3
                     [1083 frames hidden]  pandas, pyarrow, textwrap, <built-in>...