Home

Welcome to the triinv wiki!

triinv is a set of Matlab functions used to estimate the distribution of slip on faults parametrized using triangular dislocation elements (TDEs) constrained by displacement and/or stress data.

The basis of triinv is the set of algorithms relating slip on TDEs to displacement and strain throughout a homogeneous elastic half-space, as described by Meade (2007). These algorithms are implemented in Matlab as the tridisl suite, which uses the Parallel Computing Toolbox to parallelize the calculation of displacement and strain due to unit slip on each triangular dislocation element.

Various constraints can be imposed on the estimation, including Laplacian smoothing of the slip distribution, edge constraints that force the slip to be zero at fault edges, sign constraints that force unidirectional slip in the strike and/or dip directions, and Total Variation Regularization (TVR) that yields a slip distribution featuring clusters of distinct slip values rather than a spatially smooth distribution. (Requires CVX Matlab package to be installed and set up.)

triinvx.m is the master function used to estimate slip, with basic usage of:

u = triinvx(p, s, beta);

where p is a structure containing information about the triangular dislocation elements, s is a structure containing the constraining data, and beta defines the weighting of regularization applied in the slip estimation. Information about these and optional input arguments can be found on the inputs page. Slip is returned to the vector u; see the outputs page for details about its structure and other output arguments.