Hi,
I am a new user of SfePy and finite element programming in general, and I have started my own finite element project last week. After comparing SfePy with other packages like FEniCS and GetFEM++ I have decided to give SfePy a try. My problem is close to a diffusion problem. I have a fair amount of scientific Python experience in general, including some VTK.
I know I can export the computation result to VTK files in Unstructured Grid format, which is performed at the end of many examples. My question is: how can I access the final result, at each vertex, directly from the Python program itself, saving the result into say a Numpy array or equivalent data structure.
For example, the
Time Thermal Interactive example exports domain.01.vtk to domain.02.vtk, containing points, cells, and temperature at each point ("SCALARS T float 1.") Is there a way to access these numbers directly at the end of the computation, from say the Problem variable pb or State variable state0? I feel like this is probably an easy task but I could not figure out on my own just by looking into the State and Problem variables.
I need these values directly in Python as I have other part of my model that takes the temperature values at these vertices. Since the other computation will use the same mesh I would like to have the mesh and all solved values.