diff options
author | Jaron Kent-Dobias <jkentdobias@g.hmc.edu> | 2016-08-22 10:11:14 -0400 |
---|---|---|
committer | Jaron Kent-Dobias <jkentdobias@g.hmc.edu> | 2016-08-22 10:11:14 -0400 |
commit | 2bb0740b68fdb62d45adc00204b3990ca42ade77 (patch) | |
tree | a52975e3460da781467ddb70aaa8d76840d01bb4 /src/beta_scales.c | |
download | fuse_networks-2bb0740b68fdb62d45adc00204b3990ca42ade77.tar.gz fuse_networks-2bb0740b68fdb62d45adc00204b3990ca42ade77.tar.bz2 fuse_networks-2bb0740b68fdb62d45adc00204b3990ca42ade77.zip |
started repo again without all the data files gunking the works
Diffstat (limited to 'src/beta_scales.c')
-rw-r--r-- | src/beta_scales.c | 21 |
1 files changed, 21 insertions, 0 deletions
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; +} |