bug fixes

Author: tobiste

R/stress_inversion.R

@@ -57,7 +57,7 @@
 #' @examples
 #' set.seed(20250411)
 #' # Use Angelier examples
-#' par(mfrow = c(1, length(angelier1990)))
+#' par(mfrow = c(1, length(angelier1990)), mar = c(4, 1, 4, 1))
 #' invisible(lapply(angelier1990, function(x) {
 #'   # Inversion after Michael (1984)
 #'   res_michael <- slip_inversion(x, method = "michael", n_iter = 100, n = 100, res = 100)
@@ -82,11 +82,13 @@
 #'     pch = 1,
 #'     legend = c("Michael (1984)", "Angelier (1990)", "Yamaji & Sato (2006)", 
 #'       "Hansen (2013)", "WISSI"),
-#'     col = 2:6
+#'     col = 2:6, cex = 0.8
 #'   )
-#'   legend("bottomright",
+#'   legend("bottom",
 #'     pch = 1:3,
-#'     legend = c("S1", "S2", "S3")
+#'     legend = c("S1", "S2", "S3"),
+#'     horiz = TRUE,
+#'     cex = 0.8
 #'   )
 #' }))
 slip_inversion <- function(x, method = c("michael", "angelier", "hansen", "yamaji", "wissi"), ...) {

R/stress_inversion_hansen.R

@@ -310,7 +310,7 @@ slip_inversion_hansen <- function(x, flip = FALSE, type = c("9d", "6d")) {
   # Theoretically resolved shear stress on plane
   sigma_s_mean <- mean(abs(shear_stress(val_Ts[1], val_Ts[3], theta)))
 
-  SHmax <- SH(p$principal_axes[1, ], p$principal_axes[2, ], p$principal_axes[3, ], R = stress_shape$R)
+  SHmax <- SH(principal_axes[1, ], principal_axes[2, ], principal_axes[3, ], R = stress_shape$R)
 
   # shearnorm <- tau2shearnorm(Ts, x, friction = friction)
   # tendency <- tau2tendency(Ts, x, friction = friction)

R/stress_inversion_yamaji.R

@@ -171,6 +171,7 @@ yamaji_sato <- function(normals, slips, wt) {
 #' @seealso [slip_inversion_yamaji_sato_boot()]
 #'
 #' @examples
+#' set.seed(20250411)
 #' nx <- length(angelier1990)
 #' par(mfrow = c(1, nx))
 #'
@@ -266,7 +267,7 @@ slip_inversion_yamaji_sato <- function(x, weights = NULL, flip = FALSE) {
 
   # sigma_s_mean <- mean(abs(shearnorm))
 
-  SHmax <- SH(p$principal_axes[1, ], p$principal_axes[2, ], p$principal_axes[3, ], R = stress_shape$R)
+  SHmax <- SH(p$principal_axes[1, ], p$principal_axes[2, ], p$principal_axes[3, ], R = shape$R)
 
   # shearnorm <- tau2shearnorm(TR, x, friction = friction)
   # tendency <- tau2tendency(TR, x, friction = friction)

man/slip_inversion.Rd

@@ -23,6 +23,10 @@ set.seed(20250411)
 ```
 
 
+```{r load_functions, include=FALSE}
+library("structr")
+```
+
 The offset along a fault can be factorized into several components.
 
 ![Fig. 1: Graphic illustration of displacement components along a fault.