update graphics in cycle length coupling strength

This commit is contained in:
valentin.boettcher@mailbox.tu-dresden.de 2023-07-10 13:35:59 -04:00
parent 9a29fc0bef
commit b14d14cb01
No known key found for this signature in database
GPG key ID: E034E12B7AF56ACE
66 changed files with 956895 additions and 119 deletions

View file

@ -14,9 +14,12 @@ def plot_power_eff_convergence(models, steady_idx=2):
a_efficiency.set_yscale("log")
for model in models:
Ns = model.power(steady_idx=steady_idx).Ns
a_power.plot(Ns, model.power(steady_idx=steady_idx).values)
a_efficiency.plot(Ns, np.abs(model.efficiency(steady_idx=steady_idx).values))
try:
Ns = model.power(steady_idx=steady_idx).Ns
a_power.plot(Ns, model.power(steady_idx=steady_idx).values)
a_efficiency.plot(Ns, np.abs(model.efficiency(steady_idx=steady_idx).values))
except:
pass
a_power.set_xlabel("$N$")
a_power.set_ylabel("$P$")

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 2.8 MiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 2.1 MiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 3.8 MiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 43 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 40 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 3.7 MiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 265 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 37 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 2.1 MiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 39 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 88 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 39 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 39 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 43 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 37 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 40 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 38 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 42 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 44 KiB

View file

