大(27GB)文件的grep功能更快

我必须从一个文件(27GB)的文件(5MB)grep包含特定的字符串相同的字符串(和其他信息).
为了加快分析,我将27GB文件分解成1GB文件,然后应用以下脚本(在这里的一些人的帮助下).但是它不是很有效(生成一个180KB的文件需要30个小时!).

这是脚本.有没有比grep更合适的工具?还是使用grep更有效的方式?

#!/bin/bash

NR_cpuS=4
count=0


for z in `echo {a..z}` ;
do
 for x in `echo {a..z}` ;
 do
  for y in `echo {a..z}` ;
  do
   for ids in $(cat input.sam|awk '{print $1}');  
   do 
    grep $ids sample_"$z""$x""$y"|awk '{print $1" "$10" "$11}' >> output.txt &
    let count+=1
                                [[ $((count%NR_cpuS)) -eq 0 ]] && wait
   done
  done #&

解决方法

你可以尝试几件事:

1)你正在读取input.sam多次.它只需要在你的第一个循环开始之前被读取一次.将这个id保存到一个临时文件中,这个文件将被grep读取.

2)使用LC_ALL = C将您的grep命令前缀使用C语言环境而不是UTF-8.这样会加快grep.

3)使用fgrep,因为您正在搜索固定的字符串,而不是正则表达式.

4)使用-f使grep从文件中读取模式,而不是使用循环.

5)不要从多个进程写入输出文件,因为可能会导致行交错和损坏的文件.

进行这些更改后,您的脚本将成为:

awk '{print $1}' input.sam > idsFile.txt
for z in {a..z}
do
 for x in {a..z}
 do
  for y in {a..z}
  do
    LC_ALL=C fgrep -f idsFile.txt sample_"$z""$x""$y" | awk '{print $1,$10,$11}'
  done >> output.txt

另外,查看GNU Parallel是为了帮助您并行运行作业而设计的.

相关文章

Linux中的ARP防火墙主要用于防御ARP欺骗攻击,其效果取决于多...
insmod和modprobe加-f参数导致Invalid module format错误 这...
将ArchLinux安装到U盘 几个月前入门Arch的时候上网搜了不少安...
1、安装Apache。 1)执行如下命令,安装Apache服务及其扩展包...
一、先说一下用ansible批量采集机器信息的实现办法: 1、先把...
安装配置 1. 安装vsftpd 检查是否安装了vsftpd # rpm -qa | ...