Hi Anton,

Just in case: you know that SpherePack has a psd() function already?

Maybe these lines (partly stolen from test/psd.py) would do what you want:
sp0=pack.SpherePack()
sp0.fromSimulation()
pylab.semilogx(*sp0.psd(bins=30,mass=True),label='Mass PSD of (free) %d
random spheres'%len(sp0))

Bruno

On 19/07/11 16:54, [email protected] wrote:
> ------------------------------------------------------------
> revno: 2891
> committer: Anton Gladky <[email protected]>
> branch nick: yade
> timestamp: Tue 2011-07-19 16:52:47 +0200
> message:
>   Add utils.psd(). First try
> modified:
>   py/utils.py
>
>
> --
> lp:yade
> https://code.launchpad.net/~yade-dev/yade/trunk
>
> Your team Yade developers is subscribed to branch lp:yade.
> To unsubscribe from this branch go to 
> https://code.launchpad.net/~yade-dev/yade/trunk/+edit-subscription
>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-dev
> Post to     : [email protected]
> Unsubscribe : https://launchpad.net/~yade-dev
> More help   : https://help.launchpad.net/ListHelp


-- 
_______________
Bruno Chareyre
Associate Professor
ENSE³ - Grenoble INP
Lab. 3SR
BP 53 - 38041, Grenoble cedex 9 - France
Tél : +33 4 56 52 86 21
Fax : +33 4 76 82 70 43
________________

_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to     : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp

Reply via email to