← Back to team overview

yade-users team mailing list archive

Re: Contact forces

 

OK, it seems that we have the same results now!

Cheers

2010/2/12 Václav Šmilauer <eudoxos@xxxxxxxx>

>
> OK, I think I got it. How much are your numbers different? I bet they
> are all scaled down by some small value.
>
> utils.avgNumInteractions computes the average from the histogram (it
> could be computed directly like you do, but... histogram already has the
> filtering based on Aabb, so it was easier like this):
>
> def avgNumInteractions(cutoff=0.):
>        nums,counts=bodyNumInteractionsHistogram(aabbExtrema(cutoff))
>        return sum([nums[i]*counts[i] for i in
> range(len(nums))])/(1.*sum(counts))
>
> but the histogram doesn't have bin for bodies with 0 contacts.
>
> I just added that in r2020, can you check that it is correct now?
>
> Cheers, Vaclav
>
>
>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users<https://launchpad.net/%7Eyade-users>
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users<https://launchpad.net/%7Eyade-users>
> More help   : https://help.launchpad.net/ListHelp
>

References