I am trying to run a FSI simulation using FEniCS for the solid participant. My code is inspired from the perpendicular flap tutorial (Perpendicular flap | preCICE - The Coupling Library). In this example, only one field (u_n) has to be stored in checkpoints (precice.store_checkpoint(u_n, t, n) and u_cp, t_cp, n_cp = precice.retrieve_checkpoint()).
The issue / feature request that you are mentioning is unfortunately still open. At the moment you can only store a single field in the adapter. I currently do not have the time to dig deep into the implementation of the fenics-adapter and develop this feature thoroughly (including tests etc.), but you could try to start from the error that you have shown above and implement the feature on your own (A pull request would be highly appreciated! I can also help here.). The problematic line is actually this one: fenics-adapter/fenicsprecice/fenicsprecice.py at 6c8d7be8153b8bc2c573a23985dfd3836dd1a8fb · precice/fenics-adapter · GitHub. But there are some more places that have to be changed down the road. But it should not be too hard.
Another possibility: You can also just implement your own checkpointing (independent from the adapter) in your solver instead of using precice.store_checkpoing and precice.retreive_checkpoint. In the end you just need a place where you can be sure that state = (u_old, sig_old, p) is safe. At the moment the adapter offers this for you - but only for a single field.
Be careful: FEniCS sometimes updates variables in a pointer-like fashion. This was the main motivation to provide the store_checkpoint and retreive_checkpoint functions: That the user does not have to worry about these technicalities. It’s not too complicated (compared to FEM it’s pretty much a no-brainer), but these issues are sometimes hard to debug and annoying. Inside of store_checkpoint I, therefore, explicitly create a copy (user_u.copy()) to make sure FEniCS cannot mess around with the checkpointed values.
Note on preCICE v2.5.0: The fenics-adapter is not maintained anymore for preCICE v2. Is it possible for you to upgrade to preCICE v3?
I think I’ll go for the second option since I’m almost at the end of my PhD.
Unfortunately I also no longer have the time to update my code for the new release of preCICE v3, mainly for the in-house fluid solver where I had to develop my own adapter. But I’m going to talk about it in my research team for future work.
requires to modify the adapter to take into account tuples, right?
For option two, I need to call markActionFulfilled. I can find this part in the fenicsx adapter (e.g. self._interface.mark_action_fulfilled(self.action_write_iteration_checkpoint())) or dealii adapter (e.g. precice.markActionFulfilled(precice::constants::actionWriteIterationCheckpoint());) for example, but not in the fenics adapter (just after self._checkpoint = SolverState(my_u, t, n)). How this action is managed in the fenics adapter?
that’s one of the preCICE v2 / v3 breaking changes: We removed markActionFulfilled in preCICE v3. Meaning: If you are using v2, just follow the approach from the fenicsxadapter (which has not been updated to v3 yet). Do the following:
if precice.is_action.required(precice.action_write_iteration_checkpoint()):
# do your own checkpointing
precice.mark_action_fulfilled(precice.action_write_iteration_checkpoint())
(I hope the syntax is fine. I’ve become a bit rusty with preCICE v2, but the python bindings help(...) and documentation should get you on the right path)
See our porting guide. It should also help you backporting from v3 examples to v2. Not the way we intended it to be used, but if it is useful for you I’m happy
