################################################################################
# $ python3 verificar_foton_red_30092026.py
# 2026-09-30 15:14 UTC
################################################################################
=== hodge_diamante.py
H1 rango de d = 2 para k generico (ciclos: 4-2=2): True
H2 k=0.01: componente longitudinal de los ciclos / total = 1.38e-03; contenido A1 = 2.33e-16
H2 k=0.001: componente longitudinal de los ciclos / total = 9.16e-05; contenido A1 = 2.22e-16
H2 k=0.0001: componente longitudinal de los ciclos / total = 2.19e-05; contenido A1 = 1.13e-16
H3 theta_-: vector de enlaces ~ A1 (constante): True | theta_+: vector polar paralelo a k: True
H4 hexagono encontrado: True | circulacion invariante gauge: True | circulacion de un campo generico no nula: True
H5 conteo por celda: 4 enlaces = 1 (theta_-, A1, con gap) + 1 (theta_+, longitudinal, gauge) + 2 (ciclos transversos)
=== hex_curl.py
hexagonos orientados por el sitio A(0): 24
k=0.01: transversos/k^2 = 31.999644, 31.999644; |K d| = 2.2e-19; nulos 0.0e+00 1.1e-19
k=0.01: transversos/k^2 = 31.999556, 31.999556; |K d| = 3.3e-19; nulos -1.2e-19 0.0e+00
k=0.01: transversos/k^2 = 31.999526, 31.999526; |K d| = 6.2e-19; nulos -2.1e-19 4.2e-35
k=0.01: transversos/k^2 = 31.999556, 31.999556; |K d| = 3.3e-19; nulos -7.8e-19 -1.9e-19
k=0.001: transversos/k^2 = 31.999996, 31.999996; |K d| = 3.8e-21; nulos -1.9e-21 1.9e-21
k=0.001: transversos/k^2 = 31.999996, 31.999996; |K d| = 1.7e-21; nulos -5.3e-21 8.5e-22
k=0.001: transversos/k^2 = 31.999995, 31.999995; |K d| = 8.1e-21; nulos 1.9e-37 6.4e-21
k=0.001: transversos/k^2 = 31.999996, 31.999996; |K d| = 8.1e-21; nulos -2.3e-21 4.6e-21
=== hex_area.py
hexagonos distintos por el sitio A(0): 12
sum_hex (A.B)^2 por sitio = B.M.B, M=
 [[21.333333 -0.        0.      ]
 [-0.       21.333333  0.      ]
 [ 0.        0.       21.333333]]
|A|^2 de cada hexagono: [np.float64(5.333333)]
por celda: coef de |B|^2 = 7.111111111111114  =64/9? True
por hexagono (4/celda): coef |B|^2 = 1.7777777777777786 =16/9 -> (8/9) F^2
por sitio (12 hexagonos), coef de F^2 = M/2 = 10.666666666666671 (=32/3); con la normal del corpus Sum p x p (sin mitad, x4): 42.666666666666686 = 128/3, fila hexagono de la tabla de III.4.5: True
=== green_diamante.py
Vc=3.079201, prediccion 4 pi R G -> 3Vc/4 = 2.309401 (=4/sqrt3 2.309401)
R1=6.53 R2=13.06 : 4 pi (G1-G2)/(1/R1-1/R2) = 2.31545
R1=9.80 R2=19.60 : 4 pi (G1-G2)/(1/R1-1/R2) = 2.30947
R1=13.06 R2=26.13 : 4 pi (G1-G2)/(1/R1-1/R2) = 2.30540
R1=16.33 R2=32.66 : 4 pi (G1-G2)/(1/R1-1/R2) = 2.30005

