notebook.community

Edit and run

Imports and function definitions





In [1]:



    


%pylab --no-import-all inline
#%pylab --no-import-all osx
%config InlineBackend.figure_format='png' # SVG plots can get really slow if the evolutions are long



    


















    






Populating the interactive namespace from numpy and matplotlib

















In [2]:



    


import os, os.path
import h5py
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import PNEvolution
import Quaternions
from Quaternions import Quaternion, sqrtOfRotor, xHat, yHat, zHat



    











In [3]:



    


def vOfOmega(Omega):
    return Omega**(1./3.)
def OmegaOfv(v):
    return v**(3)
def norm(v) :
    return sqrt(v[0]**2+v[1]**2+v[2]**2)



    











In [4]:



    


def Frame_from_nHat_ellhat(nHat, ellHat):
    """Deduce rotor taking (x,z) onto (nHat,ellHat)
    
    Given 3-vectors `nHat` and `ellHat`, corresponding to
    the x' and z' axes of a new frame, this function
    returns the appropriate rotor taking the x and z axes
    of the original frame onto those vectors.
    
    """
    R_frame_1 = sqrtOfRotor(-Quaternion(0.,*ellHat)*zHat)
    return sqrtOfRotor(- Quaternion(0.,*nHat) * (R_frame_1*xHat*R_frame_1.conjugate())) * R_frame_1

We use cython to define a fast derivative function





In [5]:



    


%load_ext cythonmagic



    











In [6]:



    


%%cython
cimport cython
cimport numpy as np
import numpy as np
from libc.math cimport exp, sqrt, pow, log, erf

def Derivative(np.ndarray y, np.ndarray x):
    """Differentiate first argument with respect to second.
    
    This is a second-order differentiation scheme that can
    handle unequally spaced time steps.  The two arguments
    should be numpy arrays with equal first dimensions; the
    first argument can have any number of additional
    dimensions.
    
    Reference: Sundqvist and Veronis, Tellus XXII (1970), 1
    
    """
    def SQR(double a):
        return a*a
    if(x.ndim!=1 or x.shape[0]!=y.shape[0]) :
        raise TypeError("Bad shapes: x({0},{1}); y({2})".format(x.shape[0], x.ndim, [y.shape[j] for j in range(y.ndim)]))
    cdef np.ndarray z = np.empty_like(y)
    cdef int i1 = z.shape[0]-1
    cdef double hprev = x[1]-x[0]
    cdef double hnext = x[2]-x[1]
    z[0] = -((2*hprev+hnext)/(hprev*(hprev+hnext)))*y[0] + ((hnext+hprev)/(hnext*hprev))*y[1] - (hprev/(hnext*(hnext+hprev)))*y[2]
    cdef int i=1
    for i in range(1,i1):
        hnext = x[i+1]-x[i]
        z[i] = (y[i+1] - y[i-1]*SQR(hnext/hprev) - y[i]*(1-SQR(hnext/hprev))) / (hnext*(1+hnext/hprev));
        hprev = hnext
    hnext = x[i1]  -x[i1-1]
    hprev = x[i1-1]-x[i1-2]
    z[i1] = (hnext/(hprev*(hprev+hnext)))*y[i1-2] - ((hnext+hprev)/(hnext*hprev))*y[i1-1] \
        + ((hprev+2*hnext)/(hnext*(hnext+hprev)))*y[i1]
    return z

Basic parameters





In [7]:



    


StartFromTEqualsZero = True
PNOrders = [0.0, 2.0, 3.5, 6.0]
Approximants = ["T1", "T4", "T5"]

AnnexDir = '/Users/boyle/Research/Data/SimulationAnnex/Catalog/'
if not os.path.isdir(AnnexDir):
    AnnexDir = '/mnt/raid-project/nr/boyle/SimulationAnnex/Catalog/'
if not os.path.isdir(AnnexDir):
    raise IOError("Directory '{0}' does not exist".format(AnnexDir))

NRDirs = ['SKS/d10.0_q6.00_sA_0.739_0.190_-0.495_sB_-0.190_0.051_-0.227/Lev3/',
          'q1.0/SBBH/d19.0_q1.0_s0.5_0_0_s0_0_0/Lev5/', # Serguei's favorite
          'q5.0/SBBH/d15.0_q5.0_s0.5_0_0_s0/Lev5/',
          'q5.0/SBBH/d15.0_q5.0_s0_0_0.5_s0/Lev4/']
