qupled 1.1.0__tar.gz → 1.2.0__tar.gz
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.
- {qupled-1.1.0 → qupled-1.2.0}/PKG-INFO +4 -1
- {qupled-1.1.0 → qupled-1.2.0}/README.md +3 -0
- {qupled-1.1.0 → qupled-1.2.0}/dev/requirements-apt.txt +0 -1
- {qupled-1.1.0 → qupled-1.2.0}/dev/requirements-brew.txt +0 -1
- {qupled-1.1.0 → qupled-1.2.0}/docs/introduction.rst +2 -2
- {qupled-1.1.0 → qupled-1.2.0}/docs/qupled.rst +12 -24
- {qupled-1.1.0 → qupled-1.2.0}/examples/docs/solve_quantum_schemes.py +1 -0
- {qupled-1.1.0 → qupled-1.2.0}/examples/readme/create_cover.py +6 -8
- {qupled-1.1.0 → qupled-1.2.0}/examples/readme/qupled_animation_dark.svg +932 -915
- {qupled-1.1.0 → qupled-1.2.0}/examples/readme/qupled_animation_light.svg +1085 -1063
- qupled-1.2.0/examples/tests/test_examples.py +122 -0
- {qupled-1.1.0 → qupled-1.2.0}/pyproject.toml +1 -1
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/hf.py +2 -2
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/esa.hpp +6 -8
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/hf.hpp +18 -13
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/iet.hpp +19 -13
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/input.hpp +18 -53
- qupled-1.2.0/src/qupled/native/include/python_interface/inputs.hpp +10 -0
- qupled-1.2.0/src/qupled/native/include/python_interface/schemes.hpp +11 -0
- qupled-1.1.0/src/qupled/native/include/python_util.hpp → qupled-1.2.0/src/qupled/native/include/python_interface/util.hpp +2 -2
- qupled-1.2.0/src/qupled/native/include/python_interface/utilities.hpp +11 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/qstls.hpp +12 -51
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/qstlsiet.hpp +21 -53
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/qvsstls.hpp +36 -38
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/rpa.hpp +26 -34
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/stls.hpp +35 -11
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/stlsiet.hpp +21 -18
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/vector2D.hpp +6 -2
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/vector3D.hpp +2 -3
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/vsbase.hpp +37 -26
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/vsstls.hpp +34 -33
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/CMakeLists.txt +7 -7
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/esa.cpp +3 -21
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/hf.cpp +33 -29
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/iet.cpp +46 -15
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/input.cpp +3 -23
- qupled-1.2.0/src/qupled/native/src/python_interface/inputs.cpp +240 -0
- qupled-1.2.0/src/qupled/native/src/python_interface/native.cpp +39 -0
- qupled-1.2.0/src/qupled/native/src/python_interface/schemes.cpp +151 -0
- qupled-1.1.0/src/qupled/native/src/python_util.cpp → qupled-1.2.0/src/qupled/native/src/python_interface/util.cpp +1 -1
- qupled-1.2.0/src/qupled/native/src/python_interface/utilities.cpp +76 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/qstls.cpp +47 -172
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/qstlsiet.cpp +59 -167
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/qvsstls.cpp +45 -113
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/rpa.cpp +22 -55
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/stls.cpp +41 -49
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/stlsiet.cpp +27 -56
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/vector2D.cpp +6 -2
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/vector3D.cpp +0 -8
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/vsbase.cpp +68 -27
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/vsstls.cpp +24 -71
- qupled-1.2.0/src/qupled/qstls.py +68 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/qstlsiet.py +8 -13
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/qvsstls.py +1 -15
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/stls.py +9 -9
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/stlsiet.py +39 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled.egg-info/PKG-INFO +4 -1
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled.egg-info/SOURCES.txt +9 -5
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/test_esa_native.py +2 -3
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/test_hf_native.py +2 -2
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/test_qstls_iet_native.py +11 -2
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/test_qstls_native.py +11 -5
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/test_qvsstls_native.py +11 -5
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/test_rpa_native.py +3 -3
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/test_stls_iet_native.py +11 -2
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/test_stls_native.py +3 -3
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/test_vsstls_native.py +11 -2
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_hf.py +1 -1
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_qstls.py +1 -54
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_qstlsiet.py +4 -13
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_qvsstls.py +0 -9
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_stls.py +16 -16
- qupled-1.2.0/tests/test_stlsiet.py +78 -0
- qupled-1.1.0/examples/tests/test_examples.py +0 -73
- qupled-1.1.0/src/qupled/native/include/python_wrappers.hpp +0 -230
- qupled-1.1.0/src/qupled/native/src/python_modules.cpp +0 -222
- qupled-1.1.0/src/qupled/native/src/python_wrappers.cpp +0 -299
- qupled-1.1.0/src/qupled/qstls.py +0 -151
- qupled-1.1.0/tests/test_stlsiet.py +0 -45
- {qupled-1.1.0 → qupled-1.2.0}/.clang-format +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/.devcontainer/Dockerfile +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/.devcontainer/devcontainer.json +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/.github/workflows/build-and-test-base.yml +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/.github/workflows/build-and-test.yml +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/.github/workflows/formatting.yml +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/.github/workflows/release.yml +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/.gitignore +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/.readthedocs.yaml +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/LICENSE +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/MANIFEST.in +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/dev/devtool.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/dev/requirements-pip.txt +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/devtool +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/docs/_static/css/rdt_theme_python_properties.css +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/docs/conf.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/docs/contribute.rst +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/docs/examples.rst +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/docs/index.rst +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/docs/make.bat +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/docs/requirements.txt +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/examples/docs/fixed_adr.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/examples/docs/initial_guess_stls.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/examples/docs/solve_qvsstls.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/examples/docs/solve_rpa_and_esa.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/examples/docs/solve_stls.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/examples/docs/solve_stls_iet.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/examples/docs/solve_vsstls.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/setup.cfg +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/setup.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/__init__.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/database.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/esa.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/mpi.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/chemical_potential.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/database.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/dual.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/free_energy.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/internal_energy.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/logger.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/mpi_util.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/num_util.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/numerics.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/rdf.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/thermo_util.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/include/vector_util.hpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/chemical_potential.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/free_energy.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/internal_energy.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/logger.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/mpi_util.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/num_util.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/numerics.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/rdf.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/thermo_util.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/native/src/vector_util.cpp +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/output.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/rpa.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled/vsstls.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled.egg-info/dependency_links.txt +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled.egg-info/not-zip-safe +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled.egg-info/requires.txt +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/src/qupled.egg-info/top_level.txt +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/tests/native/conftest.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_database.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_esa.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_mpi.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_output.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_rpa.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/tests/test_vsstls.py +0 -0
- {qupled-1.1.0 → qupled-1.2.0}/tox.ini +0 -0
@@ -1,6 +1,6 @@
|
|
1
1
|
Metadata-Version: 2.4
|
2
2
|
Name: qupled
|
3
|
-
Version: 1.
|
3
|
+
Version: 1.2.0
|
4
4
|
Summary: qupled: a package to investigate quantum plasmas via the dielectric formalism
|
5
5
|
Author-email: Federico Lucco Castello <federico.luccocastello@gmail.com>
|
6
6
|
License-Expression: GPL-3.0-or-later
|
@@ -60,6 +60,9 @@ More detailed information on the package together with a list of examples is ava
|
|
60
60
|
Qupled has been used in the following publications:
|
61
61
|
|
62
62
|
<ol>
|
63
|
+
<li>
|
64
|
+
<a href="https://onlinelibrary.wiley.com/doi/10.1002/ctpp.70014">Tolias, P., Kalkavouras, F., Dornheim, T. & Lucco Castello, F. (2025). Dynamic Properties of the Warm Dense Uniform Electron Gas With the qSTLS Dielectric Scheme. <i>Contributions to Plasma Physics</i>, 0:e70014</a>
|
65
|
+
</li>
|
63
66
|
<li>
|
64
67
|
<a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.109.125134">Tolias, P., Lucco Castello, F., Kalkavouras, F., & Dornheim, T. (2024). Revisiting the Vashishta-Singwi dielectric scheme for the warm dense uniform electron fluid. <i>Physical Review B</i>, <i>109</i>(12)</a>
|
65
68
|
</li>
|
@@ -30,6 +30,9 @@ More detailed information on the package together with a list of examples is ava
|
|
30
30
|
Qupled has been used in the following publications:
|
31
31
|
|
32
32
|
<ol>
|
33
|
+
<li>
|
34
|
+
<a href="https://onlinelibrary.wiley.com/doi/10.1002/ctpp.70014">Tolias, P., Kalkavouras, F., Dornheim, T. & Lucco Castello, F. (2025). Dynamic Properties of the Warm Dense Uniform Electron Gas With the qSTLS Dielectric Scheme. <i>Contributions to Plasma Physics</i>, 0:e70014</a>
|
35
|
+
</li>
|
33
36
|
<li>
|
34
37
|
<a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.109.125134">Tolias, P., Lucco Castello, F., Kalkavouras, F., & Dornheim, T. (2024). Revisiting the Vashishta-Singwi dielectric scheme for the warm dense uniform electron fluid. <i>Physical Review B</i>, <i>109</i>(12)</a>
|
35
38
|
</li>
|
@@ -56,13 +56,13 @@ For linux distributions all these dependencies can be installed with
|
|
56
56
|
|
57
57
|
.. code-block:: console
|
58
58
|
|
59
|
-
sudo apt-get install -y cmake libboost-all-dev libopenmpi-dev libgsl-dev libomp-dev
|
59
|
+
sudo apt-get install -y cmake libboost-all-dev libopenmpi-dev libgsl-dev libomp-dev python3-dev libsqlite3-dev libsqlitecpp-dev
|
60
60
|
|
61
61
|
For macOS they can be installed directly from homebrew
|
62
62
|
|
63
63
|
.. code-block:: console
|
64
64
|
|
65
|
-
brew install cmake gsl libomp openmpi
|
65
|
+
brew install cmake gsl libomp openmpi boost-python3 sqlite sqlitecpp
|
66
66
|
|
67
67
|
Install with pip
|
68
68
|
~~~~~~~~~~~~~~~~
|
@@ -56,7 +56,7 @@ After the solution is completed the results are stored in an object :obj:`qupled
|
|
56
56
|
and written to the output database.
|
57
57
|
|
58
58
|
.. autoclass:: qupled.rpa.Rpa
|
59
|
-
:
|
59
|
+
:show-inheritance:
|
60
60
|
:members:
|
61
61
|
|
62
62
|
.. autoclass:: qupled.rpa.Input
|
@@ -73,7 +73,7 @@ After the solution is completed the results are stored in an object :obj:`qupled
|
|
73
73
|
and written to the output database.
|
74
74
|
|
75
75
|
.. autoclass:: qupled.stls.Stls
|
76
|
-
:
|
76
|
+
:show-inheritance:
|
77
77
|
:members:
|
78
78
|
|
79
79
|
.. autoclass:: qupled.stls.Input
|
@@ -98,7 +98,7 @@ After the solution is completed the results are stored in an object :obj:`qupled
|
|
98
98
|
and written to the output database.
|
99
99
|
|
100
100
|
.. autoclass:: qupled.stlsiet.StlsIet
|
101
|
-
:
|
101
|
+
:show-inheritance:
|
102
102
|
:members:
|
103
103
|
|
104
104
|
.. autoclass:: qupled.stlsiet.Input
|
@@ -109,6 +109,10 @@ and written to the output database.
|
|
109
109
|
:show-inheritance:
|
110
110
|
:members:
|
111
111
|
|
112
|
+
.. autoclass:: qupled.stlsiet.Guess
|
113
|
+
:members:
|
114
|
+
:exclude-members: to_native
|
115
|
+
|
112
116
|
VSStls scheme
|
113
117
|
~~~~~~~~~~~~~
|
114
118
|
|
@@ -119,7 +123,7 @@ After the solution is completed the results are stored in an object :obj:`qupled
|
|
119
123
|
and written to the output database.
|
120
124
|
|
121
125
|
.. autoclass:: qupled.vsstls.VSStls
|
122
|
-
:
|
126
|
+
:show-inheritance:
|
123
127
|
:members:
|
124
128
|
:exclude-members: get_free_energy_integrand
|
125
129
|
|
@@ -148,7 +152,7 @@ After the solution is completed the results are stored in an object :obj:`qupled
|
|
148
152
|
and written to the output database.
|
149
153
|
|
150
154
|
.. autoclass:: qupled.esa.ESA
|
151
|
-
:
|
155
|
+
:show-inheritance:
|
152
156
|
:members:
|
153
157
|
|
154
158
|
.. autoclass:: qupled.esa.Input
|
@@ -168,21 +172,13 @@ After the solution is completed the results are stored in an object :obj:`qupled
|
|
168
172
|
and written to the output database.
|
169
173
|
|
170
174
|
.. autoclass:: qupled.qstls.Qstls
|
171
|
-
:
|
175
|
+
:show-inheritance:
|
172
176
|
:members:
|
173
177
|
:exclude-members: find_fixed_adr_in_database
|
174
178
|
|
175
179
|
.. autoclass:: qupled.qstls.Input
|
176
180
|
:show-inheritance:
|
177
181
|
:members:
|
178
|
-
|
179
|
-
.. autoclass:: qupled.qstls.Result
|
180
|
-
:show-inheritance:
|
181
|
-
:members:
|
182
|
-
|
183
|
-
.. autoclass:: qupled.qstls.Guess
|
184
|
-
:members:
|
185
|
-
:exclude-members: to_native
|
186
182
|
|
187
183
|
Qstls-IET scheme
|
188
184
|
~~~~~~~~~~~~~~~~
|
@@ -195,17 +191,13 @@ After the solution is completed the results are stored in an object :obj:`qupled
|
|
195
191
|
and written to the output database.
|
196
192
|
|
197
193
|
.. autoclass:: qupled.qstlsiet.QstlsIet
|
198
|
-
:
|
194
|
+
:show-inheritance:
|
199
195
|
:members:
|
200
196
|
:exclude-members: find_fixed_adr_in_database
|
201
197
|
|
202
198
|
.. autoclass:: qupled.qstlsiet.Input
|
203
199
|
:show-inheritance:
|
204
200
|
:members:
|
205
|
-
|
206
|
-
.. autoclass:: qupled.qstlsiet.Result
|
207
|
-
:show-inheritance:
|
208
|
-
:members:
|
209
201
|
|
210
202
|
QVSStls scheme
|
211
203
|
~~~~~~~~~~~~~~~~
|
@@ -216,17 +208,13 @@ called :obj:`qupled.qvsstls.Input`. After the solution is completed the results
|
|
216
208
|
object :obj:`qupled.qvsstls.Result` and written to the output database.
|
217
209
|
|
218
210
|
.. autoclass:: qupled.qvsstls.QVSStls
|
219
|
-
:
|
211
|
+
:show-inheritance:
|
220
212
|
:members:
|
221
213
|
:exclude-members: get_free_energy_integrand
|
222
214
|
|
223
215
|
.. autoclass:: qupled.qvsstls.Input
|
224
216
|
:show-inheritance:
|
225
217
|
:members:
|
226
|
-
|
227
|
-
.. autoclass:: qupled.qvsstls.Result
|
228
|
-
:show-inheritance:
|
229
|
-
:members:
|
230
218
|
|
231
219
|
Output database
|
232
220
|
---------------
|
@@ -81,30 +81,28 @@ def create_one_svg_file(darkmode: bool, scheme: qstls.Qstls, error: np.array):
|
|
81
81
|
|
82
82
|
def solve_qstls(guess_run_id: int):
|
83
83
|
scheme = qstls.Qstls()
|
84
|
-
|
85
|
-
theta = 1.0
|
86
|
-
inputs = qstls.Input(rs, theta)
|
84
|
+
inputs = qstls.Input(coupling=15.0, degeneracy=1.0)
|
87
85
|
inputs.mixing = 0.3
|
88
86
|
inputs.resolution = 0.1
|
89
87
|
inputs.cutoff = 10
|
90
88
|
inputs.matsubara = 16
|
91
89
|
inputs.threads = 16
|
92
90
|
inputs.iterations = 0
|
93
|
-
adr_file = f"adr_fixed_theta{theta:.3f}_matsubara{inputs.matsubara}_QSTLS.bin"
|
94
91
|
inputs.guess = (
|
95
92
|
scheme.get_initial_guess(guess_run_id)
|
96
93
|
if guess_run_id is not None
|
97
94
|
else inputs.guess
|
98
95
|
)
|
99
|
-
inputs.fixed = adr_file if os.path.exists(adr_file) else inputs.fixed
|
100
96
|
scheme.compute(inputs)
|
101
97
|
return scheme
|
102
98
|
|
103
99
|
|
104
100
|
def plot_density_response(plt: plt, scheme: qstls.Qstls, settings: PlotSettings):
|
105
101
|
results = scheme.results
|
106
|
-
|
107
|
-
|
102
|
+
inputs = scheme.inputs
|
103
|
+
wvg_squared = results.wvg[:, np.newaxis] ** 2
|
104
|
+
denominator = wvg_squared + inputs.coupling * (results.idr - results.adr)
|
105
|
+
dr = results.idr * wvg_squared / denominator
|
108
106
|
plt.subplot(2, 2, 3)
|
109
107
|
parameters = np.array([0, 1, 2, 3, 4])
|
110
108
|
numParameters = parameters.size
|
@@ -124,7 +122,7 @@ def plot_density_response(plt: plt, scheme: qstls.Qstls, settings: PlotSettings)
|
|
124
122
|
plt.xlim(0, settings.xlim)
|
125
123
|
plt.xlabel("Wave-vector", fontsize=settings.labelsz)
|
126
124
|
plt.title("Density response", fontsize=settings.labelsz, fontweight="bold")
|
127
|
-
plt.legend(fontsize=settings.ticksz, loc="
|
125
|
+
plt.legend(fontsize=settings.ticksz, loc="upper right")
|
128
126
|
plt.xticks(fontsize=settings.ticksz)
|
129
127
|
plt.yticks(fontsize=settings.ticksz)
|
130
128
|
|