scan.hpp 3.9 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133
  1. /*M///////////////////////////////////////////////////////////////////////////////////////
  2. //
  3. // IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
  4. //
  5. // By downloading, copying, installing or using the software you agree to this license.
  6. // If you do not agree to this license, do not download, install,
  7. // copy or use the software.
  8. //
  9. //
  10. // License Agreement
  11. // For Open Source Computer Vision Library
  12. //
  13. // Copyright (C) 2000-2008, Intel Corporation, all rights reserved.
  14. // Copyright (C) 2009, Willow Garage Inc., all rights reserved.
  15. // Copyright (C) 2013, OpenCV Foundation, all rights reserved.
  16. // Third party copyrights are property of their respective owners.
  17. //
  18. // Redistribution and use in source and binary forms, with or without modification,
  19. // are permitted provided that the following conditions are met:
  20. //
  21. // * Redistribution's of source code must retain the above copyright notice,
  22. // this list of conditions and the following disclaimer.
  23. //
  24. // * Redistribution's in binary form must reproduce the above copyright notice,
  25. // this list of conditions and the following disclaimer in the documentation
  26. // and/or other materials provided with the distribution.
  27. //
  28. // * The name of the copyright holders may not be used to endorse or promote products
  29. // derived from this software without specific prior written permission.
  30. //
  31. // This software is provided by the copyright holders and contributors "as is" and
  32. // any express or implied warranties, including, but not limited to, the implied
  33. // warranties of merchantability and fitness for a particular purpose are disclaimed.
  34. // In no event shall the Intel Corporation or contributors be liable for any direct,
  35. // indirect, incidental, special, exemplary, or consequential damages
  36. // (including, but not limited to, procurement of substitute goods or services;
  37. // loss of use, data, or profits; or business interruption) however caused
  38. // and on any theory of liability, whether in contract, strict liability,
  39. // or tort (including negligence or otherwise) arising in any way out of
  40. // the use of this software, even if advised of the possibility of such damage.
  41. //
  42. //M*/
  43. #pragma once
  44. #ifndef OPENCV_CUDEV_WARP_SCAN_HPP
  45. #define OPENCV_CUDEV_WARP_SCAN_HPP
  46. #include "../common.hpp"
  47. #include "warp.hpp"
  48. #include "shuffle.hpp"
  49. namespace cv { namespace cudev {
  50. //! @addtogroup cudev
  51. //! @{
  52. #if __CUDACC_VER_MAJOR__ >= 9
  53. // Starting from CUDA 9.0, support for Fermi is dropped.
  54. // So CV_CUDEV_ARCH >= 300 is implied.
  55. template <typename T>
  56. __device__ T warpScanInclusive(uint mask, T data)
  57. {
  58. const uint laneId = Warp::laneId();
  59. // scan on shufl functions
  60. #pragma unroll
  61. for (int i = 1; i <= (WARP_SIZE / 2); i *= 2)
  62. {
  63. const T val = shfl_up_sync(mask, data, i);
  64. if (laneId >= i)
  65. data += val;
  66. }
  67. return data;
  68. }
  69. template <typename T>
  70. __device__ __forceinline__ T warpScanExclusive(uint mask, T data)
  71. {
  72. return warpScanInclusive(mask, data) - data;
  73. }
  74. #else // __CUDACC_VER_MAJOR__ >= 9
  75. template <typename T>
  76. __device__ T warpScanInclusive(T data, volatile T* smem, uint tid)
  77. {
  78. #if CV_CUDEV_ARCH >= 300
  79. CV_UNUSED(smem);
  80. CV_UNUSED(tid);
  81. const uint laneId = Warp::laneId();
  82. // scan on shufl functions
  83. #pragma unroll
  84. for (int i = 1; i <= (WARP_SIZE / 2); i *= 2)
  85. {
  86. const T val = shfl_up(data, i);
  87. if (laneId >= i)
  88. data += val;
  89. }
  90. return data;
  91. #else
  92. const uint laneId = Warp::laneId();
  93. smem[tid] = data;
  94. #pragma unroll
  95. for (int i = 1; i <= (WARP_SIZE / 2); i *= 2)
  96. if (laneId >= i)
  97. smem[tid] += smem[tid - i];
  98. return smem[tid];
  99. #endif
  100. }
  101. template <typename T>
  102. __device__ __forceinline__ T warpScanExclusive(T data, volatile T* smem, uint tid)
  103. {
  104. return warpScanInclusive(data, smem, tid) - data;
  105. }
  106. #endif // __CUDACC_VER_MAJOR__ >= 9
  107. //! @}
  108. }}
  109. #endif