Go backward to Parallel Solutions
Go up to Top
Go forward to Correctness

Gaussian Elimination

Program Gauss
   assign
      <[] u :: {pivot with row u}
         <|| v,j : v notequal u ::
            M[v,j] := M[v,j]-
               M[v,u]*M[u,j]/M[u,u]
               if M[u,u] notequal 0>
      || <|| j ::
            M[u,j] := M[u,j]/M[u,u], if M[u,u] notequal 0>>
   [] {exchange rows}
      <|| u,v : u notequal v ::
         <|| j :: M[u,j],M[v,j] := M[v,j],M[u,j]
            if M[u,u]=0 and M[v,v]=0 and
               (M[u,v] notequal 0 or M[v,u] notequal 0)>>
end {Gauss}


Wolfgang.Schreiner@risc.uni-linz.ac.at
Id: unity3.tex,v 1.1 1996/04/19 12:29:22 schreine Exp schreine

Prev Up Next