 		}
 	}
 	coutMaster << " (number of clusters)" << std::endl;
-	coutMaster << "n = " << std::setw(7) << n << " (number of tests)" << std::endl;
+	coutMaster << " n      = " << n << " (number of tests)" << std::endl;
 	coutMaster << "------------------------------------------------------------------" << std::endl;	
 
 

		for(Ti = 0; Ti < T_size; Ti++){
			T = T_list[Ti];

			timer1.restart();

			/* create global vector */
			TRY( VecCreate(PETSC_COMM_WORLD, &x_global) );
			TRY( VecSetSizes(x_global, K*T, PETSC_DETERMINE) );
			TRY( VecSetType(x_global, VECMPI) );
			TRY( VecSetFromOptions(x_global) );
			TRY( VecAssemblyBegin(x_global) );
			TRY( VecAssemblyEnd(x_global) );
		
			/* --- INITIALIZATION --- */

			/* prepare feasible set */
			feasibleset =  new SimplexFeasibleSet_Local(T,K);  	

			#ifdef USE_CUDA
				TRY( VecCreateSeqCUDA(PETSC_COMM_SELF, K*T, &x_local) );
				gpuErrchk(cudaDeviceSynchronize());
			#else
				TRY( VecCreateSeq(PETSC_COMM_SELF, K*T, &x_local)  );
			#endif

			/* log */
			std::ostringstream oss_print_to_log;

	TRY( PetscRandomSetFromOptions(rnd) );
	TRY( PetscRandomSetSeed(rnd,13) );

	Vec layout;
	int Tlocal;
	GeneralVector<PetscVector> *x;

	int Ti, Ki;
