// ----------------------------------------------------------------------
//
// PhysVec.cc
//
// History:
//   01-Jun-1998  WEB  Initial version
//   15-Jun-1998  WEB  Added namespace support
//   05-Mar-2002  MF   Eliminated use of RotationInterface classes
//
// ----------------------------------------------------------------------

#include "PhysicsVectors/PhysVec.h"


ZM_BEGIN_NAMESPACE( zmpv )	/*  namespace zmpv  {  */


HepLorentzRotation transposeOf (const HepLorentzRotation & lt) {
  HepRep4x4 rep_;
  HepRep4x4 lt_ = lt.rep4x4();

  rep_.xx_=lt_.xx_;  rep_.xy_=lt_.yx_;  rep_.xz_=lt_.zx_;  rep_.xt_=lt_.tx_;
  rep_.yx_=lt_.xy_;  rep_.yy_=lt_.yy_;  rep_.yz_=lt_.zy_;  rep_.yt_=lt_.ty_;
  rep_.zx_=lt_.xz_;  rep_.zy_=lt_.yz_;  rep_.zz_=lt_.zz_;  rep_.zt_=lt_.tz_;
  rep_.tx_=lt_.xt_;  rep_.ty_=lt_.yt_;  rep_.tz_=lt_.zt_;  rep_.tt_=lt_.tt_;

  return HepLorentzRotation ( rep_ );
}

HepLorentzRotation contraTransposeOf (const HepLorentzRotation & lt) {
  HepRep4x4 rep_;
  HepRep4x4 lt_ = lt.rep4x4();

  rep_.xx_=lt_.xx_;   rep_.xy_=lt_.yx_;   rep_.xz_=lt_.zx_;  rep_.xt_= -lt_.tx_;
  rep_.yx_=lt_.xy_;   rep_.yy_=lt_.yy_;   rep_.yz_=lt_.zy_;  rep_.yt_= -lt_.ty_;
  rep_.zx_=lt_.xz_;   rep_.zy_=lt_.yz_;   rep_.zz_=lt_.zz_;  rep_.zt_= -lt_.tz_;
  rep_.tx_= -lt_.xt_; rep_.ty_= -lt_.yt_; rep_.tz_= -lt_.zt_;rep_.tt_=  lt_.tt_;

  return HepLorentzRotation ( rep_ );
}


ZM_END_NAMESPACE( zmpv )	/*  }  // namespace zmpv  */