Hi all,<br><br>I have found that the denominator of the L function implemented in the method L_Func(double r, double d, double phi) of the file itkDeformableSimplexMesh3DFilter.txx in itk 3.6.0 does not correspond with the mentioned in the paper of Delingette.<br>
Instead of implementing the denominator as:<br><br>double denom = eps * (sqrt(tmpSqr) + r);<br><br>it should be:<br><br>double denom = eps * sqrt (tmpSqr) + r;<br><br>Regards,<br><br>Paula<br>