segfault in gcta::fit_bivar_reml

[cribe78]

I am a HPC support engineer, not the end user of the software, but I am reporting a bug that one of the researchers I work with has encountered. The user was encountering the segfault in gcta versions 1.94.3 and 1.95.0, I have built a debug version of gcta from source and replicated the error. I do not have permission to post the input files here, but I will share as much information as I can.

The command with results in the error is:
gcta64 --grm SCP --pheno cardiac_surveyPheno.txt --reml-bivar 4 5 --thread-num 4 --out trait4trait5genetic_corr_results

Here is a backtrace from the error:

Program terminated with signal SIGSEGV, Segmentation fault.

--Type <RET> for more, q to quit, c to continue without paging--
Missing separate debuginfos, use: dnf debuginfo-install glibc-2.34-125.el9_5.8.x86_64 libgcc-11.5.0-5.el9_5.x86_64 libgomp-11.5.0-5.el9_5.x86_64 libstdc++-11.5.0-5.el9_5.x86_64 libzstd-1.5.1-2.el9.x86_64 zlib-1.2.11-40.el9.x86_64#0  gcta::fit_bivar_reml (this=0xf8cbb0, grm_file=..., phen_file=..., qcovar_file=..., covar_file=..., keep_indi_file=..., remove_indi_file=..., sex_file=...,
    mphen=3, mphen2=4, grm_cutoff=-2, adj_grm_fac=-2, dosage_compen=-2, m_grm_flag=false, pred_rand_eff=false, est_fix_eff=false, reml_mtd=0, MaxIter=100,
    reml_priors=..., reml_priors_var=..., drop=..., no_lrt=true, prevalence=-1, prevalence2=-1, no_constrain=false, ignore_Ce=false, fixed_rg_val=...,
    bivar_no_constrain=false) at /home/chrisribe/r/78173/gcta/main/bivar_reml.cpp:172
172                     (_Asp[pos]).insertBackUncompressed(i, j) = _grm(_keep[nms1[i]], _keep[nms1[j]]);
[Current thread is 1 (Thread 0x7fb5d2ba77c0 (LWP 2105711))]
(gdb) bt
#0  gcta::fit_bivar_reml (this=0xf8cbb0, grm_file="SCP", phen_file="cardiac_surveyPheno.txt", qcovar_file="", covar_file="", keep_indi_file="",
    remove_indi_file="", sex_file="", mphen=3, mphen2=4, grm_cutoff=-2, adj_grm_fac=-2, dosage_compen=-2, m_grm_flag=false, pred_rand_eff=false,
    est_fix_eff=false, reml_mtd=0, MaxIter=100, reml_priors=std::vector of length 0, capacity 0, reml_priors_var=std::vector of length 0, capacity 0,
    drop=std::vector of length 1, capacity 1 = {...}, no_lrt=true, prevalence=-1, prevalence2=-1, no_constrain=false, ignore_Ce=false,
    fixed_rg_val=std::vector of length 0, capacity 0, bivar_no_constrain=false) at /home/chrisribe/r/78173/gcta/main/bivar_reml.cpp:172
#1  0x000000000043cd80 in option (option_num=12, option_str=0x7ffe04ed6558) at /home/chrisribe/r/78173/gcta/main/option.cpp:1461
#2  0x0000000000420b23 in main (argc=12, argv=0x7ffe04ed6558) at /home/chrisribe/r/78173/gcta/src/main.cpp:317
(gdb) p i
$1 = 60
(gdb) p j
$2 = 0
(gdb) p pos
$3 = 0