This repository provides an implementation of the multipoint flux approximation method in 2d and 3d. The core implementation is in C++, with python bindings provided for a limited set of the functionality. The primary motivation is to provide accelerated discretizations for PorePy, with suitable wrappers under implementation.
The software provided here, both the C++ code and Python wrappers, should be considered immature and be used with some caution.
Currently, the best (least bad) starting point for the software is likely the PorePy wrappers, see
here, in particular
the file test_fv.py.
See Licence.