У меня есть папка с множеством файлов (xyz1, xyz2, вплоть до xyz5025), и мне нужно запустить скрипт для каждого из них, получив xyz1.faa, xyz2.faa и т. Д. В качестве выходных данных.
Команда для одного файла:
./transeq xyz1 xyz1.faa -table 11Есть ли способ сделать это автоматически? Может быть, комбо?
Ответы:
for file in xyz*
do
./transeq "$file" "${file}.faa" -table 11
doneЭто простой forцикл, который будет перебирать каждый файл, который начинается xyzв текущем каталоге, и вызывать ./transeqпрограмму с именем файла в качестве первого аргумента, именем файла, за которым следует «.faa» в качестве второго аргумента, за которым следует «-table 11» ,
"$file".faaего немного легче вводить как часть интерактивной однострочной строки, и он безопасен, поскольку .faaне содержит метасимволов оболочки, которые необходимо заключать в кавычки.
xyz*глобус также подберет файлы .faa. Для bash запустите shopt -s extglob( ссылка ), а затем используйте, for file in xyz!(*.faa) ...чтобы исключить файлы .faa от отправки через цикл.
Если вы устанавливаете GNU Parallel, вы можете сделать это параллельно так:
parallel ./transeq {} {}.faa -table 11 ::: xyz*Если ваша программа интенсивно использует процессор, она должна немного ускориться.
Вы можете сделать что-то вроде этого в bashкомандной строке:
printf '%s\n' {1..5025} | xargs -l -I {} -t ./transeq xyz{} xyz{}.faa -table 11Мы генерируем целые числа от 1 до 5025, по одной на строку, затем подаем их один за другим в xargs, который инкапсулирует целое число в {}и затем трансплантирует его в командную строку ./transeq соответствующим образом.
Если у вас нет возможности расширения скобок, {n..m}вы можете вызвать seqутилиту для генерации этих чисел.
Или вы всегда можете эмулировать генерацию чисел с помощью:
yes | sed -n =\;5025q | xargs ...for i in {1..5025}; do ./transeq "xyz$i" "xyz$i".faa -table 11; doneэто гораздо проще думать и печатать. Если вы хотите, чтобы он печатал команды перед их выполнением, используйте set -x.
for i in {1..5025}для достижения точно такого же результата, как ваш. Вы также можете писать for ((i=1 ; i<=5025 ; i++)); do ./transeq "xyz$i" "xyz$i".faa -table 11; doneна bash, но я обычно использую {a..b}синтаксис диапазона, потому что он быстрее печатается.
Предполагая, что у вас более одного ядра, и каждый вызов может выполняться независимо от остальных, вы получите значительное ускорение при параллельных запусках.
Относительно простой способ сделать это с помощью -Pпараметра xargs- например, если у вас есть 4 ядра:
echo xyz{1..5025} | \
xargs -n 1 -P 4 -I{} /path/to/transeq xyz{} xyz{}.faa -table 11-n 1Говорит , xargsчтобы выбрать только один аргумент из списка для каждого вызова (по умолчанию он будет проходить много) , и -P 4говорит ему , чтобы породить 4 процессов одновременно - когда кто -то умирает, а новая породившие.
ИМХО, вам не нужно устанавливать GNU параллельно для этого простого случая - xargsдостаточно.
for file in xyz*; do ./transeq "$file" "${file}.faa" -table 11; done. Я пишу такие вещи все время. И если вы хотите убедиться, что имена файлов и т. Д. Расширяются так, как вы хотите, просто поставьтеechoсразу послеdoпервого раза, а затем вернитесь в историю оболочки и удалите ее во второй раз.