yade-users team mailing list archive

[Robert Caulk <question675926@xxxxxxxxxxxxxxxxxxxxx>]

Question #675926 on Yade changed:
https://answers.launchpad.net/yade/+question/675926

    Status: Open => Answered

Robert Caulk proposed the following answer:
Hello,

I guess you want something like this:

def getRadiusArray():
    radii=[]
    for b in O.bodies:
	if isinstance(b.shape,Sphere):
        	radii.append(b.shape.radius)
   
getRadiusArray()
values, base = np.histogram(data, bins=40)
cumulative = np.array(np.cumsum(values),dtype=float)
perc = cumulative/max(cumulative)*100
plt.semilogx(base[:-1], cumulative, c='blue')

Cheers,

Robert

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.