Hi folks, <div><br></div><div>I&#39;m working on the integration of a function like:</div><div><br></div><div>K(r,t) = 1/(2piDr)exp[-r^2/Dt]                       [1]</div><div><br></div><div>over Voronoi cells (r is the distance from the point at which is associated the cell). </div>
<div><br></div><div>I googled a lot and I found this two useful hints: </div><div><br></div><div><a href="http://stackoverflow.com/questions/5941113/looking-for-python-package-for-numerical-integration-over-a-tessellated-domain">http://stackoverflow.com/questions/5941113/looking-for-python-package-for-numerical-integration-over-a-tessellated-domain</a></div>
<div><br></div><div><a href="http://mathforum.org/kb/message.jspa?messageID=4963570&amp;tstart=0">http://mathforum.org/kb/message.jspa?messageID=4963570&amp;tstart=0</a></div><div><br></div><div>but I&#39;m still not able to understand how I should do this integration.</div>
<div><div>I have the function that returns the segments ( list of   [(x_star,y_start),(x_sop,y,stop)] ) necessary to construct the Voronoi cells</div><div>associated to a set of points. </div></div><div><br></div><div>Could someone suggest how to proceed ? </div>
<div><br></div><div>There&#39;s also a numerical problem connected with the integration due to the singularity in r==0.</div><div><br></div><div>Could you suggest which is a reasonable stable integration method available in scipy </div>
<div>that could handle the function [1]?</div><div><br></div><div><br></div><div>Cheers, </div><div>Eraldo </div><div><br></div><div><br></div><div> </div><div><br></div><div> </div><div><br></div><div> </div><div><br></div>