Results

Results from benchmark runs. Place tarballs containing the requested results in this folder and describe the run in the description field.

PyLith, 1 Proc, Linear Tet, 1 km Resolution, dt=0.1yr: PyLith results run on 1 processor of a Power Mac G5. Linear tetrahedral mesh at 1 km resolution. Constant time step size of 0.1 year.
PyLith, 1 Proc, Linear Hex, 1 km Resolution, dt=0.1yr: PyLith results run on 1 processor of a Power Mac G5. Linear hexahedral mesh at 1 km resolution. Constant time step size of 0.1 year.
PyLith, 1 Proc, Linear Tet, 500 m Resolution, dt=0.1yr: PyLith results run on 1 processor of a Power Mac G5. Linear tetrahedral mesh at 500 m resolution. Constant time step size of 0.1 year.
PyLith Revised Results, 500m, New BC and Split Node Input: New solution using revised BC and split node inputs. The revised BC take care of the problems of defining BC on the fault plane (or in some cases the projected fault plane). The new split node inputs no longer assume a bilinear slip distribution in the region where the fault tapers overlap, and now assumes a taper consistent with what is used for the analytical solution.
PyLith Revised Results, 500m, Altered BC for Viscoelastic Solution: New version where BC have been altered from those of previous version to make viscoelastic results consistent with those from GeoFEST. The revised BC do not pin y-component on the y=0 plane, and no BC are applied along the intersection of the fault plane (or its projection) along y=0 and z=-24.
PyLith, 1 Proc, Linear Tet, 250 m Resolution, dt=0.1yr: PyLith results run on 1 processor of an Opteron 2.4GHz Linux machine. Linear tetrahedral mesh at 250 m resolution. Constant time step size of 0.1 year.
GeoFEST/PYRAMID 1km: Parallel results using 64 processors of Intel/Linux Cluster with GeoFEST-4.5 and Pyramid-2.1.3
GeoFEST/PYRAMID 500m: Parallel results using 64 processors of Intel/Linux Cluster with GeoFEST-4.5 and Pyramid-2.1.3
GeoFEST/PYRAMID 250m: Parallel results using 128 processors of Intel/Linux Cluster with GeoFEST-4.5 and Pyramid-2.1.3
GeoFEST Linear Tet 1km resolution dt=0.1yr (updated): The taper error has been fixed.
GeoFEST Linear-Tet 500m Re-Run
GeoFEST Linear-Tet 250m Re-Run
Femlab 1 km resolution, t = 0 years: This model has ~162000 linear tetrahedral elements and is twice the size in y of the model description, since there is no symmetric boundary. This yields a resolution close to 1 km. The model and solver require about 800 MB and is solved in about 3 minutes on a 1.8 GHz AMD Opteron. An iterative solver was used, which uses the Incomplete LU preconditioner with a drop tolerance of 0.01. Decreasing this value has very little effect on the error but takes longer to solve.
Femlab 1 km resolution, t = 1 year: Viscoelastic problem requires ~3.5GB and takes about 4.5 hrs to run. Drop tolerance is 0.01.
Femlab 1 km resolution, t = 5 years: Viscoelastic problem requires ~3.5GB and takes about 4.5 hrs to run. Drop tolerance is 0.01.
Femlab 1 km resolution, t = 10 years: Viscoelastic problem requires ~3.5GB and takes about 4.5 hrs to run. Drop tolerance is 0.01.