在C lang中使用GEOS API函数会返回分段错误

Question

我试图使用我的C程序中的一些libgeos（https://github.com/libgeos/geos）函数。以下是我的代码。

#include <stdio.h>
#include <geos_c.h>

int main() {

    GEOSContextHandle_t ctx  =  GEOS_init_r();

    GEOSGeometry *inputGeom = GEOSGeomFromWKT("MULTIPOLYGON (((30 20, 45 40, 10 40, 30 20)), ((15 5, 40 10, 10 20, 5 10, 15 5)))");

    printf("%s\n", GEOSGeomToWKT(inputGeom));

    GEOS_finish_r(ctx);

    return 0;
}

但问题是我运行二进制文件时遇到了分段错误。下面是我用来编译程序的gcc命令。

gcc -g -Wall geos.c -o geos -lgeos_c

请参阅下面的gdb输出。

(gdb) run
Starting program: /home/srimal/proj/src/geos/geos 
warning: the debug information found in "/lib64/ld-2.27.so" does not match "/lib64/ld-linux-x86-64.so.2" (CRC mismatch).


Program received signal SIGSEGV, Segmentation fault.
__strlen_avx2 () at ../sysdeps/x86_64/multiarch/strlen-avx2.S:62
62      ../sysdeps/x86_64/multiarch/strlen-avx2.S: No such file or directory.

如果有人能帮助我解决这个问题，我将非常感激。

Answer 1

解决了这个问题。我应该使用initGEOS和finishGEOS函数。以下是一些改进的正确答案。

#include <stdio.h>
#include <stdlib.h>
#include <stdarg.h>

#include <geos_c.h>

void notice(const char *fmt, ...) {
    va_list ap;
    fprintf( stdout, "NOTICE: ");
    va_start (ap, fmt);
    vfprintf( stdout, fmt, ap);
    va_end(ap);
    fprintf( stdout, "\n" );
}

void log_and_exit(const char *fmt, ...) {
    va_list ap;
    fprintf( stdout, "ERROR: ");    
    va_start (ap, fmt);
    vfprintf( stdout, fmt, ap);
    va_end(ap);
    fprintf( stdout, "\n" );
    exit(1);
}

int main() {

    initGEOS(notice, log_and_exit);

    GEOSGeometry *inputGeom = GEOSGeomFromWKT("POINT (30 10)");

    if (inputGeom == NULL) {
        fprintf(stderr, "NULL GEOM");
        return 0;
    }

    printf("%s\n", GEOSGeomToWKT(inputGeom));

    finishGEOS();
    return 0;
}