| 1 | #include "mod2.h" |
|---|
| 2 | |
|---|
| 3 | #include <omalloc/omalloc.h> |
|---|
| 4 | #include <misc/auxiliary.h> |
|---|
| 5 | |
|---|
| 6 | #include "polys.h" |
|---|
| 7 | |
|---|
| 8 | |
|---|
| 9 | /* Returns TRUE if |
|---|
| 10 | * LM(p) | LM(lcm) |
|---|
| 11 | * LC(p) | LC(lcm) only if ring |
|---|
| 12 | * Exists i, j: |
|---|
| 13 | * LE(p, i) != LE(lcm, i) |
|---|
| 14 | * LE(p1, i) != LE(lcm, i) ==> LCM(p1, p) != lcm |
|---|
| 15 | * LE(p, j) != LE(lcm, j) |
|---|
| 16 | * LE(p2, j) != LE(lcm, j) ==> LCM(p2, p) != lcm |
|---|
| 17 | */ |
|---|
| 18 | BOOLEAN pCompareChain (poly p,poly p1,poly p2,poly lcm, const ring R) |
|---|
| 19 | { |
|---|
| 20 | int k, j; |
|---|
| 21 | |
|---|
| 22 | if (lcm==NULL) return FALSE; |
|---|
| 23 | |
|---|
| 24 | for (j=(R->N); j; j--) |
|---|
| 25 | if ( p_GetExp(p,j, R) > p_GetExp(lcm,j, R)) return FALSE; |
|---|
| 26 | if ( pGetComp(p) != pGetComp(lcm)) return FALSE; |
|---|
| 27 | for (j=(R->N); j; j--) |
|---|
| 28 | { |
|---|
| 29 | if (p_GetExp(p1,j, R)!=p_GetExp(lcm,j, R)) |
|---|
| 30 | { |
|---|
| 31 | if (p_GetExp(p,j, R)!=p_GetExp(lcm,j, R)) |
|---|
| 32 | { |
|---|
| 33 | for (k=(R->N); k>j; k--) |
|---|
| 34 | { |
|---|
| 35 | if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) |
|---|
| 36 | && (p_GetExp(p2,k, R)!=p_GetExp(lcm,k, R))) |
|---|
| 37 | return TRUE; |
|---|
| 38 | } |
|---|
| 39 | for (k=j-1; k; k--) |
|---|
| 40 | { |
|---|
| 41 | if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) |
|---|
| 42 | && (p_GetExp(p2,k, R)!=p_GetExp(lcm,k, R))) |
|---|
| 43 | return TRUE; |
|---|
| 44 | } |
|---|
| 45 | return FALSE; |
|---|
| 46 | } |
|---|
| 47 | } |
|---|
| 48 | else if (p_GetExp(p2,j, R)!=p_GetExp(lcm,j, R)) |
|---|
| 49 | { |
|---|
| 50 | if (p_GetExp(p,j, R)!=p_GetExp(lcm,j, R)) |
|---|
| 51 | { |
|---|
| 52 | for (k=(R->N); k>j; k--) |
|---|
| 53 | { |
|---|
| 54 | if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) |
|---|
| 55 | && (p_GetExp(p1,k, R)!=p_GetExp(lcm,k, R))) |
|---|
| 56 | return TRUE; |
|---|
| 57 | } |
|---|
| 58 | for (k=j-1; k!=0 ; k--) |
|---|
| 59 | { |
|---|
| 60 | if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) |
|---|
| 61 | && (p_GetExp(p1,k, R)!=p_GetExp(lcm,k, R))) |
|---|
| 62 | return TRUE; |
|---|
| 63 | } |
|---|
| 64 | return FALSE; |
|---|
| 65 | } |
|---|
| 66 | } |
|---|
| 67 | } |
|---|
| 68 | return FALSE; |
|---|
| 69 | } |
|---|
| 70 | #ifdef HAVE_RATGRING |
|---|
| 71 | BOOLEAN pCompareChainPart (poly p,poly p1,poly p2,poly lcm, const ring R) |
|---|
| 72 | { |
|---|
| 73 | int k, j; |
|---|
| 74 | |
|---|
| 75 | if (lcm==NULL) return FALSE; |
|---|
| 76 | |
|---|
| 77 | for (j=R->real_var_end; j>=R->real_var_start; j--) |
|---|
| 78 | if ( p_GetExp(p,j, R) > p_GetExp(lcm,j, R)) return FALSE; |
|---|
| 79 | if ( pGetComp(p) != pGetComp(lcm)) return FALSE; |
|---|
| 80 | for (j=R->real_var_end; j>=R->real_var_start; j--) |
|---|
| 81 | { |
|---|
| 82 | if (p_GetExp(p1,j, R)!=p_GetExp(lcm,j, R)) |
|---|
| 83 | { |
|---|
| 84 | if (p_GetExp(p,j, R)!=p_GetExp(lcm,j, R)) |
|---|
| 85 | { |
|---|
| 86 | for (k=(R->N); k>j; k--) |
|---|
| 87 | for (k=R->real_var_end; k>j; k--) |
|---|
| 88 | { |
|---|
| 89 | if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) |
|---|
| 90 | && (p_GetExp(p2,k, R)!=p_GetExp(lcm,k, R))) |
|---|
| 91 | return TRUE; |
|---|
| 92 | } |
|---|
| 93 | for (k=j-1; k>=R->real_var_start; k--) |
|---|
| 94 | { |
|---|
| 95 | if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) |
|---|
| 96 | && (p_GetExp(p2,k, R)!=p_GetExp(lcm,k, R))) |
|---|
| 97 | return TRUE; |
|---|
| 98 | } |
|---|
| 99 | return FALSE; |
|---|
| 100 | } |
|---|
| 101 | } |
|---|
| 102 | else if (p_GetExp(p2,j, R)!=p_GetExp(lcm,j, R)) |
|---|
| 103 | { |
|---|
| 104 | if (p_GetExp(p,j, R)!=p_GetExp(lcm,j, R)) |
|---|
| 105 | { |
|---|
| 106 | for (k=R->real_var_end; k>j; k--) |
|---|
| 107 | { |
|---|
| 108 | if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) |
|---|
| 109 | && (p_GetExp(p1,k, R)!=p_GetExp(lcm,k, R))) |
|---|
| 110 | return TRUE; |
|---|
| 111 | } |
|---|
| 112 | for (k=j-1; k>=R->real_var_start; k--) |
|---|
| 113 | { |
|---|
| 114 | if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) |
|---|
| 115 | && (p_GetExp(p1,k, R)!=p_GetExp(lcm,k, R))) |
|---|
| 116 | return TRUE; |
|---|
| 117 | } |
|---|
| 118 | return FALSE; |
|---|
| 119 | } |
|---|
| 120 | } |
|---|
| 121 | } |
|---|
| 122 | return FALSE; |
|---|
| 123 | } |
|---|
| 124 | #endif |
|---|
| 125 | |
|---|
