/contrib/ntp/scripts/stats/etf.S

https://bitbucket.org/freebsd/freebsd-head/ · Assembly · 15 lines · 15 code · 0 blank · 0 comment · 0 complexity · a4689f09592a62cff31e6d780ba9ec3b MD5 · raw file

  1. options(digits=4)
  2. file2 <- "etf_summary"
  3. etf <- scan(file1, list(day=0, sec=0, offset=0, stab=0))
  4. r <- lsfit(etf$sec, etf$offset)
  5. count<-length(etf$sec)
  6. mean<-r$coef[[1]]
  7. std<-sqrt(var(r$residuals))
  8. slope<-r$coef[[2]] * 1000
  9. cat("\n", file=file2 , append=TRUE, fill=FALSE, sep="")
  10. cat(file1, "\n", file=file2, append=TRUE, fill=FALSE, sep="")
  11. cat("etf1 ", count, ", T ", mean, " ns, R ", slope, " ps/s, std ", std, " us\n", file=file2, append=TRUE, fill=FALSE, sep="")
  12. str <- paste("eps/", file1, ".eps", sep="")
  13. postscript(str, , , , 5, pointsize=18)
  14. par(mgp=c(1, 0, 0), tck=0.03, mar=c(2, 2, 1, 1))
  15. plot(etf$sec, etf$offset, type="l", xlab=paste("MJD", etf$day, "Time (s)"), ylab="External Offset (ns)", ylim=c(-400, 400))