Skip to content
Snippets Groups Projects
Commit 2419179d authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

modify prefactors

parent 03fed056
Branches
Tags
No related merge requests found
Pipeline #
...@@ -441,13 +441,13 @@ int kspace<be, dt>::filter_calibrated_ell( ...@@ -441,13 +441,13 @@ int kspace<be, dt>::filter_calibrated_ell(
{ {
this->template low_pass<rnumber, fc>( this->template low_pass<rnumber, fc>(
a, a,
5.5 / ell); 2.8 / ell);
} }
else if (filter_type == std::string("Gauss")) else if (filter_type == std::string("Gauss"))
{ {
this->template Gauss_filter<rnumber, fc>( this->template Gauss_filter<rnumber, fc>(
a, a,
0.257*ell); 0.23*ell);
} }
else if (filter_type == std::string("ball")) else if (filter_type == std::string("ball"))
{ {
......
...@@ -38,7 +38,8 @@ def phi_b( ...@@ -38,7 +38,8 @@ def phi_b(
return phi return phi
def hat_phi_b( def hat_phi_b(
k, ell): k, ell,
prefactor = 1.0):
arg = k * ell / 2 arg = k * ell / 2
phi = (3. / arg**3) * (np.sin(arg) - arg*np.cos(arg)) phi = (3. / arg**3) * (np.sin(arg) - arg*np.cos(arg))
return phi return phi
...@@ -53,7 +54,7 @@ def phi_s( ...@@ -53,7 +54,7 @@ def phi_s(
return phi return phi
def hat_phi_s( def hat_phi_s(
k, ell, prefactor = 5.5): k, ell, prefactor = 2.8):
kc = prefactor / ell kc = prefactor / ell
bindices = np.where(np.abs(k) > kc) bindices = np.where(np.abs(k) > kc)
phi = np.ones(k.shape, k.dtype) phi = np.ones(k.shape, k.dtype)
...@@ -70,7 +71,7 @@ def phi_g( ...@@ -70,7 +71,7 @@ def phi_g(
def hat_phi_g( def hat_phi_g(
k, ell, k, ell,
prefactor = 0.257): prefactor = 0.23):
sigma = prefactor * ell sigma = prefactor * ell
return np.exp(-0.5 * (k * sigma)**2) return np.exp(-0.5 * (k * sigma)**2)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment