Documentation contents

atlas Library · Script File

nci_nilrad_roots.at

Source
atlas-scripts/nci_nilrad_roots.at (66 lines)
Definitions
7
Loads
Loaded by
lusztig_cells.at tensor_product.at

Definitions in source order

all_sums

L13all_sums(mat roots) = [vec]
given a set S of N ratvecs all of length n
returns all 2^N sums of vectors in subsets of S
caution: if S is empty this returns [ratvec:[]]
  i.e. the list of 1 ratvec of length 0
  (rather than [ratvec]:[] the empty list of ratvecs, which has length 0)
  this is appropriate since the expected length of the list is 2^N=2^0=1
otherwise =[null(n),...]

sums_of_size_k

L19sums_of_size_k(mat roots,int k) = [vec]
given a set S of N ratvecs all of length n
return list of sums in all subsets of size k

makeS

L23makeS(KGBElt x_G,Parabolic P) = mat
like makeS in induction.at, except only keep the roots in u

also defined in parabolics.at, induction.at

nci_nilrad_roots_restricted_to_H_theta

L41nci_nilrad_roots_restricted_to_H_theta (KGBElt x_L,Parabolic P) = mat
Given a theta-stable parabolic P and x_G\in KGB(G)
 return a set of nci roots of H, whose restrictions to H^{theta_x} are
 the weights of H^theta on u\cap s
 each nci root contributes itself
 each pair of complex roots for G contributes one of them
Note: our application will be to L_0\subset L\subset G
 q=l+u is the parabolic, we're interested in tensoring L_0\cap K-types
 for this command all that matters is \theta_x restricted to H
 so L_0 plays no role
Dangerous bend: this is a KGBElt for L (what is needed for a weight function to compute
tensor products) but the roots are for G, so embed_KGB(x_L,G)

also in this file at line 53

bug fix on 3/31/25
L45
old code: for root in makeS(x_G,P) do root od; however if makes(x_G,P) is empty this gives an error
it appears that we can simply use makeS(x_G,P) directly, but these scripts still
have a bug: hodge_branch_std(SL(2,R).trivial,10) gives "index 6 out of range" error
L46

nci_nilrad_roots_restricted_to_H_theta

L53nci_nilrad_roots_restricted_to_H_theta(Parabolic P) = (KGBElt x_L) [vec]
returns a function taking a KGBElt to [ratvec]

also in this file at line 41

sums_nci_nilrad_roots_wedge_k_restricted_to_H_theta

L58sums_nci_nilrad_roots_wedge_k_restricted_to_H_theta (Parabolic P,int k) = (KGBElt x_L) [vec]
first compute the previous function to give a set of n roots
run over the binom(n,k) subsets of size k of this set, and return their sums

sums_nci_nilrad_roots_restricted_to_H_theta

L64sums_nci_nilrad_roots_restricted_to_H_theta (Parabolic P) = (KGBElt x_L)[vec]
first compute the ante-previous function to give a set of n roots
run over the 2^n subsets of this set, and return their sums

Generated from atlas-scripts at commit 648506c (2026-10-13).