 	TRY( PetscRandomSetFromOptions(rnd) );
 	TRY( PetscRandomSetSeed(rnd,13) );
 
+	Vec layout;
+	int Tlocal;
 	GeneralVector<PetscVector> *x;
 
 	int Ti, Ki;

			for(j=0;j<n;j++){
				/* --- SET RANDOM VALUES TO GLOBAL VECTOR --- */
				TRY( VecSetRandom(x_global, rnd) );
		
				/* --- GET LOCAL VECTOR --- */
				TRY( VecGetLocalVector(x_global,x_local) );

				x = new GeneralVector<PetscVector>(x_local);

				/* --- COMPUTE PROJECTION --- */
				timer1.start();
					feasibleset->project(*x);
				timer1.stop();
	
				/* --- RESTORE GLOBAL VECTOR --- */
				TRY( VecRestoreLocalVector(x_global,x_local) );
			}

			/* --- PRINT INFO ABOUT TIMERS --- */
			coutMaster << "K = "<< std::setw(4) << K << ", T = " << std::setw(9) << T << ", time = " << std::setw(10) << timer1.get_value_sum() << std::endl;

			#ifdef USE_CUDA
				oss_print_to_log << "TIME_PETSCVECSEQ|";
			#else
				oss_print_to_log << "TIME_PETSCVECSEQCUDA|";
			#endif
			oss_print_to_log  << K << "|" << T << "|" << timer1.get_value_sum();



			/* destroy feasible set */
			free(feasibleset);
	
			/* destroy used vectors */
			TRY( VecDestroy(&x_global) );