InitialOmegas = [
    [1.421005247651657167e-10,  5.015494672470230332e-05, 2.903442067573783203e-02],
    [5.434176633374474299e-16, -6.815194916799825440e-06, 1.128607156589968470e-02],
    [7.764246764733414600e-16, -2.720085057750848984e-07, 1.580435785321985045e-02],
    [-5.133688651140011539e-15, 4.204483349162996808e-13, 1.571170967415409445e-02]]

NRDir = NRDirs[0]

PN evolutions and comparisons to numerical data





In [11]:



    


Header = """
# PNEvolution.EvolvePN_Q(Approximant, PNOrder, v_0, v_i, m1, m2, chi1_i, chi2_i, R_frame_i, L)
[1] = t/M
[2] = v
[3] = chi1_x
[4] = chi1_y
[5] = chi1_z
[6] = chi2_x
[7] = chi2_y
[8] = chi2_z
[9] = ellHat_x
[10] = ellHat_y
[11] = ellHat_z
[12] = Phi
[13] = nHat_x
[14] = nHat_y
[15] = nHat_z"""

def RunPNComparisons(PNOrders, Approximants, AnnexDir, NRDir, StartFromTEqualsZero, InitialOmega=None):
    MaxTFactor = 1.5
    
    # Make sure output directory exists
    if StartFromTEqualsZero:
        OutputDir = 'PNComparison/Initial/'
    else:
        OutputDir = 'PNComparison/Relaxed/'
    if not os.path.isdir(AnnexDir+NRDir+OutputDir):
        os.makedirs(AnnexDir+NRDir+OutputDir)
    
    # Read in data from Horizons.h5
    f = h5py.File(AnnexDir+NRDir+'Horizons.h5', 'r')
    t_NR = f['AhA.dir/chiInertial.dat'][:,0]
    if StartFromTEqualsZero:
        i_relaxed = 0
    else:
        t_relaxed = ! awk {'"/^relaxed-measurement-time/ {print \$3}"'} {AnnexDir+NRDir}/metadata.txt
        t_relaxed = float(t_relaxed[0])
        i_relaxed = abs(t_NR-t_relaxed).argmin()
    t_relaxed = t_NR[i_relaxed]
    t_NR = t_NR-t_NR[i_relaxed]
    chi1_NR = f['AhA.dir/chiInertial.dat'][:,1:]
    x1 = f['AhA.dir/CoordCenterInertial.dat'][:,1:]
    m1 = f['AhA.dir/ChristodoulouMass.dat'][0,1]
    m1_relaxed = f['AhA.dir/ChristodoulouMass.dat'][i_relaxed,1]
    chi2_NR = f['AhB.dir/chiInertial.dat'][:,1:]
    x2 = f['AhB.dir/CoordCenterInertial.dat'][:,1:]
    m2 = f['AhB.dir/ChristodoulouMass.dat'][0,1]
    m2_relaxed = f['AhB.dir/ChristodoulouMass.dat'][i_relaxed,1]
    f.close()
    
    # Derive orbital elements, etc.
    nHat_NR = array([(xa-xb)/norm(xa-xb) for xa,xb in zip(x1,x2)])
    nHatDot_NR = Derivative(nHat_NR, t_NR)
    ellVec_NR = array([cross(N,Ndot) for N,Ndot in zip(nHat_NR,nHatDot_NR)])
    ellHat_NR = array([o/norm(o) for o in ellVec_NR])
    v_NR = array([vOfOmega(norm(o)) for o in ellVec_NR])
    R_NR = array([Frame_from_nHat_ellhat(n, l) for n,l in zip(nHat_NR,ellHat_NR)])
    
    # Set PN parameters
    if(StartFromTEqualsZero):
        if(not InitialOmega):
            raise ValueError("Need InitialOmega input if the initial time is used")
        v0 = norm(InitialOmega)**(1./3.)
        v_i = v0
        x1_i = ! awk {'"/^initial-bh-position1|^initial-position1/ {print \$3 \$4 \$5}"'} {AnnexDir+NRDir}/metadata.txt
        x1_i = array([float(x) for x in x1_i[0].split(',')])
        x2_i = ! awk {'"/^initial-bh-position2|^initial-position2/ {print \$3 \$4 \$5}"'} {AnnexDir+NRDir}/metadata.txt
        x2_i = array([float(x) for x in x2_i[0].split(',')])
        nHat_i = x1_i-x2_i
        nHat_i = nHat_i / norm(nHat_i)
        R_frame_i = Frame_from_nHat_ellhat(nHat_i, array(InitialOmega)/norm(InitialOmega))
        chi1_i = ! awk {'"/^initial-bh-spin1|^initial-spin1/ {print \$3 \$4 \$5}"'} {AnnexDir+NRDir}/metadata.txt
        chi1_i = array([float(x) for x in chi1_i[0].split(',')])
        chi2_i = ! awk {'"/^initial-bh-spin2|^initial-spin2/ {print \$3 \$4 \$5}"'} {AnnexDir+NRDir}/metadata.txt
        chi2_i = array([float(x) for x in chi2_i[0].split(',')])
        m1 = ! awk {'"/^initial-bh-mass1|^initial-mass1/ {print \$3}"'} {AnnexDir+NRDir}/metadata.txt
        m1 = float(m1[0])
        m2 = ! awk {'"/^initial-bh-mass2|^initial-mass2/ {print \$3}"'} {AnnexDir+NRDir}/metadata.txt
        m2 = float(m2[0])
        chi1_i = chi1_i/m1**2
        chi2_i = chi2_i/m2**2
    else:
        v0 = v_NR[i_relaxed]
        v_i = v0
        chi1_i = chi1_NR[i_relaxed]
        chi2_i = chi2_NR[i_relaxed]
        R_frame_i = R_NR[i_relaxed]

    # Do the PN evolutions
    t={}
    v={}
    chi1={}
    chi2={}
    R_frame={}
    Phi={}
    ellHat={}
    nHat = {}
    L = {}
    for PNOrder in PNOrders:
        for Approximant in Approximants:
            Call = """# PNEvolution.EvolvePN_Q("Taylor{0}", {1}, {2}, {3}, {4}, {5}, {6}, {7}, {8})""".format(
                   Approximant, PNOrder, v0, v_i, m1, m2, str(chi1_i), str(chi2_i), str(R_frame_i))
            print Call
            Sym = "{1} {0}".format(Approximant,PNOrder)
            t[Sym],v[Sym],chi1[Sym],chi2[Sym],R_frame[Sym],Phi[Sym],L[Sym] = \
                PNEvolution.EvolvePN_Q("Taylor{0}".format(Approximant), PNOrder, v0, v_i, m1, m2, chi1_i, chi2_i, R_frame_i)
            ellHat[Sym] = PNEvolution.ellHat(R_frame[Sym])
            nHat[Sym] = PNEvolution.nHat(R_frame[Sym])
            t[Sym] = t[Sym]-t[Sym][0]
            savetxt(AnnexDir+NRDir+OutputDir+'PNOrbitalData_{0}_{1}.dat'.format(Approximant,PNOrder),
                    column_stack((t[Sym],v[Sym],chi1[Sym],chi2[Sym],ellHat[Sym],Phi[Sym],nHat[Sym])),
                    header=Call+Header)
    
    # Compare velocities
    fig = plt.figure()
    plot(t_NR, v_NR, label="NR data")
    for PNOrder in PNOrders:
        for Approximant in Approximants:
            Sym = "{1} {0}".format(Approximant,PNOrder)
            plot(t[Sym], v[Sym], label="Taylor{0} at {1} PN".format(Approximant,PNOrder))
    legend(loc='upper left')
    xlabel(r'$t/M$')
    ylabel(r'$v$')
    xlim(xmin=0)
    if axis()[1]>MaxTFactor*t_NR[-1]:
        xlim(xmax=MaxTFactor*t_NR[-1])
    ylim(ymax=0.7)
    title(NRDir);
    try:tight_layout(pad=0.1)
    except:pass
    savefig(AnnexDir+NRDir+OutputDir+'v.pdf');
    
    # Compare ellHat vectors
    fig = plt.figure()
    ax = fig.add_subplot(111, projection='3d')
    ax.plot(ellHat_NR[:,0], ellHat_NR[:,1], ellHat_NR[:,2])
    ax.plot([ellHat_NR[i_relaxed,0]], [ellHat_NR[i_relaxed,1]], [ellHat_NR[i_relaxed,2]], 'bo')
    for PNOrder in PNOrders:
        for Approximant in Approximants:
            Sym = "{1} {0}".format(Approximant,PNOrder)
            ax.plot(ellHat[Sym][:,0], ellHat[Sym][:,1], ellHat[Sym][:,2])
    title(r'$\hat{\ell}$ for ' + NRDir);
    try:tight_layout(pad=0.1)
    except:pass
    savefig(AnnexDir+NRDir+OutputDir+'ellHat.pdf');
    
    # Compare relative orientations
    figure()
    next(gca()._get_lines.color_cycle) # increment line-color cycle to match previous plots
    dTheta = {}
    for PNOrder in PNOrders:
        for Approximant in Approximants:
            Sym = "{1} {0}".format(Approximant,PNOrder)
            indices_NR = arange(len(t_NR))[(t_NR>=0) & (t_NR<=t[Sym][-1])]
            R_PN = array(Quaternions.Squad(R_frame[Sym], t[Sym]-t[Sym][0], t_NR[indices_NR]))
            dTheta[Sym] = Quaternions.angle(Quaternions.UnflipRotors(Quaternions.RDelta(R_PN, R_NR[indices_NR])))
            semilogy(t_NR[indices_NR], dTheta[Sym], label="Taylor{0} at {1} PN".format(Approximant,PNOrder))
    xlim(xmin=0)
    if axis()[1]>MaxTFactor*t_NR[-1]:
        xlim(xmax=MaxTFactor*t_NR[-1])
    ylim(ymin=8e-5)
    legend(loc='upper left')
    xlabel(r'$t/M$')
    ylabel(r'$\Phi_{\Delta}\ =\ 2\, \left| \log \left( R_{\mathrm{PN}}\, \bar{R}_{\mathrm{NR}} \right) \right|$');
    title(NRDir);
    try:tight_layout(pad=0.1)
    except:pass
    savefig(AnnexDir+NRDir+OutputDir+'Phases.pdf');

    return t,v,chi1,chi2,R_frame,Phi,ellHat



    











