In [1]:

    

import matplotlib.pyplot as plt
import numpy as np
import time
import pickle
from matplotlib.path import Path
from dirac_sheet import dirac_sheet
from ipywidgets import interact, interactive, fixed
import ipywidgets as widgets
from IPython import display


    

In [65]:

    

y,x = np.mgrid[slice(-10,10,.025),slice(-10,10,.025)]

u1=np.zeros((x.shape[0],y.shape[0]),dtype=np.complex)
p_0=1
d_theta=30/180.0*np.pi
for theta in np.linspace(-np.pi/2,np.pi/2,100):
    u1+=np.exp(-(theta**2/d_theta**2/2))*np.exp(1j*p_0*(np.cos(theta)*x+np.sin(theta)*y))

fig = plt.figure()
plt.ion()
plt.clf()
#plt.imshow(np.imag(u1*np.conj(np.gradient(u1)[1])-np.conj(u1)*np.gradient(u1)[1]),cmap='RdBu',interpolation='nearest', origin='lower')
#plt.colorbar()
plt.plot(np.sum(np.imag(u1*np.conj(np.gradient(u1)[1])-np.conj(u1)*np.gradient(u1)[1]),axis=0))
fig.canvas.draw()


    

In [ ]:

    




    

In [ ]: