Skip to content

Commit

Permalink
fix-up to previous merge
Browse files Browse the repository at this point in the history
  • Loading branch information
joezuntz committed Nov 22, 2024
1 parent 967c70e commit f8b776a
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions shear/point_mass/add_gammat_point_mass.py
Original file line number Diff line number Diff line change
Expand Up @@ -114,7 +114,7 @@ def execute(block, config):
nz_source = block[config['source_nz_section'], "bin_%d" % (j2 + 1)][1:]

ng_array_source_rep = np.tile(nz_source.reshape(1, len(z_source)), (len(z_lens), 1))
int_sourcez = simpson(ng_array_source_rep * (num / chi_smat), z_source)
int_sourcez = simpson(ng_array_source_rep * (num / chi_smat), x=z_source)

coeff_ints = sigcrit_inv_fac_Mpc_Msun

Expand All @@ -124,7 +124,7 @@ def execute(block, config):
# Since gamma_t is a scalar it should be same in both physical and comoving coordinates
# It is just easier to match the expressions in comoving coordinates to the ones on methods paper.

betaj1j2_pm = simpson(nz_lens * Is * (1./chi_lens**2), z_lens)
betaj1j2_pm = simpson(nz_lens * Is * (1./chi_lens**2), x=z_lens)
if (config['use_fiducial']):
config['betaj1j2'][str(j1) + '_' + str(j2)] = betaj1j2_pm
else:
Expand Down

0 comments on commit f8b776a

Please sign in to comment.