 
 		/* 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){

	
/* 6.) solve the problem with epssqrs and remember best solution */
	double epssqr, epssqr_best;
	double abserr;
	double abserr_best = std::numeric_limits<double>::max(); /* the error of best solution */

	Vec gammavector_best_Vec; /* here we store solution with best abserr value */

	}

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

	/* print solution */
	coutMaster << "--- THETA SOLUTION ---" << std::endl;
	mydata.print_thetavector(coutMaster);

	/* print timers */
	coutMaster << "--- TIMERS INFO ---" << std::endl;
	mysolver.printtimer(coutMaster);

	/* print short info */
	coutMaster << "--- FINAL SOLVER INFO ---" << std::endl;
	mysolver.printstatus(coutMaster);

	/* print absolute errors */
	coutMaster << "--- ABSOLUTE ERRORS: ---" << std::endl;
	coutMaster << "epssqr: " << print_vector(epssqr_list) << std::endl;
	coutMaster << "abserrs: " << print_array(abserrs, epssqr_list.size()) << std::endl;

	/* say bye */	
	coutMaster << "- end program" << std::endl;

	logging.end();
	Finalize();
