 	int xdim = 3; /* data dimension */
 
 	/* solution - for generating the problem */
-	int solution_K = 4;
+	int solution_K = 3;
 	int solution_xmem = 0;
 	
-	double solution_theta[12] = {
+	double solution_theta[9] = {
 		 0.0, 0.0, 0.0,		/* K=1,n=1,2,3:  mu */
 		 1.0, 0.0, 0.0,		/* K=2,n=1,2,3 */
-		 0.0, 1.0, 0.0,		/* K=3,n=1,2,3 */
-		 0.5, 0.5, 0.2		/* K=4,n=1,2,3 */
+		 0.0, 1.0, 0.0		/* K=3,n=1,2,3 */
 	};
 
-	double solution_xstart[4] = {
-		0.0, 0.0, 0.0, 0.0	/* n=1,2,3,4 */
+	double solution_xstart[3] = {
+		0.0, 0.0, 0.0	/* n=1,2,3 */
 	};
 
-	double noise_covariance[12] = {
+	double noise_covariance[9] = {
 		0.05, 0.05, 0.01,  
 		0.05, 0.05, 0.01,  
-		0.05, 0.05, 0.01,  
-		0.05,  0.05, 0.05
+		0.05, 0.05, 0.01  
 	};
 
 	/* model parameters */

	int *xmem = new int(num);

	int i;
	for(i=1; i<num; i++ ){
		K[i] = Kforall;
		epssqr[i] = epssqrforall;
		xmem[i] = 0;


	/* save results into VTK file */
	coutMaster << "--- SAVING VTK ---" << std::endl;
	example::KMeans3D::saveVTK("results/kmeans",".vtk",T,num,K,mydata.get_datavector(),mydata.get_gammavector());
	coutAll.synchronize();
	
	/* print timers */
