print "Settings for the Sun" print "M = %10.6g g" % Ms print "X = %10.3f" % X print "Y = %10.3f" % Y print u"ε = %10.3f ergs/g" % epsilon print u"μ = %10.3f" % mmw(X=X, Y=Y) print "Outer Boundary Conditions" print "Testing Outer Presure Conversion Mehtods" print " Inital Guess:" print " R = %10.6e cm" % Rs print " L = %10.6e erg/s" % Lsun start = time.clock() ri, li, Pi, Ti = outer_boundary(R=Rs, L=Lsun, M=Ms, mu=mmw(X=X, Y=Y), optable=Opacity, Piter=10) print "Iterable time taken: %g" % (time.clock() - start) start = time.clock() r, l, P, T = outer_boundary(R=Rs, L=Lsun, M=Ms, mu=mmw(X=X, Y=Y), optable=Opacity, Piter=False) print "Absolute time taken: %g" % (time.clock() - start) print "At m=M:" print " -- Algebraic -------------- | -- Iterable -------------- | -- Delta --------| " print u" R = %10.6e cm, | R = %10.6e cm | ΔR = %11.4e |" % ( r, ri, (r - ri) / r)
print "Problem #4&5:" print "Settings for the Sun" print "M = %10.6g g" % Ms print "X = %10.3f" % star.X print "Y = %10.3f" % star.Y print u"μ = %10.3f" % mmw(X=X, Y=Y) print "" print "Outer Boundary Conditions" print " Inital Guess:" print " R = %10.6e cm" % Rs print " L = %10.6e erg/s" % Lsun r, l, P, T = outer_boundary(R=Rs, L=Lsun, M=Ms, mu=mmw(X=X, Y=Y), optable=Opacity, Piter=P_Guess_Iteration) print "At m=M:" print " r = %10.6e cm" % r print " l = %10.6e erg/s" % l print " P = %10.6e Dyne/cm^2" % P print " T = %10.6e K" % T # print "Derivatives" # dr,dl,dP,dT = derivatives(np.array([Ms]),np.atleast_2d([r,l,P,T]),mu=mmw(X=X,Y=Y),optable=Opacity) # print " dr = %10.6e" % dr # print " dl = %10.6e" % dl # print " dP = %10.6e" % dP # print " dT = %10.6e" % dT # print ""
print "Stellar Model Photospheric Boundary Conditions" print "Settings for the Sun" print "M = %10.6g g" % Ms print "X = %10.3f" % X print "Y = %10.3f" % Y print u"ε = %10.3f ergs/g" % epsilon print u"μ = %10.3f" % mmw(X=X,Y=Y) print "Outer Boundary Conditions" print "Testing Outer Presure Conversion Mehtods" print " Inital Guess:" print " R = %10.6e cm" % Rs print " L = %10.6e erg/s" % Lsun start = time.clock() ri,li,Pi,Ti = outer_boundary(R=Rs,L=Lsun,M=Ms,mu=mmw(X=X,Y=Y),optable=Opacity,Piter=10) print "Iterable time taken: %g" % (time.clock() - start) start = time.clock() r,l,P,T = outer_boundary(R=Rs,L=Lsun,M=Ms,mu=mmw(X=X,Y=Y),optable=Opacity,Piter=False) print "Absolute time taken: %g" % (time.clock() - start) print "At m=M:" print " -- Algebraic -------------- | -- Iterable -------------- | -- Delta --------| " print u" R = %10.6e cm, | R = %10.6e cm | ΔR = %11.4e |" % (r,ri,(r-ri)/r) print u" L = %10.6e erg/s, | L = %10.6e erg/s | ΔL = %11.4e |" % (l,li,(l-li)/l) print u" P = %10.6e Dyne/cm^2, | P = %10.6e Dyne/cm^2 | ΔP = %11.4e |" % (P,Pi,(P-Pi)/P) print u" T = %10.6e K, | T = %10.6e K | ΔT = %11.4e |" % (T,Ti,(T-Ti)/T) Pc = 2.477e17 Tc = 1.571e7 m = 1e-30*Ms