34 template<typename Scalar> class JacobiRotation
37     typedef typename NumTraits<Scalar>::Real RealScalar;
43     JacobiRotation(const Scalar& c, const Scalar& s) : m_c(c), m_s(s) {}  in JacobiRotation()
45     Scalar& c() { return m_c; }  in c()
46     Scalar c() const { return m_c; }  in c()
47     Scalar& s() { return m_s; }  in s()
48     Scalar s() const { return m_s; }  in s()
66     bool makeJacobi(const RealScalar& x, const Scalar& y, const RealScalar& z);
68     void makeGivens(const Scalar& p, const Scalar& q, Scalar* z=0);
71     void makeGivens(const Scalar& p, const Scalar& q, Scalar* z, internal::true_type);
72     void makeGivens(const Scalar& p, const Scalar& q, Scalar* z, internal::false_type);
74     Scalar m_c, m_s;
82 template<typename Scalar>
83 bool JacobiRotation<Scalar>::makeJacobi(const RealScalar& x, const Scalar& y, const RealScalar& z)  in makeJacobi()
87   typedef typename NumTraits<Scalar>::Real RealScalar;  in makeJacobi()
91     m_c = Scalar(1);  in makeJacobi()
92     m_s = Scalar(0);  in makeJacobi()
125 template<typename Scalar>
127 inline bool JacobiRotation<Scalar>::makeJacobi(const MatrixBase<Derived>& m, Index p, Index q)  in makeJacobi()
148 template<typename Scalar>
149 void JacobiRotation<Scalar>::makeGivens(const Scalar& p, const Scalar& q, Scalar* z)  in makeGivens()
151 …makeGivens(p, q, z, typename internal::conditional<NumTraits<Scalar>::IsComplex, internal::true_ty…  in makeGivens()
156 template<typename Scalar>
157 void JacobiRotation<Scalar>::makeGivens(const Scalar& p, const Scalar& q, Scalar* r, internal::true…  in makeGivens()
163   if(q==Scalar(0))  in makeGivens()
165     m_c = numext::real(p)<0 ? Scalar(-1) : Scalar(1);  in makeGivens()
169   else if(p==Scalar(0))  in makeGivens()
181       Scalar ps = p / p1;  in makeGivens()
183       Scalar qs = q / p1;  in makeGivens()
190       m_c = Scalar(1)/u;  in makeGivens()
196       Scalar ps = p / q1;  in makeGivens()
198       Scalar qs = q / q1;  in makeGivens()
215 template<typename Scalar>
216 void JacobiRotation<Scalar>::makeGivens(const Scalar& p, const Scalar& q, Scalar* r, internal::fals…  in makeGivens()
220   if(q==Scalar(0))  in makeGivens()
222     m_c = p<Scalar(0) ? Scalar(-1) : Scalar(1);  in makeGivens()
223     m_s = Scalar(0);  in makeGivens()
226   else if(p==Scalar(0))  in makeGivens()
228     m_c = Scalar(0);  in makeGivens()
229     m_s = q<Scalar(0) ? Scalar(1) : Scalar(-1);  in makeGivens()
234     Scalar t = q/p;  in makeGivens()
235     Scalar u = sqrt(Scalar(1) + numext::abs2(t));  in makeGivens()
236     if(p<Scalar(0))  in makeGivens()
238     m_c = Scalar(1)/u;  in makeGivens()
244     Scalar t = p/q;  in makeGivens()
245     Scalar u = sqrt(Scalar(1) + numext::abs2(t));  in makeGivens()
246     if(q<Scalar(0))  in makeGivens()
248     m_s = -Scalar(1)/u;  in makeGivens()
304   typedef typename VectorX::Scalar Scalar;  in apply_rotation_in_the_plane()  typedef
306     PacketSize = packet_traits<Scalar>::size,  in apply_rotation_in_the_plane()
309   typedef typename packet_traits<Scalar>::type Packet;  in apply_rotation_in_the_plane()
316   Scalar* EIGEN_RESTRICT x = &xpr_x.derived().coeffRef(0);  in apply_rotation_in_the_plane()
317   Scalar* EIGEN_RESTRICT y = &xpr_y.derived().coeffRef(0);  in apply_rotation_in_the_plane()
344       Scalar xi = x[i];  in apply_rotation_in_the_plane()
345       Scalar yi = y[i];  in apply_rotation_in_the_plane()
350     Scalar* EIGEN_RESTRICT px = x + alignedStart;  in apply_rotation_in_the_plane()
351     Scalar* EIGEN_RESTRICT py = y + alignedStart;  in apply_rotation_in_the_plane()
392       Scalar xi = x[i];  in apply_rotation_in_the_plane()
393       Scalar yi = y[i];  in apply_rotation_in_the_plane()
409     Scalar* EIGEN_RESTRICT px = x;  in apply_rotation_in_the_plane()
410     Scalar* EIGEN_RESTRICT py = y;  in apply_rotation_in_the_plane()
427       Scalar xi = *x;  in apply_rotation_in_the_plane()
428       Scalar yi = *y;  in apply_rotation_in_the_plane()