Warning: Non-fatal error in glmnet library call: error code = -1

[Yikai-Wang]

Hi there,
I am using your code to solve logistic regression using the code

N = 80
d = 5
C = 5
pfac = scipy.ones([d+N,C])
pfac[:d,:C]=0
x = np.random.randn(N,d)
i = np.eye(N)
x = np.concatenate([x,i],axis=1)
y = np.random.randint(0,5,size=[N,])
y = y.astype('float')
fit = glmnet(x = x.copy(), y = y.copy(), family = 'multinomial', mtype = 'grouped', penalty_factor=pfac)

The above code works on the synthesis data. However, when I use it in the my own x and y. It prints out:

Warning: Non-fatal error in glmnet library call: error code =
 -1
Check results for accuracy. Partial or no results returned.

which should be the code of

glmnet_python/glmnet_python/lognet.py

Lines 225 to 227 in 645deee

	elif (jerr_r.value < 0):
	print("Warning: Non-fatal error in glmnet library call: error code = ", jerr_r.value)
	print("Check results for accuracy. Partial or no results returned.")

And send out the error of

ValueError: negative dimensions are not allowed

In

glmnet_python/glmnet_python/lognet.py

Line 275 in 645deee

dfmat[k, :] = scipy.zeros([1, lmu], dtype = scipy.float64)

The reason is that the calculated lmu=-1.
I want to locate the bug but find it was computed in

glmnet_python/glmnet_python/lognet.py

Lines 185 to 216 in 645deee

	glmlib.lognet_(
	ctypes.byref(ctypes.c_double(parm)),
	ctypes.byref(ctypes.c_int(nobs)),
	ctypes.byref(ctypes.c_int(nvars)),
	ctypes.byref(ctypes.c_int(nc)),
	x.ctypes.data_as(ctypes.POINTER(ctypes.c_double)),
	y.ctypes.data_as(ctypes.POINTER(ctypes.c_double)),
	offset.ctypes.data_as(ctypes.POINTER(ctypes.c_double)),
	jd.ctypes.data_as(ctypes.POINTER(ctypes.c_int)),
	vp.ctypes.data_as(ctypes.POINTER(ctypes.c_double)),
	cl.ctypes.data_as(ctypes.POINTER(ctypes.c_double)),
	ctypes.byref(ctypes.c_int(ne)),
	ctypes.byref(ctypes.c_int(nx)),
	ctypes.byref(ctypes.c_int(nlam)),
	ctypes.byref(ctypes.c_double(flmin)),
	ulam.ctypes.data_as(ctypes.POINTER(ctypes.c_double)),
	ctypes.byref(ctypes.c_double(thresh)),
	ctypes.byref(ctypes.c_int(isd)),
	ctypes.byref(ctypes.c_int(intr)),
	ctypes.byref(ctypes.c_int(maxit)),
	ctypes.byref(ctypes.c_int(kopt)),
	ctypes.byref(lmu_r),
	a0_r,
	ca_r,
	ia_r,
	nin_r,
	ctypes.byref(dev0_r),
	dev_r,
	alm_r,
	ctypes.byref(nlp_r),
	ctypes.byref(jerr_r)
	)

I am not familiar with Fortran, so it is hard for me to find out the error.
Could you help me?