Lexicographically higher multi-indices for high dimension with lp = 2
What would be the lexicographically highest multi-index of m=6
, n=2
, and lp=2.0
?
import minterpy as mp
mi = mp.MultiIndexSet.from_degree(6, 2, 2)
# A couple of last entries
mi.exponents[-5:,:]
# array([[0, 1, 0, 0, 1, 1],
# [1, 1, 0, 0, 1, 1],
# [0, 0, 1, 0, 1, 1],
# [1, 0, 1, 0, 1, 1],
# [0, 1, 1, 0, 1, 1]])
Why doesn't it go on until [0, 0, 0, 0, 0, 2]
?
And in fact for dimension 8:
mi = mp.MultiIndexSet.from_degree(8, 2, 2)
# A couple of last entries
mi.exponents[-5:,:]
# array([[0, 0, 1, 1, 0, 1, 1, 0],
# [0, 0, 0, 0, 1, 1, 1, 0],
# [1, 0, 0, 0, 1, 1, 1, 0],
# [0, 1, 0, 0, 1, 1, 1, 0],
# [0, 0, 1, 0, 1, 1, 1, 0]])
the highest dimension has never been activated (the last column is all zero).
Is this an expected behavior? Or perhaps there is something amiss with the way exponent matrix is allocated (i.e., the number of expected exponents) for lp other than 1.0?
PS: Dismiss this issue if this is indeed expected (and sorry!).