← Back to team overview

yade-dev team mailing list archive

Re: TesselationWrapper stuff

 

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
________________




References