yade-users team mailing list archive

[AnnaMaria <question685862@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

AnnaMaria is still having a problem:
Jan Thank you very much for your effort!!!!
I really appreciate it.
I have only changed  elatTot / nTot to elat_avg=elatTot / nTot
And then return  elat_avg # average elat so that I can plot.

def lateral():
    elatTot = 0.0 # sum of elat of all suitable particles
    nTot = 0 # number of suitable particles
    for p in O.bodies:
	x=p.state.refPos[0]
	z=p.state.refPos[2]
	d=sqrt(pow(x,2)+pow(z,2))
        if d>r1 and d<r2:
      	    xnew=p.state.pos[0]
            znew=p.state.pos[2]
            dnew=sqrt(pow(xnew,2)+pow(znew,2))
            elat=(d-dnew)/d
            elatTot += elat
            nTot += 1
	    elat_avg=elatTot / nTot
    return  elat_avg # average elat

def recorder():
	lateral()
	yade.plot.addData({'i':O.iter,'eps':abs(strainer.strain),'elateral':(elat_avg),'sigma':abs(strainer.avgStress)})
	plot.saveDataTxt(OUT)

#### Plot During Simulation
plot.plots={'elateral':('sigma')}
plot.plot()

I also changed the recorder part. But when I run the code, I get the following Error:
File "<string>", line 1, in <module>
File "UCS.py", line 87, in recorder
yade.plot.addData({'i':O.iter,'eps':abs(strainer.strain),'elateral':(elat_avg),'sigma':abs(strainer.avgStress)})
NameError: global name 'elat_avg' is not defined

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