MSTM

cricbessel(n, ds)

Compute the Riccati-Bessel function, return complex values.

cspherebessel(n, z)

Compute spherical Bessel functions jn(z) & yn(z)

eulerrotation(xp)

Euler rotation of a point (xp[1], xp[2], xp[3]).

gaussianbeamcoef(ctx, α, β, cbeam, nodr)

Regular VSWF expansion coefficients for a gaussian beam, localized approximation.

cbeam = 1/(k omega)

gentrancoef(ctx, itype, xptran, ri, nrow0, nrow1, ncol0, ncol1, iaddrow0, iaddcol0)

Calculates the vwh translation coefficients for a general translation from one origin to another.

msta1(x, mp)

Determine the starting point for backward recurrence such that the magnitude of Jn(x) at that point is about 10^(-MP).

pifunc(cb, ephi, nmax, ndim)

Vector spherical harmonic function.

TODO: ndim seems to be useless. Can it be safely deleted?

ricbessel(n, ds, ϵ, Ψ)

Compute the Ricatti-Bessel function.

References:

• Riccati-Bessel Functions at Wolfram MathWorld

sphereplanewavecoef(ctx, nsphere, neqns, nodr, nodrmax, α, β, rpos, hostsphere, numberfieldexp, rimedium)

Plane wave expansion coefficients at sphere origins. A phase shift is used.

taufunc(ctx, cb, nmax)

Normalized vector spherical harmonic function.

tranordertest(ctx, r, ri, lmax, ϵ)

Test to determine convergence of regular VSWF addition theorem for max order lmax and translation distance r w/ refractive index ri.

vcfuncuprec(ctx, m, n, k, l, wmax, vcn)

Upwards VC coefficient recurrence, used internally by vcfunc(ctx, m, n, k, l).

MSTMConfig(filename::String)

Read configurations from a .toml file.

TMaxtrix(filename::String)

Read a previously calculated T matrix.

mieoa(x, ri, nodr0, qeps, ri_medium)

Optically active Lorenz/Mie coefficients

multmiecoeffmult(mie, nrhs, idir, ain, rhs_list)

Generalized Mie coefficient multiplication.

If idir != 1, we do the transpose.