|
- function vrot=givapp(c,s,vin,k)
- % Apply a sequence of k Givens rotations, used within gmres codes
- %
- % C. T. Kelley, July 10, 1994
- %
- % 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);
- % w2=s(i)*vrot(i)+c(i)*vrot(i+1);
- w2=s(i)*vrot(i)+conj(c(i))*vrot(i+1);
- vrot(i:i+1)=[w1,w2];
- end
复制代码 |
评分
-
1
查看全部评分
-
|