In [13]:



    


for StartFromTEqualsZero in [True, False]:
    print("\nStartFromTEqualsZero={0}".format(StartFromTEqualsZero))
    for i,NRDir in enumerate(NRDirs):
        print("\nNRDir='{0}'".format(NRDir))
        sys.stdout.flush()
        t,v,chi1,chi2,R_frame,Phi,ellHat = RunPNComparisons(PNOrders, Approximants, AnnexDir, NRDir,
                                                            StartFromTEqualsZero, InitialOmegas[i])



    


















    






StartFromTEqualsZero=True

NRDir='SKS/d10.0_q6.00_sA_0.739_0.190_-0.495_sB_-0.190_0.051_-0.227/Lev3/'
# PNEvolution.EvolvePN_Q("TaylorT1", 0.0, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT4", 0.0, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT5", 0.0, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT1", 2.0, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT4", 2.0, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT5", 2.0, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT1", 3.5, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT4", 3.5, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT5", 3.5, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT1", 6.0, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT4", 6.0, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])
# PNEvolution.EvolvePN_Q("TaylorT5", 6.0, 0.307353340494, 0.307353340494, 0.857142848195, 0.142857183655, [ 0.73933902  0.19018283 -0.49528598], [-0.18960235  0.05103283 -0.22681722], [0.999999626998706, -0.000863714333432979, 7.5917872216519e-12, -1.15883520846624e-11])

