--- a/+scheme/Hypsyst3d.m	Thu Nov 03 22:03:09 2016 +0100
+++ b/+scheme/Hypsyst3d.m	Thu Nov 10 20:47:40 2016 +0100
@@ -139,13 +139,14 @@
                 side=max(length(X),length(Y));
                 cols=cols/side;
             end
-            ret=kron(ones(rows,cols),speye(side));
+            ret=cell(rows,cols);
             
             for ii=1:rows
                 for jj=1:cols
-                    ret((ii-1)*side+1:ii*side,(jj-1)*side+1:jj*side)=diag(matVec(ii,(jj-1)*side+1:jj*side));
+                    ret{ii,jj}=diag(matVec(ii,(jj-1)*side+1:jj*side));
                 end
             end
+            ret=cell2mat(ret);
         end