/*************************************************************************************

    Grid physics library, www.github.com/paboyle/Grid 

    Source file: ./lib/tensors/Tensor_logical.h

    Copyright (C) 2015

Author: Azusa Yamaguchi <ayamaguc@staffmail.ed.ac.uk>

    This program is free software; you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation; either version 2 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License along
    with this program; if not, write to the Free Software Foundation, Inc.,
    51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.

    See the full license in the file "LICENSE" in the top level distribution directory
*************************************************************************************/
/*  END LEGAL */
#pragma once

NAMESPACE_BEGIN(Grid);

#define LOGICAL_BINOP(Op)						\
  template<class v> accelerator_inline iScalar<v> operator Op (const iScalar<v>& lhs,const iScalar<v>& rhs) \
  {									\
    iScalar<v> ret;							\
    ret._internal = lhs._internal Op rhs._internal ;			\
    return ret;								\
  }									\
  template<class l> accelerator_inline iScalar<l> operator Op (const iScalar<l>& lhs,Integer rhs) \
  {									\
    typename iScalar<l>::scalar_type t; t=rhs;				\
    typename iScalar<l>::tensor_reduced srhs; srhs=t;			\
    return lhs Op srhs;							\
  }									\
  template<class l> accelerator_inline iScalar<l> operator Op (Integer lhs,const iScalar<l>& rhs) \
  {									\
    typename iScalar<l>::scalar_type t;t=lhs;				\
    typename iScalar<l>::tensor_reduced slhs;slhs=t;			\
    return slhs Op rhs;							\
  }

LOGICAL_BINOP(|);
LOGICAL_BINOP(&);
LOGICAL_BINOP(||);
LOGICAL_BINOP(&&);

template <class T>
accelerator_inline bool operator==(const iScalar<T> &t1, const iScalar<T> &t2)
{
  return (t1._internal == t2._internal);
}

template <class T, int N>
accelerator_inline bool operator==(const iVector<T, N> &t1, const iVector<T, N> &t2)
{
  bool res = true;

  for (unsigned int i = 0; i < N; ++i)
  {
    res = (res && (t1._internal[i] == t2._internal[i]));
  }

  return res;
}

template <class T, int N>
accelerator_inline bool operator==(const iMatrix<T, N> &t1, const iMatrix<T, N> &t2)
{
  bool res = true;

  for (unsigned int i = 0; i < N; ++i)
  for (unsigned int j = 0; j < N; ++j)
  {
    res = (res && (t1._internal[i][j] == t2._internal[i][j]));
  }
  
  return res;
}

NAMESPACE_END(Grid);