//------------------------------------------------- inline double energy(CpxNumTns& t) { double val=0; cpx* data = t.data(); for(int i=0; i<t.m()*t.n()*t.p(); i++) { val += norm(data[i]); } return val; }
inline float energy(CpxNumTns& m) { float val=0; cpx* data = m.data(); for(int i=0; i<m.m()*m.n()*m.p(); i++) { //val += data[i].re*data[i].re + data[i].im*data[i].im; val += norm(data[i]); } return val; }