// Lorentz transformation matrix (boost along x axis)
matrix L(4, 4);
- L(0, 0) = gamma;
- L(0, 1) = -beta*gamma;
- L(1, 0) = -beta*gamma;
- L(1, 1) = gamma;
- L(2, 2) = 1; L(3, 3) = 1;
+ L = gamma, -beta*gamma, 0, 0,
+ -beta*gamma, gamma, 0, 0,
+ 0, 0, 1, 0,
+ 0, 0, 0, 1;
// Electromagnetic field tensor
- matrix F(4, 4, lst(
- 0, -Ex, -Ey, -Ez,
+ matrix F(4, 4);
+ F = 0, -Ex, -Ey, -Ez,
Ex, 0, -Bz, By,
Ey, Bz, 0, -Bx,
- Ez, -By, Bx, 0
- ));
+ Ez, -By, Bx, 0;
// Indices
symbol s_mu("mu"), s_nu("nu"), s_rho("rho"), s_sigma("sigma");