TMJets should check that the initial set dimension is compatible

[mforets]

cc: @SebastianGuadalupe

[schillic]

Maybe we should check this during the system normalization or in solve! (for all algorithms).