idamnjanovic@1:   function y = SL_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@1:             y = SMALL.Problem.A(u,1);
idamnjanovic@1:            
idamnjanovic@1:         elseif (mode == 2)
idamnjanovic@1:     
idamnjanovic@1:             x2 = SMALL.Problem.A(x,2);
idamnjanovic@1:             y = x2(I);
idamnjanovic@1:            
idamnjanovic@1:      end
idamnjanovic@1:        
idamnjanovic@1:     end