Changeset 76b4bd in git for Singular/eigenval_ip.cc

Timestamp:

Mar 4, 2004, 11:08:24 AM (20 years ago)

Author:

Hans Schönemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

cec18c94b438d2504bc3c731cea94331af31cd83

Parents:

c4097df6269597fe0c7c72d02451429687a8b584

Message:

*hannes: eigenval->eigenval_ip, kernel/eigenval


git-svn-id: file:///usr/local/Singular/svn/trunk@7073 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• Singular/eigenval_ip.cc (moved) (moved from Singular/eigenval.cc) (5 diffs)

Singular/eigenval_ip.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR      *
 *****************************************/
-/* $Id: eigenval.cc,v 1.13 2002-05-31 09:00:12 Singular Exp $ */
+/* $Id: eigenval_ip.cc,v 1.1 2004-03-04 10:06:50 Singular Exp $ */
 /*
 * ABSTRACT: eigenvalues of constant square matrices
@@ -22 +22 @@
 #include "clapsing.h"
 #include "eigenval.h"
-
-
-matrix evSwap(matrix M,int i,int j)
-{
-  if(i==j)
-    return(M);
-
-  for(int k=1;k<=MATROWS(M);k++)
-  {
-    poly p=MATELEM(M,i,k);
-    MATELEM(M,i,k)=MATELEM(M,j,k);
-    MATELEM(M,j,k)=p;
-  }
-
-  for(int k=1;k<=MATCOLS(M);k++)
-  {
-    poly p=MATELEM(M,k,i);
-    MATELEM(M,k,i)=MATELEM(M,k,j);
-    MATELEM(M,k,j)=p;
-  }
-
-  return(M);
-}
+#include "eigenval_ip.h"
 
 
@@ -74 +52 @@
   return TRUE;
 }
-
-
-matrix evRowElim(matrix M,int i,int j,int k)
-{
-  if(MATELEM(M,i,k)==NULL||MATELEM(M,j,k)==NULL)
-    return(M);
-
-  poly p=pNSet(nDiv(pGetCoeff(MATELEM(M,i,k)),pGetCoeff(MATELEM(M,j,k))));
-  pNormalize(p);
-
-  for(int l=1;l<=MATCOLS(M);l++)
-  {
-    MATELEM(M,i,l)=pSub(MATELEM(M,i,l),ppMult_qq(p,MATELEM(M,j,l)));
-    pNormalize(MATELEM(M,i,l));
-  }
-  for(int l=1;l<=MATROWS(M);l++)
-  {
-    MATELEM(M,l,j)=pAdd(MATELEM(M,l,j),ppMult_qq(p,MATELEM(M,l,i)));
-    pNormalize(MATELEM(M,l,j));
-  }
-
-  pDelete(&p);
-
-  return(M);
-}
-
 
 BOOLEAN evRowElim(leftv res,leftv h)
@@ -134 +86 @@
 }
 
-
-matrix evColElim(matrix M,int i,int j,int k)
-{
-  if(MATELEM(M,k,i)==0||MATELEM(M,k,j)==0)
-    return(M);
-
-  poly p=pNSet(nDiv(pGetCoeff(MATELEM(M,k,i)),pGetCoeff(MATELEM(M,k,j))));
-  pNormalize(p);
-
-  for(int l=1;l<=MATROWS(M);l++)
-  {
-    MATELEM(M,l,i)=pSub(MATELEM(M,l,i),ppMult_qq(p,MATELEM(M,l,j)));
-    pNormalize(MATELEM(M,l,i));
-  }
-  for(int l=1;l<=MATCOLS(M);l++)
-  {
-    MATELEM(M,j,l)=pAdd(MATELEM(M,j,l),ppMult_qq(p,MATELEM(M,i,l)));
-    pNormalize(MATELEM(M,j,l));
-  }
-
-  pDelete(&p);
-
-  return(M);
-}
-
-
 BOOLEAN evColElim(leftv res,leftv h)
 {
@@ -192 +118 @@
   return TRUE;
 }
-
-
-matrix evHessenberg(matrix M)
-{
-  int n=MATROWS(M);
-  if(n!=MATCOLS(M))
-    return(M);
-
-  for(int k=1,j=2;k<n-1;k++,j=k+1)
-  {
-    while(j<=n&&MATELEM(M,j,k)==0)
-      j++;
-
-    if(j<=n)
-    {
-      M=evSwap(M,j,k+1);
-
-      for(int i=j+1;i<=n;i++)
-        M=evRowElim(M,i,k+1,k);
-    }
-  }
-
-  return(M);
-}
-
 
 BOOLEAN evHessenberg(leftv res,leftv h)