Skip to contents

Calculates the reduced stress tensor (without optimization).

Usage

reduced_stress(fault, method = c("michael", "angelier"))

Arguments

fault

"Fault" object where the rows are the observations, and the columns the coordinates.

method

character. The inversion algorithm, one of "michael" (the default) for a bootstrapped linear inversion after Micheal (1984), "angelier" for an iterative direct inversion after Angelier (1990) and Mostafa (2005), and "hansen" for direct inversion using the 9d parameter space after Hansen (2013).

Value

"ellipsoid" object. 3x3 matrix (i.e. 2nd rank tensor).

Examples

f <- angelier1990$TYM
reduced_stress(f)
#> Ellipsoid tensor
#>            [,1]       [,2]       [,3]
#> [1,] -0.8168187  0.1113011 -0.1570288
#> [2,]  0.1113011 -0.6797825 -0.1514454
#> [3,] -0.1570288 -0.1514454  1.4966011