GSL ERROR: function value is not finite in brent.c at line 58 errno 9

[zhoufc213]

This is a simple gwas analysis process. The author wrote a script to do gwas analysis. When running the script, this error occurred. It is worth noting that se (pve) =-nan is calculated, except for the error in line 58.Is there colleagues familiar with this software who can help to solve the problem? Is it possible that the convergence is caused by too few individuals