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


    

    




---------------------------------------------------------------------------
ImportError                               Traceback (most recent call last)
<ipython-input-2-2bd0a1a1b35a> in <module>()
      3 import matplotlib.pyplot as plt
      4 from mpl_toolkits.mplot3d import Axes3D
----> 5 import PNEvolution
      6 import Quaternions

/Users/boyle/Research/GWFrames/Code/PostNewtonian/C++/PNEvolution/__init__.py in <module>()
     27                 fp.close()
     28             return _mod
---> 29     _PNEvolution = swig_import_helper()
     30     del swig_import_helper
     31 else:

/Users/boyle/Research/GWFrames/Code/PostNewtonian/C++/PNEvolution/__init__.py in swig_import_helper()
     19             fp, pathname, description = imp.find_module('_PNEvolution', [dirname(__file__)])
     20         except ImportError:
---> 21             import _PNEvolution
     22             return _PNEvolution
     23         if fp is not None:

ImportError: No module named _PNEvolution

In [3]:

    

Approximant = "T1"
PNOrder = 3.5
v0 = 0.2
m1 = 0.75
m2 = 1-m1
chi1_i = [0.5,0.2,0.4]
chi2_i = [0.2,0.4,0.5]
R_frame_i = Quaternions.Quaternion(1,0,0,0)

t,v,chi1,chi2,R_frame,Phi,L = \
    PNEvolution.EvolvePN_Q("Taylor{0}".format(Approximant), PNOrder, v0, v0, m1, m2, chi1_i, chi2_i, R_frame_i)


    

In [4]:

    

data,vB,chi1B,chi2B = PNEvolution.WaveformModes(m1, m2, v, chi1, chi2)
data = array(data)


    

In [5]:

    

semilogy(t, abs(data[0]))


    

    
Out[5]:





[<matplotlib.lines.Line2D at 0x106ff8fd0>]