
% print initial input info

fprintf('\n............................................................. \n') 

fprintf('\nTest matrix = %3i, %40s \n',testprob,testmatname)
fprintf('\nn = %5i,  nb = %3i, maxit = %3i, maxeig = %5i\n',n,nb,maxit,maxeig)
fprintf('\nanorm = %10.5e, nonnormality = %10.5e\n',anorm0,nonnormal)
%if matbal == 1 
%   fprintf('\nBalanced matrix norm = %10.5e, nonnormality = %10.5e\n',...
%              anorm,nonnormal_b)
%end 
   
fprintf('\nSpectral transformation option = %2i \n',precond ) 
if precond >= 1, 
   fprintf('shift alpha = %10.5e \n',alpha) 
end 
if precond == 1, 
   fprintf('condition number and norm of shifted matrix = %10.5e, %10.5e\n',...
	    c,fanorm )
end 

if (fulldual == 0 & semidual == 0) & ( group == 0 & treatbd == 0 ), 
   fprintf('\n... Level 1 implementation ... \n') 
elseif fulldual == 1 & ( group == 0 & treatbd == 0 ), ,  
   fprintf('\n... Level 2 implementation ... \n') 
elseif semidual == 1 & ( group == 0 & treatbd == 0 ),  
   fprintf('\n... Leval 3 implementation ... \n') 
elseif (fulldual == 1 | semidual == 1 ) & (group == 1 | treatbd == 1), 
   fprintf('\n... Level 4 implementation ... \n') 
else
   fprintf('\n... Undefined implementation level, check logical flags ... \n')
   break; 
end

fprintf('[fulldual,semidual,group,treatbd] = [%3i,%3i,%3i,%3i] \n', ...
          fulldual,semidual,group,treatbd )
