c - 使用 OpenMP 时 snprintf 中的段错误

Question

我在某处读到 snprintf 是线程安全的；然而，当我运行我的代码时，它产生了一个段错误。

主要功能:

#pragma omp parallel 
{
    #pragma omp for private(j)
    for (i = 0; i < 8; i++) {
        for (j = 0; j < 1; j++) {
            foo(gene_seq); //gene_seq is a large char array
        }
    }
}

错误的代码行:

double foo(char *gene_seq){
    char *seq;

    /* some stuff above */
    region_length = end_pos-init_pos+1; // include \0 terminator
    seq = (char *)safe_malloc(sizeof(char) * region_length);
    snprintf(seq, region_length, "%s", gene_seq + init_pos); //SIGSEGV!!!
    /*more stuff*/
}

编辑将 snprintf 行周围的程序更新为:

#pragma omp critical
{
    snprintf(seq, region_length, "%s", gene_seq + init_pos);
}

它仍然无法执行。 gdb 为变量提供了以下值:

(gdb) p strlen(gene_seq)
$1 = 1405
(gdb) p init_pos
$2 = 683
(gdb) p region_length
$4 = 221   

当我单线程执行程序时，它运行良好，而 gdb 不是很有帮助。我的 glibc 版本是 2.15。如果需要，将提供更多信息。

Answer 1

您是否错过了 safe_malloc() 原型(prototype)并且是在 64 位上？

如果"is"，编译器假定 safe_alloc() 返回 32 位，因此如果指针的值大于 2^32，则返回一个无效指针。

打开所有警告 -Wall -Wextra -pedantic 并听听编译器告诉你什么。