Question about Poles and residues in the ...code |
- Date: 2008/07/03 20:16
- Name: Liger Chen
<b89202066@ntu.edu.tw>
Hi All :
the paper ''Poles and residues in the continued fraction representation of the Fermi-Dirac function'', is with a code for this method(one also could download this in the website).this metod is so cool,but I could not understand some sentence at this :
(begin 151th line)
*******
/* store residue */
for (i=1; i<=n; i++){
zp[i] = 1.0/zp[i];
}
for (i=0; i<n; i++) {
Rp[i+1] = -a[i*n]*a[i*n]*zp[i+1]*zp[i+1]*0.250;
}
*******
zp is the eigenvalue of A*x=y*B*x, so I recognize that the pole is the zp , instead of one over zp ; and a is the eigenector, so residues are just a*zp/4 , why Ozaki squared the numerator ?
May some one ask me , thanks.
Liger Chen
| |
This thread is locked.Only browsing is available.