Gmres算法解线性方程组

解大规模线性方程组的预条件Gmres方法，系数矩阵可以非对称正定

function vrot=givapp（csvink）
%givapp   apply a sequence of Givens rotations
%   input
%          c s      vectors of length k-1 defining rotations
%          vin       vector of length k to which rotations are applied
%          k         k-1 = number of rotations
%   output
%          tranformed vector after rotations are applied
% called by gmres_r
%   IFISS function: HCE; 15 March 2005.

% 
%  C. T. Kelley July 10 1994
%Copyright 1994 C. T. Kelley.  
%Reproduced and distributed with permission of the copyright holder.
%
% This code comes with no guarantee or warranty of any kind.
%
%  function vrot=givapp（c s vin k）
%
vrot=vin;
for i=1:k
    w1=c（i）*vrot（i）-s（i）*vrot（i+1）;        % Change on next line 6/3/97
    w2=s（i）*vrot（i）+conj（c（i））*vrot（i+1）;  % w2=s（i）*vrot（i）+c（i）*vrot（i+1）;
    vrot（i:i+1）=[w1w2];
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        863  2011-09-21 14:21  Gmres算法解线性方程组\givapp.m

     文件       5627  2005-04-23 15:07  Gmres算法解线性方程组\gmres_r.m

     文件     195246  2011-09-21 14:01  Gmres算法解线性方程组\一种预条件的再开始的GMRES算法.pdf

     文件         97  2011-09-21 14:20  Gmres算法解线性方程组\说明.txt

     目录          0  2011-09-21 14:22  Gmres算法解线性方程组

----------- ---------  ---------- -----  ----

               201833                    5