schain Commit - r1272:60a078275b06 · Jicamarca Repository

933

rangeList : List of profile ranges. Example: rangeList = ((4, 30), (32, 64), (128, 256))

933

rangeList : List of profile ranges. Example: rangeList = ((4, 30), (32, 64), (128, 256))

934

935

"""

935

"""

936

print "HOLA MUNDO CRUEL"

937

if rangeList is not None:

937

if rangeList is not None:

938

if type(rangeList[0]) not in (tuple, list):

938

if type(rangeList[0]) not in (tuple, list):

939

rangeList = [rangeList]

939

rangeList = [rangeList]

1357

def setup(self,dataOut ,lags = None,mode =None, fullBuffer= None ,pairsList = None,nAvg = 1):

1357

def setup(self,dataOut ,lags = None,mode =None, fullBuffer= None ,pairsList = None,nAvg = 1):

1358

1359

self.lags = lags

1359

self.lags = lags

1360

print self.lags

1360

self.mode = mode

1361

self.mode = mode

1361

self.fullBuffer= fullBuffer

1362

self.fullBuffer= fullBuffer

1362

self.nAvg = nAvg

1363

self.nAvg = nAvg

1364

self.pairsList = pairsList

1363

nChannels = dataOut.nChannels

1365

nChannels = dataOut.nChannels

1364

nProfiles = dataOut.nProfiles

1366

nProfiles = dataOut.nProfiles

1365

nHeights = dataOut.nHeights

1367

nHeights = dataOut.nHeights

1374

if pairsList == None:

1376

if pairsList == None:

1375

print "Pairs list selected by default (1,0)"

1377

print "Pairs list selected by default (1,0)"

1376

pairsList = [(0,1)]

1378

pairsList = [(0,1)]

1377

~~pair~~List= pairsList

1379

self.pairsList= pairsList

1378

1380

1379

if lags == None:

1381

if lags == None:

1380

if mode=='time':

1382

if mode=='time':

1390

elif mode =='height':

1392

elif mode =='height':

1391

self.tmp = numpy.zeros(len(pairsList), (len(self.lags), nProfiles),dtype='complex')

1393

self.tmp = numpy.zeros(len(pairsList), (len(self.lags), nProfiles),dtype='complex')

1392

1394

1393

print "lags", len(self.lags)

1395

print "lags", len(self.lags),self.lags

1394

print "mode",self.mode

1396

print "mode",self.mode

1395

print "nChannels", nChannels

1397

print "nChannels", nChannels

1396

print "nProfiles", nProfiles

1398

print "nProfiles", nProfiles

1397

print "nHeights" , nHeights

1399

print "nHeights" , nHeights

1398

print "pairsList", pairsList

1400

print "pairsList", pairsList,len(self.pairsList)

1399

print "fullBuffer", fullBuffer

1401

print "fullBuffer", fullBuffer

1400

#print "type(pairsList)",type(pairsList)

1402

#print "type(pairsList)",type(pairsList)

1401

print "tmp.shape",self.tmp.shape

1403

print "tmp.shape",self.tmp.shape

1408

if not self.isConfig:

1410

if not self.isConfig:

1409

self.setup(dataOut, lags = lags,mode = mode, fullBuffer= fullBuffer ,pairsList = pairsList,nAvg=nAvg)

1411

self.setup(dataOut, lags = lags,mode = mode, fullBuffer= fullBuffer ,pairsList = pairsList,nAvg=nAvg)

1410

self.isConfig = True

1412

self.isConfig = True

1411

1412

if dataOut.type == "Voltage":

1413

if dataOut.type == "Voltage":

1413

1414

data_pre = dataOut.data #data

1414

data_pre = dataOut.data #data

1415

print "size",data_pre.shape

1416

1417

# Here is the loop :D

1416

# Here is the loop :D

1418

for l in range(len(pairsList)):

1417

for l in range(len(self.pairsList)):

1419

ch0 = pairsList[l][0]

1418

ch0 = self.pairsList[l][0]

1420

ch1 = pairsList[l][1]

1419

ch1 = self.pairsList[l][1]

1421

1422

for i in range(len(self.lags)):

1420

for i in range(len(self.lags)):

1421

print i,"hola dios"

1423

idx = self.lags[i]

1422

idx = self.lags[i]

1424

if self.mode == 'time':

1423

if self.mode == 'time':

1425

acf0 = data_pre[ch0,:self.__nProfiles-idx,:]*numpy.conj(data_pre[ch1,idx:,:]) # pair,lag,height

1424

acf0 = data_pre[ch0,:self.__nProfiles-idx,:]*numpy.conj(data_pre[ch1,idx:,:]) # pair,lag,height

1430

self.tmp[l,i,:acf0.shape[0],:]= acf0

1429

self.tmp[l,i,:acf0.shape[0],:]= acf0

1431

else:

1430

else:

1432

self.tmp[l,i,:]= numpy.sum(acf0,axis=0)

