As the number of subfaults becomes large (~1000), these two script become very slow.
This is partly due to loop on subfaults.
All gmt reading should be done at once instead of trhough a loop
scipy.spatial.matrixdistance could also help to spped up the computation of the distance matrix
As the number of subfaults becomes large (~1000), these two script become very slow. This is partly due to loop on subfaults. All gmt reading should be done at once instead of trhough a loop scipy.spatial.matrixdistance could also help to spped up the computation of the distance matrix