def main() :
N = 20
K = 2
W = np.random.rand(N,N,K)
A = W.copy()
for k in range(K) :
A[:,:,k] = np.linalg.pinv(W[:,:,k])
print "This should come out to zero \n", isi.joint_ISI(W,A)
A = np.random.rand(N,N,K)
print "Who knows what this should come out to \n", isi.joint_ISI(W,A)
def main():
N = 20
K = 2
W = np.random.rand(N, N, K)
A = W.copy()
for k in range(K):
A[:, :, k] = np.linalg.pinv(W[:, :, k])
print "This should come out to zero \n", isi.joint_ISI(W, A)
A = np.random.rand(N, N, K)
print "Who knows what this should come out to \n", isi.joint_ISI(W, A)
def test1_2 (fil1, fil2, number) :
K=10
variables = loadmat(fil1)
A = np.zeros((16,16,K))
for k in range(K) :
A[:,:,k] = variables['A'][0,k]
variables = loadmat(fil2)
X = np.zeros((16,32968,K))
for k in range(K) :
X[:,:,k] = np.dot(A[:,:,k], variables['Sgt'][0,k][:,:])
fil = open("test_SCV_r0%i_FA.csv"%number, "w")
fil.write("Iteration,ISI\n")
fil2 = open("temp_ISI.txt", "w")
X_values = []
ISI_values = []
for k in range(2, K) :
W,_,_ = iva_l(X[:,:,0:k], verbose=True)
isi = joint_ISI(W,A[:,:,0:k])
print isi
fil.write("%i,%f\n"%(k,isi))
fil2.write("%i"%isi)
X_values.append(k)
ISI_values.append(isi)
fil.close()
fil2.close()
plt.plot(X_values, ISI_values)
plt.title("ISI Value by Number of total Subjects")
plt.xlabel("Number of Subjects")
plt.ylabel("ISI Value")
plt.savefig("normal_IVA_r0%i.pdf"%number)
plt.close("all")
