Merge pull request #53 from kalmarek/update/StarAlgebras-0.2

update StarAlgebras to v0.2.0

Project.toml

@@ -1,7 +1,7 @@
 name = "SymbolicWedderburn"
 uuid = "858aa9a9-4c7c-4c62-b466-2421203962a2"
 authors = ["Marek Kaluba <[email protected]>", "tweisser <[email protected]>"]
-version = "0.3.2"
+version = "0.3.3"
 
 [deps]
 Cyclotomics = "da8f5974-afbb-4dc8-91d8-516d5257c83b"
@@ -17,7 +17,7 @@ Cyclotomics = "0.3"
 GroupsCore = "0.4"
 PermutationGroups = "0.3.1"
 Primes = "0.4, 0.5"
-StarAlgebras = "0.1.5"
+StarAlgebras = "0.2"
 julia = "1.6"
 
 [extras]

src/minimal_projections.jl

@@ -1,4 +1,17 @@
+function _group_algebra(G::Group)
+    @assert isfinite(G)
+    b = StarAlgebras.Basis{UInt16}(vec(collect(G)))
+    RG = if order(Int, G) <= (typemax(UInt16) >> 1) # 32767
+        StarAlgebra(G, b, (length(b), length(b)); precompute = true)
+        # cache is about ~ 1Gb
+    else
+        StarAlgebra(G, b)
+    end
+    return RG
+end
+
 Base.parent(A::StarAlgebra{<:Group}) = A.object
+StarAlgebras.star(g::GroupElement) = inv(g)
 
 function StarAlgebras.AlgebraElement(
     χ::AbstractClassFunction,
@@ -98,7 +111,7 @@ function minimal_rank_projection(χ::Character, RG::StarAlgebra{<:Group};
     subgroups=CyclicSubgroups(parent(RG), min_order=2),
     iters=3
 )
-    degree(χ) == 1 && return one(RG, Rational{Int}), true
+    degree(χ) == 1 && return one(Rational{Int}, RG), true
 
     for H in subgroups
         µ = _small_idem(RG, H)
@@ -115,7 +128,7 @@ function minimal_rank_projection(χ::Character, RG::StarAlgebra{<:Group};
         push!(_id, deepcopy(subgroups))
     end
     @debug "Could not find minimal projection for $χ"
-    return one(RG, Rational{Int}), false
+    return one(Rational{Int}, RG), false
 end
 
 function minimal_projection_system(

src/sa_basis.jl

@@ -187,18 +187,6 @@ function _constituents_decomposition(ψ::Character, tbl::CharacterTable)
     return irr[present_irreps], multiplicities[present_irreps]
 end
 
-function _group_algebra(G::Group)
-    @assert isfinite(G)
-    b = StarAlgebras.Basis{UInt16}(vec(collect(G)))
-    RG = if order(Int, G) <= (typemax(UInt16)>>2)
-        StarAlgebra(G, b, (length(b), length(b)), precompute=true)
-        # cache is about ~ 1Gb
-    else
-        StarAlgebra(G, b)
-    end
-    return RG
-end
-
 function _symmetry_adapted_basis(
     T::Type,
     irr::AbstractVector{<:Character},