使用 echo 和 system 在 C 中运行软件

Question

我正在尝试运行一个名为 BLASTP 的生物程序，它接收两个字符串（代码中的 fasta_GWIDD 和 fasta_UNIPROT）并比较它们。我遇到的问题是在代码中使用了 echo/system。谁能建议我错过了什么？

    for(i=0;i<index1;i++)
    {
        sprintf(fasta_GWIDD,">%s\\n%s\n",fasta_name1[i],fasta_seq1[i]);
        setenv("GwiddVar", fasta_GWIDD, 1) ;
        sprintf(fasta_UNIPROT,">%s\\n%s\n",fasta_name2[i],fasta_seq2[i]);
        setenv("UniprotVar", fasta_UNIPROT, 1) ;
        system("blastp -query <(echo -e $GwiddVar) -subject<(echo -e $UniprotVar)");
    }

错误是：

sh: -c: line 0: syntax error near unexpected token `('
sh: -c: line 0: `blastp -query <(echo -e $GwiddVar) -subject<(echo -e $UniprotVar)'

Answer 1

看来shell不理解

<(echo -e $GwiddVar)

句法。请注意，该system命令可能使用与您习惯的不同的 shell（例如 csh 而不是 bash，等等）。它是您操作系统配置文件和配置文件中某处的所有内容，但我无法猜测您在那里拥有什么。

顺便提一句。我认为您应该能够system()通过以下任一方式检查命令正在使用哪个 shell ：

system("echo $SHELL")  // should simply write the path to current shell
system("ps -aux")      // look at it and find what is the parent of the PS

等等。

考虑到这在某些 shell 上是正确的：

blastp -query <(echo -e $GwiddVar) -subject<(echo -e $UniprotVar)

上面引用的语法显然只是为了将变量作为输入传递。我认为你做得太过分了。您正在使用echo -e $GwiddVar打印和捕获数据，这些数据已在您手头的变量中。你有没有尝试过像这样简单的事情：

blastp -query $GwiddVar -subject $UniprotVar

我不知道您尝试使用哪个 shell，但考虑到echo它获取了它的数据，那么它应该完全相同。

如果你担心空格，那么各种 shell 通常允许你使用引号：

blastp -query "$GwiddVar" -subject "$UniprotVar"

当然，这取决于外壳。如果您的程序使用不喜欢引号的 shell，那么您必须对其进行调整。不是针对您的外壳，而是针对system()已使用的外壳。

还有一点就是使用起来system很粗糙。当您有难以正确转义的参数时，您应该使用其他函数，例如 execve，它们能够获取真实的原始直接字符串数组并将它们作为 ARGV 直接传递给进程。使用这些，您将不需要（也不应该）在要传递的字符串中添加任何引号或转义任何空格。

sprintf(fasta_GWIDD,">%s\\n%s\n",fasta_name1[i],fasta_seq1[i]);
sprintf(fasta_UNIPROT,">%s\\n%s\n",fasta_name2[i],fasta_seq2[i]);

char** args = .....; // allocate an array of char*[5], malloc, or whatever
args[0] = "blastp";
args[1] = "-query";
args[2] = fasta_GWIDD;
args[3] = "-subject";
args[4] = fasta_UNIPROT;

int errcode = execve(4, args, null);
if( errcode ) ... // check the error (if any) and react

然而！请注意， execve 来自 exec 系列，因此它取代了您当前的进程。这就是为什么我只写一个草图而不显示整个准备运行的代码。您可能需要fork()在它之前，然后等待外循环中的孩子。

所以，我首先检查外壳和语法;)