% -*- texinfo -*-
% @deftypefn{Function File} {@var{x} =} smwsolve (@var{a}, @var{u}, @var{v}, @var{b})
% @deftypefnx{Function File} {} smwsolve (@var{solver}, @var{u}, @var{v}, @var{b})
% Solves the square system @code{(A + U*V')*X == B}, where @var{u} and @var{v} are
% matrices with several columns, using the Sherman-Morrison-Woodbury formula,
% so that a system with @var{a} as left-hand side is actually solved. This is
% especially advantageous if @var{a} is diagonal, sparse, triangular or
% positive definite.
% @var{a} can be sparse or full, the other matrices are expected to be full.
% Instead of a matrix @var{a}, a user may alternatively provide a function
% @var{solver} that performs the left division operation.
% @end deftypefn