GridClean is called in FInterp at line 143 with
_braid_GridClean(core, grids[level]);
and this deletes the state vector at that level. This seems (possibly) counter to the storage options in the code that stipulate that we preserve storage on coarse levels (if the "int storage" is the right value).
Do we want to only call GridClean if storage == -1 or storage < level ? Or perhaps put this logic inside of GridClean?
This came up in TorchBraid, where we set storage==0 and then wanted access to the coarse ua vectors after braid_Drive was finished.
What do you think Rob @rfalgout ?
GridClean is called in FInterp at line 143 with
_braid_GridClean(core, grids[level]);
and this deletes the state vector at that level. This seems (possibly) counter to the storage options in the code that stipulate that we preserve storage on coarse levels (if the "int storage" is the right value).
Do we want to only call GridClean if storage == -1 or storage < level ? Or perhaps put this logic inside of GridClean?
This came up in TorchBraid, where we set storage==0 and then wanted access to the coarse ua vectors after braid_Drive was finished.
What do you think Rob @rfalgout ?