Exemplo n.º 1
0
def bpm_3d_spheres(size, units, lam = .5, u0 = None,
                   points = None,
                   dn_inner = 0.0, rad_inner = 0.,                   
                   dn_outer = 0.05, rad_outer = 2.5,
                   return_scattering = False,
                   use_fresnel_approx = False
):
    """
    simulates the propagation of monochromativ wave of wavelength lam with initial conditions u0 along z in a media filled with coated spheres at positions points and inner radius r1 and dn1 and outer radius r2/ dn2

    size     -    the dimension of the image to be calulcated  in pixels (Nx,Ny,Nz)
    units    -    the unit lengths of each dimensions in microns
    lam      -    the wavelength
    u0       -    the initial field distribution, if u0 = None an incident  plane wave is assumed
    points   -    a list of coordinates [[x1,y1,z1],...] in real unit interval
                  [0...dx*nx] where spheres are placed.

    """
    dn_g = create_dn_buffer(size,units,points = points,
                   dn_inner = dn_inner, rad_inner = rad_inner,                   
                   dn_outer = dn_outer, rad_outer = rad_outer)

    return bpm_3d(size = size,units = units,lam = lam,
                  u0 = u0,dn = dn_g,
                  return_scattering= return_scattering,
                  use_fresnel_approx = use_fresnel_approx)
Exemplo n.º 2
0
def bpm_3d_spheres(size,
                   units,
                   lam=.5,
                   u0=None,
                   points=None,
                   dn_inner=0.0,
                   rad_inner=0.,
                   dn_outer=0.05,
                   rad_outer=2.5,
                   return_scattering=False,
                   use_fresnel_approx=False):
    """
    simulates the propagation of monochromativ wave of wavelength lam with initial conditions u0 along z in a media filled with coated spheres at positions points and inner radius r1 and dn1 and outer radius r2/ dn2

    size     -    the dimension of the image to be calulcated  in pixels (Nx,Ny,Nz)
    units    -    the unit lengths of each dimensions in microns
    lam      -    the wavelength
    u0       -    the initial field distribution, if u0 = None an incident  plane wave is assumed
    points   -    a list of coordinates [[x1,y1,z1],...] in real unit interval
                  [0...dx*nx] where spheres are placed.

    """
    dn_g = create_dn_buffer(size,
                            units,
                            points=points,
                            dn_inner=dn_inner,
                            rad_inner=rad_inner,
                            dn_outer=dn_outer,
                            rad_outer=rad_outer)

    return bpm_3d(size=size,
                  units=units,
                  lam=lam,
                  u0=u0,
                  dn=dn_g,
                  return_scattering=return_scattering,
                  use_fresnel_approx=use_fresnel_approx)
Exemplo n.º 3
0
    dz = 1.*Lz/(Nz-1)

    
    lam = .5

    x = dx*arange(Nx)
    y = dx*arange(Nx)
    z = dz*arange(Nz)

    X0,Y0 = meshgrid(x,y,indexing="ij")



    _xgrid = (X0%a)-.5*a
    _ygrid = (Y0%a)-.5*a

    R = sqrt(_xgrid**2+_ygrid**2)

    
    u0 = (jn(1,4*R)+2.e-10)/4./(1.e-10+R)+0.j
    u0 = u0.astype(complex64)

    # u0 = exp(-20*(_xgrid**2+_ygrid**2))

    # u0 *= exp(-1.5*(R0/dx/Nx*3.)**2)

    u, dn = bpm_3d((Nx,Nx,Nz), (dx,dx,dz), lam = lam ,u0 = u0,
                   use_fresnel_approx = False)

    overlaps = array([overlap(abs(u[0,...])**2,abs(u[i,...])**2) for i in range(Nz)])
Exemplo n.º 4
0
    dz = 1. * Lz / (Nz - 1)

    lam = .5

    x = dx * arange(Nx)
    y = dx * arange(Nx)
    z = dz * arange(Nz)

    X0, Y0 = meshgrid(x, y, indexing="ij")

    _xgrid = (X0 % a) - .5 * a
    _ygrid = (Y0 % a) - .5 * a

    R = sqrt(_xgrid**2 + _ygrid**2)

    u0 = (jn(1, 4 * R) + 2.e-10) / 4. / (1.e-10 + R) + 0.j
    u0 = u0.astype(complex64)

    # u0 = exp(-20*(_xgrid**2+_ygrid**2))

    # u0 *= exp(-1.5*(R0/dx/Nx*3.)**2)

    u, dn = bpm_3d((Nx, Nx, Nz), (dx, dx, dz),
                   lam=lam,
                   u0=u0,
                   use_fresnel_approx=False)

    overlaps = array(
        [overlap(abs(u[0, ...])**2,
                 abs(u[i, ...])**2) for i in range(Nz)])