Concrete cartesian product of zonotopes in standard form

[mforets]

If Z is an n-dimensional zonotope of order 2, we say that it is in standard form if its generators matrix is G = [M D] with M and D n x n matrices and D diagonal.

Create a function cartesian_product_standard(X::Vector{Zonotope}) that returns the n-fold cartesian product in standard form, assuming all the zonotopes in X are in standard form. Bonus points if it can handle zonotopes (X[i]) of different dimensions.

[schillic]

Other possibilities, ordered by effort:

• Add a new argument normal_form::Bool (or even normal_form::FullDiagonal if we want to add more normal forms later) to cartesian_product.
• Add a dispatch on Zonotopes whose generator matrix is given by a certain type (a block matrix [M D] where D is a Diagonal) to cartesian_product.
• Add a new set type ZonotopeMD <: AbstractZonotope which behaves like Zonotope but preserves the normal form.