Problem with plotting a solution of Stokes Problem

1

Hello everyone,

I am having a problem with plotting, as in I am not getting a heat map of the solution of the stokes problem, but instead some sort of vector field. As a background I am solving stokes problem on a irregular mesh. When I use dl.plot , I obtain the vector field, I found that if I used instead vtkplotter.dolfin, I get what, however after I plot it and see the plot I cannot do any more plots, I get the plot staying until I force quit it. I tried reinstalling the Fenics and vtkplotter.dolfin multiple times and the problem is still there. I was hoping someone can help me item use the dl.plot (dolfin) or the issue with the vtkplotter plot resolved

import dolfin as dl
from vtkplotter.dolfin import plot
import matplotlib.pyplot as plt
u0 = SA.solve_fwd_Stokes() %u0 here is a dolfin.cpp.la.PETScVector
fig = plt.figure(1)
p1 = dl.plot(u1f.sub(0))
plt.colorbar(p1)
plt.show()![Figure_0|666x500](upload://tKJXMXvrVbae9kU1GVSY3lSd3Th.png)


u0f = convert_to_function(u0) %dolfin.cpp.la.PETScVector to function 
plot(u0f)
![Screen Shot 2021-10-26 at 10.02.10 PM|690x424](upload://drecHIYGq9a9TyPUnTkZjptDxGx.png)

Currently I am running Fenics 2019.1 and the version on the 2020.3.1 . I am running macOS Big Sur (11.6). Thank you so much for the time and consideration! I am attaching the plots from the dl.plot and the one from vtkplotter

2

Your images have not been attached as you put them inside the code formatting.

I would suggest you use Paraview for post processing, by saving your function to either “xdmf” or “pvd”.