1431

self.tmp[l,i,:]= numpy.sum(acf0,axis=0)

1433

1434

if self.fullBuffer:

1432

if self.fullBuffer:

1435

1436

self.tmp = numpy.sum(numpy.reshape(self.tmp,(self.tmp.shape[0],self.tmp.shape[1],self.tmp.shape[2]/self.nAvg,self.nAvg,self.tmp.shape[3])),axis=3)

1433

self.tmp = numpy.sum(numpy.reshape(self.tmp,(self.tmp.shape[0],self.tmp.shape[1],self.tmp.shape[2]/self.nAvg,self.nAvg,self.tmp.shape[3])),axis=3)

1437

dataOut.nAvg = self.nAvg

1434

dataOut.nAvg = self.nAvg

1438

1435

1447

for j in range(shape[2]):

1444

for j in range(shape[2]):

1448

self.tmp[i,:,j]= self.tmp[i,:,j].real / numpy.max(numpy.abs(self.tmp[i,:,j]))

1445

self.tmp[i,:,j]= self.tmp[i,:,j].real / numpy.max(numpy.abs(self.tmp[i,:,j]))

1449

1446

1450

1451

#import matplotlib.pyplot as plt

1447

#import matplotlib.pyplot as plt

1452

#print "test",self.tmp.shape

1448

#print "test",self.tmp.shape

1453

#print self.tmp[0,0,0]

1449

#print self.tmp[0,0,0]

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages

                         rangeList            :    List of profile ranges. Example: rangeList = ((4, 30), (32, 64), (128, 256))
                     """
+                    print "HOLA MUNDO CRUEL"
                     if rangeList is not None:
                         if type(rangeList[0]) not in (tuple, list):
                             rangeList = [rangeList]
                 def setup(self,dataOut ,lags = None,mode =None, fullBuffer= None ,pairsList = None,nAvg = 1):
                 	self.lags      = lags
+                    print self.lags
                 	self.mode      = mode
                 	self.fullBuffer= fullBuffer
                 	self.nAvg      = nAvg
+                    self.pairsList  = pairsList
                 	nChannels      = dataOut.nChannels
                 	nProfiles      = dataOut.nProfiles
                 	nHeights       = dataOut.nHeights
                 	if pairsList == None:
                 	   print "Pairs list selected by default (1,0)"
                 	   pairsList = [(0,1)]
-                	     pairList= pairsList
+                	   self.pairsList= pairsList
                 	if lags == None:
                         if mode=='time':
                     elif mode =='height':
                 		self.tmp = numpy.zeros(len(pairsList), (len(self.lags), nProfiles),dtype='complex')
-                	print "lags", len(self.lags)
+                	print "lags", len(self.lags),self.lags
                 	print "mode",self.mode
                 	print "nChannels", nChannels
                 	print "nProfiles", nProfiles
                 	print "nHeights" , nHeights
-                	print "pairsList", pairsList
+                	print "pairsList", pairsList,len(self.pairsList)
                 	print "fullBuffer", fullBuffer
                 	#print "type(pairsList)",type(pairsList)
                 	print "tmp.shape",self.tmp.shape
                     if not self.isConfig:
                         self.setup(dataOut, lags = lags,mode = mode, fullBuffer= fullBuffer ,pairsList = pairsList,nAvg=nAvg)
                         self.isConfig = True
                     if dataOut.type == "Voltage":
                         data_pre = dataOut.data  #data
+                        print "size",data_pre.shape
                         #  Here is the loop :D
-                	    for l in range(len(pairsList)):
+                	    for l in range(len(self.pairsList)):
-                		ch0 = pairsList[l][0]
+                            ch0 = self.pairsList[l][0]
-                		ch1 = pairsList[l][1]
+                            ch1 = self.pairsList[l][1]
                             for i in range(len(self.lags)):
+                                print i,"hola dios"
                                 idx = self.lags[i]
                                 if self.mode == 'time':
                                     acf0 = data_pre[ch0,:self.__nProfiles-idx,:]*numpy.conj(data_pre[ch1,idx:,:])    # pair,lag,height
                                     self.tmp[l,i,:acf0.shape[0],:]= acf0
                                 else:
                                     self.tmp[l,i,:]= numpy.sum(acf0,axis=0)
                         if self.fullBuffer:
                             self.tmp = numpy.sum(numpy.reshape(self.tmp,(self.tmp.shape[0],self.tmp.shape[1],self.tmp.shape[2]/self.nAvg,self.nAvg,self.tmp.shape[3])),axis=3)
                             dataOut.nAvg = self.nAvg
                             for j in range(shape[2]):
                                 self.tmp[i,:,j]= self.tmp[i,:,j].real / numpy.max(numpy.abs(self.tmp[i,:,j]))
             	        #import matplotlib.pyplot as plt
             	        #print "test",self.tmp.shape
             	        #print self.tmp[0,0,0]