yade-users team mailing list archive

[Jan Stránský <question706106@xxxxxxxxxxxxxxxxxxxxx>]

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

Jan Stránský proposed the following answer:
I think there more options how to do it.
Below is one example.
Adjust it according to your needs - the information provided is still not very informative (e.g. we can only guess if "i" is O.iter or not, what "t1" and "t2" might be etc.)
The example blow is based on O.time

#### saved.txt
# time value
1  1.0
2  0.6
3  0.4
4  0.3
5  0.4
6  0.6
7  0.7
8  0.6
9  0.4
10 0.3
11 0.1
12 0.2
13 0.4
14 0.7
15 0.8
16 1.0
17 0.9
18 0.7
19 0.4
20 0.1
####

#### script.py
from yade import plot
import numpy

saved = numpy.loadtxt("saved.txt") # load saved data
savedTimes = saved[:,0]  # 1st comlumn
savedValues = saved[:,1] # 2nd comlumn

def plotAddData():
    time = O.time
    # simulation vlaue
    value = 0.5 + 0.5 * sin(time)
    # interpolated saved value
    savedValue = numpy.interp(time,savedTimes,savedValues)
    plot.addData(
        time = time,
        value = value,
        savedValue = savedValue,
    )

O.engines = [
    PyRunner(virtPeriod=0.2, command="plotAddData()"),
]

plot.plots = {"time":("value","savedValue")}
plot.plot()

O.stopAtTime = 20

O.dt = 1e-6
O.run()
####

Cheers
Jan

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