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)]])
