Compute geometry of bg when sd is replaced

[jhabriel]

Howdy!

I'm currently working on adaptive refinement and my first take on updating the mdg is to use the replace_subdomains_and_interfaces method. The method works as intended, but the boundary grid associated to the new_sd does not have its geometry computed.

This was a bit strange to me, since the geometry is almost immediately used to update the boundary conditions of a typical PorePy model. Is there a reason behind this design choice?

All the best!

[IvarStefansson]

Well, I guess it makes sense the boundary grids are not updated, given the method name. But I see no immediate reason why we couldn't rename it and include the boundary grid update. But the update might entail some additional mapping of stored values.

[keileg]

This is not by design. Most likely we forgot to update the functionality for grid replacement when the boundary grids were introduced. AFAIK, there are no tests where we first replace a grid and then actually use the grid for anything, so it makes sense this has been undetected until now.

[jhabriel]

@IvarStefansson: I believe the old-to-new boundary data is actually mapped. It is just that the geometry of the new boundary grid is not computed.

@keileg: I see. Yes, it does make sense.

Would you be interested in some unit tests for mdg.replace_subdomain_and_interfaces())? Something like replacing the grid and then check for expected geometry / projection operators to be in place?

[keileg]

Would you be interested in some unit tests for mdg.replace_subdomain_and_interfaces())? Something like replacing the grid and then check for expected geometry / projection operators to be in place?

We would be very grateful.

Where to put these is not immediately clear: There already are some tests that use refinement, but their focus on the update of the mortar grid, hence the placement in that particular file. Anyhow, if you are willing to write some tests, we will take care of matters such as location.

[IvarStefansson]

@IvarStefansson: I believe the old-to-new boundary data is actually mapped. It is just that the geometry of the new boundary grid is not computed.

I was actually thinking of mapping of stored variables, boundary conditions etc.