Janek Kozicki a écrit :
See line 342 in TesselationWrapper.cpp
Real sphereVol = 4.188790 * std::pow ( ( V_it->point().weight() ),1.5 );
// 4/3*PI*R³ = 4.188...*R³
why do you raise to power 3/2 instead of 3 ?
For us, a sphere is position and radius. In radical tesselation,
weighted points are (if you interpret them as representation of spheres)
positions and radius^2. A bit misleading sometimes.
If you insert a sphere using the TesselationWrapper layer, the "insert"
function take pos and rad (and rad is squared internaly).
If you use CGAL::Delaunay::insert straightaway, you have to provide pos
and rad^2.
In all cases accessing triangulation->vertex->weight() will return
squared radius.
what's wrong with numpy_boost<double,2> and its segfaults?
No clue.
Bruno
--
_______________
Bruno Chareyre
Associate Professor
Grenoble INP
Lab. 3SR
BP 53 - 38041, Grenoble cedex 9 - France
Tél : 33 4 56 52 86 21
Fax : 33 4 76 82 70 43
________________
_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help : https://help.launchpad.net/ListHelp
