Мой код в настоящее время выглядит следующим образом:
z = np.diagflat(c).dot(D).dot(x);
idxN, idxP = z<0, z>=0 # logical indexing
y1 = [-1 + np.exp(x)/(1+np.exp(x)) for x in z[idxN]]
y1 = np.array(y1) # Size (504,)
y2 = [-np.exp(-x)/(1+np.exp(-x)) for x in z[idxP]]
y2 = np.array(y2) # Size (496,)
Теперь я пытаюсь сформировать (1000,) список y
, который будет включать y1 или y2 в зависимости от того, равен ли индекс +/-,
Я предполагаю, что понимание списка было бы самым простым, но я борюсь с формой.Спасибо за ваш вклад.