tests/rtest_dot.mac

   1 (kill(all),0);
   2 0$
   3
   4 1 . 1, dotconstrules : true;
   5 1$
   6
   7 is(op(1 . x) = "."),  dotconstrules : false, dot1simp : false;
   8 true;
   9
  10 (1/3) . (2/3);
  11 2/9$
  12
  13 1.1 . 1.1;
  14 1.21$
  15
  16 bfloat_approx_equal(1.1 . 1.1b0, 1.21b0);
  17 true;
  18
  19 3.1 . 5;
  20 15.5$
  21
  22 1 . x;
  23 x$
  24
  25 x . 1$
  26 x$
  27
  28 0 . x$
  29 0$
  30
  31 x . 0$
  32 0$
  33
  34 /*--- dotident ---*/
  35
  36 x . dotident;
  37 x$
  38
  39 dotident . x;
  40 x$
  41
  42 x^^0;
  43 ''(dotident)$
  44
  45 x^^0, dotident : id;
  46 id$
  47
  48 /*--- exponent rules ----*/
  49
  50 a^^2 . a^^5;
  51 a^^7$
  52
  53 a^^2 . a^^-2;
  54 ''(dotident)$
  55
  56 (1+x).(1+x)^^-1;
  57 ''(dotident)$
  58
  59 (1 + x - x^2).(1 + x - x^2)^^-1;
  60 ''(dotident)$
  61
  62 (1 + x - x^2)^^2 .(1 + x - x^2)^^-2;
  63 ''(dotident)$
  64
  65 /*--- dotdistrib ---*/
  66
  67 a . (b + c), dotdistrib : true;
  68 a . b + a . c$
  69
  70 a . (b - c), dotdistrib : true;
  71 a . b - a . c$
  72
  73 a . (b +  c / 5), dotdistrib : true;
  74 a . b + a . c / 5$
  75
  76 /*--- dotexptsimp ---*/
  77
  78 is(equal(a . a, a^^2)), dotexptsimp : false, prederror : false;
  79 unknown$
  80
  81 is(equal(a . a, a^^2)), dotexptsimp : true;
  82 true$
  83
  84 a . a, dotexptsimp : true;
  85 a^^2$
  86
  87 a . a - a^^2;
  88 0$
  89
  90 a^^5 . a^^7;
  91 a^^12$
  92
  93 is(equal(a.a,a * a)), prederror : false;
  94 unknown$
  95
  96 is(equal(a.b, b. a)),prederror : false;
  97 unknown$
  98
  99 is(equal(a^^2, a^2)),prederror : false;
 100 unknown$
 101
 102 is(equal(a.b.a, a^^2 . b)),prederror : false;
 103 unknown$
 104
 105 /*  See SF Bug # 771061 */
 106
 107 expand((vt . a^^(-1) . u+1)^^(-2));
 108 ((vt.a^^(-1).u)^^2+2*(vt.a^^(-1).u)+1)^^(-1)$
 109
 110 (a + 1)^^(-1) . (a + 1)^^(-1), expon:2;
 111 (a^^2+2*a+1)^^(-1)$
 112
 113 (a + b)^^2, expand;
 114 a^^2 + a.b + b . a + b^^2$
 115
 116 (a + b)^^-2, expand;
 117 (a^^2 + a . b + b . a + b^^2)^^-1$
 118
 119 (a - b)^^2, expand;
 120 a^^2 - a . b - b . a + b^^2$
 121
 122 (5 * a + 7 *b)^^2, expand;
 123 25 * a^^2 + 35 * a . b + 35 * b . a + 49 * b^^2$
 124
 125 (a . b)^^2, expand, dotexptsimp : false;
 126 ''(block([dotexptsimp : false], a . b . a . b))$
 127
 128 (a . b)^^-1, expand, dotexptsimp : false;
 129 (a . b)^^-1$
 130
 131 a.b.a.b, dotexptsimp : true;
 132 (a.b)^^2$
 133
 134 /*--- See SF Bug # 629716 ---*/
 135
 136 expand((x.y)^^2 . y - x.y.x.y^^2), dotexptsimp : false;
 137 0$
 138
 139 5 * a . b - a . (5 * b);
 140 0$
 141
 142 (5 * a) . b - 5 * (a . b);
 143 0$
 144
 145 a . (5 * b) - 5 *(a . b);
 146 0$
 147
 148 /*--- declared scalars---*/
 149
 150 (declare(alpha, scalar),0);
 151 0$
 152
 153 alpha . x, dotscrules : true;
 154 alpha * x$
 155
 156 (alpha . a) . b, dotscrules : true;
 157 alpha * a . b$
 158
 159 a . (alpha . b), dotscrules : true;
 160 alpha * a . b$
 161
 162 declare(".",antisymmetric);
 163 'done$
 164
 165 b . alpha, dotscrules : true;
 166 alpha * b$
 167
 168 a . (b . alpha), dotscrules : true;
 169 alpha * a . b$
 170
 171 a . (alpha * a), dotscrules : true;
 172 0;
 173
 174 remove(".",antisymmetric);
 175 'done$
 176
 177 /*--- dotassoc ---*/
 178
 179 is(equal(a . (b . c) , (a . b) . c)), dotassoc : true, prederror : false;
 180 true$
 181
 182 is(equal(a . (b . c) , (a . b) . c)), dotassoc : false, prederror : false;
 183 unknown$
 184
 185 (declare(".", commutative),0);
 186 0$
 187
 188 is(equal(a . b, b . a));
 189 true$
 190
 191 is(equal(a . (b . (c .d)), (a.b) . (c.d)));
 192 true$
 193
 194 (remove(".",commutative),0);
 195 0$
 196
 197
 198
 199
 200
 201
 202
 203
 204
 205