#!/usr/bin/python

import sys
import getopt
import re def compare(f1,f2,o1,o2,si_line):
lines_count=0;
in1 = open(f1,"r")
in2 = open(f2,"r")
ou1 = open(o1,"w")
ou2 = open(o2,"w")
while 1:
line1 = in1.readline()
line2 = in2.readline()
if line1 and line2:
lines_count=lines_count+1
print >>ou1,re.sub(r'\n',"",line1)
ou2.write(line2)
if si_line<lines_count:
break
else:
break
f1.close()
f2.close()
o1.close()
o2.close()
return lines_count def main():
opts,args = getopt.getopt(sys.argv[1:],'h:',['fq1=','fq2=','out1=','out2=','Rlen=','Dataset='])
fq1=opts[0][1]
fq2=opts[1][1]
out1=opts[2][1]
out2=opts[3][1]
readlength = opts[4][1]
datasize = opts[5][1]
size_lines = 4*int(datasize)*1e9/(int(readlength)*2)
count = compare(fq1,fq2,out1,out2,size_lines)-1
real_size = (count/4)*2*int(readlength)/1e9
print datasize
print real_size if __name__ == "__main__":
main()

最新文章

  1. 玩转spring boot——结合AngularJs和JDBC
  2. 看看C# 6.0中那些语法糖都干了些什么(中篇)
  3. inline--行内元素
  4. 逻辑操作符“&amp;&amp;”的三层理解
  5. 一款可以下拉搜索html下拉框控件
  6. Qt ini配置文件读写
  7. 以Debug模式启动JBoss
  8. mysql备份还原数据库
  9. Linux及安全——ELF实践
  10. 解决phpMyAdmin中缺少mysqli扩展的错误
  11. 无线端通用的reset样式
  12. 如何在我们项目中利用开源的图表(js chart)
  13. iOS学习笔记--OC系列(1)
  14. Spark PySpark数据类型的转换原理—Writable Converter
  15. Ibatis代码自动生成工具——Abator安装与应用实例(图解)
  16. HTML标签CSS属性默认值汇总
  17. inception cenOS 安装
  18. java:从指定问价中读取80个字节写入指定文件中
  19. PHP下载文件的几种方案
  20. Slickflow.NET 开源工作流引擎基础介绍(二) -- 引擎组件和业务系统的集成

热门文章

  1. 题目(或游戏)流程控制器上传到GitHub
  2. P2P平台被清盘后,你会怎样捍卫自身利益?
  3. 新发布的廉价版iPhoneXR值得购买吗?
  4. Idea的类中使用实体类(有@Data注解)的Get/Set方法报错
  5. js数组,字符转换;key_value获取
  6. UVA 11584 入门DP
  7. ajax+springboot完整例子
  8. uniapp 初始化项目
  9. C语言-再论指针与数组
  10. 剑指offer【08】- 二叉树的深度(java)