#ifndef COVMATRIXSIMILARITY6_H
#define COVMATRIXSIMILARITY6_H

// $Id: Similarity6.h,v 1.2 2003/07/03 15:11:25 sachs Exp $
// David Sachs - Fermilab - July 2003
//
// Define Similarity transformation functions for 6x6 CovMatrices

// The normal intent of these functions is for P to be an orthogonal
// matrix, but the transformations are valid for any matrix P
// The matrix P is supplied as an array of Vector6

#include "CovMatrices/CovMatrix6.h"

namespace CovMatrices
{
  // Provide transform result = P * A * Pt
  void similarityPAPt( 
	const CovMatrix6& a, 	// original CovMatrix
	const Vector6 p[],	// transformation matrix (size 6)
	CovMatrix6& result);	// result CovMatrix

  // Provide transform: result = Pt * A * P
  void similarityPtAP( 
	const CovMatrix6& a, 	// original CovMatrix
	const Vector6 p[],	// transformation matrix (size 6)
	CovMatrix6& result);	// result CovMatrix
}

#endif