SEXT3
δ = 41 m · dα = 1.000000° sur l'angle entre amers · empreinte 80af6ca4e7e86001
Jouer sur les tolérances
Une erreur de carte par amer, une erreur angulaire commune. Baissez le δ d'un seul amer — celui que vous relèveriez au GPS — et regardez X rétrécir. Le calcul refait est certifié, pas approché : comptez quelques secondes par essai.
Le dα minimal pour que les mesures soient cohérentes vaut au moins 0.318°. Vous avez déclaré 1.000°, soit 3.1× ce minimum.
7 amers lus dans « SEXT3.csv » (2 colonnes : x, y)
7 mesures lues dans le second fichier (i, j, angle en degres) sur 21 paires possibles.
Les amers
7 points, a 1455 a 3938 m de l'origine.
X
environ 193 m de large — soit 20 fois moins que la portee des amers. Le cercle en pointille le situe ; le detail est ci-dessous.
Les traits pales vont de l'origine vers chaque amer. L'observateur n'est pas a l'origine : il est quelque part dans X, et c'est tout ce qu'on sait de lui.
X
23364.8 ± 5.5 m²
X∩
28440 ± 38 m²
X ⊊ X∩ : la clôture par Floyd-Warshall améliore au moins une borne, donc X est strictement plus petit.
13 sommets certifiés.
Les sommets
| # | x (m) | y (m) | contraintes |
|---|
| 1 | 38.311 | 98.795 | (1,5) ∩ (1,6) |
| 2 | 55.569 | 104.437 | (1,5) ∩ (7,5) |
| 3 | -11.048 | 70.289 | (1,6) ∩ (1,7) |
| 4 | -53.540 | 12.558 | (1,7) ∩ (3,7) |
| 5 | 112.013 | 59.617 | (2,3) ∩ (7,3) |
| 6 | 111.830 | 60.165 | (2,3) ∩ (7,4) |
| 7 | -49.963 | -77.715 | (3,7) ∩ (6,7) |
| 8 | 32.549 | -82.017 | (4,1) ∩ (6,1) |
| 9 | -4.718 | -88.933 | (4,1) ∩ (6,7) |
| 10 | 81.703 | -52.446 | (6,1) ∩ (7,1) |
| 11 | 86.233 | -45.916 | (6,3) ∩ (7,1) |
| 12 | 114.231 | 24.865 | (6,3) ∩ (7,3) |
| 13 | 111.268 | 60.720 | (7,4) ∩ (7,5) |
Chaque sommet est un enclos certifié ; la largeur du plus large vaut 1.4e-14 m.
Ce que cela a coûté
| lecture du tableur | 0.001 s |
|---|
| préparation du jeu | 0.007 s |
|---|
| calcul de X∩ et X | 4.989 s |
|---|
| aire et tracé | 1.009 s |
|---|
| total | 6.181 s |
|---|
Le calcul et l'analyse sont deux choses : la première est l'objet de l'étude, la seconde un outil pour la regarder. Leur coût ne dépend pas des mêmes choses — le calcul du jeu, l'analyse de la précision demandée.
Sur quoi
machine : Linux aarch64 -- Raspberry Pi 4 Model B Rev 1.1
4 coeur(s) a 1.50 GHz, 3.7 Go, 4 fil(s) demande(s)
calcul : intervalles fi64, mantisse hors solveur 256 bits
bati : Sep 10 2026, gcc 14.2.0
arrondi dirige disponible, socle : pthreads oui, rsgc non, flint non
gmp 6.3.0, mpfr 4.2.2, mpfi 1.5.4
sources : empreinte cbb2c1eddba1Les temps ci-dessus valent pour cette machine et pour elle seule. Ce qui se transporte d'une machine à l'autre, ce sont les rapports — pas les secondes.
Le jeu exact engendré (ce que lit le calcul)
# ==============================================================================
# Jeu SEXT3 -- engendre par preparer.py, NE PAS EDITER A LA MAIN
# 7 amers, portees 1455-3938 m depuis l'origine
# delta = 41 m (40 m de carte + 1 m d'arrondi des coordonnees)
# d_alpha = 1.000000 deg SUR L'ANGLE ENTRE AMERS
# dont 2.22e-04 arcsec de transcription, majores et absorbes
# 7 paires mesurees sur 21 possibles
# ==============================================================================
const NOM = "SEXT3"
const EMPREINTE = "80af6ca4e7e86001" # signature des donnees, recopiee dans le resultat
const AMERS = Tuple{Rational{BigInt},Rational{BigInt}}[
(1407//1, 840//1), # 1 azimut 59.17 deg
(3269//1, 2196//1), # 2 azimut 56.11 deg
(120//1, 1450//1), # 3 azimut 4.75 deg
(-2342//1, 1046//1), # 4 azimut 294.07 deg
(-3040//1, 403//1), # 5 azimut 277.54 deg
(-2368//1, -1344//1), # 6 azimut 240.42 deg
(1060//1, -3175//1), # 7 azimut 161.54 deg
]
const DELTA = Rational{BigInt}(41)
const T_D = Rational{BigInt}(722, 82733) # tan(d_alpha/2)
const MESURES = Dict{Tuple{Int,Int},Tuple{Rational{BigInt},Int}}(
(1,2) => (1362//44221, 0),
(1,7) => (-85466//67531, 0),
(2,3) => (41890//85361, 0),
(3,4) => (65723//93441, 0),
(4,5) => (14342//96949, 0),
(5,6) => (13948//41075, -1),
(6,7) => (63552//77831, 0) )
Le journal du calcul
machine : Linux aarch64 -- Raspberry Pi 4 Model B Rev 1.1
4 coeur(s) a 1.50 GHz, 3.7 Go, 4 fil(s) demande(s)
calcul : intervalles fi64, mantisse hors solveur 256 bits
bati : Sep 10 2026, gcc 14.2.0
arrondi dirige disponible, socle : pthreads oui, rsgc non, flint non
gmp 6.3.0, mpfr 4.2.2, mpfi 1.5.4
========================================
SEXT3 (7 amers)
========================================
minorant : les cycles imposent d_alpha >= 0.3184 deg (declare 1.0000)
pire cycle, 7 aretes : 4 3 2 1 7 6 5 4
defaut de fermeture 2.228601 deg pour 7.000000 deg de tolerance cumulee
exact sur TOUS les cycles : 0.3184
marge : 3.14x ce que les cycles exigent
7 amers, 14 contraintes orientees
elagage sur 10 boites : 6 paires sur 91
X∩ : 6 sommets en 0.126 s (elagage 0.007, solveur 0.108, filtres 0.009)
boite de depart reprise de X∩ : x[-71, 129] y[-104, 121]
germe repris de X∩ : 6 contrainte(s)
critere : la cloture apporte 0 amelioree(s) et 28 fabriquee(s) sur 42 -> X inclus strictement dans X∩, il faut calculer X
Floyd-Warshall certifie sur le graphe des mesures (7 amers)...
coherence de type I : OK
la cloture apporte : 0 borne(s) amelioree(s) et 28 fabriquee(s) sur 42 -> X inclus dans Xinter
--- iteration 1 : x[-71, 129] y[-104, 121]
origine des angles : secteur libre 95.7 deg
contraintes conservees : 42 sur 42
ensemble actif, tour 1 : |A| = 6, 6 sommet(s), 13 violee(s) (plafonne)
ensemble actif, tour 2 : |A| = 10, 10 sommet(s), 6 violee(s) (plafonne)
ensemble actif, tour 3 : |A| = 14, 12 sommet(s), 1 violee(s)
ensemble actif, tour 4 : |A| = 15, 13 sommet(s), 0 violee(s)
sommets certifies : 13
certificat 1 : 27 prouvee(s) redondante(s), 0 douteuse(s), 205 cellule(s)
--- iteration 2 : x[-54, 115] y[-89, 105]
origine des angles : secteur libre 96.9 deg
contraintes conservees : 42 sur 42
ensemble actif, tour 1 : |A| = 13, 13 sommet(s), 0 violee(s)
sommets certifies : 13
certificat 1 : 29 prouvee(s) redondante(s), 0 douteuse(s), 281 cellule(s)
point fixe atteint.
statut : ok
boite : x[-54.000000, 115.000000] y[-89.000000, 105.000000]
sommets : 13 contraintes actives : 13
bord : 13 arete(s), cycle unique
certificat : 1 passe(s), 13 actives sur 42 gardees, 486 cellule(s)
memoire : pic 39.6 Mo
temps : X∩ 0.126 s X 4.736 s aire 0.000 s bord 0.000 s calcul 4.985 s mural 4.985 s
solveur : direct=392 erreur=0 asymptote=0 non_gen=0 plusieurs=0 coords=0
dont decomposition 10.33 s, solve_ratparam_real 1.84 s
ecrit : /tmp/geoloc-z_hidmjw/SEXT3.resultat.json
========================================
SEXT3_courbes (7 amers) -- analyse de /tmp/geoloc-z_hidmjw/SEXT3.resultat.json
========================================
aire de X∩ : [28402, 28479] m2 (0.268 %, niveaux 2 ; 24 prouvee(s), 0 mesuree(s), 0 majoree(s))
bord de X∩ : 24 points (fleche 0.98 px a 2160, niveaux 2)
aire de X : [23359, 23371] m2 (0.0473 %, niveaux 2 ; 52 prouvee(s), 0 mesuree(s), 0 majoree(s))
bord de X : 52 points (fleche 0.23 px a 2160, niveaux 2)
temps : aire 0.928 s bord 0.074 s
ecrit : /tmp/geoloc-z_hidmjw/SEXT3.analyse.json
← un autre tableur
Sommets certifiés par arithmétique exacte sur ℤ ; aire garantie — la flèche de chaque sous-arc est majorée, pas échantillonnée. Le tracé est la polyligne du calcul lui-même, point par point.