Add block_diagonal function (#80)

* Add block_diagonal function

* Fix

* Fix format

src/dependence.jl

@@ -247,8 +247,11 @@ function BasisDependence{StaircaseDependence}(
         end
     end
     vars = MP.variables(basis)
-    full_basis =
-        MB.maxdegree_basis(MB.Monomial, vars, MP.maxdegree(basis.monomials))
+    full_basis = MB.maxdegree_basis(
+        MB.FullBasis{MB.Monomial,M}(),
+        vars,
+        MP.maxdegree(basis.monomials),
+    )
     d = StaircaseDependence[]
     # This sieve of [LLR08, Algorithm 1] is a performance improvement but not only.
     # It also ensures that the standard monomials have the "staircase structure".

src/moment_matrix.jl

@@ -112,6 +112,10 @@ struct BlockDiagonalMomentMatrix{T,B<:SA.ExplicitBasis,MT} <:
     blocks::Vector{MomentMatrix{T,B,MT}}
 end
 
+function block_diagonal(blocks::Vector{<:MomentMatrix})
+    return BlockDiagonalMomentMatrix(blocks)
+end
+
 function show_basis_indexed_blocks(io::IO, blocks)
     print(io, " with $(length(blocks)) blocks:")
     nd = ndigits(length(blocks))

test/moment_matrix.jl

@@ -30,7 +30,7 @@ And entries in a 2×2 SymMatrix{$Int}:
         @test nvariables(ν) == 2
     end
     @test_throws ArgumentError moment_matrix(moment_vector([1], [x]), [y])
-    block_ν = BlockDiagonalMomentMatrix([ν1, ν2])
+    block_ν = block_diagonal([ν1, ν2])
     @test block_ν isa BlockDiagonalMomentMatrix{Int,typeof(ν1.basis)}
     @test sprint(show, block_ν) == """
 BlockDiagonalMomentMatrix with 2 blocks: