From 2bb0740b68fdb62d45adc00204b3990ca42ade77 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Mon, 22 Aug 2016 10:11:14 -0400 Subject: started repo again without all the data files gunking the works --- src/beta_scales.c | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) create mode 100644 src/beta_scales.c (limited to 'src/beta_scales.c') diff --git a/src/beta_scales.c b/src/beta_scales.c new file mode 100644 index 0000000..8cbc5d1 --- /dev/null +++ b/src/beta_scales.c @@ -0,0 +1,21 @@ + +#include "fracture.h" + +double beta_scaling_flat(double beta, double x, double y) { return beta; } + +double beta_mag(double beta) { + double aa = -7.52579; + double bb = 9.63706; + double cc = 0.692515; + double dd = -2.47638; + + return gsl_sf_exp(aa + bb * gsl_sf_erf(cc * (gsl_sf_log(beta) - dd))); +} + +double beta_scaling_gaus(double beta, double x, double y) { + double sigma = 0.25; + double nu_f = 1.56; + double betap = beta * gsl_sf_exp((gsl_pow_2(x - 0.5) + gsl_pow_2(y - 0.5)) / + (2 * sigma * 2 * nu_f)); + return betap; +} -- cgit v1.2.3-70-g09d2