TimeSeries.retrieve() interpolation

Hello,

according to the docs, in FEniCS 2016.2.0, TimeSeries.retrieve() can be called with an optional argument, interpolate=True/False. This does not work:

ts.retrieve(p.vector(), t, interpolate=False)

TypeError: TimeSeries_retrieve() takes no keyword arguments

edit : this works if the keyword is ommited, ts.retrieve(p.vector(), t, False).
If i read all stored timestamps from a TimeSeries with vector_times(), and then retrieve the vectors at the corresponding timestamps, no interpolation should be necessary, the values should be read directly. This is not the case in the following example:

from dolfin import *
import numpy as np

set_log_level(PROGRESS)

mesh = UnitSquareMesh(8, 8)
V = FunctionSpace(mesh, FiniteElement('P', mesh.ufl_cell(), 1))
ex = Expression("t", degree=1, t=0.)
u = Function(V)

# store data
ts = TimeSeries('test_series')
times = np.arange(0, 1, 0.1)
for t in times:
    ex.t = t
    u = interpolate(ex, V)
    ts.store(u.vector(), t)

# retrieve data
ts = TimeSeries('test_series')
times = ts.vector_times()
print(times)
for t in times:
    ts.retrieve(u.vector(), t)

The times should probably be rounded when they are compared.

Regards, David

Comments (10)