8 template <
class T,
int Tensor_Dim>
class Tensor1
16 static_assert(
sizeof...(
d) ==
sizeof(
data) /
sizeof(
T),
17 "Incorrect number of Arguments. Constructor should "
18 "initialize the entire Tensor");
29 if(
N >= Tensor_Dim ||
N < 0)
32 s <<
"Bad index in Tensor1<T," << Tensor_Dim <<
">.operator(" <<
N
34 throw std::out_of_range(s.str());
42 if(
N >= Tensor_Dim ||
N < 0)
45 s <<
"Bad index in Tensor1<T," << Tensor_Dim <<
">.operator(" <<
N
46 <<
") const" << std::endl;
47 throw std::out_of_range(s.str());
58 template <
char i,
int Dim>
66 template <
char i,
int Dim>
88 return data[Current_Dim - 1] *
data[Current_Dim - 1]
96 template <
class T,
int Tensor_Dim>
101 for(
int i = 0;
i + 1 < Tensor_Dim; ++
i)
107 os << t(Tensor_Dim - 1);
113 template <
class T,
int Tensor_Dim>
118 for(
int i = 0;
i + 1 < Tensor_Dim; ++
i)
124 is >> t(Tensor_Dim - 1);
std::istream & operator>>(std::istream &is, FTensor::Tensor1< T, Tensor_Dim > &t)
std::ostream & operator<<(std::ostream &os, const FTensor::Tensor1< T, Tensor_Dim > &t)
JSON compatible output.
T l2_squared(const Number< Current_Dim > &) const
T operator()(const int N) const
T & operator()(const int N)
T l2_squared(const Number< 1 > &) const
constexpr Tensor1(U... d)
Tensor1< T, Tensor_Dim > normalize()
FTensor::Index< 'i', 3 > i
const Tensor1_Expr< const dTensor0< T, Dim, i >, typename promote< T, double >::V, Dim, i > d(const Tensor0< T * > &a, const Index< i, Dim > index, const Tensor1< int, Dim > &d_ijk, const Tensor1< double, Dim > &d_xyz)