@@ -395,9 +395,13 @@ optional< Rational< typename K::FT> > compute_normal_offset_lines_isec_timeC2 (
395395 +(l1->b ()*l0->a ()*l2->c ())
396396 -(l0->b ()*l1->a ()*l2->c ());
397397
398- FT r0 = CGAL_SS_i :: inexact_sqrt (square (l0->a ())+square (l0->b ()));
399- FT r1 = CGAL_SS_i :: inexact_sqrt (square (l1->a ())+square (l1->b ()));
400- FT r2 = CGAL_SS_i :: inexact_sqrt (square (l2->a ())+square (l2->b ()));
398+ FT sum_sq_0 = square (l0->a ()) + square (l0->b ());
399+ FT sum_sq_1 = square (l1->a ()) + square (l1->b ());
400+ FT sum_sq_2 = square (l2->a ()) + square (l2->b ());
401+
402+ FT r0 = CGAL_SS_i::inexact_sqrt (sum_sq_0);
403+ FT r1 = CGAL_SS_i::inexact_sqrt (sum_sq_1);
404+ FT r2 = CGAL_SS_i::inexact_sqrt (sum_sq_2);
401405
402406 den = r0 * (l2->b ()*l1->a () - l2->a ()*l1->b ()) +
403407 r1 * (l2->a ()*l0->b () - l2->b ()*l0->a ()) +
@@ -657,9 +661,13 @@ optional< Point_2<K> > construct_normal_offset_lines_isecC2 ( intrusive_ptr< Tri
657661
658662 if ( l0 && l1 && l2 )
659663 {
660- FT r0 = CGAL_SS_i :: inexact_sqrt (square (l0->a ())+square (l0->b ()));
661- FT r1 = CGAL_SS_i :: inexact_sqrt (square (l1->a ())+square (l1->b ()));
662- FT r2 = CGAL_SS_i :: inexact_sqrt (square (l2->a ())+square (l2->b ()));
664+ FT sum_sq_0 = square (l0->a ()) + square (l0->b ());
665+ FT sum_sq_1 = square (l1->a ()) + square (l1->b ());
666+ FT sum_sq_2 = square (l2->a ()) + square (l2->b ());
667+
668+ FT r0 = CGAL_SS_i::inexact_sqrt (sum_sq_0);
669+ FT r1 = CGAL_SS_i::inexact_sqrt (sum_sq_1);
670+ FT r2 = CGAL_SS_i::inexact_sqrt (sum_sq_2);
663671
664672 FT den = r0 * ( l2->a ()*l1->b () - l1->a ()*l2->b ()) +
665673 r1 * ( l0->a ()*l2->b () - l0->b ()*l2->a ()) +
0 commit comments