Dear all,
I'm trying to extend our analysis thorn to tackle multipatch grids built with Llama. I want to use Coordinates::volume_form in the computation of volume integrals, by setting the corresponding parameter Coordinates::store_volume_form = yes. Our standard use case relies on Coordinates::coordinate_system = "Thornburg04", for which it works like a charm. However, except for this system and the analogous Thornburg13 (I have not tested it so far), it seems the variable Coordinates::volume_form_state is not modified nor initialized, whatever the value of the parameter Coordinates::store_volume_form. That makes the parameter insufficient by itself to perform checks. Moreover, since the flag is not initialized, we also can't simply check its value, which can be polluted. I have tried to use Carpet::poison_new_timelevels = yes and CarpetLib::poison_new_memory = yes, but it didn't seem to solve the issue. Hence, I would like to ask if there's another way to check for uninitialized values, a safe way to use the variables I mentioned in a general case, or any information/subtlety I may have missed? Ultimately, I can directly check Coordinates::coordinate_system, and output a warning/error if it's not of the types mentioned above. That feels a bit unsatisfactory, although I may default to it. Thank you very much for any help or insight you can provide! Best, Jordan Nicoules
_______________________________________________ Users mailing list [email protected] http://lists.einsteintoolkit.org/mailman/listinfo/users
