From f98a45da5bd22be0da0d18f98eb594bc39bb706b Mon Sep 17 00:00:00 2001
From: Ganten-Hornby <65263129+Ganten-Hornby@users.noreply.github.com>
Date: Tue, 3 Dec 2024 15:50:23 +0800
Subject: [PATCH] optimize mem usage in latent2gene

---
 src/gsMap/latent_to_gene.py | 22 ++++++++++++----------
 1 file changed, 12 insertions(+), 10 deletions(-)

diff --git a/src/gsMap/latent_to_gene.py b/src/gsMap/latent_to_gene.py
index 07239a5..43b9ef6 100644
--- a/src/gsMap/latent_to_gene.py
+++ b/src/gsMap/latent_to_gene.py
@@ -9,7 +9,7 @@
 from scipy.stats import rankdata
 from sklearn.metrics.pairwise import cosine_similarity
 from sklearn.neighbors import NearestNeighbors
-from tqdm import tqdm
+from tqdm import tqdm, trange
 
 from gsMap.config import LatentToGeneConfig
 
@@ -153,10 +153,6 @@ def run_latent_to_gene(config: LatentToGeneConfig):
         adata.var_names = homologs.loc[adata.var_names, 'HUMAN_GENE_SYM'].values
         adata = adata[:, ~adata.var_names.duplicated()]
 
-    # Create mappings
-    n_cells = adata.n_obs
-    n_genes = adata.n_vars
-
     if config.annotation is not None:
         cell_annotations = adata.obs[config.annotation].values
         logger.info(f'Using cell annotations for {len(cell_annotations)} cells.')
@@ -204,14 +200,18 @@ def run_latent_to_gene(config: LatentToGeneConfig):
     else:
         adata_X = adata.X.tocsr()
 
-    ranks = np.zeros((n_cells, adata.n_vars), dtype=np.float32)
+    # Create mappings
+    n_cells = adata.n_obs
+    n_genes = adata.n_vars
 
+    ranks = np.zeros((n_cells, adata.n_vars), dtype=np.float16)
     for i in tqdm(range(n_cells), desc="Computing ranks per cell"):
         data = adata_X[i, :].toarray().flatten()
         ranks[i, :] = rankdata(data, method='average')
 
     if gM is None:
         gM = gmean(ranks, axis=0)
+        gM = gM.astype(np.float16)
 
     adata_X_bool = adata_X.astype(bool)
     if frac_whole is None:
@@ -225,15 +225,17 @@ def run_latent_to_gene(config: LatentToGeneConfig):
     # Normalize the ranks
     ranks /= gM
 
-    # Compute marker scores in parallel
-    logger.info('------Computing marker scores...')
     def compute_mk_score_wrapper(cell_pos):
         return compute_regional_mkscore(
             cell_pos, spatial_net_dict, coor_latent, config, cell_annotations, ranks, frac_whole, adata_X_bool
         )
 
-    mk_scores = [compute_mk_score_wrapper(cell_pos) for cell_pos in tqdm(range(n_cells), desc="Calculating marker scores")]
-    mk_score = np.vstack(mk_scores).T
+    logger.info('------Computing marker scores...')
+    mk_score = np.zeros((n_cells, n_genes), dtype=np.float16)
+    for cell_pos in trange(n_cells, desc="Calculating marker scores"):
+        mk_score[cell_pos, :] = compute_mk_score_wrapper(cell_pos)
+
+    mk_score = mk_score.T
     logger.info('Marker scores computed.')
 
     # Remove mitochondrial genes