示例#1
0
        /// <summary>
        /// Creates a molecule graph for use with jgrapht.
        /// Bond orders are not respected.
        /// </summary>
        /// <param name="molecule">the specified molecule</param>
        /// <returns>a graph representing the molecule</returns>
	    static public SimpleGraph getMoleculeGraph(IAtomContainer molecule) {
		    SimpleGraph graph = new SimpleGraph();
		    for (int i=0; i<molecule.AtomCount; i++	) {
			    IAtom atom = molecule.Atoms[i];
			    graph.addVertex(atom);
		    }
    		
		    for (int i=0; i<molecule.getBondCount(); i++	) {
			    IBond bond = molecule.Bonds[i];
    			
			    /*
			    int order = (int) bond.getOrder();
			    for (int j=0; j<order; j++) {
				    graph.addEdge(bond.getAtoms()[0], bond.getAtoms()[1]);
			    }
			    */
			    graph.addEdge(bond.getAtoms()[0], bond.getAtoms()[1]);
		    }
		    return graph;
	    }
        public virtual System.Collections.IList equivalenceClasses()
        {
            int[] weight = weightVector();

            System.Object[] cyclesArray = (System.Object[])SupportClass.ICollectionSupport.ToArray(cycles_Renamed_Field);
            System.Array.Sort(cyclesArray, new AnonymousClassComparator(this));

            System.Collections.ICollection essentialCycles = this.essentialCycles();

            bool[][] u = new bool[cyclesArray.Length][];
            for (int i = 0; i < cyclesArray.Length; i++)
            {
                u[i] = new bool[edgeList.Count];
            }

            bool[][] a = getCycleEdgeIncidenceMatrix(cyclesArray);
            bool[][] ai = inverseBinaryMatrix(a, cyclesArray.Length);

            for (int i = 0; i < cyclesArray.Length; i++)
            {
                for (int j = 0; j < cyclesArray.Length; j++)
                {
                    u[i][j] = ai[j][i];
                }
            }

            UndirectedGraph h = new SimpleGraph();
            h.addAllVertices(cycles_Renamed_Field);

            ConnectivityInspector connectivityInspector = new ConnectivityInspector(h);

            int left = 0;
            for (int right = 0; right < weight.Length; right++)
            {
                if ((right < weight.Length - 1) && (weight[right + 1] == weight[right]))
                    continue;

                // cyclesArray[left] to cyclesArray[right] have same weight

                // First test (compute pre-classes):
                // Check if there is a cycle that can replace a[i] as well as a[j] in a basis
                // This is done by finding a cycle C with <C,u[i]>=1 and <C,u[j]>=1

                for (int i = left; i <= right; i++)
                {
                    if (SupportClass.ICollectionSupport.Contains(essentialCycles, (SimpleCycle)cyclesArray[i]))
                        continue;

                    for (int j = i + 1; j <= right; j++)
                    {
                        if (SupportClass.ICollectionSupport.Contains(essentialCycles, (SimpleCycle)cyclesArray[j]))
                            continue;

                        // check if cyclesArray[i] and cyclesArray[j] are already in the same class
                        if (connectivityInspector.pathExists(cyclesArray[i], cyclesArray[j]))
                            continue;

                        bool sameClass = false;

                        AuxiliaryGraph2 auxGraph = new AuxiliaryGraph2(this, graph, edgeList, u[i], u[j]);

                        //UPGRADE_TODO: Method 'java.util.Iterator.hasNext' was converted to 'System.Collections.IEnumerator.MoveNext' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javautilIteratorhasNext'"
                        for (System.Collections.IEnumerator it = graph.vertexSet().GetEnumerator(); it.MoveNext(); )
                        {
                            //UPGRADE_TODO: Method 'java.util.Iterator.next' was converted to 'System.Collections.IEnumerator.Current' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javautilIteratornext'"
                            System.Object vertex = it.Current;

                            // check if the vertex is incident to an edge with u[edge] == 1
                            bool shouldSearchCycle = false;

                            System.Collections.ICollection incidentEdges = graph.edgesOf(vertex);

                            System.Collections.IEnumerator edgeIterator = incidentEdges.GetEnumerator();
                            //UPGRADE_TODO: Method 'java.util.Iterator.hasNext' was converted to 'System.Collections.IEnumerator.MoveNext' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javautilIteratorhasNext'"
                            while (edgeIterator.MoveNext())
                            {
                                //UPGRADE_TODO: Method 'java.util.Iterator.next' was converted to 'System.Collections.IEnumerator.Current' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javautilIteratornext'"
                                Edge edge = (Edge)edgeIterator.Current;
                                int index = getEdgeIndex(edge);
                                if (u[i][index] || u[j][index])
                                {
                                    shouldSearchCycle = true;
                                    break;
                                }
                            }

                            if (shouldSearchCycle)
                            {

                                System.Object auxVertex00 = auxGraph.auxVertex00(vertex);
                                System.Object auxVertex11 = auxGraph.auxVertex11(vertex);

                                System.Collections.IList auxPath = BFSShortestPath.findPathBetween(auxGraph, auxVertex00, auxVertex11);

                                double pathWeight = auxPath.Count;

                                if (pathWeight == weight[left])
                                {
                                    sameClass = true;
                                    break;
                                }
                            }
                        }

                        if (sameClass)
                        {
                            h.addEdge(cyclesArray[i], cyclesArray[j]);
                        }
                    }
                }

                // Second test (compute equivalence classes):
                // Check if there are two cycle Ci, Cj that can replace a[i], a[j]
                // and have a common cycle a[k] in their basis representation
                // This is done by finding a cycle a[k] with <u[k],u[i]>=1 and <u[k],u[j]>=1

                for (int i = left; i <= right; i++)
                {
                    if (SupportClass.ICollectionSupport.Contains(essentialCycles, (SimpleCycle)cyclesArray[i]))
                        continue;

                    for (int j = i + 1; j <= right; j++)
                    {
                        if (SupportClass.ICollectionSupport.Contains(essentialCycles, (SimpleCycle)cyclesArray[j]))
                            continue;

                        // check if cyclesArray[i] and cyclesArray[j] are already in the same class
                        if (connectivityInspector.pathExists(cyclesArray[i], cyclesArray[j]))
                            continue;

                        bool sameClass = false;

                        for (int k = 0; ((SimpleCycle)cyclesArray[k]).weight() < weight[left]; k++)
                        {

                            AuxiliaryGraph2 auxGraph = new AuxiliaryGraph2(this, graph, edgeList, u[i], u[k]);

                            bool shortestPathFound = false;
                            //UPGRADE_TODO: Method 'java.util.Iterator.hasNext' was converted to 'System.Collections.IEnumerator.MoveNext' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javautilIteratorhasNext'"
                            for (System.Collections.IEnumerator it = graph.vertexSet().GetEnumerator(); it.MoveNext(); )
                            {
                                //UPGRADE_TODO: Method 'java.util.Iterator.next' was converted to 'System.Collections.IEnumerator.Current' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javautilIteratornext'"
                                System.Object vertex = it.Current;

                                System.Object auxVertex00 = auxGraph.auxVertex00(vertex);
                                System.Object auxVertex11 = auxGraph.auxVertex11(vertex);

                                System.Collections.IList auxPath = BFSShortestPath.findPathBetween(auxGraph, auxVertex00, auxVertex11);

                                double pathWeight = auxPath.Count;

                                if (pathWeight == weight[left])
                                {
                                    shortestPathFound = true;
                                    break;
                                }
                            }

                            if (!shortestPathFound)
                                continue;

                            auxGraph = new AuxiliaryGraph2(this, graph, edgeList, u[j], u[k]);

                            //UPGRADE_TODO: Method 'java.util.Iterator.hasNext' was converted to 'System.Collections.IEnumerator.MoveNext' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javautilIteratorhasNext'"
                            for (System.Collections.IEnumerator it = graph.vertexSet().GetEnumerator(); it.MoveNext(); )
                            {
                                //UPGRADE_TODO: Method 'java.util.Iterator.next' was converted to 'System.Collections.IEnumerator.Current' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javautilIteratornext'"
                                System.Object vertex = it.Current;

                                System.Object auxVertex00 = auxGraph.auxVertex00(vertex);
                                System.Object auxVertex11 = auxGraph.auxVertex11(vertex);

                                System.Collections.IList auxPath = BFSShortestPath.findPathBetween(auxGraph, auxVertex00, auxVertex11);

                                double pathWeight = auxPath.Count;

                                if (pathWeight == weight[left])
                                {
                                    sameClass = true;
                                    break;
                                }
                            }

                            if (sameClass)
                                break;
                        }

                        if (sameClass)
                        {
                            h.addEdge(cyclesArray[i], cyclesArray[j]);
                        }
                    }
                }

                left = right + 1;
            }

            return connectivityInspector.connectedSets();
        }
		private System.Collections.IList lazyFindBiconnectedSets()
		{
			if (biconnectedSets_Renamed_Field == null)
			{
				biconnectedSets_Renamed_Field = new System.Collections.ArrayList();
				
                IList inspector = new ConnectivityInspector(graph).connectedSets();
                System.Collections.IEnumerator connectedSets = inspector.GetEnumerator();
				
				while (connectedSets.MoveNext())
				{
                    object obj = ((DictionaryEntry)connectedSets.Current).Value;
                    if (!(obj is CSGraphT.SupportClass.HashSetSupport))
                        continue;
					CSGraphT.SupportClass.SetSupport connectedSet = (CSGraphT.SupportClass.SetSupport)obj;
					if (connectedSet.Count == 1)
					{
						continue;
					}
					
					org._3pq.jgrapht.Graph subgraph = new Subgraph(graph, connectedSet, null);
					
					// do DFS
					
					// Stack for the DFS
					System.Collections.ArrayList vertexStack = new System.Collections.ArrayList();
					
					CSGraphT.SupportClass.SetSupport visitedVertices = new CSGraphT.SupportClass.HashSetSupport();
					IDictionary parent = new System.Collections.Hashtable();
					IList dfsVertices = new System.Collections.ArrayList();
					
					CSGraphT.SupportClass.SetSupport treeEdges = new CSGraphT.SupportClass.HashSetSupport();

                    System.Object currentVertex = subgraph.vertexSet()[0];//.ToArray()[0];

					vertexStack.Add(currentVertex);
					visitedVertices.Add(currentVertex);
					
					while (!(vertexStack.Count == 0))
					{
						currentVertex = SupportClass.StackSupport.Pop(vertexStack);
						
						System.Object parentVertex = parent[currentVertex];
						
						if (parentVertex != null)
						{
							Edge edge = subgraph.getEdge(parentVertex, currentVertex);
							
							// tree edge
							treeEdges.Add(edge);
						}
						
						visitedVertices.Add(currentVertex);
						
						dfsVertices.Add(currentVertex);
						
						System.Collections.IEnumerator edges = subgraph.edgesOf(currentVertex).GetEnumerator();
						while (edges.MoveNext())
						{
							// find a neighbour vertex of the current vertex 
							Edge edge = (Edge)edges.Current;
							
							if (!treeEdges.Contains(edge))
							{
								System.Object nextVertex = edge.oppositeVertex(currentVertex);
								
								if (!visitedVertices.Contains(nextVertex))
								{
									vertexStack.Add(nextVertex);
									
									parent[nextVertex] = currentVertex;
								}
								else
								{
									// non-tree edge
								}
							}
						}
					}
					
					// DFS is finished. Now create the auxiliary graph h
					// Add all the tree edges as vertices in h
					SimpleGraph h = new SimpleGraph();
					
					h.addAllVertices(treeEdges);
					
					visitedVertices.Clear();
					
					CSGraphT.SupportClass.SetSupport connected = new CSGraphT.SupportClass.HashSetSupport();
					
					for (System.Collections.IEnumerator it = dfsVertices.GetEnumerator(); it.MoveNext(); )
					{
						System.Object v = it.Current;
						
						visitedVertices.Add(v);
						
						// find all adjacent non-tree edges
						for (System.Collections.IEnumerator adjacentEdges = subgraph.edgesOf(v).GetEnumerator(); adjacentEdges.MoveNext();)
						{
							Edge l = (Edge)adjacentEdges.Current;
							if (!treeEdges.Contains(l))
							{
								h.addVertex(l);
								System.Object u = l.oppositeVertex(v);
								
								// we need to check if (u,v) is a back-edge
								if (!visitedVertices.Contains(u))
								{
									while (u != v)
									{
										System.Object pu = parent[u];
										Edge f = subgraph.getEdge(u, pu);
										
										h.addEdge(f, l);
										
										if (!connected.Contains(f))
										{
											connected.Add(f);
											u = pu;
										}
										else
										{
											u = v;
										}
									}
								}
							}
						}
					}
					
					ConnectivityInspector connectivityInspector = new ConnectivityInspector(h);
					
					biconnectedSets_Renamed_Field.Add(connectivityInspector.connectedSets());
				}
			}
			
			return biconnectedSets_Renamed_Field;
		}