Changeset 4fc824e in git
- Timestamp:
- Mar 26, 2001, 11:15:26 PM (22 years ago)
- Branches:
- (u'spielwiese', '8d54773d6c9e2f1d2593a28bc68b7eeab54ed529')
- Children:
- bf4997b6efb6a0d9e3332369c2de3f6e6cee60d8
- Parents:
- 355f86ae1e09b962323c4d49da60c5aa3a7eb6c5
- Location:
- Singular
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/eigenval.cc
r355f86 r4fc824e 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: eigenval.cc,v 1. 3 2001-03-22 19:10:59Singular Exp $ */4 /* $Id: eigenval.cc,v 1.4 2001-03-26 21:15:26 Singular Exp $ */ 5 5 /* 6 6 * ABSTRACT: eigenvalues of constant square matrices … … 43 43 if(MATELEM(M,i,k)==NULL||MATELEM(M,j,k)==NULL) 44 44 return M; 45 number n=nDiv( nCopy(pGetCoeff(MATELEM(M,i,k))),46 nCopy(pGetCoeff(MATELEM(M,j,k))));45 number n=nDiv(pGetCoeff(MATELEM(M,i,k)), 46 pGetCoeff(MATELEM(M,j,k))); 47 47 for(int l=1;l<=MATCOLS(M);l++) 48 48 { … … 61 61 if(MATELEM(M,k,i)==NULL||MATELEM(M,k,j)==NULL) 62 62 return M; 63 number n=nDiv( nCopy(pGetCoeff(MATELEM(M,k,i))),64 nCopy(pGetCoeff(MATELEM(M,k,j))));63 number n=nDiv(pGetCoeff(MATELEM(M,k,i)), 64 pGetCoeff(MATELEM(M,k,j))); 65 65 for(int l=1;l<=MATROWS(M);l++) 66 66 { … … 81 81 { 82 82 int j=k+1; 83 while(j<=n&&MATELEM(M,j,k)==NULL) 84 { 85 j++; 86 } 83 while(j<=n&&MATELEM(M,j,k)==NULL) j++; 87 84 if(j<=n) 88 85 { … … 94 91 } 95 92 j=k+1; 96 while(j<=n&&MATELEM(M,k,j)==NULL) 97 { 98 j++; 99 } 93 while(j<=n&&MATELEM(M,k,j)==NULL) j++; 100 94 if(j<=n) 101 95 { … … 182 176 poly t=pOne(); 183 177 pSetExp(t,1,1); 178 pSetm(t); 184 179 poly d0=pSub(MATELEM(M,j,j),t); 185 180 MATELEM(M,j,j)=NULL; … … 193 188 poly t=pOne(); 194 189 pSetExp(t,1,1); 190 pSetm(t); 195 191 d0=pSub(pMult(d0,pSub(MATELEM(M,j,j),t)), 196 192 pMult(pMult(d2,MATELEM(M,j-1,j)),MATELEM(M,j,j-1))); … … 222 218 { 223 219 pp=p; 224 p=pNSet(nNeg(nDiv( nCopy(pGetCoeff(p0)),nCopy(pGetCoeff(p1)))));220 p=pNSet(nNeg(nDiv(pGetCoeff(p0),pGetCoeff(p1)))); 225 221 pDelete(&pp); 226 222 } … … 241 237 BOOLEAN tridiag(leftv res,leftv h) 242 238 { 243 if( h!=NULL&&h->Typ()==MATRIX_CMD)239 if((h!=NULL) && (h->Typ()==MATRIX_CMD)) 244 240 { 245 241 matrix M=(matrix)h->Data(); … … 253 249 } 254 250 WerrorS("<matrix> expected"); 251 return TRUE; 255 252 } 256 253 257 254 BOOLEAN eigenval(leftv res,leftv h) 258 255 { 259 if( h!=NULL&&h->Typ()==MATRIX_CMD)256 if((h!=NULL) && (h->Typ()==MATRIX_CMD)) 260 257 { 261 258 matrix M=(matrix)h->Data(); … … 269 266 } 270 267 WerrorS("<matrix> expected"); 271 } 272 268 return TRUE; 269 } 273 270 #endif -
Singular/iparith.cc
r355f86 r4fc824e 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: iparith.cc,v 1.26 1 2001-03-26 19:30:20Singular Exp $ */4 /* $Id: iparith.cc,v 1.262 2001-03-26 21:15:24 Singular Exp $ */ 5 5 6 6 /* … … 716 716 return TRUE; 717 717 } 718 syMake(u,u->name,NULL); 718 719 // else: use next case !!! no break !!! 719 720 … … 5249 5250 while(v->next!=NULL) 5250 5251 { 5251 nok = iiInternalExport(v->next, 0, u->data);5252 nok = iiInternalExport(v->next, 0, (idhdl)u->data); 5252 5253 if(nok) return nok; 5253 5254 v = v->next; -
Singular/ipid.cc
r355f86 r4fc824e 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: ipid.cc,v 1.4 7 2001-03-26 19:30:21Singular Exp $ */4 /* $Id: ipid.cc,v 1.48 2001-03-26 21:15:23 Singular Exp $ */ 5 5 6 6 /* … … 955 955 #ifdef HAVE_NAMESPACES 956 956 *pck =namespaceroot->get(p, myynest, TRUE); // search in toplevel namespace 957 namespaceroot->push(IDPACKAGE(*pck), IDID(*pck)); 958 *h =namespaceroot->get(i, myynest); // search in toplevel namespace 959 namespaceroot->pop(); 957 if((*pck!=NULL)&&(IDTYP(*pck)==PACKAGE_CMD)) 958 { 959 namespaceroot->push(IDPACKAGE(*pck), IDID(*pck)); 960 *h =namespaceroot->get(i, myynest); // search in toplevel namespace 961 namespaceroot->pop(); 962 } 960 963 #else /* HAVE_NAMESPACES */ 961 964 #endif /* HAVE_NAMESPACES */ -
Singular/weight.cc
r355f86 r4fc824e 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: weight.cc,v 1. 19 2001-03-05 16:39:02 mschulzeExp $ */4 /* $Id: weight.cc,v 1.20 2001-03-26 21:15:26 Singular Exp $ */ 5 5 6 6 /* … … 245 245 if(iv!=NULL) 246 246 len=iv->length(); 247 for(int i=pVariables;i>len;i--) 248 s[i]= 1; 249 for(;i>0;i--) 250 s[i]= (*iv)[i-1]; 247 int i; 248 for(i=pVariables;i>len;i--) s[i]=1; 249 for(;i>0;i--) s[i]=(*iv)[i-1]; 251 250 return s; 252 251 }
Note: See TracChangeset
for help on using the changeset viewer.