Ticket #25097: error.txt

File error.txt, 5.8 KB (added by jipilab, 2 years ago)

An error in convex hull?

Line 
1sage: x = polygen(ZZ)
2sage: p = x^2 - 2
3sage: nf = NumberField(p,'a',embedding=1.41)
4sage: a = nf.gen()
5sage: starting_pts = [(323701667143080601/10148789250000*a + 922665678365641037/20297578500000, -23535799/6000*a - 32415101/6000, -45329/48*a - 3997/3, -955481/600*a - 5550041/2400, -887399/400*a - 1258953/400,
6....: -56539/200*a - 153867/400, -3263813/3000*a - 18374659/12000, -22095/16*a - 7119/4, -3410549/1500*a - 4949947/1500),
7....:  (-7453/8*a - 23983/16, 3293570551999745141081/685546755559212000*a + 1433772720973225868907/228515585186404000, -45329/48*a - 3997/3, -955481/600*a - 5550041/2400, -887399/400*a - 1258953/400, -56539/200*
8....: a - 153867/400, -3263813/3000*a - 18374659/12000, -22095/16*a - 7119/4, -3410549/1500*a - 4949947/1500),
9....:  (-7453/8*a - 23983/16, -23535799/6000*a - 32415101/6000, 1309965730774718773/29367267000000*a + 10043378812771733/305909031250, -955481/600*a - 5550041/2400, -887399/400*a - 1258953/400, -56539/200*a - 15
10....: 3867/400, -3263813/3000*a - 18374659/12000, -22095/16*a - 7119/4, -3410549/1500*a - 4949947/1500),
11....:  (-7453/8*a - 23983/16, -23535799/6000*a - 32415101/6000, -45329/48*a - 3997/3, 434077100740889820793/23831234671935000*a + 867256748852289707243/31774979562580000, -887399/400*a - 1258953/400, -56539/200*
12....: a - 153867/400, -3263813/3000*a - 18374659/12000, -22095/16*a - 7119/4, -3410549/1500*a - 4949947/1500),
13....:  (-7453/8*a - 23983/16, -23535799/6000*a - 32415101/6000, -45329/48*a - 3997/3, -955481/600*a - 5550041/2400, 16490965235981357383/1801563445260000*a + 13582081793730863569/600521148420000, -56539/200*a -
14....: 153867/400, -3263813/3000*a - 18374659/12000, -22095/16*a - 7119/4, -3410549/1500*a - 4949947/1500),
15....:  (-7453/8*a - 23983/16, -23535799/6000*a - 32415101/6000, -45329/48*a - 3997/3, -955481/600*a - 5550041/2400, -887399/400*a - 1258953/400, 2940985018828950787/24405614730000*a + 889786795015428237/54234699
16....: 40000, -3263813/3000*a - 18374659/12000, -22095/16*a - 7119/4, -3410549/1500*a - 4949947/1500),
17....:  (-7453/8*a - 23983/16, -23535799/6000*a - 32415101/6000, -45329/48*a - 3997/3, -955481/600*a - 5550041/2400, -887399/400*a - 1258953/400, -56539/200*a - 153867/400, 330003969652128971987/9753348602181000*
18....: a + 90392604728686683819/2600892960581600, -22095/16*a - 7119/4, -3410549/1500*a - 4949947/1500),
19....:  (-7453/8*a - 23983/16, -23535799/6000*a - 32415101/6000, -45329/48*a - 3997/3, -955481/600*a - 5550041/2400, -887399/400*a - 1258953/400, -56539/200*a - 153867/400, -3263813/3000*a - 18374659/12000, 19633
20....: 42831078699891/82747737000000*a + 27448751841514081/827477370000, -3410549/1500*a - 4949947/1500),
21....:  (-7453/8*a - 23983/16, -23535799/6000*a - 32415101/6000, -45329/48*a - 3997/3, -955481/600*a - 5550041/2400, -887399/400*a - 1258953/400, -56539/200*a - 153867/400, -3263813/3000*a - 18374659/12000, -2209
22....: 5/16*a - 7119/4, 689146812407956856171/59437744320336000*a + 345796586965483560513/19812581440112000),
23....:  (-7453/8*a - 23983/16, -23535799/6000*a - 32415101/6000, -45329/48*a - 3997/3, -955481/600*a - 5550041/2400, -887399/400*a - 1258953/400, -56539/200*a - 153867/400, -3263813/3000*a - 18374659/12000, -2209
24....: 5/16*a - 7119/4, -3410549/1500*a - 4949947/1500),
25....:  (294961896928506827457/28712270399494400*a + 445954595532213608383/35890337999368000, -236976304781979374547/17945168999684000*a - 730083098261334939123/35890337999368000, 1134734521329691699939/717806759
26....: 98736000*a + 2464605086307987653051/107671013998104000, -11165613476576973686531/215342027996208000*a - 8073167361124822614763/107671013998104000, 459940656484264825093/17945168999684000*a + 59823768298355
27....: 325401/1435613519974720, 2557328671787804865251/179451689996840000*a + 3538533472423884658117/179451689996840000, -3263813/3000*a - 18374659/12000, 8625430985444658216919/143561351997472000*a + 15219674536
28....: 42113815467/17945168999684000, -5804910538578789232111/430684055992416000*a - 1854566448959339106641/107671013998104000),
29....:  (-23515003542967663451/71009257920000*a - 66021389080738078379/142018515840000, 37290193825664968747/71009257920000*a + 158453024240177007223/213027773760000, -12545264994837002809/23669752640000*a - 1598
30....: 87734377789172183/213027773760000, 159490830379271832983/71009257920000*a + 112624106226759633881/35504628960000, -65966046062823533579/71009257920000*a - 23504722866986792051/17752314480000, -159704280184
31....: 721404183/355046289600000*a - 37591499557356466427/59174381600000, 17320215516885696427/71009257920000*a + 1506049413178283461/4438078620000, -77042490328726964549/28403703168000*a - 34054562548554137183/8
32....: 876157240000, 7762465280221581577/23669752640000*a + 21779952039122239353/47339505280000),
33....:  (-1043331582411795643/5704586385000*a - 9342806296403266933/36509352864000, 20837210017088452831/45636691080000*a + 177827269372280787889/273820146480000, -21059323214510059151/45636691080000*a - 17967140
34....: 7031438306969/273820146480000, 26528091937986083147/15212230360000*a + 112767960679459799867/45636691080000, -9308330035308535393/9127338216000*a - 16644181784649694813/11409172770000, -7972145712126556824
35....: 1/228183455400000*a - 37653757022179728439/76061151800000, -3263813/3000*a - 18374659/12000, -272439483578244786299/182546764320000*a - 96346445454224815529/45636691080000, 256261633215088277/1426146596250
36....: *a + 3055867161280206487/12169784288000)]
37sage: len(starting_pts)
3813
39sage: the_poly = Polyhedron(vertices=starting_pts,backend='normaliz')
40sage: the_poly.n_vertices()
4112
42sage: l_v = [v.vector() for v in the_poly.vertices()]
43sage: for i in range(len(starting_pts)):
44....:     print the_poly.contains(starting_pts[i])
45....:     if starting_pts[i] not in l_v:
46....:         print i
47....:         
48True
490
50False
511
52True
532
54False
553
56False
574
58True
595
60False
616
62True
637
64False
658
66True
679
68False
6910
70False
7111
72False
7312