17 m_ye = std::cos(qkl)* dxyz[1] + std::sin(qkl)*dir[1]/qkc ;
18 m_yae = -qk*std::sin(qkl)* dxyz[1] + std::cos(qkl)*dir[1]/dir[2] ;
19 m_xe = std::cosh(qkl)* dxyz[0] + std::sinh(qkl)*dir[0]/qkc ;
20 m_xae = qk*std::sinh(qkl)* dxyz[0] + std::cosh(qkl)*dir[0]/dir[2] ;