/// <summary>Initializes this Metric object from metric matrix (called by constructor) </summary> /// <param name="m">The NxN metric matrix</param> private void init(DotNetMatrix.GeneralMatrix m) { if (!Util.IsSymmetric(m)) { throw new Exception("The metric matrix must be symmetric"); } m_matrix = m.Copy(); // System.Console.WriteLine("m_matrix: " + Util.ToString(m_matrix)); // compute eigen value decomposition m_eig = new DotNetMatrix.EigenvalueDecomposition(m_matrix); // System.Console.WriteLine("m_eig: " + Util.ToString(m_eig.GetV())); m_invEigMatrix = m_eig.GetV().Transpose(); m_eigenMetric = m_eig.RealEigenvalues; // { // DotNetMatrix.GeneralMatrix D = Util.Diagonal(m_eigenMetric); // DotNetMatrix.GeneralMatrix tmp = m_eig.GetV().Multiply(D).Multiply(m_invEigMatrix); // System.Console.WriteLine("input_matrix = " + Util.ToString(tmp)); // } m_isDiagonal = Util.IsDiagonal(m_matrix); if (!m_isDiagonal) { m_isEuclidean = m_isAntiEuclidean = false; } else { m_isEuclidean = m_isAntiEuclidean = true; for (int i = 0; i < m.RowDimension; i++) { if (m_matrix.GetElement(i, i) != 1.0) { m_isEuclidean = false; } if (m_matrix.GetElement(i, i) != -1.0) { m_isAntiEuclidean = false; } } } }