Opened 6 years ago
Closed 6 years ago
#14853 closed defect (fixed)
RDF vertices of a graph are transformed into consecutive integers
Reported by: | tmonteil | Owned by: | jason, ncohen, rlm |
---|---|---|---|
Priority: | major | Milestone: | sage-5.12 |
Component: | graph theory | Keywords: | |
Cc: | ncohen | Merged in: | sage-5.12.beta3 |
Authors: | Robert Miller | Reviewers: | Nathann Cohen |
Report Upstream: | N/A | Work issues: | |
Branch: | Commit: | ||
Dependencies: | Stopgaps: |
Description
As reported in the comments of this ask answer:
sage: A=Set([RDF.random_element(min=0,max=10) for k in range(10)]) ; A {6.42320967152, 1.77698693175, 2.95922392964, 9.50745089775, 4.60546838289, 3.67300191731, 5.21254750195, 5.90579400282, 7.55309974188, 0.442799267782} sage: G = Graph() sage: G.add_vertices(A) sage: G.vertices() [0, 1, 2, 3, 4, 5, 6, 7, 9]
Attachments (1)
Change History (11)
comment:1 Changed 6 years ago by
comment:2 Changed 6 years ago by
I suspect the error is in sage/graphs/base/c_graph.pyx
:
cdef int get_vertex(object u, dict vertex_ints, dict vertex_labels, CGraph G) except ? -2: """ Returns an int representing the arbitrary hashable vertex u (whether or not u is actually in the graph), or -1 if a new association must be made for u to be a vertex. TESTS: We check that the bug described in #8406 is gone:: sage: G = Graph() sage: R.<a> = GF(3**3) sage: S.<x> = R[] sage: G.add_vertex(a**2) sage: G.add_vertex(x) sage: G.vertices() [a^2, x] And that the bug described in #9610 is gone:: sage: n = 20 sage: k = 3 sage: g = DiGraph() sage: g.add_edges( (i,Mod(i+j,n)) for i in range(n) for j in range(1,k+1) ) sage: g.vertices() [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19] sage: g.strongly_connected_components() [[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]] """ cdef int u_int if u in vertex_ints: return vertex_ints[u] try: u_int = u except StandardError: return -1 if u_int < 0 or u_int >= G.active_vertices.size or u_int in vertex_labels: return -1 return u_int
This "try/except" has apparently been written in the hope that if u
is anything other than an integer, there will be an exception; but that reasoning fails because rationals and reals are silently converted into integers. Am I right? I'm unable to import methods from cython files into the sage console, so I can't test my conjectures...
comment:3 Changed 6 years ago by
See also #14967
comment:4 Changed 6 years ago by
- Status changed from new to needs_review
comment:5 Changed 6 years ago by
for patchbot:
Apply trac_14853.patch
comment:6 Changed 6 years ago by
- Status changed from needs_review to needs_work
OKayyyyyyyyyyy... Could you replace the references to the track ticket with :trac:14967
and :trac:14853
?
:-)`
Short of this I think it's good to go :-)
Nathann
comment:7 Changed 6 years ago by
OOps. :trac:`14967`
and :trac:`14853`
sorry !
Nathann
Changed 6 years ago by
comment:9 Changed 6 years ago by
- Reviewers set to Nathann Cohen
- Status changed from needs_review to positive_review
Well.. All tests pass :-)
Thaaaaaaaaaaaaaaanks !
Nathann
comment:10 Changed 6 years ago by
- Merged in set to sage-5.12.beta3
- Resolution set to fixed
- Status changed from positive_review to closed
Unfortunately this is not limited to RDF:
My attempts at debugging this stop here:
What is that mysterious
check_vertex
method and where can I find it?EDIT: Oh, I found it. I have been grepping for "def check_vertex" but it is "cdef int check_vertex". Let me see.