TransitionListener 2.0.0__py3-none-any.whl

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (95) hide show
  1. transitionlistener/__init__.py +21 -0
  2. transitionlistener/bubbledynamics.py +2250 -0
  3. transitionlistener/bubbledynamics_fixedstep.py +2929 -0
  4. transitionlistener/check_installation.py +309 -0
  5. transitionlistener/ckm.py +91 -0
  6. transitionlistener/colors.py +288 -0
  7. transitionlistener/config.py +185 -0
  8. transitionlistener/constants.py +64 -0
  9. transitionlistener/counterterms/__init__.py +230 -0
  10. transitionlistener/counterterms/cw.py +298 -0
  11. transitionlistener/counterterms/tensors.py +343 -0
  12. transitionlistener/counterterms/twohdm.py +520 -0
  13. transitionlistener/counterterms/twohdm_solver.py +374 -0
  14. transitionlistener/errors.py +315 -0
  15. transitionlistener/examples/bp_et.yaml +64 -0
  16. transitionlistener/examples/bp_lisa.yaml +70 -0
  17. transitionlistener/examples/bp_muares.yaml +69 -0
  18. transitionlistener/examples/bp_pta.yaml +64 -0
  19. transitionlistener/examples/example_dark_photon.yaml +73 -0
  20. transitionlistener/examples/example_grid.yaml +55 -0
  21. transitionlistener/examples/example_line.yaml +51 -0
  22. transitionlistener/examples/example_point.yaml +73 -0
  23. transitionlistener/examples/example_random_conformal.yaml +42 -0
  24. transitionlistener/examples/example_ultranest.yaml +39 -0
  25. transitionlistener/finiteT.py +468 -0
  26. transitionlistener/generic_potential.py +1521 -0
  27. transitionlistener/gridplots.py +592 -0
  28. transitionlistener/gwfopt.py +604 -0
  29. transitionlistener/helper_functions.py +843 -0
  30. transitionlistener/hydrodynamics.py +647 -0
  31. transitionlistener/interface/__init__.py +95 -0
  32. transitionlistener/interface/check.py +51 -0
  33. transitionlistener/interface/cli.py +211 -0
  34. transitionlistener/interface/config_loader.py +264 -0
  35. transitionlistener/interface/logging_utils.py +30 -0
  36. transitionlistener/interface/output_schema.py +296 -0
  37. transitionlistener/interface/pipeline.py +300 -0
  38. transitionlistener/interface/samplers.py +565 -0
  39. transitionlistener/interface/scans.py +468 -0
  40. transitionlistener/interface/single_point.py +581 -0
  41. transitionlistener/interface/state.py +61 -0
  42. transitionlistener/lineplots.py +597 -0
  43. transitionlistener/logo/TL-logo_large.png +0 -0
  44. transitionlistener/logo/TL-logo_small.png +0 -0
  45. transitionlistener/models/TL_2HDM.py +512 -0
  46. transitionlistener/models/TL_conformal_dark_u1.py +202 -0
  47. transitionlistener/models/TL_dark_U1.py +198 -0
  48. transitionlistener/models/TL_dark_U1_g_parameterization.py +198 -0
  49. transitionlistener/models/TL_dark_flipflop.py +274 -0
  50. transitionlistener/models/templatePotential.py +260 -0
  51. transitionlistener/nucleation.py +412 -0
  52. transitionlistener/observability.py +790 -0
  53. transitionlistener/particles.py +436 -0
  54. transitionlistener/pathDeformation.py +1462 -0
  55. transitionlistener/percolation_adaptive_gridbuilders.py +1385 -0
  56. transitionlistener/percolation_adaptive_rate.py +1049 -0
  57. transitionlistener/percolation_adaptivestepsize.py +1632 -0
  58. transitionlistener/percolation_fixedstepsize.py +618 -0
  59. transitionlistener/phases.py +1499 -0
  60. transitionlistener/plot_settings.py +121 -0
  61. transitionlistener/plots.py +1781 -0
  62. transitionlistener/runtime_options.py +278 -0
  63. transitionlistener/tab_data/30f_fs{hd}_100fDMGP_ceffyl/bandwidths.npy +0 -0
  64. transitionlistener/tab_data/30f_fs{hd}_100fDMGP_ceffyl/density.npy +0 -0
  65. transitionlistener/tab_data/30f_fs{hd}_100fDMGP_ceffyl/freqs.npy +0 -0
  66. transitionlistener/tab_data/30f_fs{hd}_100fDMGP_ceffyl/log.txt +16 -0
  67. transitionlistener/tab_data/30f_fs{hd}_100fDMGP_ceffyl/log10rhogrid.npy +0 -0
  68. transitionlistener/tab_data/30f_fs{hd}_100fDMGP_ceffyl/log10rholabels.txt +30 -0
  69. transitionlistener/tab_data/30f_fs{hd}_100fDMGP_ceffyl/pulsar_list.txt +1 -0
  70. transitionlistener/tab_data/HLV-O2-PLI.csv +82 -0
  71. transitionlistener/tab_data/HLV-O2-noise.csv +83 -0
  72. transitionlistener/tab_data/HLVK-design-PLI.csv +61 -0
  73. transitionlistener/tab_data/HLVK-design-noise.csv +83 -0
  74. transitionlistener/tab_data/NG15_PLI_hc.csv +31 -0
  75. transitionlistener/tab_data/NG15_noise_hc.csv +56 -0
  76. transitionlistener/tab_data/PLI-dataset.csv +327 -0
  77. transitionlistener/tab_data/energydensities.dat +1001 -0
  78. transitionlistener/tab_data/finiteT_b.dat.txt +1999 -0
  79. transitionlistener/tab_data/finiteT_f.dat.txt +1000 -0
  80. transitionlistener/tab_data/geff_SM.dat +10009 -0
  81. transitionlistener/tab_data/muAres_f_Hz_PLI_h2OmegaGW.csv +131 -0
  82. transitionlistener/tab_data/muAres_f_Hz_h2Omega_noise.csv +99 -0
  83. transitionlistener/tab_data/noise-dataset.csv +327 -0
  84. transitionlistener/tab_data/numberdensities.dat +1001 -0
  85. transitionlistener/tab_data/pressures.dat +1001 -0
  86. transitionlistener/thermodynamics.py +686 -0
  87. transitionlistener/transitionObservables.py +1199 -0
  88. transitionlistener/transitionObservables_fixedstep.py +976 -0
  89. transitionlistener/transitions.py +1000 -0
  90. transitionlistener/tunneling1D.py +910 -0
  91. transitionlistener-2.0.0.dist-info/METADATA +817 -0
  92. transitionlistener-2.0.0.dist-info/RECORD +95 -0
  93. transitionlistener-2.0.0.dist-info/WHEEL +4 -0
  94. transitionlistener-2.0.0.dist-info/entry_points.txt +2 -0
  95. transitionlistener-2.0.0.dist-info/licenses/LICENSE.txt +674 -0
