Hi, I am currently learning FEniCS (but using 2019.1.0, the legacy release) and have been looking at the advection-diffusion reaction equations example, which uses a TimeSeries to retrieve the velocities previously solved at each time step in the navier stokes solution, like this:
timeseries_u2.retrieve(u2.vector(), t)
Running problems like these in parallel, e.g. mpirun -np 4 … .py gives an error at the TimeSeries retrieval, as shown below:
*** Error: Unable to successfully call PETSc function 'VecAXPY'.
*** Reason: PETSc error code is: 75 (Arguments are incompatible).
*** Where: This error was encountered inside /home/conda/feedstock_root/build_artifacts/fenics-pkgs_1662986696266/work/dolfin/dolfin/la/PETScVector.cpp.
So I was wondering if anyone had an idea about an alternative/similar object to work with parallel retrieval of past solutions. Is there something similar to TimeSeries that I can use for parallel running of the system of advection-diffusion equations? Or should I look into other parallel array storing methods?
Thanks very much for reading.