Ei = electric field
Bi = magnetic field
Pi = polarization field
Mi = magnetization field
ϵ0 = permittivity of free space
ϵr,ij = relative permittivity of material
ϵij = permittivity of material
χe,ij = electric susceptibility
μ0 = permeability of free space
μr,ij = relative permeability of material
μij = permeability of material
χm,ij = magnetic susceptibility

Di=ϵ0Ei+Pi
Pi=ϵ0χe,ijEj
Di=ϵ0Ei+ϵ0χe,ijEj
Di=ϵ0(δij+χe,ij)Ej
Di=ϵ0ϵr,ijEj
Di=ϵijEj

Bi=μ0Hi+Mi
Mi=μ0χm,ijHj
Bi=μ0Hi+μ0χm,ijHj
Bi=μ0(δij+χm,ij)Hj
Bi=μ0μr,ijHj
Bi=μijHj


Ei=(ϵ1)ijDj
Hi=(μ1)ijBj

Di,i=ρfree ... where ρfree = free charge (densitized?)
Bi,i=0

Di,t=ϵijkHk,jJi,free
Bi,t=ϵijkEk,j

Di,t=ϵijk(μ1B)k,jJi,free
Bi,t=ϵijk(ϵ1D)k,j

Di,t=ϵijk((μ1)klBl),jJi,free
Bi,t=ϵijk((ϵ1)klDl),j

Di,tγimϵmjk(μ1)klBl,j=ϵijk(μ1)kl,jBlJi,free
Bi,t+γimϵmjk(ϵ1)klDl,j=ϵijk(ϵ1)kl,jDl

Di,t1γγimϵ¯mjk(μ1)klBl,j=ϵijk(μ1)kl,jBlJi,free
Bi,t+1γγimϵ¯mjk(ϵ1)klDl,j=ϵijk(ϵ1)kl,jDl

(Di,tBi,t)+1γ(γim00γim)(0ϵ¯mjkϵ¯mjk0)((μ1)kl00(ϵ1)kl)(Bl,jDl,j)=(ϵijk(μ1)kl,jBlJi,freeϵijk(ϵ1)kl,jDl)