Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] Problem calculating the volume of alpha_shape_3_cells

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] Problem calculating the volume of alpha_shape_3_cells


Chronological Thread 
  • From: Acecross <>
  • To:
  • Subject: Re: [cgal-discuss] Problem calculating the volume of alpha_shape_3_cells
  • Date: Thu, 14 Dec 2017 01:56:56 -0700 (MST)
  • Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None ; spf=Fail ; spf=Pass
  • Ironport-phdr: 9a23:kGcxPhTFYwZl9MmpZXy4yS+ZO9psv+yvbD5Q0YIujvd0So/mwa68bRON2/xhgRfzUJnB7Loc0qyK6/mmATRIyK3CmUhKSIZLWR4BhJdetC0bK+nBN3fGKuX3ZTcxBsVIWQwt1Xi6NU9IBJS2PAWK8TW94jEIBxrwKxd+KPjrFY7OlcS30P2594HObwlSizexfa5+IA+qoQnNq8IbnZZsJqEtxxXTv3BGYf5WxWRmJVKSmxbz+MK994N9/ipTpvws6ddOXb31cKokQ7NYCi8mM30u683wqRbDVwqP6WACXWgQjxFFHhLK7BD+Xpf2ryv6qu9w0zSUMMHqUbw5Xymp4qF2QxHqlSgHLSY0/m/XhMJuj6xUohyhqQFwzIHIb4+YL+Z+c6HBcd8GWWZNQthdWipcCY28dYsPCO8BMP5EoobhuVQBsQG+BQ63BOP01zRFgWL73akg3OQ8EgHJxhYgEMgSsHnPotT6KrwdXfqow6bW1DjOae5d1zn66IjNaB8hoPeMUKp2ccXL0kYgDAfFgUuKpYP5ODOV0+IAuHWY4ep4Te+jlm8qpx9rrjSxxsohhJPFip8bx13K7yl0wpg5KNulQ0Bhe9GkCoFftySCOot2XMwiR2ZotT48yrAJo5K7ejIFx488yB7bdfOLaZSH4hXmVOqJIDd4gmxqeK6nihqv8EWs1PfwWte03VtEtCZJj9fBu34X2xHR9MSLUv598V2g2TaL2QDT8OZEIUUsmKfGL54u2KMwlpkJvkTZBCD2hF77jaGIeUo/4OSn9vrobq38ppCAL490lh3+MqM2l8OjDuQ3KAwOU3GG9uS9z73s4VD5QK5RjvAtianYsJXaJdwBqaKjAg9V1Jwj6xelADu83tQYhypPEFUQcx2Oi83lOkrFPevjJfa5mVWl1jlxlN7cObi0K57AJWLEi6ypKZx54lNfwRA60fhA6pgSBrxXc6G7YVP4qNGNVkxxCAez2euyUNg=

Works perfectly fine. Thank you very much. New code:
double get_alpha_3_volume(const Alpha_shape3 &A,
Alpha_shape3::Classification_type c) {
std::vector<Alpha_shape3::Cell_handle> cells;
A.get_alpha_shape_cells(std::back_inserter(cells), c);
double volume = 0;
for (std::vector<Alpha_shape3::Cell_handle>::iterator it =
cells.begin();
it != cells.end(); it++) {
Point3 p1 = (*it)->vertex(0)->point();
Point3 p2 = (*it)->vertex(1)->point();
Point3 p3 = (*it)->vertex(2)->point();
Point3 p4 = (*it)->vertex(3)->point();
Alpha_shape3::Tetrahedron t(p1,p2,p3,p4);
volume += t.volume();
}
...



--
Sent from: http://cgal-discuss.949826.n4.nabble.com/



Archive powered by MHonArc 2.6.18.

Top of Page