Differences
This shows you the differences between two versions of the page.
Next revision | Previous revision | ||
documentation:language_reference:objects:responsefunction:functions:invertenergy [2024/09/16 12:06] – created Maurits W. Haverkort | documentation:language_reference:objects:responsefunction:functions:invertenergy [2024/12/22 20:29] (current) – Maurits W. Haverkort | ||
---|---|---|---|
Line 3: | Line 3: | ||
### | ### | ||
- | alligned paragraph text | + | InvertEnergy(G) takes response function $G(\omega)$ as an argument and returns $-G(-\omega)^*$. As we only store the poles and their residues we can do this by negating the energy of the poles without the need to change the size of the residue |
### | ### | ||
Line 9: | Line 9: | ||
### | ### | ||
- | description text | + | This example initializes a response function in the "list of poles" representation and inverts its energy. |
### | ### | ||
==== Input ==== | ==== Input ==== | ||
<code Quanty Example.Quanty> | <code Quanty Example.Quanty> | ||
- | -- some example code | + | NEDOS = 10 |
+ | HalfBW = 1 | ||
+ | dE =HalfBW/ | ||
+ | |||
+ | a = {0} | ||
+ | b = {} | ||
+ | |||
+ | for i=1,NEDOS do | ||
+ | a[#a+1] = (i-0.5) * dE -- energy axis from 0 to 1 | ||
+ | b[#b+1] = 0.5 / NEDOS -- flat density of states | ||
+ | end | ||
+ | |||
+ | G0 = ResponseFunction.New( {a, | ||
+ | G0_inv = ResponseFunction.InvertEnergy(G0) | ||
+ | |||
+ | print(G0) | ||
+ | print(G0_inv) | ||
</ | </ | ||
==== Result ==== | ==== Result ==== | ||
<file Quanty_Output> | <file Quanty_Output> | ||
- | text produced as output | + | { { 0 , 0.05 , 0.15 , 0.25 , 0.35 , 0.45 , 0.55 , 0.65 , 0.75 , 0.85 , 0.95 } , |
+ | { 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 } , | ||
+ | mu = 0 , | ||
+ | name = G0 , | ||
+ | type = ListOfPoles } | ||
+ | { { 0 , -0.05 , -0.15 , -0.25 , -0.35 , -0.45 , -0.55 , -0.65 , -0.75 , -0.85 , -0.95 } , | ||
+ | { 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 , 0.05 } , | ||
+ | mu = 0 , | ||
+ | name = G0 , | ||
+ | type = ListOfPoles } | ||
</ | </ | ||
===== Table of contents ===== | ===== Table of contents ===== | ||
- | {{indexmenu> | + | {{indexmenu> |