The below animations demonstrate how an LU decomposition with partial pivoting proceeds using PLAPACK (frames are time-scaled).

The animations show a 3D plot of a 20x20 matrix fed to PLAPACK from Mathematica through Greg Morrow's PMI (PLAPACK Mathpackage Interface) engine. This engine allows the user to carry out calculations on a parallel machine through the Mathematica (or HiQ/LabView/MATLAB) interface. The animations themselves was created by John Gunnels.

It should be noted that Mathematica was not used for any calculations, but was used to generate the initial matrix and render the graphics.