NRDir='q1.0/SBBH/d19.0_q1.0_s0.5_0_0_s0_0_0/Lev5/'
# PNEvolution.EvolvePN_Q("TaylorT1", 0.0, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 0.0, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 0.0, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT1", 2.0, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 2.0, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 2.0, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT1", 3.5, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 3.5, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 3.5, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT1", 6.0, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 6.0, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 6.0, 0.224309478846, 0.224309478846, 0.499999999868, 0.500000000048, [  5.00000000e-01  -6.53288000e-11  -7.53680000e-12], [  5.53092600e-10   3.55121000e-11  -1.28481970e-08], [0.999999954419313, 0.000301929415611176, 3.35220331961847e-30, 0])

NRDir='q5.0/SBBH/d15.0_q5.0_s0.5_0_0_s0/Lev5/'
# PNEvolution.EvolvePN_Q("TaylorT1", 0.0, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 0.0, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 0.0, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT1", 2.0, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 2.0, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 2.0, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT1", 3.5, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 3.5, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 3.5, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT1", 6.0, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 6.0, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 6.0, 0.250952938249, 0.250952938249, 0.83333333343, 0.166666666625, [  5.00000000e-01  -2.55540000e-12   3.10540000e-12], [  8.94044600e-10  -4.54100000e-13   3.32535540e-09], [0.999999999962973, 8.60549059541362e-06, 3.86073112423494e-30, 0])

