public Vertex(Tags Data, double x, double y, double z) { this.x = x; this.y = y; this.z = z; isin = false; wasMin = false; this.Data = Data; }
public Clusters(double[,] A, Tags[] texts, Tags[] words) { r = 1; V = new List<Vertex>(); textTitles = texts; C = new List<Cluster>(); Texts = new Cluster(); double[] W = new double[A.GetLength(0)]; double[,] U = new double[A.GetLength(0), A.GetLength(1)]; double[,] VT = new double[A.GetLength(1), A.GetLength(1)]; alglib.rmatrixsvd(A, A.GetLength(0), A.GetLength(1), 1, 1, 0, out W, out U, out VT); #region /************************************************************************* Singular value decomposition of a rectangular matrix. The algorithm calculates the singular value decomposition of a matrix of size MxN: A = U * S * V^T The algorithm finds the singular values and, optionally, matrices U and V^T. The algorithm can find both first min(M,N) columns of matrix U and rows of matrix V^T (singular vectors), and matrices U and V^T wholly (of sizes MxM and NxN respectively). Take into account that the subroutine does not return matrix V but V^T. Input parameters: A - matrix to be decomposed. Array whose indexes range within [0..M-1, 0..N-1]. M - number of rows in matrix A. N - number of columns in matrix A. UNeeded - 0, 1 or 2. See the description of the parameter U. VTNeeded - 0, 1 or 2. See the description of the parameter VT. AdditionalMemory - If the parameter: * equals 0, the algorithm doesn’t use additional memory (lower requirements, lower performance). * equals 1, the algorithm uses additional memory of size min(M,N)*min(M,N) of real numbers. It often speeds up the algorithm. * equals 2, the algorithm uses additional memory of size M*min(M,N) of real numbers. It allows to get a maximum performance. The recommended value of the parameter is 2. Output parameters: W - contains singular values in descending order. U - if UNeeded=0, U isn't changed, the left singular vectors are not calculated. if Uneeded=1, U contains left singular vectors (first min(M,N) columns of matrix U). Array whose indexes range within [0..M-1, 0..Min(M,N)-1]. if UNeeded=2, U contains matrix U wholly. Array whose indexes range within [0..M-1, 0..M-1]. VT - if VTNeeded=0, VT isn’t changed, the right singular vectors are not calculated. if VTNeeded=1, VT contains right singular vectors (first min(M,N) rows of matrix V^T). Array whose indexes range within [0..min(M,N)-1, 0..N-1]. if VTNeeded=2, VT contains matrix V^T wholly. Array whose indexes range within [0..N-1, 0..N-1]. *************************************************************************/ #endregion //Добавление вершин из матриц for (int i = 0; i < U.GetLength(0); i++) { V.Add(new Vertex(words[i], U[i, 0], U[i, 1], U[i, 2])); } for (int i = 0; i < VT.GetLength(1); i++) { V.Add(new Vertex(texts[i], VT[0, i], VT[1, i], VT[2, i])); } for (int i = 0; i < texts.Length; i++) { for (int j = 0; j < V.Count; j++) { if (texts[i].GetTag == V[j].Data.GetTag) { Texts.Data.Add(V[j]); } } } }