@ -39,7 +39,7 @@ z
</g>
<g id="PolyCollection_1">
<defs>
<path id="m346fc0bb85" d="M 83.77363 -201.368344
<path id="m522b71bcad" d="M 83.77363 -201.368344
L 83.77363 -201.368344
L 84.021784 -201.368344
L 84.269938 -201.368344
@ -6050,13 +6050,13 @@ L 83.77363 -201.368344
z
" style="stroke: #bfdef4; stroke-opacity: 0.5; stroke-width: 0.24"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#m346fc0bb85" x="0" y="296.39952" style="fill: #bfdef4; fill-opacity: 0.5; stroke: #bfdef4; stroke-opacity: 0.5; stroke-width: 0.24"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#m522b71bcad" x="0" y="296.39952" style="fill: #bfdef4; fill-opacity: 0.5; stroke: #bfdef4; stroke-opacity: 0.5; stroke-width: 0.24"/>
</g>
</g>
<g id="PolyCollection_2">
<defs>
<path id="m4bff993999" d="M 83.77363 -201.368344
<path id="m579150b4a2" d="M 83.77363 -201.368344
L 83.77363 -201.368344
L 84.021784 -201.368344
L 84.269938 -201.368344
@ -12067,13 +12067,13 @@ L 83.77363 -201.368344
z
" style="stroke: #ffdfc3; stroke-opacity: 0.5; stroke-width: 0.24"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#m4bff993999" x="0" y="296.39952" style="fill: #ffdfc3; fill-opacity: 0.5; stroke: #ffdfc3; stroke-opacity: 0.5; stroke-width: 0.24"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#m579150b4a2" x="0" y="296.39952" style="fill: #ffdfc3; fill-opacity: 0.5; stroke: #ffdfc3; stroke-opacity: 0.5; stroke-width: 0.24"/>
</g>
</g>
<g id="PolyCollection_3">
<defs>
<path id="me032711082" d="M 83.77363 -201.368344
<path id="mb9d217397f" d="M 83.77363 -201.368344
L 83.77363 -201.368344
L 84.021784 -201.368344
L 84.269938 -201.368344
@ -18084,13 +18084,13 @@ L 83.77363 -201.368344
z
" style="stroke: #c3eec3; stroke-opacity: 0.5; stroke-width: 0.24"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#me032711082" x="0" y="296.39952" style="fill: #c3eec3; fill-opacity: 0.5; stroke: #c3eec3; stroke-opacity: 0.5; stroke-width: 0.24"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#mb9d217397f" x="0" y="296.39952" style="fill: #c3eec3; fill-opacity: 0.5; stroke: #c3eec3; stroke-opacity: 0.5; stroke-width: 0.24"/>
</g>
</g>
<g id="PolyCollection_4">
<defs>
<path id="mcfc6e5a0a3" d="M 83.77363 -201.368344
<path id="mf0f7ee4703" d="M 83.77363 -201.368344
L 83.77363 -201.368344
L 84.021784 -201.368344
L 84.269938 -201.368344
@ -24101,13 +24101,13 @@ L 83.77363 -201.368344
z
" style="stroke: #f5c9c9; stroke-opacity: 0.5; stroke-width: 0.24"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#mcfc6e5a0a3" x="0" y="296.39952" style="fill: #f5c9c9; fill-opacity: 0.5; stroke: #f5c9c9; stroke-opacity: 0.5; stroke-width: 0.24"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#mf0f7ee4703" x="0" y="296.39952" style="fill: #f5c9c9; fill-opacity: 0.5; stroke: #f5c9c9; stroke-opacity: 0.5; stroke-width: 0.24"/>
</g>
</g>
<g id="PolyCollection_5">
<defs>
<path id="m23245554f8" d="M 83.77363 -201.368344
<path id="m2cd9cdc7f5" d="M 83.77363 -201.368344
L 83.77363 -201.368344
L 84.021784 -201.368344
L 84.269938 -201.368344
@ -30118,13 +30118,13 @@ L 83.77363 -201.368344
z
" style="stroke: #e4d9ee; stroke-opacity: 0.5; stroke-width: 0.24"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#m23245554f8" x="0" y="296.39952" style="fill: #e4d9ee; fill-opacity: 0.5; stroke: #e4d9ee; stroke-opacity: 0.5; stroke-width: 0.24"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#m2cd9cdc7f5" x="0" y="296.39952" style="fill: #e4d9ee; fill-opacity: 0.5; stroke: #e4d9ee; stroke-opacity: 0.5; stroke-width: 0.24"/>
</g>
</g>
<g id="PolyCollection_6">
<defs>
<path id="me1442dc4ab" d="M 83.77363 -201.368344
<path id="m974683f3ad" d="M 83.77363 -201.368344
L 83.77363 -201.368344
L 84.021784 -201.368344
L 84.269938 -201.368344
@ -36135,20 +36135,20 @@ L 83.77363 -201.368344
z
" style="stroke: #e5d3cf; stroke-opacity: 0.5; stroke-width: 0.24"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#me1442dc4ab" x="0" y="296.39952" style="fill: #e5d3cf; fill-opacity: 0.5; stroke: #e5d3cf; stroke-opacity: 0.5; stroke-width: 0.24"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#m974683f3ad" x="0" y="296.39952" style="fill: #e5d3cf; fill-opacity: 0.5; stroke: #e5d3cf; stroke-opacity: 0.5; stroke-width: 0.24"/>
</g>
</g>
<g id="matplotlib.axis_1">
<g id="xtick_1">
<g id="line2d_1">
<defs>
<path id="me3ad2d4970" d="M 0 0
<path id="m38fac81e65" d="M 0 0
L 0 3.5
" style="stroke: #000000; stroke-width: 0.8"/>
</defs>
<g>
<use xlink:href="#me3ad2d4970" x="83.77363" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m38fac81e65" x="83.77363" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_1">
@ -36184,7 +36184,7 @@ z
<g id="xtick_2">
<g id="line2d_2">
<g>
<use xlink:href="#me3ad2d4970" x="183.002314" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m38fac81e65" x="183.002314" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_2">
@ -36222,7 +36222,7 @@ z
<g id="xtick_3">
<g id="line2d_3">
<g>
<use xlink:href="#me3ad2d4970" x="282.230999" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m38fac81e65" x="282.230999" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_3">
@ -36264,7 +36264,7 @@ z
<g id="xtick_4">
<g id="line2d_4">
<g>
<use xlink:href="#me3ad2d4970" x="381.459683" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m38fac81e65" x="381.459683" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_4">
@ -36303,7 +36303,7 @@ z
<g id="xtick_5">
<g id="line2d_5">
<g>
<use xlink:href="#me3ad2d4970" x="480.688368" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m38fac81e65" x="480.688368" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_5">
@ -36347,7 +36347,7 @@ z
<g id="xtick_6">
<g id="line2d_6">
<g>
<use xlink:href="#me3ad2d4970" x="579.917052" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m38fac81e65" x="579.917052" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_6">
@ -36382,7 +36382,7 @@ z
<g id="xtick_7">
<g id="line2d_7">
<g>
<use xlink:href="#me3ad2d4970" x="679.145737" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m38fac81e65" x="679.145737" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_7">
@ -36397,7 +36397,7 @@ z
<g id="xtick_8">
<g id="line2d_8">
<g>
<use xlink:href="#me3ad2d4970" x="778.374421" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m38fac81e65" x="778.374421" y="260.159296" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_8">
@ -36440,12 +36440,12 @@ z
<g id="ytick_1">
<g id="line2d_9">
<defs>
<path id="mf1fad36efc" d="M 0 0
<path id="m3baec71acd" d="M 0 0
L -3.5 0
" style="stroke: #000000; stroke-width: 0.8"/>
</defs>
<g>
<use xlink:href="#mf1fad36efc" x="46.562873" y="236.317349" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m3baec71acd" x="46.562873" y="236.317349" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_10">
@ -36507,7 +36507,7 @@ z
<g id="ytick_2">
<g id="line2d_10">
<g>
<use xlink:href="#mf1fad36efc" x="46.562873" y="189.221958" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m3baec71acd" x="46.562873" y="189.221958" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_11">
@ -36523,7 +36523,7 @@ z
<g id="ytick_3">
<g id="line2d_11">
<g>
<use xlink:href="#mf1fad36efc" x="46.562873" y="142.126567" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m3baec71acd" x="46.562873" y="142.126567" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_12">
@ -36539,7 +36539,7 @@ z
<g id="ytick_4">
<g id="line2d_12">
<g>
<use xlink:href="#mf1fad36efc" x="46.562873" y="95.031176" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m3baec71acd" x="46.562873" y="95.031176" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_13">
@ -36554,7 +36554,7 @@ z
<g id="ytick_5">
<g id="line2d_13">
<g>
<use xlink:href="#mf1fad36efc" x="46.562873" y="47.935785" style="stroke: #000000; stroke-width: 0.8"/>
<use xlink:href="#m3baec71acd" x="46.562873" y="47.935785" style="stroke: #000000; stroke-width: 0.8"/>
</g>
</g>
<g id="text_14">
@ -36748,13 +36748,13 @@ z
<g id="LineCollection_1">
<path d="M 83.77363 95.031176
L 83.77363 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #1f77b4; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #1f77b4; stroke-width: 0.5"/>
<path d="M 331.845341 81.928022
L 331.845341 81.835339
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #1f77b4; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #1f77b4; stroke-width: 0.5"/>
<path d="M 579.917052 81.740477
L 579.917052 81.645779
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #1f77b4; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #1f77b4; stroke-width: 0.5"/>
</g>
<g id="line2d_14">
<path d="M 83.77363 95.031176
@ -36996,18 +36996,18 @@ L 824.018292 81.633874
L 826.747991 81.713147
L 827.988763 81.714066
L 827.988763 81.714066
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #7fbee9; stroke-width: 0.5; stroke-linecap: square"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #7fbee9; stroke-width: 0.5; stroke-linecap: square"/>
</g>
<g id="LineCollection_2">
<path d="M 83.77363 95.031176
L 83.77363 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #ff7f0e; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #ff7f0e; stroke-width: 0.5"/>
<path d="M 331.845341 70.386012
L 331.845341 70.360613
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #ff7f0e; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #ff7f0e; stroke-width: 0.5"/>
<path d="M 579.917052 44.846719
L 579.917052 44.810727
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #ff7f0e; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #ff7f0e; stroke-width: 0.5"/>
</g>
<g id="line2d_15">
<path d="M 83.77363 95.031176
@ -37082,18 +37082,18 @@ L 819.055204 18.880396
L 824.762756 19.534531
L 827.988763 19.554958
L 827.988763 19.554958
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #ffbf86; stroke-width: 0.5; stroke-linecap: square"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #ffbf86; stroke-width: 0.5; stroke-linecap: square"/>
</g>
<g id="LineCollection_3">
<path d="M 83.77363 95.031176
L 83.77363 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #2ca02c; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #2ca02c; stroke-width: 0.5"/>
<path d="M 331.845341 95.031176
L 331.845341 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #2ca02c; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #2ca02c; stroke-width: 0.5"/>
<path d="M 579.917052 95.031176
L 579.917052 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #2ca02c; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #2ca02c; stroke-width: 0.5"/>
</g>
<g id="line2d_16">
<path d="M 83.77363 95.031176
@ -37171,18 +37171,18 @@ L 825.259065 95.143759
L 826.747991 95.036918
L 827.988763 95.031176
L 827.988763 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #87de87; stroke-width: 0.5; stroke-linecap: square"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #87de87; stroke-width: 0.5; stroke-linecap: square"/>
</g>
<g id="LineCollection_4">
<path d="M 83.77363 95.031176
L 83.77363 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #d62728; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #d62728; stroke-width: 0.5"/>
<path d="M 331.845341 164.05691
L 331.845341 163.904228
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #d62728; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #d62728; stroke-width: 0.5"/>
<path d="M 579.917052 206.666396
L 579.917052 206.447888
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #d62728; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #d62728; stroke-width: 0.5"/>
</g>
<g id="line2d_17">
<path d="M 83.77363 95.031176
@ -37266,18 +37266,18 @@ L 698.700306 248.494473
L 702.670776 248.524605
L 827.988763 248.524388
L 827.988763 248.524388
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #eb9293; stroke-width: 0.5; stroke-linecap: square"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #eb9293; stroke-width: 0.5; stroke-linecap: square"/>
</g>
<g id="LineCollection_5">
<path d="M 83.77363 95.031176
L 83.77363 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #9467bd; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #9467bd; stroke-width: 0.5"/>
<path d="M 331.845341 95.031176
L 331.845341 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #9467bd; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #9467bd; stroke-width: 0.5"/>
<path d="M 579.917052 95.031176
L 579.917052 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #9467bd; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #9467bd; stroke-width: 0.5"/>
</g>
<g id="line2d_18">
<path d="M 83.77363 95.031176
@ -37417,18 +37417,18 @@ L 702.174468 95.032643
L 715.078498 95.031176
L 827.988763 95.031176
L 827.988763 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #c9b3de; stroke-width: 0.5; stroke-linecap: square"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #c9b3de; stroke-width: 0.5; stroke-linecap: square"/>
</g>
<g id="LineCollection_6">
<path d="M 83.77363 95.031176
L 83.77363 95.031176
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #8c564b; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #8c564b; stroke-width: 0.5"/>
<path d="M 331.845341 126.048672
L 331.845341 126.017988
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #8c564b; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #8c564b; stroke-width: 0.5"/>
<path d="M 579.917052 143.011874
L 579.917052 142.964945
" clip-path="url(#pd926ee2045)" style="fill: none; stroke: #8c564b; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke: #8c564b; stroke-width: 0.5"/>
</g>
<g id="line2d_19">
<path d="M 83.77363 95.031176
@ -37598,11 +37598,11 @@ L 825.507219 159.601042
L 827.740609 159.543131
L 827.988763 159.543114
L 827.988763 159.543114
" clip-path="url(#pd926ee2045)" style="fill: none; stroke-dasharray: 1.85,0.8; stroke-dashoffset: 0; stroke: #cca79f; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke-dasharray: 1.85,0.8; stroke-dashoffset: 0; stroke: #cca79f; stroke-width: 0.5"/>
</g>
<g id="line2d_20">
<defs>
<path id="mce91a6f3f8" d="M 0 1
<path id="meb3df71cc8" d="M 0 1
C 0.265203 1 0.51958 0.894634 0.707107 0.707107
C 0.894634 0.51958 1 0.265203 1 0
C 1 -0.265203 0.894634 -0.51958 0.707107 -0.707107
@ -37614,15 +37614,15 @@ C -0.51958 0.894634 -0.265203 1 0 1
z
"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#mce91a6f3f8" x="83.77363" y="95.031176" style="fill: #1f77b4"/>
<use xlink:href="#mce91a6f3f8" x="331.845341" y="81.881681" style="fill: #1f77b4"/>
<use xlink:href="#mce91a6f3f8" x="579.917052" y="81.693128" style="fill: #1f77b4"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#meb3df71cc8" x="83.77363" y="95.031176" style="fill: #1f77b4"/>
<use xlink:href="#meb3df71cc8" x="331.845341" y="81.881681" style="fill: #1f77b4"/>
<use xlink:href="#meb3df71cc8" x="579.917052" y="81.693128" style="fill: #1f77b4"/>
</g>
</g>
<g id="line2d_21">
<defs>
<path id="m7032e5276b" d="M 0 1
<path id="mcd122ec835" d="M 0 1
C 0.265203 1 0.51958 0.894634 0.707107 0.707107
C 0.894634 0.51958 1 0.265203 1 0
C 1 -0.265203 0.894634 -0.51958 0.707107 -0.707107
@ -37634,15 +37634,15 @@ C -0.51958 0.894634 -0.265203 1 0 1
z
"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#m7032e5276b" x="83.77363" y="95.031176" style="fill: #ff7f0e"/>
<use xlink:href="#m7032e5276b" x="331.845341" y="70.373313" style="fill: #ff7f0e"/>
<use xlink:href="#m7032e5276b" x="579.917052" y="44.828723" style="fill: #ff7f0e"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#mcd122ec835" x="83.77363" y="95.031176" style="fill: #ff7f0e"/>
<use xlink:href="#mcd122ec835" x="331.845341" y="70.373313" style="fill: #ff7f0e"/>
<use xlink:href="#mcd122ec835" x="579.917052" y="44.828723" style="fill: #ff7f0e"/>
</g>
</g>
<g id="line2d_22">
<defs>
<path id="m9430796c91" d="M 0 1
<path id="m65faea4cff" d="M 0 1
C 0.265203 1 0.51958 0.894634 0.707107 0.707107
C 0.894634 0.51958 1 0.265203 1 0
C 1 -0.265203 0.894634 -0.51958 0.707107 -0.707107
@ -37654,15 +37654,15 @@ C -0.51958 0.894634 -0.265203 1 0 1
z
"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#m9430796c91" x="83.77363" y="95.031176" style="fill: #2ca02c"/>
<use xlink:href="#m9430796c91" x="331.845341" y="95.031176" style="fill: #2ca02c"/>
<use xlink:href="#m9430796c91" x="579.917052" y="95.031176" style="fill: #2ca02c"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#m65faea4cff" x="83.77363" y="95.031176" style="fill: #2ca02c"/>
<use xlink:href="#m65faea4cff" x="331.845341" y="95.031176" style="fill: #2ca02c"/>
<use xlink:href="#m65faea4cff" x="579.917052" y="95.031176" style="fill: #2ca02c"/>
</g>
</g>
<g id="line2d_23">
<defs>
<path id="mabe4a44535" d="M 0 1
<path id="m6451126cae" d="M 0 1
C 0.265203 1 0.51958 0.894634 0.707107 0.707107
C 0.894634 0.51958 1 0.265203 1 0
C 1 -0.265203 0.894634 -0.51958 0.707107 -0.707107
@ -37674,15 +37674,15 @@ C -0.51958 0.894634 -0.265203 1 0 1
z
"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#mabe4a44535" x="83.77363" y="95.031176" style="fill: #d62728"/>
<use xlink:href="#mabe4a44535" x="331.845341" y="163.980569" style="fill: #d62728"/>
<use xlink:href="#mabe4a44535" x="579.917052" y="206.557142" style="fill: #d62728"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#m6451126cae" x="83.77363" y="95.031176" style="fill: #d62728"/>
<use xlink:href="#m6451126cae" x="331.845341" y="163.980569" style="fill: #d62728"/>
<use xlink:href="#m6451126cae" x="579.917052" y="206.557142" style="fill: #d62728"/>
</g>
</g>
<g id="line2d_24">
<defs>
<path id="mfa1085208e" d="M 0 1
<path id="mcb579a4d54" d="M 0 1
C 0.265203 1 0.51958 0.894634 0.707107 0.707107
C 0.894634 0.51958 1 0.265203 1 0
C 1 -0.265203 0.894634 -0.51958 0.707107 -0.707107
@ -37694,19 +37694,19 @@ C -0.51958 0.894634 -0.265203 1 0 1
z
"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#mfa1085208e" x="83.77363" y="95.031176" style="fill: #9467bd"/>
<use xlink:href="#mfa1085208e" x="331.845341" y="95.031176" style="fill: #9467bd"/>
<use xlink:href="#mfa1085208e" x="579.917052" y="95.031176" style="fill: #9467bd"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#mcb579a4d54" x="83.77363" y="95.031176" style="fill: #9467bd"/>
<use xlink:href="#mcb579a4d54" x="331.845341" y="95.031176" style="fill: #9467bd"/>
<use xlink:href="#mcb579a4d54" x="579.917052" y="95.031176" style="fill: #9467bd"/>
</g>
</g>
<g id="line2d_25">
<path d="M 83.77363 95.031176
L 331.845341 126.03333
L 579.917052 142.988409
" clip-path="url(#pd926ee2045)" style="fill: none; stroke-dasharray: 1.85,0.8; stroke-dashoffset: 0; stroke: #8c564b; stroke-width: 0.5"/>
" clip-path="url(#p8dfadca05e)" style="fill: none; stroke-dasharray: 1.85,0.8; stroke-dashoffset: 0; stroke: #8c564b; stroke-width: 0.5"/>
<defs>
<path id="mddbdc8c8c4" d="M 0 1
<path id="mc4d3c3b078" d="M 0 1
C 0.265203 1 0.51958 0.894634 0.707107 0.707107
C 0.894634 0.51958 1 0.265203 1 0
C 1 -0.265203 0.894634 -0.51958 0.707107 -0.707107
@ -37718,10 +37718,10 @@ C -0.51958 0.894634 -0.265203 1 0 1
z
"/>
</defs>
<g clip-path="url(#pd926ee2045)">
<use xlink:href="#mddbdc8c8c4" x="83.77363" y="95.031176" style="fill: #8c564b"/>
<use xlink:href="#mddbdc8c8c4" x="331.845341" y="126.03333" style="fill: #8c564b"/>
<use xlink:href="#mddbdc8c8c4" x="579.917052" y="142.988409" style="fill: #8c564b"/>
<g clip-path="url(#p8dfadca05e)">
<use xlink:href="#mc4d3c3b078" x="83.77363" y="95.031176" style="fill: #8c564b"/>
<use xlink:href="#mc4d3c3b078" x="331.845341" y="126.03333" style="fill: #8c564b"/>
<use xlink:href="#mc4d3c3b078" x="579.917052" y="142.988409" style="fill: #8c564b"/>
</g>
</g>
<g id="patch_3">
@ -37766,7 +37766,7 @@ L 63.662873 174.835821
<g id="line2d_26"/>
<g id="line2d_27">
<g>
<use xlink:href="#mce91a6f3f8" x="63.662873" y="179.335821" style="fill: #1f77b4"/>
<use xlink:href="#meb3df71cc8" x="63.662873" y="179.335821" style="fill: #1f77b4"/>
</g>
</g>
<g id="text_16">
@ -37900,7 +37900,7 @@ L 63.662873 188.302197
<g id="line2d_28"/>
<g id="line2d_29">
<g>
<use xlink:href="#m7032e5276b" x="63.662873" y="192.802197" style="fill: #ff7f0e"/>
<use xlink:href="#mcd122ec835" x="63.662873" y="192.802197" style="fill: #ff7f0e"/>
</g>
</g>
<g id="text_17">
@ -38074,7 +38074,7 @@ L 63.662873 201.768574
<g id="line2d_30"/>
<g id="line2d_31">
<g>
<use xlink:href="#m9430796c91" x="63.662873" y="206.268574" style="fill: #2ca02c"/>
<use xlink:href="#m65faea4cff" x="63.662873" y="206.268574" style="fill: #2ca02c"/>
</g>
</g>
<g id="text_18">
@ -38127,7 +38127,7 @@ L 63.662873 215.23495
<g id="line2d_32"/>
<g id="line2d_33">
<g>
<use xlink:href="#mabe4a44535" x="63.662873" y="219.73495" style="fill: #d62728"/>
<use xlink:href="#m6451126cae" x="63.662873" y="219.73495" style="fill: #d62728"/>
</g>
</g>
<g id="text_19">
@ -38203,7 +38203,7 @@ L 63.662873 228.701326
<g id="line2d_34"/>
<g id="line2d_35">
<g>
<use xlink:href="#mfa1085208e" x="63.662873" y="233.201326" style="fill: #9467bd"/>
<use xlink:href="#mcb579a4d54" x="63.662873" y="233.201326" style="fill: #9467bd"/>
</g>
</g>
<g id="text_20">
@ -38232,7 +38232,7 @@ L 72.662873 246.667702
</g>
<g id="line2d_37">
<g>
<use xlink:href="#mddbdc8c8c4" x="63.662873" y="246.667702" style="fill: #8c564b"/>
<use xlink:href="#mc4d3c3b078" x="63.662873" y="246.667702" style="fill: #8c564b"/>
</g>
</g>
<g id="text_21">
@ -38248,7 +38248,7 @@ L 72.662873 246.667702
</g>
</g>
<defs>
<clipPath id="pd926ee2045">
<clipPath id="p8dfadca05e">
<rect x="46.562873" y="7.2" width="818.636647" height="252.959296"/>
</clipPath>
</defs>

Before

Width:  |  Height:  |  Size: 968 KiB

After

Width:  |  Height:  |  Size: 968 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 44 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 42 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 43 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 81 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 40 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 816 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 814 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 50 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 812 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 398 KiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 3.1 MiB

File diff suppressed because it is too large Load diff

After

Width:  |  Height:  |  Size: 809 KiB

View file

@ -1,4 +1,4 @@
#+PROPERTY: header-args :session otto_cycle_length :kernel python :pandoc no :async yes :tangle tangle/coupling_strength.py
#+PROPERTY: header-args :session otto_cycle_length :kernel python :pandoc no :async yes :tangle tangle/coupling_strength.py :noweb yes
Here we scan the cycle length and coupling strength.
@ -33,6 +33,7 @@ We take the same baseline as in [[id:c06111fd-d719-433d-a316-c163f6e1d384][cycle
But we vary the cycle speed while keeping a fixed proportion
coupling-change/cycle time.
#+NAME: make-model
#+begin_src jupyter-python :results none
def make_model(Θ, δ):
(p_H, p_L) = ot.timings(.06, .06)
@ -93,18 +94,33 @@ coupling-change/cycle time.
#+RESULTS:
:RESULTS:
| <Figure | size | 1200x400 | with | 1 | Axes> | <AxesSubplot: | xlabel= | $\tau$ | ylabel= | Energy | > |
[[file:./.ob-jupyter/dcb84273573f02eab9be5cb49c8855b61ef57f84.svg]]
: (pid=gcs_server) [2023-05-16 18:24:06,597 E 443560 443560] (gcs_server) gcs_server.cc:283: Failed to get the resource load: GrpcUnavailable: RPC Error message: Connection reset by peer; RPC Error details:
: (pid=gcs_server) [2023-05-16 18:24:06,597 E 443560 443560] (gcs_server) gcs_server.cc:283: Failed to get the resource load: GrpcUnavailable: RPC Error message: Connection reset by peer; RPC Error details:
: (pid=gcs_server) [2023-05-16 18:24:06,597 E 443560 443560] (gcs_server) gcs_server.cc:283: Failed to get the resource load: GrpcUnavailable: RPC Error message: Connection reset by peer; RPC Error details:
: (pid=gcs_server) [2023-05-16 18:24:06,597 E 443560 443560] (gcs_server) gcs_server.cc:283: Failed to get the resource load: GrpcUnavailable: RPC Error message: Connection reset by peer; RPC Error details:
: (pid=gcs_server) [2023-05-16 18:25:13,636 E 443560 443560] (gcs_server) gcs_server.cc:283: Failed to get the resource load: GrpcUnavailable: RPC Error message: Connection reset by peer; RPC Error details:
: (pid=gcs_server) [2023-05-16 18:25:13,636 E 443560 443560] (gcs_server) gcs_server.cc:283: Failed to get the resource load: GrpcUnavailable: RPC Error message: Connection reset by peer; RPC Error details:
: (pid=gcs_server) [2023-05-16 18:25:13,636 E 443560 443560] (gcs_server) gcs_server.cc:283: Failed to get the resource load: GrpcUnavailable: RPC Error message: Connection reset by peer; RPC Error details:
: (pid=gcs_server) [2023-05-16 18:25:13,637 E 443560 443560] (gcs_server) gcs_server.cc:283: Failed to get the resource load: GrpcUnavailable: RPC Error message: Connection reset by peer; RPC Error details:
: (pid=gcs_server) [2023-05-16 18:25:13,637 E 443560 443560] (gcs_server) gcs_server.cc:283: Failed to get the resource load: GrpcUnavailable: RPC Error message: Connection reset by peer; RPC Error details:
[[file:./.ob-jupyter/164cf59963539ffe267c06050906eacc5c84a306.svg]]
:END:
#+begin_src jupyter-python
for i in range(len(Θs)):
fig, ax = plt.subplots()
for model in models[len(δs) * i :len(δs) * (i+1)]:
pu.plot_with_σ(model.t, model.system_energy(), ax=ax)
#+end_src
#+RESULTS:
:RESULTS:
[[file:./.ob-jupyter/e260500c9f6e3ee5bffaa8e55362d765f315dbd9.svg]]
[[file:./.ob-jupyter/b689998cdfdd74781109afb6a523502318876d06.svg]]
[[file:./.ob-jupyter/88fcef6f0d71dbae4e9a3e3cbaa67b1a3db1363c.svg]]
[[file:./.ob-jupyter/73c5fc49ba925c3d94b2261ba9acb883c7d63471.svg]]
:END:
#+begin_src jupyter-python
for model in models:
plt.plot(model.t, abs(model.total_energy_from_power().value - model.total_energy().value))
#+end_src
#+RESULTS:
[[file:./.ob-jupyter/c273a43c171394f3692125872ff6c55120ae3433.svg]]
#+begin_src jupyter-python
#[model.efficiency(steady_idx=2).value * 100 for model in models][10]
@ -129,7 +145,7 @@ coupling-change/cycle time.
#+RESULTS:
:RESULTS:
| <Figure | size | 1200x400 | with | 2 | Axes> | (<AxesSubplot: xlabel= $N$ ylabel= $P$ > <AxesSubplot: xlabel= $N$ ylabel= $\eta$ >) |
[[file:./.ob-jupyter/b83fcdfb93a1d0a389b5cd9da2636a6608fa8984.svg]]
[[file:./.ob-jupyter/9153bc446bd9633e64c8e85412f9927be1726228.svg]]
:END:
@ -237,11 +253,11 @@ coupling-change/cycle time.
#+RESULTS:
:RESULTS:
[[file:./.ob-jupyter/ee03f1655e72e9136c7065f01d0b7ba04b650017.svg]]
[[file:./.ob-jupyter/f1a1d3e17092ad84d554061c23b8309c424d0de1.svg]]
[[file:./.ob-jupyter/06a89fcc4c201fe7e9c39443a887282b76e54249.svg]]
[[file:./.ob-jupyter/b03717337099eef176d08e0082ceeb6305a8acae.svg]]
[[file:./.ob-jupyter/2c31e3832d8aaecafdf50e0f29694c30c77591c3.svg]]
[[file:./.ob-jupyter/2f0b0c42049806b930471968f418c6e09a12bb14.svg]]
[[file:./.ob-jupyter/37e8a4ed48d03590f17e8fa905d06a079ef9489d.svg]]
[[file:./.ob-jupyter/58030ee1b1a60e263a7ecff46e5e0978c0ebf080.svg]]
[[file:./.ob-jupyter/3e5f7ce49106a6dccaec67c139a3960c3b5977a3.svg]]
[[file:./.ob-jupyter/0f0039092a2d1bba7e2ee7161fd53005eea9c3d3.svg]]
:END:
#+begin_src jupyter-python
@ -275,13 +291,52 @@ coupling-change/cycle time.
#+end_src
#+RESULTS:
\[[file:./.ob-jupyter/1b3959a046911e0ebcbce3e6d6f5df7f5867ffbb.svg]]
[[file:./.ob-jupyter/40d7d6c96595d9d2028f7c9d6b5bb9e1646e07b4.svg]]
#+begin_src jupyter-python
* Weak Coupling Limit
#+begin_src jupyter-python :tangle tangle/weak_coupling_model.py
<<boilerplate>>
<<make-model>>
#+end_src
#+RESULTS:
#+begin_src jupyter-python :results none :tangle tangle/weak_coupling_model.py
weak_coupling_model = make_model(50, .01)
weak_coupling_model.num_cycles = 10
weak_coupling_model.k_max = 2
#+end_src
** Integrate
#+begin_src jupyter-python :tangle tangle/weak_coupling_model.py
ot.integrate_online(weak_coupling_model, 10, analyze_kwargs=dict(every=10_000))
#+end_src
#+RESULTS:
: [INFO hops.core.integration 137953] Choosing the nonlinear integrator.
: [INFO root 137953] Starting analysis process.
: [INFO root 137953] Started analysis process with pid 142086.
: [INFO hops.core.hierarchy_data 137953] Creating the streaming fifo at: /home/hiro/Documents/Projects/UNI/master/eflow_paper/python/otto_motor/subprojects/cycle_length_coupling_strength/results_84524c4b2c9c8428030cfd17a749cdee70253c3d428b2416de22620c4b230cbd.fifo
: [INFO hops.core.integration 137953] Using 16 integrators.
: [INFO hops.core.integration 137953] Some 10 trajectories have to be integrated.
: [INFO hops.core.integration 137953] Using 66 hierarchy states.
: 100% 10/10 [01:56<00:00, 11.60s/it]
:
** Analysis
#+begin_src jupyter-python
ot.plot_energy(weak_coupling_model)
#+end_src
#+RESULTS:
:RESULTS:
| <Figure | size | 1200x400 | with | 1 | Axes> | <AxesSubplot: | xlabel= | $\tau$ | ylabel= | Energy | > |
[[file:./.ob-jupyter/cee1630d28ef76e85c16f224fc9c6516c3680e3f.svg]]
:END:
* Things to Look At
- power and efficiency
- magnitude of interaction energies

View file

@ -56,6 +56,14 @@ ot.integrate_online_multi(models, 50_000, increment=10_000, analyze_kwargs=dict(
ot.plot_energy(models[5])
for i in range(len(Θs)):
fig, ax = plt.subplots()
for model in models[len(δs) * i :len(δs) * (i+1)]:
pu.plot_with_σ(model.t, model.system_energy(), ax=ax)
for model in models:
plt.plot(model.t, abs(model.total_energy_from_power().value - model.total_energy().value))
#[model.efficiency(steady_idx=2).value * 100 for model in models][10]
models[10].strobe, models[1].strobe
@ -191,4 +199,4 @@ ot.plot_3d_heatmap(
a_efficiency.set_zlabel(r"$\sigma_\eta/|\eta|$")
fs.export_fig("coupling_speed_scan_power_efficiency_uncertainty")
ot.plot_energy(weak_coupling_model)

View file

@ -0,0 +1,51 @@
import figsaver as fs
import plot_utils as pu
from hiro_models.one_qubit_model import StocProcTolerances
from hiro_models.otto_cycle import OttoEngine
import hiro_models.model_auxiliary as aux
import numpy as np
import qutip as qt
import utilities as ut
import stocproc
import matplotlib.pyplot as plt
import otto_utilities as ot
import ray
ray.shutdown()
#ray.init(address='auto')
ray.init()
from hops.util.logging_setup import logging_setup
import logging
logging_setup(logging.INFO)
plt.rcParams['figure.figsize'] = (12,4)
def make_model(Θ, δ):
(p_H, p_L) = ot.timings(.06, .06)
return OttoEngine(
δ=[δ, δ],
ω_c=[1, 1],
ψ_0=qt.basis([2], [1]),
description=f"Classic Cycle",
k_max=4,
bcf_terms=[5] * 2,
truncation_scheme="simplex",
driving_process_tolerances=[StocProcTolerances(1e-3, 1e-3)] * 2,
thermal_process_tolerances=[StocProcTolerances(1e-3, 1e-3)] * 2,
T=[0.5, 4],
therm_methods=["tanhsinh", "tanhsinh"],
Δ=1,
num_cycles=3,
Θ=Θ,
dt=0.001,
timings_H=p_H,
timings_L=p_L,
streaming_mode=True,
shift_to_resonance=(False, False),
L_shift=(0, 0),
)
weak_coupling_model = make_model(50, .01)
weak_coupling_model.num_cycles = 10
weak_coupling_model.k_max = 2
ot.integrate_online(weak_coupling_model, 10, analyze_kwargs=dict(every=10_000))