Working:

- implemented object-oriented interface:
  - structure
  - material
  - incident field

- "core" module:
  - setup side-by-side matrix (using non-retarded "123" propagator)
  - LU inversion
  - kwargs based processing of incident field parameters using generalized propagator at each wavelength

- Added testscript which includes comparison and relative error checking with old pyGDM implementation