myprojects-hg-reborn

Fork

rm(list=ls())
x1<-read.table("eigenvect001",header=FALSE)
x1<-x1[ ,1]

x2<-read.table("eigenvect002",header=FALSE)
x2<-x2[ ,1]

x3<-read.table("eigenvect003",header=FALSE)
x3<-x3[ ,1]

x4<-read.table("eigenvect004",header=FALSE)
x4<-x4[ ,1]

x5<-read.table("eigenvect005",header=FALSE)
x5<-x5[ ,1]

t<-seq(1:length(x1))

eigval<-read.table("real--eigval",header=FALSE)

eigval<-eigval[ ,1]

eigval<-sort(eigval)



pdf("spectrum.pdf")
plot(t,eigval,xlab="eigenvalue number",ylab="eigenvalue")

dev.off()

ground<-read.table("ground-state",header=FALSE)
ground<-ground[ ,1]

t<-seq(1:length(ground))
pdf("ground.pdf")
plot(t,ground,xlab="site number",ylab="ground state")

dev.off()

imag<-read.table("imag--eigval")

imag<-imag[ ,1]

flag=0
for ( i in 1:length(imag))
{
if (imag[i] != 0)
{
flag=1
}
}

if (flag==1) {
print("Warning: there are some complex eigenvalues");
}
if (flag !=1)
{
print("All the eigenvectors are real");
}
pdf("omega-imag-vs-omega-real.pdf")
plot(eigval, imag,xlab="Re(omega)",ylab="Im(omega)")
dev.off()
# xtest<-read.table("eigenvect3rd",header=FALSE)
# xtest<-xtest[ ,1]