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 ? what's wrong with numpy_boost<double,2> and its segfaults? -- Janek Kozicki http://janek.kozicki.pl/ | _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
