songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Use cross product instead of Plane_Plane intersection line

This commit is contained in:
Giles Bathgate 2022-04-02 18:10:06 +01:00
parent 3d44e3f14d
commit 9cdf8ba469
1 changed files with 2 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
Nef_S2/include/CGAL/Nef_S2/Sphere_circle.h
View File

@ -173,11 +173,9 @@ Sphere_point<R> intersection(const Sphere_circle<R>& c1,
|c1| and |c2|. \precond |c1 != c2| as sets.}*/ |c1| and |c2|. \precond |c1 != c2| as sets.}*/
{ {
CGAL_assertion(!equal_as_sets(c1,c2)); CGAL_assertion(!equal_as_sets(c1,c2));
typename R::Line_3 lres;
CGAL_NEF_TRACEN("circle_intersection "<<c1<<" "<<c2); CGAL_NEF_TRACEN("circle_intersection "<<c1<<" "<<c2);
CGAL::Object o = CGAL::intersection(c1.plane(),c2.plane()); return R().construct_cross_product_vector_3_object()(
if ( !CGAL::assign(lres,o) ) CGAL_error(); c1.orthogonal_vector(), c2.orthogonal_vector());
return CGAL::ORIGIN + lres.direction().vector();
} }