  We develop a correspondence between arbitrary tensors and matrices based on
the use of Kronecker products and associated identities.
  Utilizing the rules of matrix differentiation we derive the vacuum Einstein
field equations as a differential-matrix equation. This formulation may
facilitate their efficient use in numerical relativistic models.
