Changeset 15d67db in git


Ignore:
Timestamp:
Feb 20, 2002, 6:32:12 PM (21 years ago)
Author:
Mathias Schulze <mschulze@…>
Branches:
(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'f875bbaccd0831e36aaed09ff6adeb3eb45aeb94')
Children:
2ca72f72ba5811479c336c3a36c9d120db4edfc7
Parents:
a7c2dbcae69f5b83fdb67594276913303071d020
Message:
*** empty log message ***


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

Legend:

Unmodified
Added
Removed
  • Singular/eigenval.cc

    ra7c2db r15d67db  
    22*  Computer Algebra System SINGULAR      *
    33*****************************************/
    4 /* $Id: eigenval.cc,v 1.10 2002-02-20 15:50:34 mschulze Exp $ */
     4/* $Id: eigenval.cc,v 1.11 2002-02-20 17:32:12 mschulze Exp $ */
    55/*
    66* ABSTRACT: eigenvalues of constant square matrices
     
    8282
    8383  poly p=pNSet(nDiv(pGetCoeff(MATELEM(M,i,k)),pGetCoeff(MATELEM(M,j,k))));
     84  pNormalize(p);
    8485
    8586  for(int l=1;l<=MATCOLS(M);l++)
     87  {
    8688    MATELEM(M,i,l)=pSub(MATELEM(M,i,l),ppMult_qq(p,MATELEM(M,j,l)));
    87 
     89    pNormalize(MATELEM(M,i,l));
     90  }
    8891  for(int l=1;l<=MATROWS(M);l++)
     92  {
    8993    MATELEM(M,l,j)=pAdd(MATELEM(M,l,j),ppMult_qq(p,MATELEM(M,l,i)));
     94    pNormalize(MATELEM(M,l,j));
     95  }
    9096
    9197  pDelete(&p);
     
    135141
    136142  poly p=pNSet(nDiv(pGetCoeff(MATELEM(M,k,i)),pGetCoeff(MATELEM(M,k,j))));
     143  pNormalize(p);
    137144
    138145  for(int l=1;l<=MATROWS(M);l++)
     146  {
    139147    MATELEM(M,l,i)=pSub(MATELEM(M,l,i),ppMult_qq(p,MATELEM(M,l,j)));
    140 
     148    pNormalize(MATELEM(M,l,i));
     149  }
    141150  for(int l=1;l<=MATCOLS(M);l++)
     151  {
    142152    MATELEM(M,j,l)=pAdd(MATELEM(M,j,l),ppMult_qq(p,MATELEM(M,i,l)));
     153    pNormalize(MATELEM(M,j,l));
     154  }
    143155
    144156  pDelete(&p);
     
    284296            e->m[k]=pNSet(nDiv(e1,pGetCoeff(pNext(e0->m[i]))));
    285297          nDelete(&e1);
     298          pNormalize(e->m[k]);
    286299          (*m)[k]=(*m0)[i];
    287300          k++;
     
    290303        {
    291304          e->m[k]=e0->m[i];
     305          pNormalize(e->m[k]);
    292306          e0->m[i]=NULL;
    293307          (*m)[k]=(*m0)[i];
  • Singular/gms.cc

    ra7c2db r15d67db  
    22*  Computer Algebra System SINGULAR      *
    33*****************************************/
    4 /* $Id: gms.cc,v 1.6 2002-02-20 15:50:34 mschulze Exp $ */
     4/* $Id: gms.cc,v 1.7 2002-02-20 17:32:12 mschulze Exp $ */
    55/*
    66* ABSTRACT: Gauss-Manin system normal form
Note: See TracChangeset for help on using the changeset viewer.