#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;
}