ivec.h

Go to the documentation of this file.

/*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 Digital Holographic Library
//
// Openholo library is free software;
// you can redistribute it and/or modify it under the terms of the BSD 2-Clause license.
//
// Copyright (C) 2017-2024, Korea Electronics Technology Institute. All rights reserved.
// E-mail : contact.openholo@gmail.com
// Web : http://www.openholo.org
//
// Redistribution and use in source and binary forms, with or without modification,
// are permitted provided that the following conditions are met:
//
//  1. Redistribution's of source code must retain the above copyright notice,
//     this list of conditions and the following disclaimer.
//
//  2. 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.
//
// 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 copyright holder 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.
//
// This software contains opensource software released under GNU Generic Public License,
// NVDIA Software License Agreement, or CUDA supplement to Software License Agreement.
// Check whether software you use contains licensed software.
//
//M*/

#ifndef __ivec_h
#define __ivec_h
//|
//| ivec : n-dimensional ivector
//|

#include <stdio.h>
#ifndef OPH_DLL
#ifdef _WIN32
#define OPH_DLL __declspec(dllexport)
#else
#define OPH_DLL __attribute__((visibility("default")))
#endif
#endif

namespace oph {

    struct OPH_DLL ivec2 {
        int v[2];
        static const int n;

        inline ivec2()
        {
            v[0] = v[1] = 0;
        };

        inline ivec2(int a)
        {
            v[0] = a;  v[1] = a;
        }

        inline ivec2(int v_1, int v_2)
        {
            v[0] = v_1;  v[1] = v_2;
        }

        inline ivec2(const ivec2& a)
        {
            v[0] = a[0];  v[1] = a[1];
        }

        inline ivec2& operator=(const ivec2& a)
        {
            v[0] = a[0];  v[1] = a[1];
            return *this;
        }


        inline int& operator[] (int i) { return v[i]; }
        inline const int& operator[] (int i) const { return v[i]; }
        inline int& operator() (int i) { return v[i % 2]; }
        inline const int& operator() (int i) const { return v[i % 2]; }
    };




    //| binary op : componentwise


