This thread is locked.Only browsing is available.
Top Page > Browsing
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
e
Page: [1]