 
 		/* compute absolute error of computed solution */
 		abserr = mydata.compute_abserr_reconstructed(solution);
-		
-		coutMaster << " - abserr = " << abserr << "\n";
-		mysolver.printtimer(coutMaster);
-		mysolver.printstatus(coutMaster);	
+		abserrs[depth] = abserr;
+		
+		coutMaster << " - abserr = " << abserr << std::endl;
+//		mysolver.printtimer(coutMaster);
+//		mysolver.printstatus(coutMaster);	
 	
 		/* if this solution is better then previous, then store it */
 		if(abserr < abserr_best){

	}

/* 8.) store best solution */
	coutMaster << "--- SAVING OUTPUT ---\n";
	oss << filename_out;
	mydata.saveSignal1D(oss.str(),false);
	oss.str("");

	/* print solution */
	coutMaster << "--- THETA SOLUTION ---\n";
	mydata.print_thetavector(coutMaster);

	/* print timers */
	coutMaster << "--- TIMERS INFO ---\n";
	mysolver.printtimer(coutMaster);

	/* print short info */
	coutMaster << "--- FINAL SOLVER INFO ---\n";
	mysolver.printstatus(coutMaster);

	/* say bye */	
	coutMaster << "- end program\n";

	logging.end();
	Finalize();

	Vec thetavector_best_Vec; /* here we store solution with best abserr value */
	TRYCXX( VecDuplicate(mydata.get_thetavector()->get_vector(),&thetavector_best_Vec) );
	
	/* go throught given list of epssqr */
	for(int depth = 0; depth < epssqr_list.size();depth++){
		epssqr = epssqr_list[depth];
		coutMaster << "--- SOLVING THE PROBLEM with epssqr = " << epssqr << " ---\n";
