home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 2.ddi / MUTOOLS2.DI$ / CGIVENS.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  463 b   |  26 lines
  1. % function [c,s] = cgivens(x,y)
  2. %
  3. %  produces a Givens rotation such that
  4. %            | c        s |   | x |     | z |
  5. %            |            | * |   |  =  |   |
  6. %            |-conj(s)  c |   | y |     | 0 |
  7. %
  8. %  
  9. function [c,s] = cgivens(x,y)
  10. if abs(x) == 0.0
  11.   c = 0; 
  12.   s = 1; 
  13.   return 
  14. end
  15. if abs(y) == 0.0
  16.   c=1.0; 
  17.   s=0.0; 
  18.   return
  19. end
  20. nrm = norm([x y]);
  21. k1 = x/abs(x);
  22. c = abs(x)/nrm;
  23. s = k1*y'/nrm;
  24. %
  25. % Copyright MUSYN INC 1991,  All Rights Reserved
  26.