diff --git a/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/flake.lock b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/flake.lock index 2d978fc..258136d 100644 --- a/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/flake.lock +++ b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/flake.lock @@ -32,11 +32,11 @@ }, "nixpkgs": { "locked": { - "lastModified": 1649012074, - "narHash": "sha256-vQUROcJ3FfT3GTB/nJrXwVvjuq8WfK0ImN+RUgDVN1c=", + "lastModified": 1651007983, + "narHash": "sha256-GNay7yDPtLcRcKCNHldug85AhAvBpTtPEJWSSDYBw8U=", "owner": "NixOS", "repo": "nixpkgs", - "rev": "bc4b9eef3ce3d5a90d8693e8367c9cbfc9fc1e13", + "rev": "e10da1c7f542515b609f8dfbcf788f3d85b14936", "type": "github" }, "original": { diff --git a/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/flake.nix b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/flake.nix index 2249534..eb3b862 100644 --- a/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/flake.nix +++ b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/flake.nix @@ -16,7 +16,7 @@ export PYTHONPATH=/home/hiro/src/two_qubit_model/:$PYTHONPATH export PYTHONPATH=/home/hiro/src/hops/:$PYTHONPATH export PYTHONPATH=/home/hiro/src/hopsflow/:$PYTHONPATH -# export PYTHONPATH=/home/hiro/src/stocproc/:$PYTHONPATH + export PYTHONPATH=/home/hiro/src/stocproc/:$PYTHONPATH ''; }); noPackage = true; diff --git a/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/poetry.lock b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/poetry.lock index 875b249..946679f 100644 --- a/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/poetry.lock +++ b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/poetry.lock @@ -356,7 +356,7 @@ testing = ["covdefaults (>=1.2.0)", "coverage (>=4)", "pytest (>=4)", "pytest-co [[package]] name = "fonttools" -version = "4.33.2" +version = "4.33.3" description = "Tools to manipulate font files" category = "main" optional = false @@ -469,7 +469,7 @@ typer = "^0.4.0" type = "git" url = "git@gitlab.hrz.tu-chemnitz.de:s8896854--tu-dresden.de/hops.git" reference = "time_dep_coupling" -resolved_reference = "02b50c653956165d3a7bf9336ae9f229c7306503" +resolved_reference = "8eecf5a2969a2ffc931c6e6c24f9cd20feaecd6f" [[package]] name = "hopsflow" @@ -665,7 +665,7 @@ qtconsole = "*" [[package]] name = "jupyter-client" -version = "7.2.2" +version = "7.3.0" description = "Jupyter protocol implementation and client libraries" category = "main" optional = false @@ -793,6 +793,17 @@ python-versions = ">=3.5" [package.dependencies] traitlets = "*" +[[package]] +name = "memory-profiler" +version = "0.60.0" +description = "A module for monitoring memory usage of a python program" +category = "main" +optional = false +python-versions = ">=3.4" + +[package.dependencies] +psutil = "*" + [[package]] name = "mistune" version = "0.8.4" @@ -1402,7 +1413,7 @@ tests = ["pytest", "typeguard", "pygments", "littleutils", "cython"] [[package]] name = "stocproc" -version = "1.0.5" +version = "1.0.6" description = "Generate continuous time stationary stochastic processes from a given auto correlation function." category = "main" optional = false @@ -1419,7 +1430,7 @@ scipy = "^1.7.3" type = "git" url = "https://github.com/vale981/stocproc" reference = "master" -resolved_reference = "20dee534d76db6f9b9656250dd594db36c5542aa" +resolved_reference = "9e32e3d6ae1740f19e1d0905e46389a2f20e7a06" [[package]] name = "terminado" @@ -1599,7 +1610,7 @@ notebook = ">=4.4.1" [metadata] lock-version = "1.1" python-versions = ">=3.9,<3.11" -content-hash = "a5b3a21338348e199692c8e61c50d50eede495e99132d5f9ac7126f75d016146" +content-hash = "3461211ab361262c64343efb4c8331db458ceda1f339f423c44da8ea884afb2d" [metadata.files] aiosignal = [ @@ -1799,8 +1810,8 @@ filelock = [ {file = "filelock-3.6.0.tar.gz", hash = "sha256:9cd540a9352e432c7246a48fe4e8712b10acb1df2ad1f30e8c070b82ae1fed85"}, ] fonttools = [ - {file = "fonttools-4.33.2-py3-none-any.whl", hash = "sha256:b4da40696829845ea8c1cb33ce51c552179754cbee7ab4e8b96a6bcf421f437a"}, - {file = "fonttools-4.33.2.zip", hash = "sha256:696fe922a271584c3ec8325ba31d4001a4fd6c4953b22900b767f1cb53ce1044"}, + {file = "fonttools-4.33.3-py3-none-any.whl", hash = "sha256:f829c579a8678fa939a1d9e9894d01941db869de44390adb49ce67055a06cc2a"}, + {file = "fonttools-4.33.3.zip", hash = "sha256:c0fdcfa8ceebd7c1b2021240bd46ef77aa8e7408cf10434be55df52384865f8e"}, ] frozenlist = [ {file = "frozenlist-1.3.0-cp310-cp310-macosx_10_9_universal2.whl", hash = "sha256:d2257aaba9660f78c7b1d8fea963b68f3feffb1a9d5d05a18401ca9eb3e8d0a3"}, @@ -1975,8 +1986,8 @@ jupyter = [ {file = "jupyter-1.0.0.zip", hash = "sha256:3e1f86076bbb7c8c207829390305a2b1fe836d471ed54be66a3b8c41e7f46cc7"}, ] jupyter-client = [ - {file = "jupyter_client-7.2.2-py3-none-any.whl", hash = "sha256:44045448eadc12493d819d965eb1dc9d10d1927698adbb9b14eb9a3a4a45ba53"}, - {file = "jupyter_client-7.2.2.tar.gz", hash = "sha256:8fdbad344a8baa6a413d86d25bbf87ce21cb2b4aa5a8e0413863b9754eb8eb8a"}, + {file = "jupyter_client-7.3.0-py3-none-any.whl", hash = "sha256:671dd2d90d03f41716b09627a4eb06bb37875f92bf6563cc2ce4fe71c61c5cda"}, + {file = "jupyter_client-7.3.0.tar.gz", hash = "sha256:3bcc8e08a294d0fa9406e48cfe17e11ef0efdb7c504fe8cc335128e3ef8f3dac"}, ] jupyter-console = [ {file = "jupyter_console-6.4.3-py3-none-any.whl", hash = "sha256:e630bcb682c0088dda45688ad7c2424d4a825c8acf494cb036ced03ed0424841"}, @@ -2125,6 +2136,9 @@ matplotlib-inline = [ {file = "matplotlib-inline-0.1.3.tar.gz", hash = "sha256:a04bfba22e0d1395479f866853ec1ee28eea1485c1d69a6faf00dc3e24ff34ee"}, {file = "matplotlib_inline-0.1.3-py3-none-any.whl", hash = "sha256:aed605ba3b72462d64d475a21a9296f400a19c4f74a31b59103d2a99ffd5aa5c"}, ] +memory-profiler = [ + {file = "memory_profiler-0.60.0.tar.gz", hash = "sha256:6a12869511d6cebcb29b71ba26985675a58e16e06b3c523b49f67c5497a33d1c"}, +] mistune = [ {file = "mistune-0.8.4-py2.py3-none-any.whl", hash = "sha256:88a1051873018da288eee8538d476dffe1262495144b33ecb586c4ab266bb8d4"}, {file = "mistune-0.8.4.tar.gz", hash = "sha256:59a3429db53c50b5c6bcc8a07f8848cb00d7dc8bdb431a4ab41920d201d4756e"}, diff --git a/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/profiling.org b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/profiling.org new file mode 100644 index 0000000..7e58846 --- /dev/null +++ b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/profiling.org @@ -0,0 +1,1162 @@ +#+PROPERTY: header-args :session 09_one_qubit_otto_prof :kernel python :pandoc no :async yes :tangle no +This is a clone of an old version of [[file:first_experiments.org]] +wherein I try to find the bottleneck in HOPS. + +* Boilerplate +#+begin_src jupyter-python :results none + import figsaver as fs + from hiro_models.one_qubit_model import QubitModelMutliBath, StocProcTolerances + import hiro_models.model_auxiliary as aux + import numpy as np + import qutip as qt +#+end_src + +Init ray and silence stocproc. +#+begin_src jupyter-python + import ray + ray.shutdown() + #ray.init(address='auto') + ray.init() +#+end_src + +#+RESULTS: +: RayContext(dashboard_url='', python_version='3.9.12', ray_version='1.12.0', ray_commit='f18fc31c7562990955556899090f8e8656b48d2d', address_info={'node_ip_address': '141.30.17.221', 'raylet_ip_address': '141.30.17.221', 'redis_address': None, 'object_store_address': '/tmp/ray/session_2022-04-28_15-07-59_498396_50472/sockets/plasma_store', 'raylet_socket_name': '/tmp/ray/session_2022-04-28_15-07-59_498396_50472/sockets/raylet', 'webui_url': '', 'session_dir': '/tmp/ray/session_2022-04-28_15-07-59_498396_50472', 'metrics_export_port': 62867, 'gcs_address': '141.30.17.221:60066', 'address': '141.30.17.221:60066', 'node_id': '723796a8b4111e83bbf569c493d3caf14eb7e578a2ef834a9c150003'}) + +#+begin_src jupyter-python :results none + from hops.util.logging_setup import logging_setup + import logging + logging_setup(logging.INFO) +#+end_src + +* Cycle +#+begin_src jupyter-python :results none + from hops.util.dynamic_matrix import SmoothStep, Periodic, ConstantMatrix, ScaleTime, Shift, Piecewise +#+end_src + +Now we build a basic otto cycle. +Let's lay down some basic functionality. + +#+begin_src jupyter-python + H_op = (1/2 * (qt.sigmaz() + qt.identity(2))).full() + L_op = (1/2 * qt.sigmax()).full() + print(H_op, L_op, sep="\n") +#+end_src + +#+RESULTS: +: [[1.+0.j 0.+0.j] +: [0.+0.j 0.+0.j]] +: [[0. +0.j 0.5+0.j] +: [0.5+0.j 0. +0.j]] + +#+begin_src jupyter-python :results none + comp_ratio = .1 + t_exp_end = .1 + t_hot_modulation = .1 + t_hot_end = .5 + t_comp_end = .6 + t_cold_modulation = .1 + t_cold_end = 1 + t = np.linspace(0, 100, 1000) +#+end_src + +** Expansion +#+begin_src jupyter-python + %matplotlib inline + H_exp = SmoothStep(comp_ratio * H_op, 0, t_exp_end, 3) + H_op * (1 - comp_ratio) + plt.plot(t, H_exp.operator_norm(t)) +#+end_src + +#+RESULTS: +:RESULTS: +| | +[[file:./.ob-jupyter/e830a0d011a3a50eedfd9290cca26345db6694b4.svg]] +:END: + +** Hot Thermalization +#+begin_src jupyter-python + L_hot = SmoothStep(L_op, t_exp_end, t_exp_end + t_hot_modulation, 3) - SmoothStep( + L_op, t_hot_end - t_hot_modulation, t_hot_end, 3 + ) + plt.plot(t, L_hot.operator_norm(t)) +#+end_src + +#+RESULTS: +:RESULTS: +| | +[[file:./.ob-jupyter/1053a13025f11b7e28b671c0a1c610f9cd739b30.svg]] +:END: + +** Compression +#+begin_src jupyter-python + H_comp = ConstantMatrix(H_op) - SmoothStep(comp_ratio * H_op, t_hot_end, t_comp_end, 3) + plt.plot(t, H_comp.operator_norm(t)) +#+end_src + +#+RESULTS: +:RESULTS: +| | +[[file:./.ob-jupyter/2c38f9764b31aae0b6d675acac18596ead77d991.svg]] +:END: + +** Cold Thermalization +#+begin_src jupyter-python + L_cold = SmoothStep(L_op, t_comp_end, t_comp_end + t_cold_modulation, 3) - SmoothStep( + L_op, t_cold_end - t_cold_modulation, t_cold_end, 3 + ) + plt.plot(t, L_cold.operator_norm(t)) +#+end_src + +#+RESULTS: +:RESULTS: +| | +[[file:./.ob-jupyter/4490857649e12acaac00a4d4b0970fa9786ea89d.svg]] +:END: + +** Full Cycle +Now we turn the system around after each fill-cycle. +#+begin_src jupyter-python + H_cyc = ScaleTime(Periodic(Piecewise([H_exp, H_comp], [0, t_exp_end, 1]), 1), 0.01 * 2) + L = [ScaleTime(Periodic(L_i, 1), 0.01 * 2) for L_i in (L_cold, L_hot)] + plt.plot(t, H_cyc.operator_norm(t)) + plt.plot(t, L[0].operator_norm(t)) + plt.plot(t, L[1].operator_norm(t)) +#+end_src + +#+RESULTS: +:RESULTS: +| | +[[file:./.ob-jupyter/914d0769bab50fe57b96b80ba867e3b9224c9ae4.svg]] +:END: + +* Model +#+begin_src jupyter-python :results none + model = QubitModelMutliBath( + δ=[.01, .01], + ω_c=[2] * 2, + t=np.linspace(0, 800, 10), + ψ_0=(qt.basis([2], [0]) + qt.basis([2], [1])) / np.sqrt(2), + description=f"A simple otto cycle.", + k_max=1, + bcf_terms=[4, 4], + truncation_scheme="simplex", + driving_process_tolerances=[StocProcTolerances(1e-3, 1e-3)] * 2, + thermal_process_tolerances=[StocProcTolerances(1e-3, 1e-3)] * 2, + T = [0, 4], + L = L, + H = H_cyc, + ) +#+end_src + + +#+begin_src jupyter-python + config = model.hops_config +#+end_src + +#+RESULTS: +#+begin_example + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + [INFO stocproc.stocproc 50472] non neg freq only + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + [INFO stocproc.method_ft 50472] get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + [INFO stocproc.method_ft 50472] acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + [INFO stocproc.method_ft 50472] acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + [INFO stocproc.method_ft 50472] acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + [INFO stocproc.method_ft 50472] acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + [INFO stocproc.method_ft 50472] acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + [INFO stocproc.method_ft 50472] acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + [INFO stocproc.method_ft 50472] acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + [INFO stocproc.method_ft 50472] requires dt < 1.543e-02 + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + [INFO stocproc.method_ft 50472] get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + [INFO stocproc.method_ft 50472] increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + [INFO stocproc.method_ft 50472] increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + [INFO stocproc.method_ft 50472] increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + [INFO stocproc.method_ft 50472] return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + [INFO stocproc.method_ft 50472] requires dx < 9.706e-02 + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + [INFO stocproc.method_ft 50472] increase N to match dt_new*(N-1) < t_max + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + [INFO stocproc.stocproc 50472] Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + [INFO stocproc.stocproc 50472] Number of Nodes : 131072 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + [INFO stocproc.stocproc 50472] yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + [INFO stocproc.stocproc 50472] yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + [INFO stocproc.stocproc 50472] yields t_max : 1.448e+03 + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + [INFO stocproc.stocproc 50472] non neg freq only + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - get_dt_for_accurate_interpolation, please wait ... + [INFO stocproc.method_ft 50472] get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + [INFO stocproc.method_ft 50472] acc interp N 33 dt 9.88e-01 -> diff 5.01e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + [INFO stocproc.method_ft 50472] acc interp N 65 dt 4.94e-01 -> diff 3.02e-01 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + [INFO stocproc.method_ft 50472] acc interp N 129 dt 2.47e-01 -> diff 9.52e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + [INFO stocproc.method_ft 50472] acc interp N 257 dt 1.23e-01 -> diff 1.98e-02 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + [INFO stocproc.method_ft 50472] acc interp N 513 dt 6.17e-02 -> diff 4.37e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + [INFO stocproc.method_ft 50472] acc interp N 1025 dt 3.09e-02 -> diff 1.06e-03 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + [INFO stocproc.method_ft 50472] acc interp N 2049 dt 1.54e-02 -> diff 2.62e-04 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + stocproc.method_ft - INFO - requires dt < 1.543e-02 + [INFO stocproc.method_ft 50472] requires dt < 1.543e-02 + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - get_N_a_b_for_accurate_fourier_integral, please wait ... + [INFO stocproc.method_ft 50472] get_N_a_b_for_accurate_fourier_integral, please wait ... + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 32 yields: interval [0.00e+00,1.46e+01] diff 2.07e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-03 N 32 yields: interval [0.00e+00,1.98e+01] diff 3.82e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 64 yields: interval [0.00e+00,1.46e+01] diff 2.06e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-04 N 32 yields: interval [0.00e+00,2.48e+01] diff 5.99e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-03 N 64 yields: interval [0.00e+00,1.98e+01] diff 9.64e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 128 yields: interval [0.00e+00,1.46e+01] diff 2.22e-02 + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + [INFO stocproc.method_ft 50472] increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-05 N 32 yields: interval [0.00e+00,2.98e+01] diff 8.56e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-04 N 64 yields: interval [0.00e+00,2.48e+01] diff 1.52e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-03 N 128 yields: interval [0.00e+00,1.98e+01] diff 2.42e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 256 yields: interval [0.00e+00,1.46e+01] diff 2.26e-02 + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + [INFO stocproc.method_ft 50472] increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-06 N 32 yields: interval [0.00e+00,3.47e+01] diff 1.15e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-05 N 64 yields: interval [0.00e+00,2.98e+01] diff 2.18e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-04 N 128 yields: interval [0.00e+00,2.48e+01] diff 3.81e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-03 N 256 yields: interval [0.00e+00,1.98e+01] diff 1.95e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-02 N 512 yields: interval [0.00e+00,1.46e+01] diff 2.27e-02 + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - increasing N while shrinking the interval does lower the error -> try next level + [INFO stocproc.method_ft 50472] increasing N while shrinking the interval does lower the error -> try next level + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-07 N 32 yields: interval [0.00e+00,3.96e+01] diff 1.48e-01 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-06 N 64 yields: interval [0.00e+00,3.47e+01] diff 2.96e-02 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-05 N 128 yields: interval [0.00e+00,2.98e+01] diff 5.49e-03 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + [INFO stocproc.method_ft 50472] J_w_min:1.00e-04 N 256 yields: interval [0.00e+00,2.48e+01] diff 9.54e-04 + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - return, cause tol of 0.001 was reached + [INFO stocproc.method_ft 50472] return, cause tol of 0.001 was reached + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + stocproc.method_ft - INFO - requires dx < 9.706e-02 + [INFO stocproc.method_ft 50472] requires dx < 9.706e-02 + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + stocproc.method_ft - INFO - increase N to match dt_new*(N-1) < t_max + [INFO stocproc.method_ft 50472] increase N to match dt_new*(N-1) < t_max + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + [INFO stocproc.stocproc 50472] Fourier Integral Boundaries: [0.000e+00, 5.689e+02] + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + stocproc.stocproc - INFO - Number of Nodes : 131072 + [INFO stocproc.stocproc 50472] Number of Nodes : 131072 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dx : 4.341e-03 + [INFO stocproc.stocproc 50472] yields dx : 4.341e-03 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields dt : 1.104e-02 + [INFO stocproc.stocproc 50472] yields dt : 1.104e-02 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + stocproc.stocproc - INFO - yields t_max : 1.448e+03 + [INFO stocproc.stocproc 50472] yields t_max : 1.448e+03 + /home/hiro/src/hops/hops/util/bcf.py:280: UserWarning: this implementation uses mpmath to evaluate the zeta_function! for a better performance consider the 'OhmEnv' package + warnings.warn( + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + stocproc.stocproc - INFO - non neg freq only + [INFO stocproc.stocproc 50472] non neg freq only + stocproc.stocproc - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.stocproc - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.stocproc - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.stocproc - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.stocproc - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.stocproc - INFO - get_dt_for_accurate_interpolation, please wait ... + stocproc.stocproc - INFO - get_dt_for_accurate_interpolation, please wait ... + [INFO stocproc.stocproc 50472] get_dt_for_accurate_interpolation, please wait ... + stocproc.method_ft - INFO - acc interp N 33 dt 3.96e+01 -> diff 3.99e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 3.96e+01 -> diff 3.99e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 3.96e+01 -> diff 3.99e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 3.96e+01 -> diff 3.99e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 3.96e+01 -> diff 3.99e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 3.96e+01 -> diff 3.99e-01 + stocproc.method_ft - INFO - acc interp N 33 dt 3.96e+01 -> diff 3.99e-01 + [INFO stocproc.method_ft 50472] acc interp N 33 dt 3.96e+01 -> diff 3.99e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 1.98e+01 -> diff 3.98e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 1.98e+01 -> diff 3.98e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 1.98e+01 -> diff 3.98e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 1.98e+01 -> diff 3.98e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 1.98e+01 -> diff 3.98e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 1.98e+01 -> diff 3.98e-01 + stocproc.method_ft - INFO - acc interp N 65 dt 1.98e+01 -> diff 3.98e-01 + [INFO stocproc.method_ft 50472] acc interp N 65 dt 1.98e+01 -> diff 3.98e-01 + stocproc.method_ft - INFO - acc interp N 129 dt 9.89e+00 -> diff 3.95e-01 + stocproc.method_ft - INFO - acc interp N 129 dt 9.89e+00 -> diff 3.95e-01 + stocproc.method_ft - INFO - acc interp N 129 dt 9.89e+00 -> diff 3.95e-01 + stocproc.method_ft - INFO - acc interp N 129 dt 9.89e+00 -> diff 3.95e-01 + stocproc.method_ft - INFO - acc interp N 129 dt 9.89e+00 -> diff 3.95e-01 + stocproc.method_ft - INFO - acc interp N 129 dt 9.89e+00 -> diff 3.95e-01 + stocproc.method_ft - INFO - acc interp N 129 dt 9.89e+00 -> diff 3.95e-01 + [INFO stocproc.method_ft 50472] acc interp N 129 dt 9.89e+00 -> diff 3.95e-01 + stocproc.method_ft - INFO - acc interp N 257 dt 4.95e+00 -> diff 3.84e-01 + stocproc.method_ft - INFO - acc interp N 257 dt 4.95e+00 -> diff 3.84e-01 + stocproc.method_ft - INFO - acc interp N 257 dt 4.95e+00 -> diff 3.84e-01 + stocproc.method_ft - INFO - acc interp N 257 dt 4.95e+00 -> diff 3.84e-01 + stocproc.method_ft - INFO - acc interp N 257 dt 4.95e+00 -> diff 3.84e-01 + stocproc.method_ft - INFO - acc interp N 257 dt 4.95e+00 -> diff 3.84e-01 + stocproc.method_ft - INFO - acc interp N 257 dt 4.95e+00 -> diff 3.84e-01 + [INFO stocproc.method_ft 50472] acc interp N 257 dt 4.95e+00 -> diff 3.84e-01 + stocproc.method_ft - INFO - acc interp N 513 dt 2.47e+00 -> diff 3.43e-01 + stocproc.method_ft - INFO - acc interp N 513 dt 2.47e+00 -> diff 3.43e-01 + stocproc.method_ft - INFO - acc interp N 513 dt 2.47e+00 -> diff 3.43e-01 + stocproc.method_ft - INFO - acc interp N 513 dt 2.47e+00 -> diff 3.43e-01 + stocproc.method_ft - INFO - acc interp N 513 dt 2.47e+00 -> diff 3.43e-01 + stocproc.method_ft - INFO - acc interp N 513 dt 2.47e+00 -> diff 3.43e-01 + stocproc.method_ft - INFO - acc interp N 513 dt 2.47e+00 -> diff 3.43e-01 + [INFO stocproc.method_ft 50472] acc interp N 513 dt 2.47e+00 -> diff 3.43e-01 + stocproc.method_ft - INFO - acc interp N 1025 dt 1.24e+00 -> diff 2.40e-01 + stocproc.method_ft - INFO - acc interp N 1025 dt 1.24e+00 -> diff 2.40e-01 + stocproc.method_ft - INFO - acc interp N 1025 dt 1.24e+00 -> diff 2.40e-01 + stocproc.method_ft - INFO - acc interp N 1025 dt 1.24e+00 -> diff 2.40e-01 + stocproc.method_ft - INFO - acc interp N 1025 dt 1.24e+00 -> diff 2.40e-01 + stocproc.method_ft - INFO - acc interp N 1025 dt 1.24e+00 -> diff 2.40e-01 + stocproc.method_ft - INFO - acc interp N 1025 dt 1.24e+00 -> diff 2.40e-01 + [INFO stocproc.method_ft 50472] acc interp N 1025 dt 1.24e+00 -> diff 2.40e-01 + stocproc.method_ft - INFO - acc interp N 2049 dt 6.18e-01 -> diff 1.06e-01 + stocproc.method_ft - INFO - acc interp N 2049 dt 6.18e-01 -> diff 1.06e-01 + stocproc.method_ft - INFO - acc interp N 2049 dt 6.18e-01 -> diff 1.06e-01 + stocproc.method_ft - INFO - acc interp N 2049 dt 6.18e-01 -> diff 1.06e-01 + stocproc.method_ft - INFO - acc interp N 2049 dt 6.18e-01 -> diff 1.06e-01 + stocproc.method_ft - INFO - acc interp N 2049 dt 6.18e-01 -> diff 1.06e-01 + stocproc.method_ft - INFO - acc interp N 2049 dt 6.18e-01 -> diff 1.06e-01 + [INFO stocproc.method_ft 50472] acc interp N 2049 dt 6.18e-01 -> diff 1.06e-01 + stocproc.method_ft - INFO - acc interp N 4097 dt 3.09e-01 -> diff 2.75e-02 + stocproc.method_ft - INFO - acc interp N 4097 dt 3.09e-01 -> diff 2.75e-02 + stocproc.method_ft - INFO - acc interp N 4097 dt 3.09e-01 -> diff 2.75e-02 + stocproc.method_ft - INFO - acc interp N 4097 dt 3.09e-01 -> diff 2.75e-02 + stocproc.method_ft - INFO - acc interp N 4097 dt 3.09e-01 -> diff 2.75e-02 + stocproc.method_ft - INFO - acc interp N 4097 dt 3.09e-01 -> diff 2.75e-02 + stocproc.method_ft - INFO - acc interp N 4097 dt 3.09e-01 -> diff 2.75e-02 + [INFO stocproc.method_ft 50472] acc interp N 4097 dt 3.09e-01 -> diff 2.75e-02 + stocproc.method_ft - INFO - acc interp N 8193 dt 1.55e-01 -> diff 5.79e-03 + stocproc.method_ft - INFO - acc interp N 8193 dt 1.55e-01 -> diff 5.79e-03 + stocproc.method_ft - INFO - acc interp N 8193 dt 1.55e-01 -> diff 5.79e-03 + stocproc.method_ft - INFO - acc interp N 8193 dt 1.55e-01 -> diff 5.79e-03 + stocproc.method_ft - INFO - acc interp N 8193 dt 1.55e-01 -> diff 5.79e-03 + stocproc.method_ft - INFO - acc interp N 8193 dt 1.55e-01 -> diff 5.79e-03 + stocproc.method_ft - INFO - acc interp N 8193 dt 1.55e-01 -> diff 5.79e-03 + [INFO stocproc.method_ft 50472] acc interp N 8193 dt 1.55e-01 -> diff 5.79e-03 + stocproc.method_ft - INFO - acc interp N 16385 dt 7.73e-02 -> diff 1.35e-03 + stocproc.method_ft - INFO - acc interp N 16385 dt 7.73e-02 -> diff 1.35e-03 + stocproc.method_ft - INFO - acc interp N 16385 dt 7.73e-02 -> diff 1.35e-03 + stocproc.method_ft - INFO - acc interp N 16385 dt 7.73e-02 -> diff 1.35e-03 + stocproc.method_ft - INFO - acc interp N 16385 dt 7.73e-02 -> diff 1.35e-03 + stocproc.method_ft - INFO - acc interp N 16385 dt 7.73e-02 -> diff 1.35e-03 + stocproc.method_ft - INFO - acc interp N 16385 dt 7.73e-02 -> diff 1.35e-03 + [INFO stocproc.method_ft 50472] acc interp N 16385 dt 7.73e-02 -> diff 1.35e-03 + stocproc.method_ft - INFO - acc interp N 32769 dt 3.86e-02 -> diff 3.31e-04 + stocproc.method_ft - INFO - acc interp N 32769 dt 3.86e-02 -> diff 3.31e-04 + stocproc.method_ft - INFO - acc interp N 32769 dt 3.86e-02 -> diff 3.31e-04 + stocproc.method_ft - INFO - acc interp N 32769 dt 3.86e-02 -> diff 3.31e-04 + stocproc.method_ft - INFO - acc interp N 32769 dt 3.86e-02 -> diff 3.31e-04 + stocproc.method_ft - INFO - acc interp N 32769 dt 3.86e-02 -> diff 3.31e-04 + stocproc.method_ft - INFO - acc interp N 32769 dt 3.86e-02 -> diff 3.31e-04 + [INFO stocproc.method_ft 50472] acc interp N 32769 dt 3.86e-02 -> diff 3.31e-04 + stocproc.stocproc - INFO - requires dt < 3.864e-02 + stocproc.stocproc - INFO - requires dt < 3.864e-02 + stocproc.stocproc - INFO - requires dt < 3.864e-02 + stocproc.stocproc - INFO - requires dt < 3.864e-02 + stocproc.stocproc - INFO - requires dt < 3.864e-02 + stocproc.stocproc - INFO - requires dt < 3.864e-02 + stocproc.stocproc - INFO - requires dt < 3.864e-02 + [INFO stocproc.stocproc 50472] requires dt < 3.864e-02 + stocproc.stocproc - INFO - yields N = 20707 (time domain) + stocproc.stocproc - INFO - yields N = 20707 (time domain) + stocproc.stocproc - INFO - yields N = 20707 (time domain) + stocproc.stocproc - INFO - yields N = 20707 (time domain) + stocproc.stocproc - INFO - yields N = 20707 (time domain) + stocproc.stocproc - INFO - yields N = 20707 (time domain) + stocproc.stocproc - INFO - yields N = 20707 (time domain) + [INFO stocproc.stocproc 50472] yields N = 20707 (time domain) + stocproc.stocproc - INFO - find accurate discretisation in frequency domain + stocproc.stocproc - INFO - find accurate discretisation in frequency domain + stocproc.stocproc - INFO - find accurate discretisation in frequency domain + stocproc.stocproc - INFO - find accurate discretisation in frequency domain + stocproc.stocproc - INFO - find accurate discretisation in frequency domain + stocproc.stocproc - INFO - find accurate discretisation in frequency domain + stocproc.stocproc - INFO - find accurate discretisation in frequency domain + [INFO stocproc.stocproc 50472] find accurate discretisation in frequency domain + stocproc.stocproc - INFO - wmax:16.002043640915172 + stocproc.stocproc - INFO - wmax:16.002043640915172 + stocproc.stocproc - INFO - wmax:16.002043640915172 + stocproc.stocproc - INFO - wmax:16.002043640915172 + stocproc.stocproc - INFO - wmax:16.002043640915172 + stocproc.stocproc - INFO - wmax:16.002043640915172 + stocproc.stocproc - INFO - wmax:16.002043640915172 + [INFO stocproc.stocproc 50472] wmax:16.002043640915172 + n:32 d:0.805723129398869 tol:0.001 + n:64 d:0.6376491842373901 tol:0.001 + n:128 d:0.4119929516989666 tol:0.001 + n:256 d:0.26896884232518714 tol:0.001 + n:512 d:0.1711755293296233 tol:0.001 + n:1024 d:0.14524918216344024 tol:0.001 + n:2048 d:0.07837654392414552 tol:0.001 + n:4096 d:0.030531290938563486 tol:0.001 + n:8192 d:0.009187541105285848 tol:0.001 + stocproc.stocproc - INFO - perform numeric check of entire time axis [0,800.0027421234619] N:20707 + stocproc.stocproc - INFO - perform numeric check of entire time axis [0,800.0027421234619] N:20707 + stocproc.stocproc - INFO - perform numeric check of entire time axis [0,800.0027421234619] N:20707 + stocproc.stocproc - INFO - perform numeric check of entire time axis [0,800.0027421234619] N:20707 + stocproc.stocproc - INFO - perform numeric check of entire time axis [0,800.0027421234619] N:20707 + stocproc.stocproc - INFO - perform numeric check of entire time axis [0,800.0027421234619] N:20707 + stocproc.stocproc - INFO - perform numeric check of entire time axis [0,800.0027421234619] N:20707 + [INFO stocproc.stocproc 50472] perform numeric check of entire time axis [0,800.0027421234619] N:20707 + n:16384 d:2.2746881237172653e-05 tol:0.001 + done! +#+end_example + +Let's look at the size of the thermal process. +#+begin_src jupyter-python :results none + import sys + proc = model.thermal_process(1) +#+end_src + + +#+begin_src jupyter-python + sys.getsizeof(proc.omega_k) * 2 / 1024 ** 2 +#+end_src + +#+RESULTS: +: 0.250213623046875 + +Ok it gets sigificantly larger with longer time, but not too bad.. + +#+begin_src jupyter-python + from hops.core.integration import HOPSSupervisor + sup = HOPSSupervisor(config, 1, data_path=f".profiling_{np.random.random(6)}") +#+end_src + +#+RESULTS: +: [INFO hops.core.integration 50472] Choosing the nonlinear integrator. + +#+begin_src jupyter-python + sup.integrate_single_process() +#+end_src + +#+RESULTS: +:RESULTS: +: [INFO hops.core.integration 50472] Some 1 trajectories have to be integrated. +: [INFO hops.core.integration 50472] Using 9 hierarchy states. +: /home/hiro/src/stocproc/stocproc/stocproc.py:790: ComplexWarning: Casting complex values to real discards the imaginary part +: import pdb +: 0% 0/1 [01:21() + ----> 1 sup.integrate_single_process() + + File ~/src/hops/hops/core/integration.py:1199, in HOPSSupervisor.integrate_single_process(self, clear_pd) +  1186 integrator = self.actor( +  1187 t, +  1188 self.params, +  1189 data.result_filter, +  1190 ) +  1192 for seed, index in tqdm( +  1193 indices, +  1194 disable=self._hide_progress, +  (...) +  1197 mininterval=1, +  1198 ): + -> 1199 _, _, (incomplete, psi0, aux_states, stoc_proc) = integrator.integrate( +  1200  seed, index +  1201  ) +  1203 data.new_samples( +  1204 idx=index, +  1205 incomplete=incomplete, +  (...) +  1211 rng_seed=seed, +  1212 ) + + File ~/src/hops/hops/core/integration.py:718, in HOPSActor.integrate(self, seed, id) +  715 log.debug(f"Integrating trajectory {id} with seed {seed}.") +  716 self.update_random_numbers(seed) + --> 718 result = self.result_filter(self.call_solver()) +  719 return ( +  720 seed, +  721 id, +  722 result, +  723 ) + + File ~/src/hops/hops/core/integration.py:688, in HOPSActor.call_solver(self) +  683 def call_solver(self): +  684 """Call the solver and return the trajectory. +  685 +  686  Can be overwritten to customize the solver. +  687  """ + --> 688 solver_return = scipy.integrate.solve_ivp( +  689  self._rhs_wrapper, +  690  (0, self._t_points[-1]), +  691  self.initial_state(), +  692  t_eval=self._t_points, +  693  **self.params.IntP.solver_args, +  694  ) +  696 if solver_return.status < 0: +  697 raise RuntimeError(solver_return.message) + + File /nix/store/4p9d610cn1yqsfb9hizm1ychc9dz9kdx-python3-3.9.12-env/lib/python3.9/site-packages/scipy/integrate/_ivp/ivp.py:580, in solve_ivp(fun, t_span, y0, method, t_eval, dense_output, events, vectorized, args, **options) +  578 status = None +  579 while status is None: + --> 580 message = solver.step() +  582 if solver.status == 'finished': +  583 status = 0 + + File /nix/store/4p9d610cn1yqsfb9hizm1ychc9dz9kdx-python3-3.9.12-env/lib/python3.9/site-packages/scipy/integrate/_ivp/base.py:181, in OdeSolver.step(self) +  179 else: +  180 t = self.t + --> 181 success, message = self._step_impl() +  183 if not success: +  184 self.status = 'failed' + + File /nix/store/4p9d610cn1yqsfb9hizm1ychc9dz9kdx-python3-3.9.12-env/lib/python3.9/site-packages/scipy/integrate/_ivp/rk.py:144, in RungeKutta._step_impl(self) +  141 h = t_new - t +  142 h_abs = np.abs(h) + --> 144 y_new, f_new = rk_step(self.fun, t, y, self.f, h, self.A, +  145  self.B, self.C, self.K) +  146 scale = atol + np.maximum(np.abs(y), np.abs(y_new)) * rtol +  147 error_norm = self._estimate_error_norm(self.K, h, scale) + + File /nix/store/4p9d610cn1yqsfb9hizm1ychc9dz9kdx-python3-3.9.12-env/lib/python3.9/site-packages/scipy/integrate/_ivp/rk.py:67, in rk_step(fun, t, y, f, h, A, B, C, K) +  64 K[s] = fun(t + c * h, y + dy) +  66 y_new = y + h * np.dot(K[:-1].T, B) + ---> 67 f_new = fun(t + h, y_new) +  69 K[-1] = f_new +  71 return y_new, f_new + + File /nix/store/4p9d610cn1yqsfb9hizm1ychc9dz9kdx-python3-3.9.12-env/lib/python3.9/site-packages/scipy/integrate/_ivp/base.py:138, in OdeSolver.__init__..fun(t, y) +  136 def fun(t, y): +  137 self.nfev += 1 + --> 138 return self.fun_single(t, y) + + File /nix/store/4p9d610cn1yqsfb9hizm1ychc9dz9kdx-python3-3.9.12-env/lib/python3.9/site-packages/scipy/integrate/_ivp/base.py:20, in check_arguments..fun_wrapped(t, y) +  19 def fun_wrapped(t, y): + ---> 20 return np.asarray(fun(t, y), dtype=dtype) + + File ~/src/hops/hops/core/integration.py:681, in HOPSActor._rhs_wrapper(self, t, V_psi) +  678 if (t % 20) < 0.1: +  679 gc.collect() + --> 681 return self.rhs(t, V_psi) + + File ~/src/hops/hops/core/integration.py:1051, in NonLinearHOPSActor.rhs(self, t, V_psi_eta_lambda) +  1048 psi_zero = self.psi_zero(V_psi) +  1050 ddt_eta_lambda = self.ddt_eta_lambda(t, psi_zero, eta_lambda) + -> 1051 ddt_psi = self.generic_rhs( +  1052  t, +  1053  V_psi, +  1054  eta_det, +  1055 ) +  1057 # include normalization by subtracting the amount of change in the +  1058 # direction of the psi_sys-vector: ddt_psi -> ddt_psi - dn * V_psi +  1059 dn = 0 + + File ~/src/hops/hops/core/integration.py:647, in HOPSActor.generic_rhs(self, t, V_psi, eta_det) +  640 # filling the data here ... +  641 K_mat = self.K(t, psi, eta_stoc, eta_det) +  643 ddt_psi_view[:, :] = ( +  644 self._o_vec * psi +  645 + psi.dot(K_mat.T) +  646 + sum( + --> 647 M_up.dot(psi).dot(B.T) for M_up, B in zip(self._m_ups, self.Bs(t, psi)) +  648 ) +  649 + sum( +  650 M_down.dot(psi).dot(C.T) +  651 for M_down, C in zip(self._m_downs, self.Cs(t, psi)) +  652 ) +  653 ) +  654 # ... leads to the correct linear alignment for the result vector +  656 return ddt_psi + + File ~/src/hops/hops/core/integration.py:951, in NonLinearHOPSActor.Bs(self, t, psi) +  949 psi_zero = psi[0, :] +  950 exp_vals = np.zeros(self.minus_L_dagger.shape, dtype=np.complex128) + --> 951 for i, exp_val in enumerate(self.coupling_expectation(t, psi_zero)): +  952 np.fill_diagonal(exp_vals[i], exp_val) +  954 return self.minus_L_dagger(t) + exp_vals + + File ~/src/hops/hops/core/integration.py:752, in HOPSActor.coupling_expectation(self, t, psi_zero) +  738 r""" +  739 Calculate the expectation values of :math:`L^{†,(n)}` in the +  740 current HOPS state. +  (...) +  748 :param psi_zero: The zeroth order HOPS state. +  749 """ +  751 psi_conj = psi_zero.conj() + --> 752 return np.dot(self.L_dagger(t).dot(psi_zero), psi_conj) / np.dot( +  753 psi_conj, psi_zero +  754 ) + + KeyboardInterrupt: +#+end_example +:END: + +#+begin_src jupyter-python :results none + import cProfile + + profiling_file = f'.profiling_data/profile_stats_{model.hexhash}' + cProfile.run('sup.integrate_single_process(True)', profiling_file) +#+end_src diff --git a/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/pyproject.toml b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/pyproject.toml index b820d1c..d692fd7 100644 --- a/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/pyproject.toml +++ b/python/energy_flow_proper/09_dynamic_two_bath_one_qubit/pyproject.toml @@ -13,6 +13,7 @@ matplotlib = "^3.5.0" jupyter = "^1.0.0" qutip = "^4.6.2" lmfit = "^1.0.3" +memory-profiler = "^0.60.0" [tool.poetry.dev-dependencies] black = "^21.12b0"