fixed bug that division was forced as numpy even though x could be cs…

…r_matrix using np.divide instead of / made the code very slow (seemingly arresting) here.
theislab · Aug 29, 2019 · be9234f · be9234f
1 parent ddcd0f0
commit be9234f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/batchglm/models/base_glm/utils.py b/batchglm/models/base_glm/utils.py
@@ -137,7 +137,7 @@ def closedform_glm_scale(
     :return: tuple (groupwise_scales, logphi, rmsd)
     """
     if size_factors is not None:
-        x = np.divide(x, size_factors)
+        x = x / size_factors
 
     # to circumvent nonlocal error
     provided_groupwise_means = groupwise_means