I want to use the RBF module in bitpit to achieve mesh morphing. And I modified "/home/hzy/bitpit-bitpit-1.9.0/test/integration_tests/RBF/test_RBF_0004.cpp" and preliminary achieved RBF interpolation of displacement field. But when the number of control points is larger than 10000, the solver is very slow, about 20 min. Can we use MPI paralle to speed up the solving procedure. perhaps split the control points into much small sub sets or just use MPI to solve the large linear algebra system? Does anybody can provide an examples? Thanks!