# A simple, scalable, and robust implicit Material Point Method

In this dataset, all of the data and scripts required to generate the figures from the paper "A simple, scalable, and robust implicit Material Point Method" may be found.
Additionally, when visualisations are generated the whole set of converged load-steps are included as VTKs that may be opened in paraview.

The data is grouped into logical folders roughly correlating with the structure of the paper:
 - Column/ Contains the data to generate the stress convergence plot, and additional visualisation figure from Section 4.1
 - Beam/ Contains the data to generate the load-displacement, aggregation visualisation, residual-over-loadstep figures. Additionally a script that generates Table 1 (total iterations to solve) is included with the relevant data.
 - Elasto-plastic collapse/ Contains the data to generate the three elasto-plastic collapse example figures, this includes the visualisation of stresses, and the two residual-over-iteration plots comparing the damping strategy and fictitious mass-matrix strategy.
 - The Performance/ folder has all the data for Section 5, including strong and weak scaling tests executed on Hamilton 8 (Durham HPC resource), and representitive examples on a workstation.

There are no dependencies on the python scripts, apart from those available on pip (matplotlib, numpy, pandas, vtk).
