double precision

python/spherelib/fitsfunc2.py

@@ -133,7 +133,7 @@ def write_covmat (filename, covmat, nbmodes=None):
     c = 0
     for i in range(ndet):
         for j in range(i+1):
-            cols.append(pyfits.Column(name='cross %d' % c, format='E', array=array(covmat[i,j,:])))
+            cols.append(pyfits.Column(name='cross %d' % c, format='D', array=array(covmat[i,j,:])))
             c = c+1
             
     cols  = pyfits.ColDefs(cols)
@@ -186,7 +186,7 @@ def write_table (filename, mat, colname=None):
 
     cols = []
     for c in range(ncols):
-        cols.append(pyfits.Column(name=colname[c], format='E', array=array(mat[:,c])))
+        cols.append(pyfits.Column(name=colname[c], format='D', array=array(mat[:,c])))
 
     cols  = pyfits.ColDefs(cols)
     tbhdu = pyfits.new_table(cols)

python/spherelib/myIO.py

@@ -161,7 +161,7 @@ def write_map (file, map):
         except pio.pioError,val: 
             return str(val)
     else:
-        hp.write_map(file, map)
+        hp.write_map(file, map,dtype=np.float64 )
 
 def read_alm (file):
     """ Read Healpix alm from file