mirror of
https://github.com/paboyle/Grid.git
synced 2024-11-10 15:55:37 +00:00
38 lines
837 B
C
38 lines
837 B
C
|
#ifndef GRID_TENSOR_UNARY_H
|
||
|
|
#define GRID_TENSOR_UNARY_H
|
||
|
|
namespace Grid {
|
||
|
|
|
||
|
|
#define UNARY_REAL(func)\
|
||
|
|
template<class obj> inline auto func(const iScalar<obj> &z) -> iScalar<obj>\
|
||
|
|
{\
|
||
|
|
iScalar<obj> ret;\
|
||
|
|
ret._internal = func( (z._internal));\
|
||
|
|
return ret;\
|
||
|
|
}\
|
||
|
|
template<class obj,int N> inline auto func(const iVector<obj,N> &z) -> iVector<obj,N>\
|
||
|
|
{\
|
||
|
|
iVector<obj,N> ret;\
|
||
|
|
for(int c1=0;c1<N;c1++){\
|
||
|
|
ret._internal[c1] = func( (z._internal[c1]));\
|
||
|
|
}\
|
||
|
|
return ret;\
|
||
|
|
}\
|
||
|
|
template<class obj,int N> inline auto func(const iMatrix<obj,N> &z) -> iMatrix<obj,N>\
|
||
|
|
{\
|
||
|
|
iMatrix<obj,N> ret;\
|
||
|
|
for(int c1=0;c1<N;c1++){\
|
||
|
|
for(int c2=0;c2<N;c2++){\
|
||
|
|
ret._internal[c1][c2] = func( (z._internal[c1][c2]));\
|
||
|
|
}}\
|
||
|
|
return ret;\
|
||
|
|
}
|
||
|
|
|
||
|
|
UNARY_REAL(sqrt);
|
||
|
|
UNARY_REAL(rsqrt);
|
||
|
|
UNARY_REAL(sin);
|
||
|
|
UNARY_REAL(cos);
|
||
|
|
|
||
|
|
|
||
|
|
}
|
||
|
|
#endif
|