Wavelet-based mesh refinement

[zgxn]

Motivation

Add a mesh refinement strategy that is based on multiresolution analysis of arbitrary input data. The method is described in (Ozgen-Xian et al., 2020).

New Files

• wavelet_analysis.py: Utilities to carry out a Discrete Haar-Wavelet Transform, which is the core of the meshing strategy
• examples/WaveletExample.ipynb: Jupyter notebook that illustrates the usage of the added methodology

Modified Files

• dem_class.py: Added the function build_wavelet_triplane(), which is the high-level function to call from an external script
• meshing.py: Added the function _wavelet_surface_mesh(), which calls the methods from wavelet_analysis.py to generate the multiresolution surface mesh

Testing

Run the Jupyter notebook examples/WaveletExample.ipynb to verify that everything works as intended.

References

Özgen-Xian, I., Kesserwani, G., Caviedes-Voullième, D., Molins, S., Xu, Z., Dwivedi, D., Moulton, J. D., Steefel, C. I. (2020) Wavelet-based local mesh refinement for rainfall-runoff simulations, Journal of Hydroinformatics, 22, 1059-1077. doi: 10.2166/hydro.2020.198

[daniellivingston]

This is phenomenal. Thanks Ilhan!