31   typedef typename R::RT RT;  in equal_xy()  typedef
34   const RT& phw = p.hw();  in equal_xy()
35   const RT& qhw = q.hw();  in equal_xy()
55 template <class RT>
57 compare_power_distanceH2(const RT& phx, const RT& phy, const RT& phw,  in compare_power_distanceH2()
58                          const Quotient<RT>& pwt,  in compare_power_distanceH2()
59                          const RT& qhx, const RT& qhy, const RT& qhw,  in compare_power_distanceH2()
60                          const Quotient<RT>& qwt,  in compare_power_distanceH2()
61                          const RT& rhx, const RT& rhy, const RT& rhw)  in compare_power_distanceH2()
64   RT dphx = rhx * phw - phx * rhw;  in compare_power_distanceH2()
65   RT dphy = rhy * phw - phy * rhw;  in compare_power_distanceH2()
66   RT dqhx = rhx * qhw - qhx * rhw;  in compare_power_distanceH2()
67   RT dqhy = rhy * qhw - qhy * rhw;  in compare_power_distanceH2()
68   RT dphw = CGAL_NTS square(phw);  in compare_power_distanceH2()
69   RT dqhw = CGAL_NTS square(qhw);  in compare_power_distanceH2()
70   RT drhw = CGAL_NTS square(rhw);  in compare_power_distanceH2()
72   RT npwt = pwt.numerator();  in compare_power_distanceH2()
73   RT dpwt = pwt.denominator();  in compare_power_distanceH2()
74   RT nqwt = qwt.numerator();  in compare_power_distanceH2()
75   RT dqwt = qwt.denominator();  in compare_power_distanceH2()
78   RT dh1 = (CGAL_NTS square(dphx) + CGAL_NTS square(dphy))*dpwt - npwt * dphw * drhw;  in compare_power_distanceH2()
79   RT dh2 = (CGAL_NTS square(dqhx) + CGAL_NTS square(dqhy))*dqwt - nqwt * dqhw * drhw;  in compare_power_distanceH2()
84 template <class RT>
86 power_testH2( const RT &phx, const RT &phy, const RT &phw, const Quotient<RT> &pwt,  in power_testH2()
87               const RT &qhx, const RT &qhy, const RT &qhw, const Quotient<RT> &qwt,  in power_testH2()
88               const RT &rhx, const RT &rhy, const RT &rhw, const Quotient<RT> &rwt,  in power_testH2()
89               const RT &thx, const RT &thy, const RT &thw, const Quotient<RT> &twt)  in power_testH2()
91     RT npwt = pwt.numerator();  in power_testH2()
92     RT dpwt = pwt.denominator();  in power_testH2()
93     RT nqwt = qwt.numerator();  in power_testH2()
94     RT dqwt = qwt.denominator();  in power_testH2()
95     RT nrwt = rwt.numerator();  in power_testH2()
96     RT drwt = rwt.denominator();  in power_testH2()
97     RT ntwt = twt.numerator();  in power_testH2()
98     RT dtwt = twt.denominator();  in power_testH2()
100     RT dphx = phx*phw;  in power_testH2()
101     RT dphy = phy*phw;  in power_testH2()
102     RT dphw = CGAL_NTS square(phw);  in power_testH2()
103     RT dpz = (CGAL_NTS square(phx) + CGAL_NTS square(phy))*dpwt  - npwt*dphw;  in power_testH2()
105     RT dqhx = qhx*qhw;  in power_testH2()
106     RT dqhy = qhy*qhw;  in power_testH2()
107     RT dqhw = CGAL_NTS square(qhw);  in power_testH2()
108     RT dqz = (CGAL_NTS square(qhx) + CGAL_NTS square(qhy))*dqwt - nqwt*dqhw;  in power_testH2()
110     RT drhx = rhx*rhw;  in power_testH2()
111     RT drhy = rhy*rhw;  in power_testH2()
112     RT drhw = CGAL_NTS square(rhw);  in power_testH2()
113     RT drz = (CGAL_NTS square(rhx) + CGAL_NTS square(rhy)) *drwt - nrwt*drhw;  in power_testH2()
115     RT dthx = thx*thw;  in power_testH2()
116     RT dthy = thy*thw;  in power_testH2()
117     RT dthw = CGAL_NTS square(thw);  in power_testH2()
118     RT dtz = (CGAL_NTS square(thx) + CGAL_NTS square(thy))*dtwt - ntwt*dthw;  in power_testH2()
131 template <class RT>
133 power_testH2( const RT &phx, const RT &phy, const RT &phw, const Quotient<RT> &pwt,  in power_testH2()
134               const RT &qhx, const RT &qhy, const RT &qhw, const Quotient<RT> &qwt,  in power_testH2()
135               const RT &thx, const RT &thy, const RT &thw, const Quotient<RT> &twt)  in power_testH2()
137     RT npwt = pwt.numerator();  in power_testH2()
138     RT dpwt = pwt.denominator();  in power_testH2()
139     RT nqwt = qwt.numerator();  in power_testH2()
140     RT dqwt = qwt.denominator();  in power_testH2()
141     RT ntwt = twt.numerator();  in power_testH2()
142     RT dtwt = twt.denominator();  in power_testH2()
146     RT pa, qa, ta;  in power_testH2()
161     RT dphw = CGAL_NTS square(phw);  in power_testH2()
162     RT dpz = (CGAL_NTS square(phx) + CGAL_NTS square(phy))*dpwt - npwt*dphw;  in power_testH2()
164     RT dqhw = CGAL_NTS square(qhw);  in power_testH2()
165     RT dqz = (CGAL_NTS square(qhx) + CGAL_NTS square(qhy))*dqwt - nqwt*dqhw;  in power_testH2()
167     RT dthw = CGAL_NTS square(thw);  in power_testH2()
168     RT dtz = (CGAL_NTS square(thx) + CGAL_NTS square(thy))*dtwt - ntwt*dthw;  in power_testH2()