FINDPIV.M
上传用户:sfyaiting
上传日期:2009-10-25
资源大小:320k
文件大小:0k
源码类别:

GPS编程

开发平台:

Matlab

  1. function [k,p] = findpiv(A,k,p,tol)
  2. %FINDPIV Used by PLU to find a pivot for Gaussian elimination.
  3. % [r,p] = FINDPIV(A(k:m,p:n),k,p,tol) finds the first element in
  4. % the specified submatrix which is larger than tol in absolute value.
  5. % It returns indices r and p so that A(r,p) is the pivot.
  6. [m,n] = size(A);
  7. r = find(abs(A(:))>tol);
  8. if isempty(r)
  9.    return
  10. end
  11. r = r(1);
  12. j = fix((r-1)/m)+1;
  13. p = p+j-1;
  14. k = k+r-(j-1)*m-1;