[julia-users] Segfault when passing large arrays to Julia from C

[David Smith]

The following code works as posted, but when I increase N to 10000000, I 
get a segfault. I don't see why that would break anything, as the size_t 
type should be large enough to handle that.

Can someone point me in the right direction?

Thanks!


#include <julia.h> 
#include <stdio.h> 

int main() 
{ 
    size_t N = 1000000; 
    jl_init(NULL); 
    JL_SET_STACK_BASE; 
    jl_value_t* array_type = jl_apply_array_type( jl_float64_type, 1 ); 
    jl_array_t* x = jl_alloc_array_1d(array_type, N); 
    jl_array_t* y = jl_alloc_array_1d(array_type, N); 
    JL_GC_PUSH2(&x, &y); 
    double* xData = jl_array_data(x); 
    double* yData = jl_array_data(y); 
    for (size_t i=0; i<jl_array_len(x); i++){ 
        xData[i] = i; 
        yData[i] = i; 
    } 
    jl_eval_string("myfft(x,y) = fft(complex(x,y))"); 
    jl_function_t *func = jl_get_function(jl_current_module, "myfft"); 
    jl_value_t* jlret = jl_call2(func, (jl_value_t*) x, (jl_value_t*)y); 
    double *ret = jl_unbox_voidpointer(jlret); 
    for(size_t i=0; i<10; i++) 
        printf("(%f,%f) = %f + %fi\n", xData[i], yData[i], ret[2*i], 
ret[2*i+1]); 
    JL_GC_POP(); 
    return 0; 
}