 	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";


		/* compute absolute error of computed solution */
		abserr = mydata.compute_abserr_reconstructed(solution);
		
		coutMaster << " - abserr = " << abserr << "\n";
		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; /* actual error */
	double abserrs[epssqr_list.size()]; /* all errors */
	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 */
