Skip to content
Snippets Groups Projects
Commit e487ee0b authored by Pumpe, Daniel (dpumpe)'s avatar Pumpe, Daniel (dpumpe)
Browse files

update & fixes

parent a6f6ef81
Branches
Tags
2 merge requests!100Docstring operators,!92Fix fft side effect
Pipeline #
...@@ -83,10 +83,10 @@ if __name__ == "__main__": ...@@ -83,10 +83,10 @@ if __name__ == "__main__":
a_s = sigma_s ** 2. * lambda_s * total_volume a_s = sigma_s ** 2. * lambda_s * total_volume
# creation of spaces # creation of spaces
x1 = RGSpace(npix, dtype=np.float64, distances=total_volume / npix, x1 = RGSpace(npix, distances=total_volume / npix,
zerocenter=False) zerocenter=False)
k1 = RGRGTransformation.get_codomain(x1) k1 = RGRGTransformation.get_codomain(x1)
p1 = PowerSpace(harmonic_domain=k1, log=False, dtype=np.float64) p1 = PowerSpace(harmonic_domain=k1, log=False)
# creating Power Operator with given spectrum # creating Power Operator with given spectrum
spec = (lambda k: a_s / (1 + (k / k_0) ** 2) ** 2) spec = (lambda k: a_s / (1 + (k / k_0) ** 2) ** 2)
...@@ -97,8 +97,8 @@ if __name__ == "__main__": ...@@ -97,8 +97,8 @@ if __name__ == "__main__":
Fft_op = FFTOperator(domain=x1, target=k1, Fft_op = FFTOperator(domain=x1, target=k1,
domain_dtype=np.float64, domain_dtype=np.float64,
target_dtype=np.complex128) target_dtype=np.complex128)
R_op = ResponseOperator(x1, sigma=length_convolution, R_op = ResponseOperator(x1, sigma=[length_convolution],
exposure=exposure) exposure=[exposure])
# drawing a random field # drawing a random field
sk = p_field.power_synthesize(real_signal=True, mean=0.) sk = p_field.power_synthesize(real_signal=True, mean=0.)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment