/*M///////////////////////////////////////////////////////////////////////////////////////
//
//  IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
//
//  By downloading, copying, installing or using the software you agree to this license.
//  If you do not agree to this license, do not download, install,
//  copy or use the software.
//
//
//                           License Agreement
//                For Open Source Computer Vision Library
//
// Copyright (C) 2000-2008, Intel Corporation, all rights reserved.
// Copyright (C) 2009, Willow Garage Inc., all rights reserved.
// Third party copyrights are property of their respective owners.
//
// Redistribution and use in source and binary forms, with or without modification,
// are permitted provided that the following conditions are met:
//
//   * Redistribution's of source code must retain the above copyright notice,
//     this list of conditions and the following disclaimer.
//
//   * Redistribution's in binary form must reproduce the above copyright notice,
//     this list of conditions and the following disclaimer in the documentation
//     and/or other materials provided with the distribution.
//
//   * The name of the copyright holders may not be used to endorse or promote products
//     derived from this software without specific prior written permission.
//
// This software is provided by the copyright holders and contributors "as is" and
// any express or implied warranties, including, but not limited to, the implied
// warranties of merchantability and fitness for a particular purpose are disclaimed.
// In no event shall the Intel Corporation or contributors be liable for any direct,
// indirect, incidental, special, exemplary, or consequential damages
// (including, but not limited to, procurement of substitute goods or services;
// loss of use, data, or profits; or business interruption) however caused
// and on any theory of liability, whether in contract, strict liability,
// or tort (including negligence or otherwise) arising in any way out of
// the use of this software, even if advised of the possibility of such damage.
//
//M*/

#ifndef __OPENCV_GPU_REDUCTION_DETAIL_HPP__
#define __OPENCV_GPU_REDUCTION_DETAIL_HPP__

namespace cv { namespace gpu { namespace device
{
    namespace utility_detail
    {
        ///////////////////////////////////////////////////////////////////////////////
        // Reductor

