static void TestPotential(Func <LorentzVectorVariableU, LorentzVectorU> createPotential, EuclideanVector3 v)
{
LorentzVectorVariableU coordinates = new LorentzVectorVariableU("t", "x", "y", "z");
LorentzVectorU potential = createPotential(coordinates);
LorentzVectorOperatorL del = coordinates.Del;
ElectromagneticField field0 = new ElectromagneticField(potential, del);
TestField(field0);
ElectromagneticField field1 = field0.Transform(v);
TestField(field1);
TestFieldTransforms(field0.ElectricField, field0.MagneticField, field1.ElectricField, field1.MagneticField, v);
CompareInvariantScalars(field0, field1);
ElectromagneticField field2 = field1.Transform(-v);
TestField(field2);
TestFieldTransforms(field1.ElectricField, field1.MagneticField, field2.ElectricField, field2.MagneticField, -v);
CompareInvariantScalars(field1, field2);
Assert.AreEqual(field0.Potential, field2.Potential, "A dual transform");
Assert.AreEqual(field0.Del.ToString(), field2.Del.ToString(), "Del dual transform");
Assert.AreEqual(field0.ElectricField, field2.ElectricField, "E dual transform");
Assert.AreEqual(field0.MagneticField, field2.MagneticField, "B dual transform");
Assert.AreEqual(field0.CurrentDensity, field2.CurrentDensity, "J dual transform");
}
static void TestLorentzTransform(EuclideanVector3 v)
{
LorentzVectorVariableU x0 = new LorentzVectorVariableU("t", "x", "y", "z");
LorentzVectorU x1 = x0.Transform(v);
LorentzVectorU x2 = x1.Transform(-v);
Assert.AreEqual(x0, x2);
Assert.AreEqual(LorentzMatrixUL.One, LorentzTransform.Matrix(v) * LorentzTransform.Matrix(-v));
}
