Open cortner opened 2 years ago
To reproduce:
RnYlm = ACE.Utils.RnYlm_1pbasis(; maxdeg = 6 ) Zk = ACE.Categorical1pBasis([:a, ]; varsym = :z, idxsym = :k, label = "Zk") B1p = RnYlm * Zk Bsel = ACE.SimpleSparseBasis(3, 6); basis = ACE.SymmetricBasis(ACE.EuclideanMatrix(), B1p, Bsel) length(basis)
To reproduce: