From 43fa9dfdaf8b4a254d1dd59ccb0c759b3600f724 Mon Sep 17 00:00:00 2001 From: bbauvin <baptiste.bauvin@centrale-marseille.fr> Date: Thu, 15 Sep 2016 09:52:17 -0400 Subject: [PATCH] Trying to fix segmentation fault --- Code/MonoMutliViewClassifiers/Multiview/GetMultiviewDb.py | 1 + Code/MonoMutliViewClassifiers/utils/Dataset.py | 6 ++++++ 2 files changed, 7 insertions(+) diff --git a/Code/MonoMutliViewClassifiers/Multiview/GetMultiviewDb.py b/Code/MonoMutliViewClassifiers/Multiview/GetMultiviewDb.py index 7c5a3386..2164b88a 100644 --- a/Code/MonoMutliViewClassifiers/Multiview/GetMultiviewDb.py +++ b/Code/MonoMutliViewClassifiers/Multiview/GetMultiviewDb.py @@ -430,6 +430,7 @@ def getAdjacenceMatrix(sortedRNASeq, k=1): # indices.append(patient[i-1]+patient[i]*nbGenes) # data.append(True) indptr[patientIndex+1] = pointer + mat = sparse.csr_matrix((data, indices, indptr), shape=(sortedRNASeq.shape[0], sortedRNASeq.shape[1]*sortedRNASeq.shape[1]), dtype=bool) return mat diff --git a/Code/MonoMutliViewClassifiers/utils/Dataset.py b/Code/MonoMutliViewClassifiers/utils/Dataset.py index 93db8a31..67b660a0 100644 --- a/Code/MonoMutliViewClassifiers/utils/Dataset.py +++ b/Code/MonoMutliViewClassifiers/utils/Dataset.py @@ -8,6 +8,12 @@ def getV(DATASET, viewIndex, usedIndices=None): if not DATASET.get("View"+str(viewIndex)).attrs["sparse"]: return DATASET.get("View"+str(viewIndex))[usedIndices, :] else: + sparse_mat = sparse.csr_matrix((DATASET.get("View"+str(viewIndex)).get("data").value, + DATASET.get("View"+str(viewIndex)).get("indices").value, + DATASET.get("View"+str(viewIndex)).get("indptr").value), + shape=DATASET.get("View"+str(viewIndex)).attrs["shape"])[usedIndices,:] + print sparse_mat.shape + print sparse_mat.indptr return sparse.csr_matrix((DATASET.get("View"+str(viewIndex)).get("data").value, DATASET.get("View"+str(viewIndex)).get("indices").value, DATASET.get("View"+str(viewIndex)).get("indptr").value), -- GitLab