+ int dim = d + d*d + d*d*d;
+
+//printf("X: \n");
+//for (int kk=0; kk<d*n; kk++) printf("%f ",X[kk]);
+//printf("\n");
+//printf("Y: \n");
+//for (int kk=0; kk<n; kk++) printf("%f ",Y[kk]);
+//printf("\n");
+//printf("M: \n");
+//for (int kk=0; kk<dim; kk++) printf("%f ",M[kk]);
+//printf("\n");
+
+ // (Re)Initialize W:
+ for (int j=0; j<dim; j++)
+ {
+ for (int k=0; k<dim; k++)
+ W[j*dim+k] = 0.0;
+ }
+