Author: eudoxos
Date: 2008-10-27 22:41:21 +0100 (Mon, 27 Oct 2008)
New Revision: 1561
Modified:
trunk/lib/smoothing/WeightedAverage2d.hpp
Log:
1. Fix algebraic typo in max distance check.
Modified: trunk/lib/smoothing/WeightedAverage2d.hpp
===================================================================
--- trunk/lib/smoothing/WeightedAverage2d.hpp 2008-10-27 20:26:39 UTC (rev
1560)
+++ trunk/lib/smoothing/WeightedAverage2d.hpp 2008-10-27 21:41:21 UTC (rev
1561)
@@ -119,8 +119,8 @@
virtual Real getWeight(const Vector2r& meanPt, const T& e){
Vector2r pos=getPosition(e);
Real rSq=pow(meanPt[0]-pos[0],2)+pow(meanPt[1]-pos[1],2);
- if(rSq>relThreshold*stDev*stDev) return 0.;
- return (1./sqrt(2*M_PI*stDev*stDev))*exp(-rSq/(2*stDev*stDev));
+ if(rSq>pow(relThreshold*stDev,2)) return 0.; // discard points
further than relThreshold*stDev, by default 3*stDev
+ return (1./(stDev*sqrt(2*M_PI)))*exp(-rSq/(2*stDev*stDev));
}
virtual vector<Vector2i> filterCells(const Vector2r& refPt){return
WeightedAverage<T,Tvalue>::grid->circleFilter(refPt,stDev*relThreshold);}
};
_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to : [EMAIL PROTECTED]
Unsubscribe : https://launchpad.net/~yade-dev
More help : https://help.launchpad.net/ListHelp
_______________________________________________
yade-dev mailing list
yade-dev@lists.berlios.de
https://lists.berlios.de/mailman/listinfo/yade-dev
