/*

* Simd Library (http://ermig1979.github.io/Simd).

*

* Copyright (c) 2011-2017 Yermalayeu Ihar.

*

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

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

* in the Software without restriction, including without limitation the rights

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

* copies of the Software, and to permit persons to whom the Software is

* furnished to do so, subject to the following conditions:

*

* The above copyright notice and this permission notice shall be included in

* all copies or substantial portions of the Software.

*

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

* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,

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

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

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

* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE

* SOFTWARE.

*/

#include "Simd/SimdMemory.h"

#include "Simd/SimdStore.h"

#include "Simd/SimdCompare.h"

namespace Simd

{

#ifdef SIMD_AVX2_ENABLE    

    namespace Avx2

    {

        template <bool align> void LbpEstimate(const uint8_t * src, ptrdiff_t stride, uint8_t * dst)

        {

            __m256i threshold = Load<false>((__m256i*)src);

            __m256i lbp = _mm256_setzero_si256();

            lbp = _mm256_or_si256(lbp, _mm256_and_si256(GreaterOrEqual8u(Load<align>((__m256i*)(src - 1 - stride)), threshold), K8_01));

            lbp = _mm256_or_si256(lbp, _mm256_and_si256(GreaterOrEqual8u(Load<false>((__m256i*)(src - stride)), threshold), K8_02));

            lbp = _mm256_or_si256(lbp, _mm256_and_si256(GreaterOrEqual8u(Load<false>((__m256i*)(src + 1 - stride)), threshold), K8_04));

            lbp = _mm256_or_si256(lbp, _mm256_and_si256(GreaterOrEqual8u(Load<false>((__m256i*)(src + 1)), threshold), K8_08));

            lbp = _mm256_or_si256(lbp, _mm256_and_si256(GreaterOrEqual8u(Load<false>((__m256i*)(src + 1 + stride)), threshold), K8_10));

            lbp = _mm256_or_si256(lbp, _mm256_and_si256(GreaterOrEqual8u(Load<false>((__m256i*)(src + stride)), threshold), K8_20));

            lbp = _mm256_or_si256(lbp, _mm256_and_si256(GreaterOrEqual8u(Load<align>((__m256i*)(src - 1 + stride)), threshold), K8_40));

            lbp = _mm256_or_si256(lbp, _mm256_and_si256(GreaterOrEqual8u(Load<align>((__m256i*)(src - 1)), threshold), K8_80));

            Store<false>((__m256i*)dst, lbp);

        }

Unexecuted instantiation: void Simd::Avx2::LbpEstimate<true>(unsigned char const*, long, unsigned char*)

Unexecuted instantiation: void Simd::Avx2::LbpEstimate<false>(unsigned char const*, long, unsigned char*)

        template <bool align> void LbpEstimate(

            const uint8_t * src, size_t srcStride, size_t width, size_t height, uint8_t * dst, size_t dstStride)

        {

            assert(width >= A + 2);

            if (align)

                assert(Aligned(src) && Aligned(srcStride) && Aligned(dst) && Aligned(dstStride));

            size_t alignedWidth = AlignLo(width - 2, A) + 1;

            memset(dst, 0, width);

            src += srcStride;

            dst += dstStride;

            for (size_t row = 2; row < height; ++row)

            {

                dst[0] = 0;

                for (size_t col = 1; col < alignedWidth; col += A)

                    LbpEstimate<align>(src + col, srcStride, dst + col);

                if (alignedWidth != width - 1)

                    LbpEstimate<false>(src + width - 1 - A, srcStride, dst + width - 1 - A);

                dst[width - 1] = 0;

                src += srcStride;

                dst += dstStride;

            }

            memset(dst, 0, width);

        }

Unexecuted instantiation: void Simd::Avx2::LbpEstimate<true>(unsigned char const*, unsigned long, unsigned long, unsigned long, unsigned char*, unsigned long)

Unexecuted instantiation: void Simd::Avx2::LbpEstimate<false>(unsigned char const*, unsigned long, unsigned long, unsigned long, unsigned char*, unsigned long)

        void LbpEstimate(const uint8_t * src, size_t srcStride, size_t width, size_t height, uint8_t * dst, size_t dstStride)

        {

            if (Aligned(src) && Aligned(srcStride) && Aligned(dst) && Aligned(dstStride))

                LbpEstimate<true>(src, srcStride, width, height, dst, dstStride);

            else

                LbpEstimate<false>(src, srcStride, width, height, dst, dstStride);

        }

    }

#endif// SIMD_AVX2_ENABLE

}