LinearAlgebra`BLAS`
LinearAlgebra`BLAS`

ROTG

ROTG[a,b,c,s]

computes a Givens rotation {c,s} for given scalars a and b.

Details

  • To use ROTG, you first need to load the BLAS Package using Needs["LinearAlgebra`BLAS`"].
  • The following arguments must be given:
  • ainput/output symbolscalar; the symbol value is modified in place
    binput/output symbolscalar; the symbol value is modified in place
    coutput symbolscalar; the symbol value is modified in place
    soutput symbolscalar; the symbol value is modified in place
  • ROTG[a,b,c,s] computes {c,s} such that {{c,s},{-Conjugate[s],c}.{a,b}=={r,0}, where the signed radius r is given by rSign[b]Norm[{a,b}] if a==0 or a,b with Abs[a]Abs[b], and by rSign[a]Norm[{a,b}] otherwise. In the degenerate case a==b==0, by convention c==1 and s==0.
  • The scalars c and s can be viewed as the cosine and sine, respectively, of the angle of rotation and satisfy the conditions c and c2+Abs[s]21.
  • The scalar a is set to r. The scalar b is not modified if either a or b is complex. Otherwise, b is set to the value z given by:
  • sif Abs[a]>Abs[b]
    1/cif 0<Abs[a]Abs[b]
    1if 0==Abs[a]<Abs[b]
    0if 0==Abs[a]==Abs[b]
  • When given, the value z can be used to reconstruct c and s as follows:
  • c1/z sif Abs[z]>1
    c0s1if z1
    cszif 0<Abs[z]<1
    c1s0if z==0

Examples

open all close all

Basic Examples  (1)

Load the BLAS package:

In[1]:=
Click for copyable input

Compute a Givens rotation for given scalars:

In[2]:=
Click for copyable input
Out[2]=
In[3]:=
Click for copyable input
Out[3]=

Scope  (3)

Properties & Relations  (1)

Possible Issues  (2)