mirror of
https://github.com/paboyle/Grid.git
synced 2024-11-14 17:55:38 +00:00
93 lines
2.6 KiB
C++
93 lines
2.6 KiB
C++
/*************************************************************************************
|
|
|
|
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 */
|
|
#ifndef GRID_TENSOR_LOGICAL_H
|
|
#define GRID_TENSOR_LOGICAL_H
|
|
|
|
namespace Grid {
|
|
|
|
#define LOGICAL_BINOP(Op)\
|
|
template<class v> strong_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> strong_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> strong_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>
|
|
strong_inline bool operator==(const iScalar<T> &t1, const iScalar<T> &t2)
|
|
{
|
|
return (t1._internal == t2._internal);
|
|
}
|
|
|
|
template <class T, int N>
|
|
strong_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>
|
|
strong_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;
|
|
}
|
|
|
|
}
|
|
#endif
|