Script File
test_unitarity.at
test list [(Param p,bool answer)], test p for unitarity and compare result with answer prints: parameter, computed unitarity, claimed unitarity, number of reducibility points, number of deformations
Definitions in source order
uflag 2 overloads
L18
uflag (int predicted) = stringL21
uflag (bool flag) = stringas_predicted
L23
as_predicted (int predicted, bool unitary) = boolpassedflag 2 overloads
L25
passedflag(int predicted, bool unitary) = stringL27
passedflag (bool flag) = stringtest 2 overloads
L34
test ([(Param,int)] parameters,bool verbose) = (bool,[(Param,KTypePol,int,bool)])test list of parameters for same group for unitarity/non-unitarity parameters=[(param p, int predicted)] predicted takes values 1,0, or -1 which respectively mean parameter is predicted to be unitary / no prediction / predicted non-unitary
L65
test ([(Param,int)] parameters)default quiet
also in this file at line 70; also defined in all_finite_order.at, nilpotent_centralizer.at, sommers.at
print_test
L67
print_test ([(Param,int)] parameters) = voidalso in this file at line 72
test
L70
test([Param] parameters) = (bool,[(Param,KTypePol,int,bool)])also in this file at line 34, line 65; also defined in all_finite_order.at, nilpotent_centralizer.at, sommers.at
print_test
L72
print_test([Param] parameters) = voidalso in this file at line 67
test_one
L74
test_one (Param param,int unitary) = (bool,[(Param,KTypePol,int,bool)])print_test_one
L76
print_test_one (Param param,int unitary) = voidtest_one_unitary
L78
test_one_unitary (Param param) = (bool,[(Param,KTypePol,int,bool)])test_spherical_unipotent
L85
test_spherical_unipotent (RealForm G) = (bool,[(Param,KTypePol,int,bool)])G should be absolutely simple, SL(n,R), Sp(2n,R), S(n,n) or SO(n+1,n) or simple complex, see spherical_unipotent_representations in representations.at
print_test_spherical_unipotent
L87
print_test_spherical_unipotent(RealForm G) = voidtest_Aq
L91
test_Aq ([Param] B) = (bool,[(Param,KTypePol,int,bool)])also in this file at line 132
print_test_Aq
L104
print_test_Aq ([Param] B) = voidalso in this file at line 136
print_test_Upq_Aq 2 overloads
L117
print_test_Upq_Aq ([Param] B) = voidspecial test in U(p,q) case
L128
print_test_Upq_Aq (int p, int q, ratvec gamma) = voidtest_Aq
L132
test_Aq (RealForm G) = (bool,[(Param,KTypePol,int,bool)])also in this file at line 91
print_test_Aq
L136
print_test_Aq (RealForm G) = voidalso in this file at line 104
test_all_real_induced_one_dimensional
L141
test_all_real_induced_one_dimensional(RealForm G) = (bool,[(Param,KTypePol,int,bool)])print_test_all_real_induced_one_dimensional
L144
print_test_all_real_induced_one_dimensional(RealForm G) = voidGenerated from atlas-scripts at commit 7e1b958 (2026-09-17).
Commented-out code, lines 147–185 (38 lines)
comment out fixed case tests, though maybe useful for atlas sanity testing { a series of tests if increasing length } { very fast tests } set unitary_if (bool b) = int: if b then 1 else -1 fi set success ((bool,[Param,KTypePol,int,bool]) (b,)) = bool: b set test1()=test(spherical_unipotent_representations(Sp(4,R))).success set test2()=test(spherical_unipotent_representations(split_form(G2))).success set test3()=test(spherical_unipotent_representations(Sp(6,R))).success set test4()= let G2=quasisplit_form(inner_class(adjoint("G2"),"e")) then p=trivial(G2) then (B,t)=block(p) in test (for p@i in B do (p,unitary_if(i<5 or i=t)) od,true).success { slightly longer, up to a few minutes } set test5()= test(spherical_unipotent_representations (SO(5,4))).success set test6()= let p=trivial(Sp(4,R)) then (B,t)=block(p) in test(for p@i in B do (p,unitary_if(i<7 or i=t)) od,true).success { first 49 of 59 spherical unitary parameters for F4 } set test7()= bool: test(for i:49 do minimal_spherical_principal_series (split_form(F4),F4_spherical_unitary[i]) od).success { first 100 of spherical unitary for E7 requires more memory} set test8()= bool: test(for i:100 do minimal_spherical_principal_series(split_form(E7),E7_spherical_unitary[i]) od).success { trivial of F4 takes up to two hours } set test9()=bool: prints("Testing trivial of F4"); is_unitary(trivial(split_form(F4)))