idamnjanovic@8:   function y = SL_A(A, mode, m, n, x, I, dim)
idamnjanovic@1:   % Ivan Damnjanovic 2009
idamnjanovic@1:   % This is auxilary function to allow implicit matrices from SPARCO 
idamnjanovic@1:   % to be used with SparsLab solvers
idamnjanovic@1:   
idamnjanovic@1:   global SMALL
idamnjanovic@1:     if (mode == 1)
idamnjanovic@1:     
idamnjanovic@1:             u = zeros(dim, 1);
idamnjanovic@1:             u(I) = x;
idamnjanovic@8:             y = A(u,1);
idamnjanovic@1:            
idamnjanovic@1:         elseif (mode == 2)
idamnjanovic@1:     
idamnjanovic@8:             x2 = A(x,2);
idamnjanovic@1:             y = x2(I);
idamnjanovic@1:            
idamnjanovic@1:      end
idamnjanovic@1:        
idamnjanovic@1:     end