    inline ivec2 operator + (const ivec2& a, const ivec2& b)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = a[i] + b[i]; }
        return c;
    }

    inline ivec2 operator + (int a, const ivec2& b)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = a + b[i]; }
        return c;
    }

    inline ivec2 operator + (const ivec2& a, int b)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = a[i] + b; }
        return c;
    }



    inline ivec2 operator - (const ivec2& a, const ivec2& b)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = a[i] - b[i]; }
        return c;
    }

    inline ivec2 operator - (int a, const ivec2& b)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = a - b[i]; }
        return c;
    }

    inline ivec2 operator - (const ivec2& a, int b)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = a[i] - b; }
        return c;
    }



    inline ivec2 operator * (const ivec2& a, const ivec2& b)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = a[i] * b[i]; }
        return c;
    }

    inline ivec2 operator * (int a, const ivec2& b)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = a * b[i]; }
        return c;
    }

    inline ivec2 operator * (const ivec2& a, int b)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = a[i] * b; }
        return c;
    }


    //
    //inline ivec2 operator / (const ivec2& a, const ivec2& b)
    //{
    //    ivec2 c;
    //    for(int i = 0; i < 2;++i) { c[i] = a[i] / b[i]; }
    //    return c;
    //}
    //
    //inline ivec2 operator / (int a, const ivec2& b)
    //{
    //    ivec2 c;
    //    for(int i = 0; i < 2;++i) { c[i] = a / b[i]; }
    //    return c;
    //}
    //
    //inline ivec2 operator / (const ivec2& a, int b)
    //{
    //    ivec2 c;
    //    for(int i = 0; i < 2;++i) { c[i] = a[i] / b; }
    //    return c;
    //}



    //| cumulative op : componentwise


    inline ivec2 operator += (ivec2& a, const ivec2& b)
    {
        return a = (a + b);
    }

    inline ivec2 operator += (ivec2& a, int b)
    {
        return a = (a + b);
    }



    inline ivec2 operator -= (ivec2& a, const ivec2& b)
    {
        return a = (a - b);
    }

    inline ivec2 operator -= (ivec2& a, int b)
    {
        return a = (a - b);
    }



    inline ivec2 operator *= (ivec2& a, const ivec2& b)
    {
        return a = (a * b);
    }

    inline ivec2 operator *= (ivec2& a, int b)
    {
        return a = (a * b);
    }



    //inline ivec2 operator /= (ivec2& a, const ivec2& b)
    //{
    //    return a = (a / b);
    //}
    //
    //inline ivec2 operator /= (ivec2& a, int b)
    //{
    //    return a = (a / b);
    //}



    //| logical op : componentwise


    inline int operator == (const ivec2& a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] == b[i]); }
        return c;
    }

    inline int operator == (int a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a == b[i]); }
        return c;
    }

    inline int operator == (const ivec2& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] == b); }
        return c;
    }



    inline int operator < (const ivec2& a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] < b[i]); }
        return c;
    }

    inline int operator < (int a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a < b[i]); }
        return c;
    }

    inline int operator < (const ivec2& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] < b); }
        return c;
    }



    inline int operator <= (const ivec2& a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] <= b[i]); }
        return c;
    }

    inline int operator <= (int a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a <= b[i]); }
        return c;
    }

    inline int operator <= (const ivec2& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] <= b); }
        return c;
    }



    inline int operator > (const ivec2& a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] > b[i]); }
        return c;
    }

    inline int operator > (int a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a > b[i]); }
        return c;
    }

    inline int operator > (const ivec2& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] > b); }
        return c;
    }




    inline int operator >= (const ivec2& a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] >= b[i]); }
        return c;
    }
    inline int operator >= (int a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a >= b[i]); }
        return c;
    }

    inline int operator >= (const ivec2& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] >= b); }
        return c;
    }

    inline int operator != (const ivec2& a, const ivec2& b)
    {
        int c = 1;
        for (int i = 0; i < 2; ++i) { c = c & (a[i] != b[i]); }
        return c;
    }

    //| unary op : componentwise
    inline ivec2 operator - (const ivec2& a)
    {
        ivec2 c;
        for (int i = 0; i < 2; ++i) { c[i] = -a[i]; }
        return c;
    }

    struct OPH_DLL ivec3 {
        int v[3];
        static const int n;

        inline ivec3()
        {
            v[0] = v[1] = v[2] = 0;
        };

        inline ivec3(int a)
        {
            v[0] = v[1] = v[2] = a;
        }

        inline ivec3(int v_1, int v_2, int v_3)
        {
            v[0] = v_1;  v[1] = v_2;  v[2] = v_3;
        }

        inline ivec3(const ivec3& a)
        {
            v[0] = a[0];  v[1] = a[1];  v[2] = a[2];
        }

        inline ivec3& operator=(const ivec3& a)
        {
            v[0] = a[0];  v[1] = a[1];  v[2] = a[2];
            return *this;
        }


        inline int& operator[] (int i) { return v[i]; }
        inline const int& operator[] (int i) const { return v[i]; }
        inline int& operator() (int i) { return v[i % 3]; }
        inline const int& operator() (int i) const { return v[i % 3]; }
    };
    //| binary op : componentwise

    inline ivec3 operator + (const ivec3& a, const ivec3& b)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = a[i] + b[i]; }
        return c;
    }

    inline ivec3 operator + (int a, const ivec3& b)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = a + b[i]; }
        return c;
    }

    inline ivec3 operator + (const ivec3& a, int b)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = a[i] + b; }
        return c;
    }



    inline ivec3 operator - (const ivec3& a, const ivec3& b)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = a[i] - b[i]; }
        return c;
    }

    inline ivec3 operator - (int a, const ivec3& b)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = a - b[i]; }
        return c;
    }

    inline ivec3 operator - (const ivec3& a, int b)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = a[i] - b; }
        return c;
    }



    inline ivec3 operator * (const ivec3& a, const ivec3& b)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = a[i] * b[i]; }
        return c;
    }

    inline ivec3 operator * (int a, const ivec3& b)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = a * b[i]; }
        return c;
    }

    inline ivec3 operator * (const ivec3& a, int b)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = a[i] * b; }
        return c;
    }

    //
    //
    //inline ivec3 operator / (const ivec3& a, const ivec3& b)
    //{
    //    ivec3 c;
    //    for(int i = 0; i < 3;++i) { c[i] = a[i] / b[i]; }
    //    return c;
    //}
    //
    //inline ivec3 operator / (int a, const ivec3& b)
    //{
    //    ivec3 c;
    //    for(int i = 0; i < 3;++i) { c[i] = a / b[i]; }
    //    return c;
    //}
    //
    //inline ivec3 operator / (const ivec3& a, int b)
    //{
    //    ivec3 c;
    //    for(int i = 0; i < 3;++i) { c[i] = a[i] / b; }
    //    return c;
    //}



    //| cumulative op : componentwise


    inline ivec3 operator += (ivec3& a, const ivec3& b)
    {
        return a = (a + b);
    }

    inline ivec3 operator += (ivec3& a, int b)
    {
        return a = (a + b);
    }



    inline ivec3 operator -= (ivec3& a, const ivec3& b)
    {
        return a = (a - b);
    }

    inline ivec3 operator -= (ivec3& a, int b)
    {
        return a = (a - b);
    }



    inline ivec3 operator *= (ivec3& a, const ivec3& b)
    {
        return a = (a * b);
    }

    inline ivec3 operator *= (ivec3& a, int b)
    {
        return a = (a * b);
    }


    //
    //inline ivec3 operator /= (ivec3& a, const ivec3& b)
    //{
    //    return a = (a / b);
    //}
    //
    //inline ivec3 operator /= (ivec3& a, int b)
    //{
    //    return a = (a / b);
    //}



    //| logical op : componentwise


    inline int operator == (const ivec3& a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] == b[i]); }
        return c;
    }

    inline int operator == (int a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a == b[i]); }
        return c;
    }

    inline int operator == (const ivec3& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] == b); }
        return c;
    }



    inline int operator < (const ivec3& a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] < b[i]); }
        return c;
    }

    inline int operator < (int a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a < b[i]); }
        return c;
    }

    inline int operator < (const ivec3& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] < b); }
        return c;
    }



    inline int operator <= (const ivec3& a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] <= b[i]); }
        return c;
    }

    inline int operator <= (int a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a <= b[i]); }
        return c;
    }

    inline int operator <= (const ivec3& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] <= b); }
        return c;
    }



    inline int operator > (const ivec3& a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] > b[i]); }
        return c;
    }

    inline int operator > (int a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a > b[i]); }
        return c;
    }

    inline int operator > (const ivec3& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] > b); }
        return c;
    }



    inline int operator >= (const ivec3& a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] >= b[i]); }
        return c;
    }

    inline int operator >= (int a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a >= b[i]); }
        return c;
    }

    inline int operator >= (const ivec3& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] >= b); }
        return c;
    }

    inline int operator != (const ivec3& a, const ivec3& b)
    {
        int c = 1;
        for (int i = 0; i < 3; ++i) { c = c & (a[i] != b[i]); }
        return c;
    }

    //| unary op : componentwise
    inline ivec3 operator - (const ivec3& a)
    {
        ivec3 c;
        for (int i = 0; i < 3; ++i) { c[i] = -a[i]; }
        return c;
    }

    struct OPH_DLL ivec4 {
        int v[4];
        static const int n;

        inline ivec4()
        {
            v[0] = v[1] = v[2] = v[3] = 0;
        };

        inline ivec4(int a)
        {
            v[0] = a;  v[1] = a;  v[2] = a;  v[3] = a;
        }

        inline ivec4(int v_1, int v_2, int v_3, int v_4)
        {
            v[0] = v_1;  v[1] = v_2;  v[2] = v_3;  v[3] = v_4;
        }

        inline ivec4(const ivec4& a)
        {
            v[0] = a[0];  v[1] = a[1];  v[2] = a[2];  v[3] = a[3];
        }

        inline ivec4& operator=(const ivec4& a)
        {
            v[0] = a[0];  v[1] = a[1];  v[2] = a[2];  v[3] = a[3];
            return *this;
        }


        inline int& operator[] (int i) { return v[i]; }
        inline const int& operator[] (int i) const { return v[i]; }
        inline int& operator() (int i) { return v[i % 4]; }
        inline const int& operator() (int i) const { return v[i % 4]; }
    };




    //| binary op : componentwise


    inline ivec4 operator + (const ivec4& a, const ivec4& b)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = a[i] + b[i]; }
        return c;
    }

    inline ivec4 operator + (int a, const ivec4& b)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = a + b[i]; }
        return c;
    }

    inline ivec4 operator + (const ivec4& a, int b)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = a[i] + b; }
        return c;
    }



    inline ivec4 operator - (const ivec4& a, const ivec4& b)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = a[i] - b[i]; }
        return c;
    }

    inline ivec4 operator - (int a, const ivec4& b)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = a - b[i]; }
        return c;
    }

    inline ivec4 operator - (const ivec4& a, int b)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = a[i] - b; }
        return c;
    }



    inline ivec4 operator * (const ivec4& a, const ivec4& b)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = a[i] * b[i]; }
        return c;
    }

    inline ivec4 operator * (int a, const ivec4& b)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = a * b[i]; }
        return c;
    }

    inline ivec4 operator * (const ivec4& a, int b)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = a[i] * b; }
        return c;
    }


    //
    //inline ivec4 operator / (const ivec4& a, const ivec4& b)
    //{
    //    ivec4 c;
    //    for(int i = 0; i < 4;++i) { c[i] = a[i] / b[i]; }
    //    return c;
    //}
    //
    //inline ivec4 operator / (int a, const ivec4& b)
    //{
    //    ivec4 c;
    //    for(int i = 0; i < 4;++i) { c[i] = a / b[i]; }
    //    return c;
    //}
    //
    //inline ivec4 operator / (const ivec4& a, int b)
    //{
    //    ivec4 c;
    //    for(int i = 0; i < 4;++i) { c[i] = a[i] / b; }
    //    return c;
    //}



    //| cumulative op : componentwise


    inline ivec4 operator += (ivec4& a, const ivec4& b)
    {
        return a = (a + b);
    }

    inline ivec4 operator += (ivec4& a, int b)
    {
        return a = (a + b);
    }



    inline ivec4 operator -= (ivec4& a, const ivec4& b)
    {
        return a = (a - b);
    }

    inline ivec4 operator -= (ivec4& a, int b)
    {
        return a = (a - b);
    }



    inline ivec4 operator *= (ivec4& a, const ivec4& b)
    {
        return a = (a * b);
    }

    inline ivec4 operator *= (ivec4& a, int b)
    {
        return a = (a * b);
    }



    //inline ivec4 operator /= (ivec4& a, const ivec4& b)
    //{
    //    return a = (a / b);
    //}
    //
    //inline ivec4 operator /= (ivec4& a, int b)
    //{
    //    return a = (a / b);
    //}



    //| logical op : componentwise


    inline int operator == (const ivec4& a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] == b[i]); }
        return c;
    }

    inline int operator == (int a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a == b[i]); }
        return c;
    }

    inline int operator == (const ivec4& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] == b); }
        return c;
    }



    inline int operator < (const ivec4& a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] < b[i]); }
        return c;
    }

    inline int operator < (int a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a < b[i]); }
        return c;
    }

    inline int operator < (const ivec4& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] < b); }
        return c;
    }



    inline int operator <= (const ivec4& a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] <= b[i]); }
        return c;
    }

    inline int operator <= (int a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a <= b[i]); }
        return c;
    }

    inline int operator <= (const ivec4& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] <= b); }
        return c;
    }



    inline int operator > (const ivec4& a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] > b[i]); }
        return c;
    }

    inline int operator > (int a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a > b[i]); }
        return c;
    }

    inline int operator > (const ivec4& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] > b); }
        return c;
    }



    inline int operator >= (const ivec4& a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] >= b[i]); }
        return c;
    }

    inline int operator != (const ivec4& a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] != b[i]); }
        return c;
    }

    inline int operator >= (int a, const ivec4& b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a >= b[i]); }
        return c;
    }

    inline int operator >= (const ivec4& a, int b)
    {
        int c = 1;
        for (int i = 0; i < 4; ++i) { c = c & (a[i] >= b); }
        return c;
    }



    //| unary op : componentwise
    inline ivec4 operator - (const ivec4& a)
    {
        ivec4 c;
        for (int i = 0; i < 4; ++i) { c[i] = -a[i]; }
        return c;
    }

}; //namespace oph

#endif // !__ivec_h