Ejemplo n.º 1
0
 def set_const(self, consts):
     
     if self.consts != None:
         self.consts[:len(consts)]   =   consts        
     else:
         self.consts                 =   consts
     
     
     if self.system in ['spiral_w_wave_sqr', 'spiral_w_wave_str'] \
                  and self.rps == None:
         self.initialize_prime_maps()
     
     if self.system in ['spiral_w_wave_sqr', 'spiral_w_wave_str']:
         
         self.u()
         self.du()
         
         rdev, phidev            =   [1,1] 
         
         while not phidev < 1e-10 or not rdev < 1e-10: #1e-6 < phidev or 1e-6 < rdev:
             
             self.curve_Lphi     =   curve_Lphi(self, test = False)
             
             self.phips          =   phips(self)
             
             self.u()
             self.du()
             
             rdev                =   test_prime_maps(self)[0]
             
             self.curve_Lr       =   curve_Lr(self)
             self.rps            =   rps(self) 
             
             self.u()
             self.du()
             
             phidev              =   test_prime_maps(self)[1]
             
             #print rdev, phidev
     
     self.u()  
     self.du()
     self.ddu()        
     self.dA()
     
     '''
Ejemplo n.º 2
0
 def initialize_prime_maps(self):
     
     store_opt               =   self.opt
         
     self.opt                =   [1, 0, self.opt[2]] 
     n                       =   0
     rdev, phidev            =   1, 1
     while not phidev < 1e-10 or not rdev < 1e-10: #1e-6 < phidev or 1e-6 < rdev:
         
         self.u()
         self.du()
         
         if n > 0:
             phidev          =   test_prime_maps(self)[1]
             #path_integral_r(self)
             
         self.curve_Lphi     =   curve_Lphi(self, test = False)
         self.phips          =   phips(self)
         #plot_surf(self)
        
         if n == 0:
             self.opt        =   [1, 3, self.opt[2]] 
         else:
             rdev    =   test_prime_maps(self)[0]
             #print phidev, rdev
             
        
         self.u()
         self.du()
         
         self.curve_Lr       =   curve_Lr(self)
         self.rps            =   rps(self) 
         
         
         n                  +=   1
         self.opt            =   store_opt