tests/testthat/test-Layer2Derivatives.R

In ak-git/RSM: Resistivity Sounding Measurements

test_that("dRdrho1N", {
  0.01 -> dx
  20 -> lmm
  9 -> rho1
  layer2Ohms(rho1 = rho1, rho2 = 1, hmm = lmm * 0.2, smm = lmm * 0.5, lmm = lmm) -> R1
  layer2Ohms(rho1 = rho1 + dx, rho2 = 1, hmm = lmm * 0.2, smm = lmm * 0.5, lmm = lmm) -> R2
  ((R2 - R1) / dx) * (rho1 / R1) -> der
  expect_equal(round(dRdrho1N(-0.8, 0.2, 0.5), 3), round(der, 3))
})

test_that("dRdrho2N", {
  0.01 -> dx
  20 -> lmm
  1 -> rho2
  layer2Ohms(rho1 = 9, rho2 = rho2, hmm = lmm * 0.2, smm = lmm * 0.5, lmm = lmm) -> R1
  layer2Ohms(rho1 = 9, rho2 = rho2 + dx, hmm = lmm * 0.2, smm = lmm * 0.5, lmm = lmm) -> R2
  ((R2 - R1) / dx) * (rho2 / R1) -> der
  expect_equal(round(dRdrho2N(-0.8, 0.2, 0.5), 3), round(der, 3))
})