def test_raw_fields(): sp = Sphere(r=.5, n=1.6, center=(10, 10, 5)) wavelen = .66 index = 1.33 pol = to_vector((0, 1)) sch = detector_grid(3, .1) wavevec=2*np.pi/(wavelen/index) pos = Mie._transform_to_desired_coordinates( sch, (10, 10, 5), wavevec=wavevec) rf = Mie()._raw_fields( pos, sp, medium_wavevec=wavevec, medium_index=index, illum_polarization=pol) assert_allclose(rf, [[(0.0015606995428858754-0.0019143174710834162j), (-0.0003949071974815011-0.0024154494284017187j), (-0.002044525390662322-0.001302770747742109j), (-0.0003949071974815009-0.002415449428401719j), (-0.002055824337886397-0.0012853546864338861j), (-0.00230285180386436+0.000678693819245102j), (-0.0020445253906623225-0.0013027707477421095j), (-0.0023028518038643603+0.0006786938192451026j), (-0.0010011090105680883+0.0021552249454706712j)], [(-0.0010507058414478587+0.0036584360153097306j), (0.0020621595919700776+0.003210547679920805j), (0.0037794246074692407+0.000585690417403587j), (0.0020542215584045407+0.0031619947065620246j), (0.0037426710578253295+0.000527040269055415j), (0.002871631795307833-0.002470099566862354j), (0.0036968090916832948+0.0005330478443315597j), (0.002824872178181336-0.0024563186266035124j), (2.261564613123139e-06-0.003751168280253104j)], [(0.0010724312167657794+0.0039152445632936j), (0.003651474601303447+0.0017688083711547462j), (0.003740131549224567-0.001566271371618957j), (0.0036883581831347947+0.0017866751223785315j), (0.0037648739662344477-0.001614943488355339j), (0.0012643679510138835-0.003894481935619062j), (0.003816460764514863-0.0015982360934887314j), (0.0012772696647997395-0.0039342215472070105j), (-0.0021320123934202356-0.0035427449839031066j)]])