Segfault calling into Julia from C/MPI

I have written a Julia package and compiled a shared object library to call one @ccallable function. When I test it with a serial C program, it works fine. When I call it from a more complex program using MPI (even when run with -np 1), it segfaults in libc.so.6.

Thread 1 "compute-density" received signal SIGSEGV, Segmentation fault.
0x0000000000000000 in ?? ()

I can use gdb to find where this is happening, but I get no information from the Julia part; the last meaningful information I get is where the C program called the Julia function.

Is there anything special that I need to do building the library so that it can be called from MPI? Is there any way to enable gdb to give me Julia information (line numbers, variable values, etc.)?

I found the cause, I had failed to include the init_julia() call at the start of the C program. It would still be nice know how to use gdb in a Julia library running in C.