[Users] Plotting all refinement levels using Kuibit

José Ferreira jpmferreira at ua.pt
Tue Oct 8 12:05:24 CDT 2024


Dear toolkit community,


I’m struggling to find a way to plot a 1D quantity in Kuibit without 
merging the different meshes originated from the mesh refinement.

As a workaround, I have made a Python script that reads the file of a 
grid function and plots every point along the axis and its corresponding 
value. I can share this script if somebody is interested, but it is not 
very good as it breaks once you use checkpoints and has hardcoded paths.

However, it gets the job done, as you can see in an example plot below, 
where it is hopefully obvious that there is an increase in the number of 
points as z approaches 0.

Obviously, a solution that uses Kuibit would be much better. However, I 
cannot come up with such a solution and so I’m writing this e-mail.


 From the documentation, the object that represents the simulation data 
(including the mesh refinement levels), is the |HierarchicalGridData|, 
and I know how to go from a grid function to this object.

For instance, to get the |HierarchicalGridData| that represents the real 
part of a scalar field along the x-axis at t=0, we do

|sd = SimDir("<output>") gf = sd.gf phi1=gf.x["phi1"].get_time(0) |

All I have to do now is plot this somehow. As it turns out, this is 
impossible, because according to the documentation

    We cannot plot directly this object, because it is a complicated
    object. To plot it, we have to merge the refinement levels to a
    single |UniformGridData|. We can do this with
    |refinement_levels_merged()|.

But, all of the information regarding the mesh refinement levels is 
here, so it should be possible to extract this information and plot it, 
right?

Looking at the methods and variables in |HierarchicalGridData| I came up 
with nothing.

Does anybody have a solution to this?


Thank you in advance,

Best Regards,

José Ferreira


​
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20241008/098dcead/attachment.htm>


More information about the Users mailing list