35     Scalar ax = internal::abs(v.coeff(i));  in lapackNorm()
38       ssq += internal::abs2(ax/scale);  in lapackNorm()
42       ssq = Scalar(1) + ssq * internal::abs2(scale/ax);  in lapackNorm()
46   return scale * internal::sqrt(ssq);  in lapackNorm()
76   return internal::sqrt(v(0));  in divacNorm()
80 Packet4f internal::plt(const Packet4f& a, Packet4f& b) { return _mm_cmplt_ps(a,b); }  in plt()
81 Packet2d internal::plt(const Packet2d& a, Packet2d& b) { return _mm_cmplt_pd(a,b); }  in plt()
83 Packet4f internal::pandnot(const Packet4f& a, Packet4f& b) { return _mm_andnot_ps(a,b); }  in pandnot()
84 Packet2d internal::pandnot(const Packet2d& a, Packet2d& b) { return _mm_andnot_pd(a,b); }  in pandnot()
129     relerr  = internal::sqrt(eps);      // tolerance for neglecting asml  in pblueNorm()
135   typedef typename internal::packet_traits<Scalar>::type Packet;  in pblueNorm()
136   const int ps = internal::packet_traits<Scalar>::size;  in pblueNorm()
137   Packet pasml = internal::pset1(Scalar(0));  in pblueNorm()
138   Packet pamed = internal::pset1(Scalar(0));  in pblueNorm()
139   Packet pabig = internal::pset1(Scalar(0));  in pblueNorm()
140   Packet ps2m = internal::pset1(s2m);  in pblueNorm()
141   Packet ps1m = internal::pset1(s1m);  in pblueNorm()
142   Packet pb2  = internal::pset1(b2);  in pblueNorm()
143   Packet pb1  = internal::pset1(b1);  in pblueNorm()
146     Packet ax = internal::pabs(v.template packet<Aligned>(j));  in pblueNorm()
147     Packet ax_s2m = internal::pmul(ax,ps2m);  in pblueNorm()
148     Packet ax_s1m = internal::pmul(ax,ps1m);  in pblueNorm()
149     Packet maskBig = internal::plt(pb2,ax);  in pblueNorm()
150     Packet maskSml = internal::plt(ax,pb1);  in pblueNorm()
164     pabig = internal::padd(pabig, internal::pand(maskBig, internal::pmul(ax_s2m,ax_s2m)));  in pblueNorm()
165     pasml = internal::padd(pasml, internal::pand(maskSml, internal::pmul(ax_s1m,ax_s1m)));  in pblueNorm()
166 …pamed = internal::padd(pamed, internal::pandnot(internal::pmul(ax,ax),internal::pand(maskSml,maskB…  in pblueNorm()
168   Scalar abig = internal::predux(pabig);  in pblueNorm()
169   Scalar asml = internal::predux(pasml);  in pblueNorm()
170   Scalar amed = internal::predux(pamed);  in pblueNorm()
173     abig = internal::sqrt(abig);  in pblueNorm()
182       amed = internal::sqrt(amed);  in pblueNorm()
194       abig = internal::sqrt(amed);  in pblueNorm()
195       amed = internal::sqrt(asml) / s1m;  in pblueNorm()
199       return internal::sqrt(asml)/s1m;  in pblueNorm()
204     return internal::sqrt(amed);  in pblueNorm()
211     return abig * internal::sqrt(Scalar(1) + internal::abs2(asml/abig));  in pblueNorm()
237   double yf = basef * internal::abs(internal::random<double>());  in check_accuracy()
238   double yd = based * internal::abs(internal::random<double>());  in check_accuracy()
242 …std::cout << "reference\t" << internal::sqrt(double(s))*yf << "\t" << internal::sqrt(double(s))*yd…  in check_accuracy()
258 …vf[i] = internal::abs(internal::random<double>()) * std::pow(double(10), internal::random<int>(ef0…  in check_accuracy_var()
259 …vd[i] = internal::abs(internal::random<double>()) * std::pow(double(10), internal::random<int>(ed0…  in check_accuracy_var()
276   double y = 1.1345743233455785456788e12 * internal::random<double>();  in main()