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()
'''
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