Hi Robert,
thank you for the help.
The interplation was working according to the examples from tests.
This is a part of the code from my linear elastic problem, the result extraction (displacements, stress and strain) and the additional interpolation of certain points at the end:
# field and field variables
field = Field.from_args('displacement', numpy.float64, 3, omega,
approx_order=1)
u = FieldVariable('u', 'unknown', field, order=0)
v = FieldVariable('v', 'test', field, primary_var_name='u')
...
# run simulation
vec = pb.solve()
# postprocessing and saving results
output_dict = vec.create_output_dict(fill_value=None,
extend=True,
linearization=None) # u als key
# struct array displacements, strain, stress
output_dict = post_process(output_dict, pb, state, extend=True)
nodes_displacements = output_dict['u'].data
strain_tensors = output_dict['cauchy_strain'].data
stress_tensors = output_dict['cauchy_stress'].data
...
p1 = numpy.array([1120, 53.375, 60.804])
coors = numpy.array([p1])
interpolations = u.evaluate_at(coors)
My question is: Is it possible to run the point interpolation after everything is finished and closed just by loading the saved result data. What I mean is, I run the FE simulation in a first step and I save
the results (displacements, stress, strain). Later, I would like to load the saved data and run the point interpolation. As far as I understand the code, I would need the
FieldVariable and the Field, but also the results data of the problem. Could you give some advice what would be necessary for saving to reload the required data structures for the interpolation?
Thank you,
Kathrin