 function x = precLU(PRE, rhs) 
%% function rhs = precLU(PRE, rhs) 
%% ILU or other LU-type preconditioning operation
%% PRE = struct for preconditioner
%%-------------------------------------------------
L = PRE.L;
U = PRE.U; 
x = U \ (L \ rhs);

