Changeset 80ad5e in git
- Timestamp:
- Apr 7, 1998, 8:13:57 PM (25 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '00e2e9c41af3fde1273eb3633f4c0c7c3db2579d')
- Children:
- 97454d23ed1d5477f8790a863d5d481deacedbcc
- Parents:
- d86e64ab96a8a8f265a4d1b473857a697f643276
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/hdegree.cc
rd86e64 r80ad5e 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: hdegree.cc,v 1.1 0 1997-12-15 22:46:23 obachmanExp $ */4 /* $Id: hdegree.cc,v 1.11 1998-04-07 18:13:57 Singular Exp $ */ 5 5 /* 6 6 * ABSTRACT - dimension, multiplicity, HC, kbase … … 73 73 Exponent_t mc; 74 74 hexist = hInit(S, Q, &hNexist); 75 if ( !hNexist)75 if (hNexist==0) 76 76 return pVariables; 77 77 hwork = (scfmon)Alloc(hNexist * sizeof(scmon)); … … 79 79 hpure = (scmon)Alloc((1 + (pVariables * pVariables)) * sizeof(Exponent_t)); 80 80 mc = hisModule; 81 if ( !mc)81 if (mc==0) 82 82 { 83 83 hrad = hexist; … … 90 90 loop 91 91 { 92 if (mc )92 if (mc!=0) 93 93 hComp(hexist, hNexist, mc, hrad, &hNrad); 94 if (hNrad )94 if (hNrad!=0) 95 95 { 96 96 hNvar = pVariables; 97 97 hRadical(hrad, &hNrad, hNvar); 98 98 hSupp(hrad, hNrad, hvar, &hNvar); 99 if (hNvar )99 if (hNvar!=0) 100 100 { 101 101 memset(hpure, 0, (pVariables + 1) * sizeof(Exponent_t)); … … 119 119 Free((ADDRESS)hwork, hNexist * sizeof(scmon)); 120 120 hDelete(hexist, hNexist); 121 if (hisModule )121 if (hisModule!=0) 122 122 Free((ADDRESS)hrad, hNexist * sizeof(scmon)); 123 123 return pVariables - hCo; … … 139 139 { 140 140 hCo = dn; 141 for (iv=pVariables; iv ; iv--)142 { 143 if (pure[iv] )141 for (iv=pVariables; iv>0; iv--) 142 { 143 if (pure[iv]!=0) 144 144 hInd[iv] = 0; 145 145 else 146 146 hInd[iv] = 1; 147 147 } 148 if (Nrad )148 if (Nrad!=0) 149 149 { 150 150 pn = *rad; … … 153 153 { 154 154 x = var[iv]; 155 if (pn[x] )155 if (pn[x]!=0) 156 156 { 157 157 hInd[x] = 0; … … 169 169 while(pure[var[iv]]) iv--; 170 170 hStepR(rad, Nrad, var, iv, &rad0); 171 if (rad0 )171 if (rad0!=0) 172 172 { 173 173 iv--; … … 195 195 { 196 196 hCo = Npure + 1; 197 for (x=pVariables; x ; x--)198 { 199 if (pure[x] )197 for (x=pVariables; x>0; x--) 198 { 199 if (pure[x]!=0) 200 200 hInd[x] = 0; 201 201 else … … 265 265 Free((ADDRESS)hwork, hNexist * sizeof(scmon)); 266 266 hDelete(hexist, hNexist); 267 if (hisModule )267 if (hisModule!=0) 268 268 Free((ADDRESS)hrad, hNexist * sizeof(scmon)); 269 269 return Set; … … 303 303 for (iv=pVariables; iv!=0 ; iv--) 304 304 { 305 if (pure[iv] )305 if (pure[iv]!=0) 306 306 (*Set)[iv-1] = 0; 307 307 else … … 331 331 { 332 332 x = var[iv]; 333 if (pn[x] )333 if (pn[x]!=0) 334 334 { 335 335 pure[x] = 1; … … 350 350 loop 351 351 { 352 if( !pure[var[iv]])352 if(pure[var[iv]]==0) 353 353 { 354 354 if(hNotZero(rad, Nrad, var, iv)) … … 360 360 } 361 361 iv--; 362 if ( !iv)362 if (iv==0) 363 363 return; 364 364 } 365 365 } 366 while(pure[var[iv]] ) iv--;366 while(pure[var[iv]]!=0) iv--; 367 367 hStepR(rad, Nrad, var, iv, &rad0); 368 368 iv--; … … 486 486 for (iv=pVariables; iv; iv--) 487 487 { 488 if (pure[iv] )488 if (pure[iv]!=0) 489 489 (*Set)[iv-1] = 0; 490 490 else … … 506 506 if (dn > hCo) 507 507 { 508 if ( !Nrad)508 if (Nrad==0) 509 509 hCheckIndep(pure); 510 510 else 511 511 { 512 512 pn = *rad; 513 for (iv = Nvar; iv ; iv--)513 for (iv = Nvar; iv>0; iv--) 514 514 { 515 515 x = var[iv]; 516 if (pn[x] )516 if (pn[x]!=0) 517 517 { 518 518 pure[x] = 1; … … 526 526 } 527 527 iv = Nvar; 528 while(pure[var[iv]] ) iv--;528 while(pure[var[iv]]!=0) iv--; 529 529 hStepR(rad, Nrad, var, iv, &rad0); 530 530 iv--; … … 574 574 hRadical(hrad, &hNrad, hNvar); 575 575 hSupp(hrad, hNrad, hvar, &hNvar); 576 if (hNvar )576 if (hNvar!=0) 577 577 { 578 578 hCo = hNvar; … … 582 582 hDimSolve(hpure, hNpure, hrad, hNrad, hvar, hNvar); 583 583 } 584 if ( hCo&& (hCo < pVariables))584 if ((hCo!=0) && (hCo < pVariables)) 585 585 { 586 586 hIndMult(hpure, hNpure, hrad, hNrad, hvar, hNvar); … … 649 649 scmon pn; 650 650 scfmon sn; 651 if ( !iv)651 if (iv==0) 652 652 return pure[var[1]]; 653 else if ( !Nstc)653 else if (Nstc==0) 654 654 { 655 655 sum = 1; 656 for (i = Nvar; i ; i--)656 for (i = Nvar; i>0; i--) 657 657 sum *= pure[var[i]]; 658 658 return sum; … … 695 695 for (i = 1; i <= pVariables; i++) 696 696 { 697 if (pure[i] )697 if (pure[i]!=0) 698 698 { 699 699 i0++; … … 723 723 if (dn == hCo) 724 724 { 725 if ( !Nrad)725 if (Nrad==0) 726 726 hProject(pure, hsel); 727 727 else 728 728 { 729 729 pn = *rad; 730 for (iv = Nvar; iv ; iv--)730 for (iv = Nvar; iv>0; iv--) 731 731 { 732 732 x = var[iv]; 733 if (pn[x] )733 if (pn[x]!=0) 734 734 { 735 735 pure[x] = 1; … … 750 750 loop 751 751 { 752 if( !pure[var[iv]])752 if(pure[var[iv]]==0) 753 753 { 754 754 if(hNotZero(rad, Nrad, var, iv)) … … 760 760 } 761 761 iv--; 762 if ( !iv)762 if (iv==0) 763 763 return; 764 764 } … … 793 793 Exponent_t mc; 794 794 hexist = hInit(S, Q, &hNexist); 795 if ( !hNexist)795 if (hNexist==0) 796 796 { 797 797 hCo = 0; … … 806 806 mc = hisModule; 807 807 hrad = (scfmon)Alloc(hNexist * sizeof(scmon)); 808 if ( !mc)808 if (mc==0) 809 809 { 810 810 memcpy(hrad, hexist, hNexist * sizeof(scmon)); … … 820 820 loop 821 821 { 822 if (mc )822 if (mc!=0) 823 823 { 824 824 hComp(hexist, hNexist, mc, hrad, &hNrad); … … 826 826 memcpy(hstc, hrad, hNrad * sizeof(scmon)); 827 827 } 828 if (hNrad )828 if (hNrad!=0) 829 829 { 830 830 hNvar = pVariables; 831 831 hRadical(hrad, &hNrad, hNvar); 832 832 hSupp(hrad, hNrad, hvar, &hNvar); 833 if (hNvar )833 if (hNvar!=0) 834 834 { 835 835 hCo = hNvar; … … 850 850 hMu = 0; 851 851 } 852 if ( hNvar&& (hCo == di))853 { 854 if ( di&& (di < pVariables))852 if ((hNvar!=0) && (hCo == di)) 853 { 854 if ((di!=0) && (di < pVariables)) 855 855 hDimMult(hpure, hNpure, hrad, hNrad, hvar, hNvar); 856 else if ( !di)856 else if (di==0) 857 857 hMu++; 858 858 else
Note: See TracChangeset
for help on using the changeset viewer.