sgn reader

from ctypes import *
import mmap

class SgnHeader(Structure):
    _pack_ = 1
    _fields_ = [ ("ID", c_char * 64),
                 ("headLen", c_ulong),
                 ("RemarkLen", c_ulong),
                 ("usPackingType", c_ushort),
                 ("Squared", c_ushort),
                 ("SampleSize", c_ushort),
                 ("FreqDiskr", c_double),
                 ("dBand", c_double),
                 ("dCF", c_double),
                 ("tSamplesOrigin", c_ulonglong) ]

class Ipp16sc(Structure):
    _pack_ = 1
    _fields_ = [ ("re", c_short), ("im", c_short) ]

headerLen = sizeof( SgnHeader )
sampleLen = sizeof( Ipp16sc )

def doit():
    global headerLen
    global sampleLen

    f = open( "d:/InComing/Analiz/file/0_tone.sgn", "r+b" )
    m = mmap.mmap( f.fileno(), 0 )
    header = SgnHeader.from_buffer_copy( m[:headerLen] )
    samplesCount = ( len( m ) - headerLen ) / sampleLen
    s = Ipp16sc * samplesCount
    samples = s.from_buffer_copy( m[headerLen:] )
    for i in samples[:8192]:
        i.re * i.re + i.im * i.im
    m.close()
    f.close()

from timeit import Timer
t = Timer( "doit()", "from __main__ import doit" )
print t.timeit( number = 10 )