from conductivity import getBoundary
from solveBC import sweep, findrho
from printStatus import printRatio
import numpy as np
import pylab as pl
from fig import fig, saveFig, printText
sigma0s=[]
taus=[]
ms=[]
a2s=np.logspace(-2,5,40)
for a2i in range(len(a2s)):
printRatio(a2i,len(a2s))
a2=a2s[a2i]
Tr=2 #use T = 2 T_c
bbs,sols=sweep(lambda x,y:getBoundary(x,y,[0,a2]),[1e-6],oscN=1)
rho=-np.array(bbs)[:,:,1,1].real
rhoc=min([min(r) for r in rho])#rho, in units used at Tc
rho=rhoc/Tr**2
zh=1
T=3./(zh*4.*np.pi)
Tc=T*np.sqrt(rho/rhoc)
bb,osc=getBoundary(1,0,[0,a2])
assert osc==0
hphi=rho/(-bb[1][1])
ws=[1e-3]#np.logspace(-3,3,200)
sigmas=[]
from conductivity import getBoundary
from solveBC import sweep, findrho
from printStatus import printRatio
import numpy as np
import pylab as pl
from fig import fig, saveFig, printText
a2=1e0
bbs,sols=sweep(lambda x,y:getBoundary(x,y,[0,a2]),[1e-6],oscN=1)
rho=-np.array(bbs)[:,:,1,1].real
rhoc=min([min(r) for r in rho])#rho, in units used at Tc
Tr=2 #use T = 2 T_c
rho=rhoc/Tr**2
zh=1
T=3./(zh*4.*np.pi)
Tc=T*np.sqrt(rho/rhoc)
bb,osc=getBoundary(1,0,[0,a2])
guess=rho/(-bb[1][1])
from scipy.optimize import fsolve
hphi=fsolve(lambda hphi:rho+getBoundary(hphi,0,[0,a2])[0][1][1], guess)
assert osc==0
ws=np.logspace(-3,3,200)
sigmas=[]
for i in range(len(ws)):
printRatio(i,len(ws))