@@ -0,0 +1,278 @@
1
+ """Optional runtime overrides for scripts and benchmark runners.
2
+
3
+ The core configuration classes live in :mod:`transitionlistener.config`. This
4
+ module keeps input-dictionary parsing and named precision bundles out of that
5
+ defaults-only file.
6
+
7
+ Part of TransitionListener v2.0
8
+ Documentation: https://tasillo.de/TransitionListener/
9
+
10
+ Authors:
11
+ Jonas Matuszak <jonas.matuszak@kit.edu>
12
+ Carlo Tasillo <carlo.tasillo@ific.uv.es>
13
+ """
14
+
15
+ import copy
16
+ from collections.abc import Mapping
17
+
18
+ from transitionlistener.config import TracingConf
19
+ from transitionlistener import errors
20
+
21
+
22
+ ALGORITHM_MODES = ("adaptive_step_size", "fixed_step_size")
23
+
24
+
25
+ PRECISION_OVERRIDE_KEYS = (
26
+ "precision_trace_field_accuracy",
27
+ "precision_trace_temp_accuracy",
28
+ "precision_trace_dtstart",
29
+ "precision_trace_tjump",
30
+ "precision_trace_dtmin",
31
+ "precision_trace_local_n",
32
+ "precision_trace_local_edge",
33
+ "precision_deform_converge_0",
34
+ "precision_deform_fRatioConv",
35
+ )
36
+
37
+ PERCOLATION_OVERRIDE_KEYS = (
38
+ "percolation_integral_method",
39
+ "percolation_time_temperature_mode",
40
+ "percolation_n_action",
41
+ "percolation_n_action_min",
42
+ "percolation_n_action_increment",
43
+ "percolation_n_action_max",
44
+ "percolation_maxit",
45
+ "percolation_max_action_temperatures",
46
+ "percolation_acc_tperc",
47
+ "percolation_acc_tfinal",
48
+ "percolation_acc_rh",
49
+ "percolation_weak_threshold",
50
+ "percolation_jitter_GH4_threshold",
51
+ "percolation_action_jitter_tunneltight_rescue",
52
+ "percolation_n_jitter_save",
53
+ )
54
+
55
+ GW_OVERRIDE_KEYS = (
56
+ "gw_wall_velocity",
57
+ )
58
+
59
+ PRECISION_MODE_NAMES = ("default", "robust", "xtrace", "tunneltight", "benchmark")
60
+
61
+
62
+ def normalize_precision_mode(precision_mode: str | None) -> str:
63
+ """Return a validated precision-mode label, defaulting ``None`` to ``"default"``."""
64
+ if precision_mode is None:
65
+ return "default"
66
+ mode = str(precision_mode)
67
+ if mode not in PRECISION_MODE_NAMES:
68
+ raise ValueError(f"Unknown precision_mode {precision_mode!r}.")
69
+ return mode
70
+
71
+
72
+ def precision_mode_profile(mode: str) -> dict:
73
+ """Materialise the named tracing/tunneling override bundle for ``mode``."""
74
+ default_conf = TracingConf()
75
+ default_trace = default_conf.tracing_args
76
+ default_tunnel = default_conf.tunneling_params["deformation_deform_params"]
77
+ profile = {
78
+ "precision_trace_field_accuracy": default_conf.tracing_field_accuracy,
79
+ "precision_trace_temp_accuracy": default_conf.tracing_temp_accuracy,
80
+ "precision_trace_dtstart": default_trace["dtstart"],
81
+ "precision_trace_tjump": default_trace["tjump"],
82
+ "precision_trace_dtmin": default_trace["dtmin"],
83
+ "precision_trace_local_n": default_trace["local_min_args"]["n"],
84
+ "precision_trace_local_edge": default_trace["local_min_args"]["edge"],
85
+ "precision_deform_converge_0": default_tunnel["converge_0"],
86
+ "precision_deform_fRatioConv": default_tunnel["fRatioConv"],
87
+ }
88
+
89
+ if mode == "robust":
90
+ profile.update(
91
+ precision_trace_field_accuracy=3e-4,
92
+ precision_trace_temp_accuracy=3e-4,
93
+ precision_trace_dtstart=5e-5,
94
+ precision_trace_tjump=5e-6,
95
+ precision_trace_dtmin=3e-7,
96
+ precision_trace_local_n=200,
97
+ precision_trace_local_edge=0.03,
98
+ )
99
+ if mode in {"xtrace", "benchmark"}:
100
+ profile.update(
101
+ precision_trace_field_accuracy=1e-4,
102
+ precision_trace_temp_accuracy=1e-4,
103
+ precision_trace_dtstart=2e-5,
104
+ precision_trace_tjump=2e-6,
105
+ precision_trace_dtmin=1e-7,
106
+ precision_trace_local_n=300,
107
+ precision_trace_local_edge=0.02,
108
+ )
109
+ if mode in {"tunneltight", "benchmark"}:
110
+ profile.update(
111
+ precision_deform_converge_0=1.0,
112
+ precision_deform_fRatioConv=5e-3,
113
+ )
114
+ return profile
115
+
116
+
117
+ def apply_precision_overrides(tracing_conf, overrides: dict[str, float | int | None]) -> None:
118
+ """Apply individual tracing and path-deformation overrides in place."""
119
+ if overrides["precision_trace_field_accuracy"] is not None:
120
+ tracing_conf.tracing_field_accuracy = float(overrides["precision_trace_field_accuracy"])
121
+ if overrides["precision_trace_temp_accuracy"] is not None:
122
+ tracing_conf.tracing_temp_accuracy = float(overrides["precision_trace_temp_accuracy"])
123
+
124
+ tracing = copy.deepcopy(tracing_conf.tracing_args)
125
+ for key, trace_key in (
126
+ ("precision_trace_dtstart", "dtstart"),
127
+ ("precision_trace_tjump", "tjump"),
128
+ ("precision_trace_dtmin", "dtmin"),
129
+ ):
130
+ if overrides[key] is not None:
131
+ tracing[trace_key] = float(overrides[key])
132
+
133
+ local = dict(tracing.get("local_min_args", {}))
134
+ if overrides["precision_trace_local_n"] is not None:
135
+ local["n"] = int(round(float(overrides["precision_trace_local_n"])))
136
+ if overrides["precision_trace_local_edge"] is not None:
137
+ local["edge"] = float(overrides["precision_trace_local_edge"])
138
+ tracing["local_min_args"] = local
139
+ tracing_conf.tracing_args = tracing
140
+
141
+ tunneling = copy.deepcopy(tracing_conf.tunneling_params)
142
+ deform = dict(tunneling.get("deformation_deform_params", {}))
143
+ if overrides["precision_deform_converge_0"] is not None:
144
+ deform["converge_0"] = float(overrides["precision_deform_converge_0"])
145
+ if overrides["precision_deform_fRatioConv"] is not None:
146
+ deform["fRatioConv"] = float(overrides["precision_deform_fRatioConv"])
147
+ tunneling["deformation_deform_params"] = deform
148
+ tracing_conf.tunneling_params = tunneling
149
+
150
+
151
+ def apply_precision_mode(tracing_conf, precision_mode: str | None) -> str:
152
+ """Resolve ``precision_mode`` and apply its predefined override profile."""
153
+ mode = normalize_precision_mode(precision_mode)
154
+ if mode != "default":
155
+ apply_precision_overrides(tracing_conf, precision_mode_profile(mode))
156
+ return mode
157
+
158
+
159
+ def _as_int(value, minimum: int) -> int:
160
+ """Cast ``value`` to ``int`` while enforcing a lower bound."""
161
+ return max(int(round(float(value))), int(minimum))
162
+
163
+
164
+ def _as_float(value, minimum: float | None = None) -> float:
165
+ """Cast ``value`` to ``float`` while optionally enforcing a lower bound."""
166
+ number = float(value)
167
+ if minimum is not None:
168
+ number = max(number, float(minimum))
169
+ return number
170
+
171
+
172
+ def apply_percolation_overrides(percolation_conf, overrides: dict[str, object]) -> None:
173
+ """Apply runtime percolation controls parsed from model or scan inputs."""
174
+ overrides = {key: overrides.get(key) for key in PERCOLATION_OVERRIDE_KEYS}
175
+
176
+ if overrides["percolation_integral_method"] is not None:
177
+ method = str(overrides["percolation_integral_method"])
178
+ if method not in {"ode", "double_integral"}:
179
+ raise ValueError("percolation_integral_method must be 'ode' or 'double_integral'.")
180
+ percolation_conf.integral_method = method
181
+ if overrides["percolation_time_temperature_mode"] is not None:
182
+ mode = str(overrides["percolation_time_temperature_mode"])
183
+ if mode not in {"sound_speed", "bag"}:
184
+ raise ValueError("percolation_time_temperature_mode must be 'sound_speed' or 'bag'.")
185
+ percolation_conf.time_temperature_mode = mode
186
+
187
+ int_minimums = {
188
+ "percolation_n_action": ("n_action", 2),
189
+ "percolation_n_action_min": ("n_action_min", 2),
190
+ "percolation_n_action_increment": ("n_action_increment", 1),
191
+ "percolation_n_action_max": ("n_action_max", 2),
192
+ "percolation_maxit": ("maxit", 1),
193
+ "percolation_max_action_temperatures": ("max_action_temperatures", 1),
194
+ "percolation_n_jitter_save": ("n_jitter_save", 0),
195
+ }
196
+ for key, (attr, minimum) in int_minimums.items():
197
+ if overrides[key] is not None:
198
+ setattr(percolation_conf, attr, _as_int(overrides[key], minimum))
199
+
200
+ float_minimums = {
201
+ "percolation_acc_tperc": ("acc_tperc", 0.0),
202
+ "percolation_acc_tfinal": ("acc_tfinal", 0.0),
203
+ "percolation_acc_rh": ("acc_rh", 0.0),
204
+ "percolation_weak_threshold": ("weak_threshold", 0.0),
205
+ "percolation_jitter_GH4_threshold": ("jitter_GH4_threshold", 0.0),
206
+ }
207
+ for key, (attr, minimum) in float_minimums.items():
208
+ if overrides[key] is not None:
209
+ setattr(percolation_conf, attr, _as_float(overrides[key], minimum=minimum))
210
+
211
+ bool_overrides = {
212
+ "percolation_action_jitter_tunneltight_rescue": "action_jitter_tunneltight_rescue",
213
+ }
214
+ for key, attr in bool_overrides.items():
215
+ if overrides[key] is not None:
216
+ setattr(percolation_conf, attr, bool(overrides[key]))
217
+
218
+ percolation_conf.n_action_max = max(
219
+ int(percolation_conf.n_action_max),
220
+ int(percolation_conf.n_action_min),
221
+ )
222
+
223
+
224
+ def apply_gw_overrides(gw_conf, overrides: dict[str, object]) -> None:
225
+ """Apply runtime GW-observable overrides consumed from scan inputs."""
226
+ if overrides.get("gw_wall_velocity") is not None:
227
+ gw_conf.wall_velocity = overrides["gw_wall_velocity"]
228
+
229
+
230
+ def apply_input_overrides(potential, args: tuple, dargs: dict) -> tuple[tuple, dict]:
231
+ """Pop shared runtime-override keys from the model input dictionary.
232
+
233
+ Consumes ``precision_mode``, ``percolation_algorithm_mode`` and the keyed
234
+ override dicts; leaves the remaining entries untouched for model-specific
235
+ parameter validation.
236
+ """
237
+ args = list(args)
238
+ dargs = dict(dargs)
239
+ inputparam_dict = None
240
+
241
+ if args and isinstance(args[0], Mapping):
242
+ inputparam_dict = dict(args[0])
243
+ args[0] = inputparam_dict
244
+ else:
245
+ for key in ("inputparam_dict", "input_dict"):
246
+ value = dargs.get(key)
247
+ if isinstance(value, Mapping):
248
+ inputparam_dict = dict(value)
249
+ dargs[key] = inputparam_dict
250
+ break
251
+
252
+ if inputparam_dict is None:
253
+ return tuple(args), dargs
254
+
255
+ precision_mode = inputparam_dict.pop("precision_mode", None)
256
+ algorithm_mode = inputparam_dict.pop("percolation_algorithm_mode", None)
257
+ precision_overrides = {key: inputparam_dict.pop(key, None) for key in PRECISION_OVERRIDE_KEYS}
258
+ percolation_overrides = {key: inputparam_dict.pop(key, None) for key in PERCOLATION_OVERRIDE_KEYS}
259
+ gw_overrides = {key: inputparam_dict.pop(key, None) for key in GW_OVERRIDE_KEYS}
260
+
261
+ try:
262
+ apply_precision_mode(potential.config.tracingConf, precision_mode)
263
+ except ValueError as exc:
264
+ raise errors.InitPotentialError(str(exc)) from exc
265
+
266
+ if algorithm_mode is not None:
267
+ algorithm_mode = str(algorithm_mode)
268
+ if algorithm_mode not in ALGORITHM_MODES:
269
+ raise errors.InitPotentialError(
270
+ f"percolation_algorithm_mode must be one of {ALGORITHM_MODES}, "
271
+ f"got {algorithm_mode!r}."
272
+ )
273
+ potential.config.percolationConf.algorithm_mode = algorithm_mode
274
+
275
+ apply_precision_overrides(potential.config.tracingConf, precision_overrides)
276
+ apply_percolation_overrides(potential.config.percolationConf, percolation_overrides)
277
+ apply_gw_overrides(potential.config.gwConf, gw_overrides)
278
+ return tuple(args), dargs
@@ -0,0 +1,16 @@
1
+ Date created:time.struct_time(tm_year=2023, tm_mon=5, tm_mday=15, tm_hour=16, tm_min=52, tm_sec=35, tm_wday=0, tm_yday=135, tm_isdst=1)
2
+ la-forge cores: ['/data/taylor_group/william_lamb/analyses/NANOGrav/NG15_v1p1/dm_gp/ng15_dmgp/fshd_dmgp/chain.core']
3
+ KDE function: FFTKDE
4
+ KDE kernel choice: epanechnikov
5
+ Bandwidth choice: <function sj at 0x2b8ea42b2040>
6
+ Bandwidth values: [0.01382799 0.0076616 0.00881017 0.01775468 0.03435934 0.0352577
7
+ 0.03740746 0.02612554 0.03060202 0.0316155 0.03392853 0.02787687
8
+ 0.0277001 0.03093307 0.03309433 0.05967801 0.03425261 0.02855942
9
+ 0.03035076 0.03150779 0.028027 0.03046142 0.03247571 0.02890307
10
+ 0.03327375 0.01868055 0.0296922 0.03165958 0.03124277 0.0287242 ]
11
+ log10rho grid: [-9.99999894 -9.99952853 -9.99905811 ... -5.29728896 -5.29681855
12
+ -5.29634814]
13
+ Thin chain for KDE? False
14
+ Thin chain for bandwidth? False
15
+ Bootstrap? False
16
+ Number of bootstrapped samples: None
@@ -0,0 +1,30 @@
1
+ gwb_log10_rho_0
2
+ gwb_log10_rho_1
3
+ gwb_log10_rho_2
4
+ gwb_log10_rho_3
5
+ gwb_log10_rho_4
6
+ gwb_log10_rho_5
7
+ gwb_log10_rho_6
8
+ gwb_log10_rho_7
9
+ gwb_log10_rho_8
10
+ gwb_log10_rho_9
11
+ gwb_log10_rho_10
12
+ gwb_log10_rho_11
13
+ gwb_log10_rho_12
14
+ gwb_log10_rho_13
15
+ gwb_log10_rho_14
16
+ gwb_log10_rho_15
17
+ gwb_log10_rho_16
18
+ gwb_log10_rho_17
19
+ gwb_log10_rho_18
20
+ gwb_log10_rho_19
21
+ gwb_log10_rho_20
22
+ gwb_log10_rho_21
23
+ gwb_log10_rho_22
24
+ gwb_log10_rho_23
25
+ gwb_log10_rho_24
26
+ gwb_log10_rho_25
27
+ gwb_log10_rho_26
28
+ gwb_log10_rho_27
29
+ gwb_log10_rho_28
30
+ gwb_log10_rho_29
@@ -0,0 +1,82 @@
1
+ 9.40112e+0;2.76753e-2
2
+ 9.40112e+0;4.18354e-2
3
+ 9.40112e+0;6.18801e-2
4
+ 9.53419e+0;9.32950e-3
5
+ 9.53419e+0;1.67830e-2
6
+ 9.53419e+0;8.57495e-2
7
+ 9.53419e+0;1.21436e-1
8
+ 9.53419e+0;2.13754e-1
9
+ 9.66915e+0;3.28484e-3
10
+ 9.66915e+0;5.41684e-3
11
+ 9.66915e+0;3.92986e-1
12
+ 9.66915e+0;1.16579e+0
13
+ 9.66915e+0;1.80101e+0
14
+ 9.80602e+0;1.37271e-4
15
+ 9.80602e+0;2.46938e-4
16
+ 9.80602e+0;3.42184e-4
17
+ 9.80602e+0;7.32536e-4
18
+ 9.80602e+0;1.08352e-3
19
+ 9.80602e+0;1.90722e-3
20
+ 9.80602e+0;5.68784e-1
21
+ 9.94483e+0;8.69448e-5
22
+ 1.00856e+1;5.87821e-5
23
+ 1.03732e+1;2.16172e-5
24
+ 1.03732e+1;3.41303e-5
25
+ 1.09731e+1;1.28279e-5
26
+ 1.11284e+1;5.25926e-6
27
+ 1.11284e+1;8.67272e-6
28
+ 1.17720e+1;2.73912e-6
29
+ 1.26291e+1;1.66119e-6
30
+ 1.28079e+1;9.03593e-7
31
+ 1.41322e+1;4.40909e-7
32
+ 1.49495e+1;2.61640e-7
33
+ 1.64952e+1;1.30475e-7
34
+ 1.79467e+1;6.94503e-8
35
+ 1.92534e+1;4.39918e-8
36
+ 2.06552e+1;2.78657e-8
37
+ 2.24728e+1;1.80394e-8
38
+ 2.41090e+1;1.09403e-8
39
+ 2.66017e+1;7.23827e-9
40
+ 2.97677e+1;5.00191e-9
41
+ 3.28456e+1;3.93821e-9
42
+ 3.67547e+1;3.85405e-9
43
+ 4.05550e+1;4.11441e-9
44
+ 4.47482e+1;4.79157e-9
45
+ 5.15015e+1;6.22150e-9
46
+ 6.01131e+1;8.43744e-9
47
+ 6.44899e+1;1.04881e-8
48
+ 7.74194e+1;1.76798e-8
49
+ 8.78598e+1;2.67302e-8
50
+ 1.10017e+2;5.84983e-8
51
+ 1.41691e+2;1.36658e-7
52
+ 1.67724e+2;3.12345e-7
53
+ 2.10023e+2;6.83559e-7
54
+ 2.48611e+2;1.49584e-6
55
+ 2.90182e+2;2.93605e-6
56
+ 3.20185e+2;5.05754e-6
57
+ 3.89820e+2;1.13112e-5
58
+ 4.36214e+2;2.17229e-5
59
+ 4.95039e+2;3.82433e-5
60
+ 5.61798e+2;6.58792e-5
61
+ 6.74432e+2;1.60727e-4
62
+ 7.13438e+2;2.59363e-4
63
+ 7.87204e+2;4.46771e-4
64
+ 8.93362e+2;7.21013e-4
65
+ 9.71973e+2;1.24197e-3
66
+ 1.08765e+3;2.09340e-3
67
+ 1.15056e+3;3.37808e-3
68
+ 1.26952e+3;5.33417e-3
69
+ 1.38123e+3;8.79722e-3
70
+ 1.48180e+3;1.35920e-2
71
+ 1.68162e+3;2.44548e-2
72
+ 1.77888e+3;4.03299e-2
73
+ 1.99059e+3;6.36843e-2
74
+ 2.07633e+3;9.21761e-2
75
+ 2.25903e+3;1.45549e-1
76
+ 2.38968e+3;2.20035e-1
77
+ 2.56367e+3;3.11637e-1
78
+ 2.71194e+3;4.51069e-1
79
+ 2.86879e+3;6.52885e-1
80
+ 3.03470e+3;9.04776e-1
81
+ 3.16541e+3;1.22685e+0
82
+ 3.30174e+3;1.81478e+0
@@ -0,0 +1,83 @@
1
+ 9.75860341854095; 1.6969370409948528
2
+ 10.051183547053066; 0.15285381473707907
3
+ 10.83272053313637; 0.05526813235374402
4
+ 11.141491225276706; 0.03647712975797426
5
+ 11.783860105428879; 0.019699297540412455
6
+ 12.345732680414734; 0.014168763342078727
7
+ 12.574508391289518; 0.018076543993131453
8
+ 13.413631374423405; 0.026614874390942678
9
+ 13.541889986582527; 0.018870485383875316
10
+ 13.798818393836958; 0.01318922160245812
11
+ 14.060330903812087; 0.009486376819685782
12
+ 14.063236818468965; 0.007122775349846694
13
+ 15.009303029403045; 0.0051970035062758714
14
+ 16.16650360439359; 0.00437605059909336
15
+ 16.16834118294355; 0.003737956270639411
16
+ 16.628508576095143; 0.002612584286564685
17
+ 16.943998971284135; 0.0018260236771676065
18
+ 18.253191295287216; 0.001240216160631325
19
+ 19.123926730364506; 0.0008668290955299814
20
+ 20.410858957034385; 0.0006145993979623295
21
+ 21.38385915330521; 0.0004484312742271195
22
+ 21.988829926422422; 0.00039417845883177343
23
+ 23.89998634780728; 0.00043576343009657585
24
+ 25.03595256196393; 0.00038304321567839383
25
+ 25.04035064706042; 0.00030023707379011027
26
+ 25.753556630743166; 0.00020391787663604123
27
+ 27.741706904536795; 0.0001509321016006575
28
+ 29.332682859076666; 0.00012174279711467096
29
+ 32.78681677781824; 0.00010549155967322482
30
+ 35.97090913439195; 0.00010549155967322482
31
+ 39.82716831724237; 0.00012349967990084068
32
+ 43.28561079989749; 0.0001509321016006575
33
+ 45.33079590844245; 0.00019255955574170756
34
+ 47.91366364420047; 0.00025280964791396564
35
+ 49.71100530490953; 0.0003565619477810621
36
+ 53.52652057144452; 0.00046812686148363234
37
+ 56.55824682550625; 0.0009582776823856468
38
+ 56.57168956629359; 0.0006892433433153593
39
+ 58.67986486813849; 0.001351551888015081
40
+ 59.19811796353287; 0.002612584286564685
41
+ 59.20974091903139; 0.0019899480645211177
42
+ 61.425109614051756; 0.003192905741105303
43
+ 61.998349936779164; 0.0031027085702644646
44
+ 68.03478473065273; 0.0022638348855076127
45
+ 71.2743690502429; 0.0017744398901962187
46
+ 71.94918996375488; 0.0014312745483124989
47
+ 76.07230345994304; 0.0012225730831311066
48
+ 83.46610348143551; 0.0011059027454667885
49
+ 84.22763024587383; 0.0014312745483124989
50
+ 87.39711762255274; 0.0017243133061688968
51
+ 93.24761289343643; 0.0019337335376626208
52
+ 97.66049674291227; 0.002231630004132432
53
+ 100.39539750171778; 0.0028882062543547486
54
+ 105.1204742725507; 0.004701052661057711
55
+ 105.13676815841596; 0.0037918991007943347
56
+ 107.07396736958466; 0.005582976386919792
57
+ 108.05312013400646; 0.007021448112302839
58
+ 110.06452341353835; 0.007987846249262798
59
+ 113.14561949176104; 0.010487168970081005
60
+ 119.61717141135327; 0.01033798051344707
61
+ 127.6509025868699; 0.008704925013934384
62
+ 133.69466301077483; 0.009762150174252833
63
+ 138.72702921257823; 0.011593542518657064
64
+ 142.61049040783; 0.015221054031453877
65
+ 149.33937966376334; 0.02116233188494443
66
+ 152.0988925367621; 0.029004124340753366
67
+ 169.9639618607662; 0.03647712975797426
68
+ 181.3284908717801; 0.0452229617292691
69
+ 182.99234674756642; 0.05448189849116219
70
+ 189.86461918944724; 0.07256101639512741
71
+ 198.83751377041122; 0.0912565943835976
72
+ 206.32190515459232; 0.10837645274956616
73
+ 228.4406327742625; 0.1268770436688991
74
+ 234.85005499106734; 0.15285381473707907
75
+ 239.16987455007657; 0.2349378385561308
76
+ 239.20694636377695; 0.18950236106547225
77
+ 277.2532124732199; 0.49490961722250765
78
+ 277.2990525052609; 0.3935183543935018
79
+ 277.319109902453; 0.35596483721102606
80
+ 279.8059392902582; 0.5711516760796542
81
+ 304.09695723447396; 0.7183107324694664
82
+ 318.35325638133065; 1.4916355437224058
83
+ 318.43879326532186; 1.0277234386310978
@@ -0,0 +1,61 @@
1
+ 9.56020e+0;3.87251e-1
2
+ 9.56020e+0;2.11932e-1
3
+ 9.56020e+0;6.64881e-2
4
+ 9.56020e+0;3.63872e-2
5
+ 9.56020e+0;1.99137e-2
6
+ 9.56020e+0;6.10422e-3
7
+ 9.56020e+0;3.34068e-3
8
+ 9.56020e+0;1.87115e-3
9
+ 9.56020e+0;1.04805e-3
10
+ 9.56020e+0;3.13899e-4
11
+ 9.56020e+0;2.48610e-7
12
+ 9.56020e+0;9.22712e-8
13
+ 9.61494e+0;2.97431e-6
14
+ 9.77849e+0;1.02026e-6
15
+ 9.99978e+0;3.89062e-8
16
+ 1.01508e+1;1.07831e+0
17
+ 1.01508e+1;4.71896e-5
18
+ 1.03556e+1;1.01580e-4
19
+ 1.04596e+1;1.08833e-1
20
+ 1.04596e+1;1.05270e-2
21
+ 1.04596e+1;1.98628e-8
22
+ 1.06175e+1;1.06215e-8
23
+ 1.06494e+1;2.14217e-5
24
+ 1.06974e+1;5.77932e-9
25
+ 1.11057e+1;1.32926e-5
26
+ 1.13814e+1;2.40826e-9
27
+ 1.34943e+1;1.15669e-9
28
+ 1.76723e+1;6.36927e-10
29
+ 2.42080e+1;4.76775e-10
30
+ 3.36615e+1;5.59846e-10
31
+ 4.47490e+1;8.85307e-10
32
+ 5.51943e+1;1.41128e-9
33
+ 6.70653e+1;2.28600e-9
34
+ 8.27196e+1;3.95039e-9
35
+ 9.92858e+1;6.85834e-9
36
+ 1.16954e+2;1.21830e-8
37
+ 1.36655e+2;2.05107e-8
38
+ 1.57561e+2;3.52169e-8
39
+ 1.93371e+2;6.90160e-8
40
+ 2.16642e+2;1.15242e-7
41
+ 2.57064e+2;1.93475e-7
42
+ 3.04648e+2;3.33691e-7
43
+ 3.48634e+2;5.64050e-7
44
+ 3.98971e+2;9.69783e-7
45
+ 4.49785e+2;1.58749e-6
46
+ 4.99528e+2;2.44279e-6
47
+ 5.59411e+2;4.76119e-6
48
+ 6.74432e+2;9.42940e-6
49
+ 7.84994e+2;1.95847e-5
50
+ 9.29223e+2;4.51537e-5
51
+ 1.23780e+3;1.58081e-4
52
+ 1.54128e+3;5.25240e-4
53
+ 1.91916e+3;1.41634e-3
54
+ 2.38968e+3;5.08917e-3
55
+ 2.97557e+3;1.73564e-2
56
+ 3.34848e+3;3.42141e-2
57
+ 3.41096e+3;6.54647e-2
58
+ 3.68738e+3;1.36555e-1
59
+ 3.92103e+3;2.54245e-1
60
+ 4.08818e+3;4.01693e-1
61
+ 4.36557e+3;6.78335e-1
@@ -0,0 +1,83 @@
1
+ 9.316977590875046; 1.6255414978000047
2
+ 9.721668002532182; 0.0008544977508843073
3
+ 9.81614604423656; 0.0004886875002645911
4
+ 9.819392243936585; 0.00030896510416232495
5
+ 9.822132059936079; 0.00020984586347074256
6
+ 9.825075676778685; 0.0001384988865193285
7
+ 10.011071706948115; 0.00010251150294695928
8
+ 10.105960278877467; 0.00008151015163634493
9
+ 10.48982395186064; 0.00006120129000492882
10
+ 10.788152040579694; 0.00004401916324837686
11
+ 11.197349502307846; 0.00003550613821494273
12
+ 12.286430579516738; 0.000029472040698232353
13
+ 14.519191862457241; 0.000024115397279013237
14
+ 16.229266284005394; 0.000020305970273593715
15
+ 19.177749848026874; 0.000017595361085895403
16
+ 23.29893993029483; 0.00001709830545113642
17
+ 26.772963935602785; 0.000018106866357476684
18
+ 31.926303195834087; 0.000019732341980632513
19
+ 35.67774338351009; 0.000023434155556068795
20
+ 39.86998949210312; 0.00002783043427713925
21
+ 43.73699290596427; 0.00003258127697964509
22
+ 45.800192406487696; 0.00004595253257506219
23
+ 48.41180636679338; 0.00005697020669645024
24
+ 50.700768708188654; 0.00006962473936955147
25
+ 53.0952125074767; 0.00009140967207947027
26
+ 55.08977094308996; 0.00012001090740643299
27
+ 59.32180523826773; 0.00014458190762233448
28
+ 64.06164312501953; 0.0002028687624940692
29
+ 64.9127138006202; 0.0002493560959704688
30
+ 68.46518363864827; 0.0002881000652791707
31
+ 77.21031937616974; 0.00025455431942535383
32
+ 81.43701040714382; 0.0002881000652791707
33
+ 85.89252531748875; 0.0003398029966783695
34
+ 90.58911090646768; 0.0004176687805904151
35
+ 98.12424348308745; 0.0005350049326133014
36
+ 100.76075223656943; 0.0006853044594713297
37
+ 103.46656207628307; 0.0008961275117209026
38
+ 109.11759792063536; 0.0011962351083890696
39
+ 115.08412569477808; 0.001470352127274545
40
+ 116.61650963905149; 0.0017342240396098494
41
+ 124.66044436380521; 0.001807283002329019
42
+ 135.0456780990428; 0.0019627640188868136
43
+ 144.34148360231288; 0.0024628240093637495
44
+ 150.19787361743724; 0.003570442442638777
45
+ 150.22022328302634; 0.0029048066263445693
46
+ 154.23814298912765; 0.004388608563675243
47
+ 162.66217155552715; 0.005858326602163601
48
+ 162.68153476464917; 0.004966949358917984
49
+ 169.29211391535247; 0.006630350175633007
50
+ 173.85635851433418; 0.007504112767506133
51
+ 183.38734438625065; 0.00766054790435064
52
+ 198.63545155825318; 0.01022601832085543
53
+ 198.65909700544466; 0.00867007229066023
54
+ 206.6947193870408; 0.014825017822291073
55
+ 206.72547586678448; 0.012061197091876545
56
+ 220.89615747672934; 0.021940390995292013
57
+ 220.93231438791926; 0.017485540924620364
58
+ 249.04847157826603; 0.027530218201513416
59
+ 249.08553035893544; 0.02239777325674107
60
+ 252.35337138599272; 0.03454418448538926
61
+ 262.6390345527702; 0.0390964955052541
62
+ 277.02893482508296; 0.04159290682665848
63
+ 292.18986213202027; 0.04805544907684919
64
+ 300.0407367155544; 0.061555719483554025
65
+ 308.13006397573014; 0.06966767188061272
66
+ 312.24685718061846; 0.07723847716738635
67
+ 320.6604996083827; 0.08923948791940949
68
+ 338.21427828243816; 0.10099967014152711
69
+ 356.73430367493853; 0.11197533239675822
70
+ 371.18604919690813; 0.1763002555859836
71
+ 371.24680598493495; 0.14050366443930912
72
+ 396.7836929492636; 0.1875574782174132
73
+ 412.92546625554525; 0.23534212775692878
74
+ 429.7111242555967; 0.3077414022850713
75
+ 429.7494881041687; 0.2719086618108586
76
+ 447.21905228103327; 0.3555570508207372
77
+ 465.43338902610554; 0.4193659284385621
78
+ 471.63079156244333; 0.4946260565756724
79
+ 490.8247633597129; 0.6079695115041858
80
+ 524.6104542504156; 0.762863948740444
81
+ 545.9849115450015; 0.8813947482450022
82
+ 568.1795126227466; 1.1525418690518685
83
+ 583.4372958553371; 1.5071031028520836
@@ -0,0 +1,31 @@
1
+ f [Hz], h_c [1]
2
+ 1.9870807315792565e-7, 1.7068314890414264e-13
3
+ 1.7008574091943525e-7, 1.4231275542136064e-13
4
+ 1.3823067594099873e-7, 1.1378329644281288e-13
5
+ 1.1307209550114642e-7, 8.846419531043188e-14
6
+ 9.615975570312793e-8, 7.172574170222732e-14
7
+ 8.230870057503779e-8, 5.98037240455262e-14
8
+ 7.183590878076275e-8, 4.9863345081342994e-14
9
+ 6.269565373169187e-8, 4.335638732674537e-14
10
+ 5.471838616033257e-8, 3.717508180198244e-14
11
+ 4.389760621733326e-8, 2.972258769382155e-14
12
+ 3.567609928270252e-8, 2.40987287587688e-14
13
+ 2.78886176288316e-8, 1.9267653065498593e-14
14
+ 2.356407600360698e-8, 1.6753305014381723e-14
15
+ 1.940063332419105e-8, 1.4772194095210424e-14
16
+ 1.5869634211670776e-8, 1.2844484617167135e-14
17
+ 1.2567402012848246e-8, 1.132559931636509e-14
18
+ 9.634999158060278e-9, 1.012694695054068e-14
19
+ 7.1978024024725115e-9, 9.84765602602091e-15
20
+ 6.042391536892428e-9, 9.84765602602091e-15
21
+ 4.785060736139321e-9, 1.0269548899696928e-14
22
+ 3.939607424522497e-9, 1.1013250775461587e-14
23
+ 3.1810829150682037e-9, 1.2316808418440966e-14
24
+ 2.6704468127047355e-9, 1.3968625575042362e-14
25
+ 2.271024623624146e-9, 1.5841969268888387e-14
26
+ 1.943901420896981e-9, 1.7966548603379178e-14
27
+ 1.5901029609541676e-9, 2.247139182194254e-14
28
+ 1.3176656438745895e-9, 2.850152647120018e-14
29
+ 1.1205806866124665e-9, 3.5152855091568895e-14
30
+ 9.780002772804193e-10, 4.275434478314857e-14
31
+ 8.0520130590184e-10, 5.6550558691613935e-14