Tensor_logical.h

Go to the documentation of this file.

/*************************************************************************************
 
    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);