%
% All vectors are assumed to be *column* vectors.
%
+ % The rather verbose name of this function was chosen to avoid
+ % conflicts with other implementations.
+ %
n = length(x0);
M = eye(n);
% The standard CGM is equivalent to the preconditioned CGM is you
% use the identity matrix as your preconditioner.
- [x, k] = preconditioned_conjugate_gradient_method(A,
- M,
- b,
- x0,
- tolerance,
+ [x, k] = preconditioned_conjugate_gradient_method(A, ...
+ M, ...
+ b, ...
+ x0, ...
+ tolerance, ...
max_iterations);
end