I have a 2D finite difference scheme to implement. I use C++ but do not want to write my own 2D grid as I would like to reuse the code for many other problems. Hence, I want to use a well-established code that provides all the functionalities that I may need. Deall ii fits this purpose but since deal is meant for FEM, I want to know if it is advisable to also use it for FDM.
Your thoughts are highly welcomed.