        template <int n> struct WarpReductor
        {
            template <typename T, typename Op> static __device__ __forceinline__ void reduce(volatile T* data, T& partial_reduction, int tid, const Op& op)
            {
                if (tid < n)
                    data[tid] = partial_reduction;
                if (n > 32) __syncthreads();

                if (n > 32)
                {
                    if (tid < n - 32)
                        data[tid] = partial_reduction = op(partial_reduction, data[tid + 32]);
                    if (tid < 16)
                    {
                        data[tid] = partial_reduction = op(partial_reduction, data[tid + 16]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  8]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  4]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  2]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  1]);
                    }
                }
                else if (n > 16)
                {
                    if (tid < n - 16)
                        data[tid] = partial_reduction = op(partial_reduction, data[tid + 16]);
                    if (tid < 8)
                    {
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  8]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  4]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  2]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  1]);
                    }
                }
                else if (n > 8)
                {
                    if (tid < n - 8)
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  8]);
                    if (tid < 4)
                    {
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  4]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  2]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  1]);
                    }
                }
                else if (n > 4)
                {
                    if (tid < n - 4)
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  4]);
                    if (tid < 2)
                    {
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  2]);
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  1]);
                    }
                }
                else if (n > 2)
                {
                    if (tid < n - 2)
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  2]);
                    if (tid < 2)
                    {
                        data[tid] = partial_reduction = op(partial_reduction, data[tid +  1]);
                    }
                }
            }
        };
        template <> struct WarpReductor<64>
        {
            template <typename T, typename Op> static __device__ void reduce(volatile T* data, T& partial_reduction, int tid, const Op& op)
            {
                data[tid] = partial_reduction;
                __syncthreads();

                if (tid < 32)
                {
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 32]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 16]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 8 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 4 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 2 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 1 ]);
                }
            }
        };
        template <> struct WarpReductor<32>
        {
            template <typename T, typename Op> static __device__ void reduce(volatile T* data, T& partial_reduction, int tid, const Op& op)
            {
                data[tid] = partial_reduction;

                if (tid < 16)
                {
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 16]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 8 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 4 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 2 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 1 ]);
                }
            }
        };
        template <> struct WarpReductor<16>
        {
            template <typename T, typename Op> static __device__ void reduce(volatile T* data, T& partial_reduction, int tid, const Op& op)
            {
                data[tid] = partial_reduction;

                if (tid < 8)
                {
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 8 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 4 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 2 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 1 ]);
                }
            }
        };
        template <> struct WarpReductor<8>
        {
            template <typename T, typename Op> static __device__ void reduce(volatile T* data, T& partial_reduction, int tid, const Op& op)
            {
                data[tid] = partial_reduction;

                if (tid < 4)
                {
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 4 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 2 ]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 1 ]);
                }
            }
        };

        template <bool warp> struct ReductionDispatcher;
        template <> struct ReductionDispatcher<true>
        {
            template <int n, typename T, typename Op> static __device__ void reduce(volatile T* data, T& partial_reduction, int tid, const Op& op)
            {
                WarpReductor<n>::reduce(data, partial_reduction, tid, op);
            }
        };
        template <> struct ReductionDispatcher<false>
        {
            template <int n, typename T, typename Op> static __device__ void reduce(volatile T* data, T& partial_reduction, int tid, const Op& op)
            {
                if (tid < n)
                    data[tid] = partial_reduction;
                __syncthreads();


                if (n == 512) { if (tid < 256) { data[tid] = partial_reduction = op(partial_reduction, data[tid + 256]); } __syncthreads(); }
                if (n >= 256) { if (tid < 128) { data[tid] = partial_reduction = op(partial_reduction, data[tid + 128]); } __syncthreads(); }
                if (n >= 128) { if (tid <  64) { data[tid] = partial_reduction = op(partial_reduction, data[tid +  64]); } __syncthreads(); }

                if (tid < 32)
                {
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 32]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid + 16]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid +  8]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid +  4]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid +  2]);
                    data[tid] = partial_reduction = op(partial_reduction, data[tid +  1]);
                }
            }
        };

        ///////////////////////////////////////////////////////////////////////////////
        // PredValWarpReductor

        template <int n> struct PredValWarpReductor;
        template <> struct PredValWarpReductor<64>
        {
            template <typename T, typename V, typename Pred>
            static __device__ void reduce(T& myData, V& myVal, volatile T* sdata, V* sval, int tid, const Pred& pred)
            {
                if (tid < 32)
                {
                    myData = sdata[tid];
                    myVal = sval[tid];

                    T reg = sdata[tid + 32];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 32];
                    }

                    reg = sdata[tid + 16];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 16];
                    }

                    reg = sdata[tid + 8];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 8];
                    }

                    reg = sdata[tid + 4];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 4];
                    }

                    reg = sdata[tid + 2];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 2];
                    }

                    reg = sdata[tid + 1];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 1];
                    }
                }
            }
        };
        template <> struct PredValWarpReductor<32>
        {
            template <typename T, typename V, typename Pred>
            static __device__ void reduce(T& myData, V& myVal, volatile T* sdata, V* sval, int tid, const Pred& pred)
            {
                if (tid < 16)
                {
                    myData = sdata[tid];
                    myVal = sval[tid];

                    T reg = sdata[tid + 16];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 16];
                    }

                    reg = sdata[tid + 8];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 8];
                    }

                    reg = sdata[tid + 4];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 4];
                    }

                    reg = sdata[tid + 2];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 2];
                    }

                    reg = sdata[tid + 1];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 1];
                    }
                }
            }
        };

        template <> struct PredValWarpReductor<16>
        {
            template <typename T, typename V, typename Pred>
            static __device__ void reduce(T& myData, V& myVal, volatile T* sdata, V* sval, int tid, const Pred& pred)
            {
                if (tid < 8)
                {
                    myData = sdata[tid];
                    myVal = sval[tid];

                    T reg = reg = sdata[tid + 8];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 8];
                    }

                    reg = sdata[tid + 4];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 4];
                    }

                    reg = sdata[tid + 2];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 2];
                    }

                    reg = sdata[tid + 1];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 1];
                    }
                }
            }
        };
        template <> struct PredValWarpReductor<8>
        {
            template <typename T, typename V, typename Pred>
            static __device__ void reduce(T& myData, V& myVal, volatile T* sdata, V* sval, int tid, const Pred& pred)
            {
                if (tid < 4)
                {
                    myData = sdata[tid];
                    myVal = sval[tid];

                    T reg = reg = sdata[tid + 4];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 4];
                    }

                    reg = sdata[tid + 2];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 2];
                    }

                    reg = sdata[tid + 1];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 1];
                    }
                }
            }
        };

        template <bool warp> struct PredValReductionDispatcher;
        template <> struct PredValReductionDispatcher<true>
        {
            template <int n, typename T, typename V, typename Pred> static __device__ void reduce(T& myData, V& myVal, volatile T* sdata, V* sval, int tid, const Pred& pred)
            {
                PredValWarpReductor<n>::reduce(myData, myVal, sdata, sval, tid, pred);
            }
        };
        template <> struct PredValReductionDispatcher<false>
        {
            template <int n, typename T, typename V, typename Pred> static __device__ void reduce(T& myData, V& myVal, volatile T* sdata, V* sval, int tid, const Pred& pred)
            {
                myData = sdata[tid];
                myVal = sval[tid];

                if (n >= 512 && tid < 256)
                {
                    T reg = sdata[tid + 256];

                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 256];
                    }
                    __syncthreads();
                }
                if (n >= 256 && tid < 128)
                {
                    T reg = sdata[tid + 128];

                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 128];
                    }
                    __syncthreads();
                }
                if (n >= 128 && tid < 64)
                {
                    T reg = sdata[tid + 64];

                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval[tid] = myVal = sval[tid + 64];
                    }
                    __syncthreads();
                }

                if (tid < 32)
                {
                    if (n >= 64)
                    {
                        T reg = sdata[tid + 32];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval[tid] = myVal = sval[tid + 32];
                        }
                    }
                    if (n >= 32)
                    {
                        T reg = sdata[tid + 16];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval[tid] = myVal = sval[tid + 16];
                        }
                    }
                    if (n >= 16)
                    {
                        T reg = sdata[tid + 8];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval[tid] = myVal = sval[tid + 8];
                        }
                    }
                    if (n >= 8)
                    {
                        T reg = sdata[tid + 4];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval[tid] = myVal = sval[tid + 4];
                        }
                    }
                    if (n >= 4)
                    {
                        T reg = sdata[tid + 2];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval[tid] = myVal = sval[tid + 2];
                        }
                    }
                    if (n >= 2)
                    {
                        T reg = sdata[tid + 1];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval[tid] = myVal = sval[tid + 1];
                        }
                    }
                }
            }
        };

        ///////////////////////////////////////////////////////////////////////////////
        // PredVal2WarpReductor

        template <int n> struct PredVal2WarpReductor;
        template <> struct PredVal2WarpReductor<64>
        {
            template <typename T, typename V1, typename V2, typename Pred>
            static __device__ void reduce(T& myData, V1& myVal1, V2& myVal2, volatile T* sdata, V1* sval1, V2* sval2, int tid, const Pred& pred)
            {
                if (tid < 32)
                {
                    myData = sdata[tid];
                    myVal1 = sval1[tid];
                    myVal2 = sval2[tid];

                    T reg = sdata[tid + 32];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 32];
                        sval2[tid] = myVal2 = sval2[tid + 32];
                    }

                    reg = sdata[tid + 16];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 16];
                        sval2[tid] = myVal2 = sval2[tid + 16];
                    }

                    reg = sdata[tid + 8];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 8];
                        sval2[tid] = myVal2 = sval2[tid + 8];
                    }

                    reg = sdata[tid + 4];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 4];
                        sval2[tid] = myVal2 = sval2[tid + 4];
                    }

                    reg = sdata[tid + 2];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 2];
                        sval2[tid] = myVal2 = sval2[tid + 2];
                    }

                    reg = sdata[tid + 1];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 1];
                        sval2[tid] = myVal2 = sval2[tid + 1];
                    }
                }
            }
        };
        template <> struct PredVal2WarpReductor<32>
        {
            template <typename T, typename V1, typename V2, typename Pred>
            static __device__ void reduce(T& myData, V1& myVal1, V2& myVal2, volatile T* sdata, V1* sval1, V2* sval2, int tid, const Pred& pred)
            {
                if (tid < 16)
                {
                    myData = sdata[tid];
                    myVal1 = sval1[tid];
                    myVal2 = sval2[tid];

                    T reg = sdata[tid + 16];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 16];
                        sval2[tid] = myVal2 = sval2[tid + 16];
                    }

                    reg = sdata[tid + 8];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 8];
                        sval2[tid] = myVal2 = sval2[tid + 8];
                    }

                    reg = sdata[tid + 4];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 4];
                        sval2[tid] = myVal2 = sval2[tid + 4];
                    }

                    reg = sdata[tid + 2];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 2];
                        sval2[tid] = myVal2 = sval2[tid + 2];
                    }

                    reg = sdata[tid + 1];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 1];
                        sval2[tid] = myVal2 = sval2[tid + 1];
                    }
                }
            }
        };

        template <> struct PredVal2WarpReductor<16>
        {
            template <typename T, typename V1, typename V2, typename Pred>
            static __device__ void reduce(T& myData, V1& myVal1, V2& myVal2, volatile T* sdata, V1* sval1, V2* sval2, int tid, const Pred& pred)
            {
                if (tid < 8)
                {
                    myData = sdata[tid];
                    myVal1 = sval1[tid];
                    myVal2 = sval2[tid];

                    T reg = reg = sdata[tid + 8];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 8];
                        sval2[tid] = myVal2 = sval2[tid + 8];
                    }

                    reg = sdata[tid + 4];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 4];
                        sval2[tid] = myVal2 = sval2[tid + 4];
                    }

                    reg = sdata[tid + 2];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 2];
                        sval2[tid] = myVal2 = sval2[tid + 2];
                    }

                    reg = sdata[tid + 1];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 1];
                        sval2[tid] = myVal2 = sval2[tid + 1];
                    }
                }
            }
        };
        template <> struct PredVal2WarpReductor<8>
        {
            template <typename T, typename V1, typename V2, typename Pred>
            static __device__ void reduce(T& myData, V1& myVal1, V2& myVal2, volatile T* sdata, V1* sval1, V2* sval2, int tid, const Pred& pred)
            {
                if (tid < 4)
                {
                    myData = sdata[tid];
                    myVal1 = sval1[tid];
                    myVal2 = sval2[tid];

                    T reg = reg = sdata[tid + 4];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 4];
                        sval2[tid] = myVal2 = sval2[tid + 4];
                    }

                    reg = sdata[tid + 2];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 2];
                        sval2[tid] = myVal2 = sval2[tid + 2];
                    }

                    reg = sdata[tid + 1];
                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 1];
                        sval2[tid] = myVal2 = sval2[tid + 1];
                    }
                }
            }
        };

        template <bool warp> struct PredVal2ReductionDispatcher;
        template <> struct PredVal2ReductionDispatcher<true>
        {
            template <int n, typename T, typename V1, typename V2, typename Pred>
            static __device__ void reduce(T& myData, V1& myVal1, V2& myVal2, volatile T* sdata, V1* sval1, V2* sval2, int tid, const Pred& pred)
            {
                PredVal2WarpReductor<n>::reduce(myData, myVal1, myVal2, sdata, sval1, sval2, tid, pred);
            }
        };
        template <> struct PredVal2ReductionDispatcher<false>
        {
            template <int n, typename T, typename V1, typename V2, typename Pred>
            static __device__ void reduce(T& myData, V1& myVal1, V2& myVal2, volatile T* sdata, V1* sval1, V2* sval2, int tid, const Pred& pred)
            {
                myData = sdata[tid];
                myVal1 = sval1[tid];
                myVal2 = sval2[tid];

                if (n >= 512 && tid < 256)
                {
                    T reg = sdata[tid + 256];

                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 256];
                        sval2[tid] = myVal2 = sval2[tid + 256];
                    }
                    __syncthreads();
                }
                if (n >= 256 && tid < 128)
                {
                    T reg = sdata[tid + 128];

                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 128];
                        sval2[tid] = myVal2 = sval2[tid + 128];
                    }
                    __syncthreads();
                }
                if (n >= 128 && tid < 64)
                {
                    T reg = sdata[tid + 64];

                    if (pred(reg, myData))
                    {
                        sdata[tid] = myData = reg;
                        sval1[tid] = myVal1 = sval1[tid + 64];
                        sval2[tid] = myVal2 = sval2[tid + 64];
                    }
                    __syncthreads();
                }

                if (tid < 32)
                {
                    if (n >= 64)
                    {
                        T reg = sdata[tid + 32];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval1[tid] = myVal1 = sval1[tid + 32];
                            sval2[tid] = myVal2 = sval2[tid + 32];
                        }
                    }
                    if (n >= 32)
                    {
                        T reg = sdata[tid + 16];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval1[tid] = myVal1 = sval1[tid + 16];
                            sval2[tid] = myVal2 = sval2[tid + 16];
                        }
                    }
                    if (n >= 16)
                    {
                        T reg = sdata[tid + 8];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval1[tid] = myVal1 = sval1[tid + 8];
                            sval2[tid] = myVal2 = sval2[tid + 8];
                        }
                    }
                    if (n >= 8)
                    {
                        T reg = sdata[tid + 4];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval1[tid] = myVal1 = sval1[tid + 4];
                            sval2[tid] = myVal2 = sval2[tid + 4];
                        }
                    }
                    if (n >= 4)
                    {
                        T reg = sdata[tid + 2];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval1[tid] = myVal1 = sval1[tid + 2];
                            sval2[tid] = myVal2 = sval2[tid + 2];
                        }
                    }
                    if (n >= 2)
                    {
                        T reg = sdata[tid + 1];

                        if (pred(reg, myData))
                        {
                            sdata[tid] = myData = reg;
                            sval1[tid] = myVal1 = sval1[tid + 1];
                            sval2[tid] = myVal2 = sval2[tid + 1];
                        }
                    }
                }
            }
        };
    } // namespace utility_detail
}}} // namespace cv { namespace gpu { namespace device

#endif // __OPENCV_GPU_REDUCTION_DETAIL_HPP__