changed ldstore LD to 32bit

RobinM-code · Feb 19, 2021 · f11cbf2 · f11cbf2
1 parent b8613f3
commit f11cbf2
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/ldstore/bcor.py b/ldstore/bcor.py
@@ -72,15 +72,15 @@ def readCorr( self, snps ):
 
         if len( snps ) == 0:
             self.__fh.seek( self.__corr_block_offset )
-            corr = np.identity( self.__nSNPs )
+            corr = np.identity( self.__nSNPs, dtype=np.float32)
             for snp_x in range( self.__nSNPs - 1 ):
                 for snp_y in range( snp_x + 1, self.__nSNPs ):
                     corr[ snp_y, snp_x ] = self.__readCorrPair( snp_x, snp_y, False )
             return np.tril( corr ) + np.triu( corr.T, 1 )
         else:
             areSNPsIncluded( snps, self.__nSNPs, self.__fname )
 
-            corr = np.zeros( [ self.__nSNPs, len( snps ) ] )
+            corr = np.zeros( [ self.__nSNPs, len( snps ) ], dtype=np.float32)
             for snp_x, index in zip( snps, list(range( len( snps ))) ):
                 for snp_y in range( self.__nSNPs ):
                     if snp_x == snp_y: