Skip to content

Commit 152f208

Browse files
committed
committor problems, dependency versions in setup.py
1 parent b7cd47d commit 152f208

File tree

5 files changed

+122
-118
lines changed

5 files changed

+122
-118
lines changed

csnanalysis/csn.py

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -238,13 +238,13 @@ def calc_committors(self,basins,labels=None,basin_labels=None,add_basins=False,t
238238
else:
239239
# use trimmed transition matrix
240240
comm = committor(self.transmat,basins,tol=tol,maxstep=maxstep)
241-
full_comm = np.zeros(self.nnodes,dtype=float64)
241+
full_comm = np.zeros((self.transmat.shape[0],len(basins)),dtype=float)
242242
for i,ind in enumerate(self.trim_indices):
243243
full_comm[ind] = comm[i]
244244

245245
if labels is None:
246246
labels = ['p' + str(i) for i in range(len(basins))]
247-
for i,b in enumerate(basins):
247+
for i in range(len(basins)):
248248
self.add_attr(labels[i], full_comm[:,i])
249249

250250
if add_basins:

csnanalysis/matrix.py

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -137,7 +137,8 @@ def committor(transmat,basins,tol=1e-6,maxstep=20):
137137
"""
138138

139139
# make sink_matrix
140-
sink_mat = make_sink(transmat,list(np.array(basins).flatten()))
140+
flat_basin_list = [state for b in basins for state in b]
141+
sink_mat = make_sink(transmat,flat_basin_list)
141142

142143
sink_results = trans_mult_iter(sink_mat,tol,maxstep)
143144

0 commit comments

Comments
 (0)