NRDir='q5.0/SBBH/d15.0_q5.0_s0_0_0.5_s0/Lev4/'
# PNEvolution.EvolvePN_Q("TaylorT1", 0.0, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 0.0, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 0.0, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT1", 2.0, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 2.0, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 2.0, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT1", 3.5, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 3.5, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 3.5, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT1", 6.0, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT4", 6.0, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])
# PNEvolution.EvolvePN_Q("TaylorT5", 6.0, 0.250461598777, 0.250461598777, 0.833333333366, 0.166666666645, [ -2.22000000e-14  -2.40000000e-15   5.00000000e-01], [ -1.53720000e-12   7.86700000e-12   1.61732820e-09], [1, -1.33800949621651e-11, 2.92478856639308e-35, 0])

StartFromTEqualsZero=False

NRDir='SKS/d10.0_q6.00_sA_0.739_0.190_-0.495_sB_-0.190_0.051_-0.227/Lev3/'
# PNEvolution.EvolvePN_Q("TaylorT1", 0.0, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT4", 0.0, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT5", 0.0, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT1", 2.0, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT4", 2.0, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT5", 2.0, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT1", 3.5, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT4", 3.5, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT5", 3.5, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT1", 6.0, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT4", 6.0, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])
# PNEvolution.EvolvePN_Q("TaylorT5", 6.0, 0.308748166624, 0.308748166624, 0.857142505489, 0.142856965002, [ 0.703729    0.31316407 -0.4838715 ], [-0.17048677 -0.07178834 -0.23619506], [0.778267178555199, 0.0417393311497245, 0.0871005525223144, -0.62046073265696])

NRDir='q1.0/SBBH/d19.0_q1.0_s0.5_0_0_s0_0_0/Lev5/'
# PNEvolution.EvolvePN_Q("TaylorT1", 0.0, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT4", 0.0, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT5", 0.0, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT1", 2.0, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT4", 2.0, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT5", 2.0, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT1", 3.5, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT4", 3.5, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT5", 3.5, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT1", 6.0, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT4", 6.0, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])
# PNEvolution.EvolvePN_Q("TaylorT5", 6.0, 0.224611820203, 0.224611820203, 0.499999994797, 0.499999999644, [  4.96599000e-01   5.27498000e-02  -2.91492000e-04], [  5.67841000e-08   4.21627000e-08   9.17758000e-08], [0.377515261975278, 0.00301177004669471, -0.00383586362773967, 0.925990519588161])

NRDir='q5.0/SBBH/d15.0_q5.0_s0.5_0_0_s0/Lev5/'
# PNEvolution.EvolvePN_Q("TaylorT1", 0.0, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT4", 0.0, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT5", 0.0, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT1", 2.0, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT4", 2.0, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT5", 2.0, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT1", 3.5, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT4", 3.5, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT5", 3.5, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT1", 6.0, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT4", 6.0, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])
# PNEvolution.EvolvePN_Q("TaylorT5", 6.0, 0.251324042787, 0.251324042787, 0.833333341231, 0.166666661219, [  4.98408000e-01   2.70544000e-02   1.01614000e-04], [  1.54123000e-09  -4.07429000e-09   3.20314000e-07], [0.305724813170255, 0.00622613400136661, -0.0112041247863449, 0.952033634623894])

NRDir='q5.0/SBBH/d15.0_q5.0_s0_0_0.5_s0/Lev4/'










    






/usr/local/lib/python2.7/site-packages/matplotlib/pyplot.py:412: RuntimeWarning: More than 20 figures have been opened. Figures created through the pyplot interface (`matplotlib.pyplot.figure`) are retained until explicitly closed and may consume too much memory. (To control this warning, see the rcParam `figure.max_num_figures`).
  max_open_warning, RuntimeWarning)










    






# PNEvolution.EvolvePN_Q("TaylorT1", 0.0, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT4", 0.0, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT5", 0.0, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT1", 2.0, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT4", 2.0, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT5", 2.0, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT1", 3.5, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT4", 3.5, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT5", 3.5, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT1", 6.0, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT4", 6.0, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])
# PNEvolution.EvolvePN_Q("TaylorT5", 6.0, 0.251011466869, 0.251011466869, 0.833333329077, 0.166666663946, [  1.08187000e-12  -3.58928000e-13   4.99552000e-01], [ -2.29015000e-12   2.46625000e-13   2.87615000e-07], [0.0853022740117124, 3.00207171099692e-12, -7.52248469623431e-12, 0.996355118443437])










    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    

















    
























In [ ]:

Content source: moble/PostNewtonian

Similar notebooks:

notebook.community | gallery | about