print(filmMaterial.getTensor())
"""
matrix([[ 4.5625, -0.5625, 0.7955],
[-0.5625, 4.5625, -0.7955],
[ 0.7955, -0.7955, 5.125 ]])
"""
air = Berreman4x4.IsotropicNonDispersiveMaterial(n_i)
front = Berreman4x4.IsotropicHalfSpace(air) # Front half-space
Kx = front.get_Kx_from_Phi(Phi_i)
print("\nValue of Kx: {:.4f}".format(Kx))
"""
Kx = 0.9397
"""
film = Berreman4x4.HomogeneousLayer(filmMaterial, d)
Delta = film.getDeltaMatrix(Kx)
print("\nDelta matrix:")
print(Delta)
"""
matrix([[-0.1459, 0.1459, 0. , 0.8277],
[ 0. , 0. , -1. , 0. ],
[ 0.439 , -3.556 , 0. , -0.1459],
[ 4.439 , -0.439 , 0. , -0.1459]])
"""
q, Psi = scipy.linalg.eig(Delta)
print("\nEigenvalues of the Delta matrix (eq 6.64, p.241):")
print(numpy.real(q))
""" (eq 6.64)
real(q) = array([-2.174 , -1.7655, 1.7655, 1.8822])
imag(q) = array([ 0., 0., 0., 0.])
