The position of the data points in 3D. The size of the vector should equal to the number of points time three.
The validity of each axis
The weight of this dataset
The position of the data poinrs in 3D
The data, one for each point
The line of sight of the satellite
The weight of this dataset
An array of Vectord, each entry of the array represents the inverted displ of a surface.
Allow to do slip inversion using multiple type of dataset. This is still a work in progress (optimization, new dataset, ...)
This is almost based on
Maerten, F., Resor, P., Pollard, D., & Maerten, L. (2005). Inverting for slip on three-dimensional fault surfaces using angular dislocations. Bulletin of the Seismological Society of America, 95(5), 1654-1665.
except that we use an iterative approach to solve for the unknown burger components on fault surfaces. Consequently, the regularization operator (aka, Tikhonov regularization) is also simply incorporated into the iterative process. The FNNLS (Fast Non Negativity Least Squares) algorithm is simply replaced by the notion of DIC (Displacement Inequality Constraint).
The iterative approach is similar to the Block Gauss-Seidel solver (Forward with '
seidel
' parameter), except that we use a Block Least-Squares formulation, with an optimized building process for the underlaying matrices and vectors.Separate the notion of data (gps, insar, ...) from the solver itself.
Another example