text: import os
  import numpy as np 
  
  datafolder="../QMCPACK/288/288x115x69x69/split" 
  ebs=[i*1e-4 for i in range(1,10)]+[i*1e-3 for i in range(1,10)]+[i*1e-2 for i in range(1,11)]
  
  cr=np.zeros((29,289),dtype=np.float32)
  psnr=np.zeros((29,289),dtype=np.float32)
  ctime=np.zeros((29,289),dtype=np.float32)
  dtime=np.zeros((29,289),dtype=np.float32)
  idxrange=range(288)
  for i,eb in enumerate(ebs):
      cr[i+1][0]=eb
      psnr[i+1][0]=eb
      ctime[i+1][0]=eb
      dtime[i+1][0]=eb
      for j,idx in enumerate(idxrange):
          filename="%d.dat" % idx
          cr[0][j+1]=idx
          psnr[0][j+1]=idx
          ctime[0][j+1]=idx
          dtime[0][j+1]=idx
          filepath=os.path.join(datafolder,filename)
          comm="sz_interp %s -3 115 69 69 %f>tempqmcpackinterp.txt" % (filepath,eb)
          os.system(comm)
          
          with open("tempqmcpackinterp.txt","r") as f:
              lines=f.read().splitlines()
              print(lines)
              p=eval(lines[-3].split(',')[0].split('=')[-1])
              
              r=eval(lines[-3].split(',')[-1].split('=')[-1])
              ct=eval(lines[-2].split('=')[-1])
              dt=eval(lines[-1].split('=')[-1])
              
              
          cr[i+1][j+1]=r
          psnr[i+1][j+1]=p
          ctime[i+1][j+1]=ct
          dtime[i+1][j+1]=dt
          os.system("rm -f tempqmcpackinterp.txt")
  
  
  np.savetxt("interp_qmcpack288_cr.txt",cr,delimiter='\t')
  np.savetxt("interp_qmcpack288_psnr.txt",psnr,delimiter='\t')
  np.savetxt("interp_qmcpack288_ctime.txt",ctime,delimiter='\t')
  np.savetxt("interp_qmcpack288_dtime.txt",dtime,delimiter='\t')
  