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 (reset (dotident), 0);
  37 0;
  38
  39 x . dotident;
  40 x$
  41
  42 dotident . x;
  43 x$
  44
  45 x^^0;
  46 ''(dotident)$
  47
  48 (kill(id), dotident: id);
  49 id;
  50
  51 x^^0;
  52 id$
  53
  54 /*--- exponent rules ----*/
  55
  56 a^^2 . a^^5;
  57 a^^7$
  58
  59 a^^2 . a^^-2;
  60 ''(dotident)$
  61
  62 (1+x).(1+x)^^-1;
  63 ''(dotident)$
  64
  65 (1 + x - x^2).(1 + x - x^2)^^-1;
  66 ''(dotident)$
  67
  68 (1 + x - x^2)^^2 .(1 + x - x^2)^^-2;
  69 ''(dotident)$
  70
  71 a^^-2 . a^^3;
  72 a;
  73
  74 a^^-3 . a^^5;
  75 a^^2;
  76
  77 /*--- dotdistrib ---*/
  78
  79 a . (b + c), dotdistrib : true;
  80 a . b + a . c$
  81
  82 a . (b - c), dotdistrib : true;
  83 a . b - a . c$
  84
  85 a . (b +  c / 5), dotdistrib : true;
  86 a . b + a . c / 5$
  87
  88 /*--- dotexptsimp ---*/
  89
  90 is(equal(a . a, a^^2)), dotexptsimp : false, prederror : false;
  91 unknown$
  92
  93 is(equal(a . a, a^^2)), dotexptsimp : true;
  94 true$
  95
  96 a . a, dotexptsimp : true;
  97 a^^2$
  98
  99 a . a - a^^2;
 100 0$
 101
 102 a^^5 . a^^7;
 103 a^^12$
 104
 105 is(equal(a.a,a * a)), prederror : false;
 106 unknown$
 107
 108 is(equal(a.b, b. a)),prederror : false;
 109 unknown$
 110
 111 is(equal(a^^2, a^2)),prederror : false;
 112 unknown$
 113
 114 is(equal(a.b.a, a^^2 . b)),prederror : false;
 115 unknown$
 116
 117 /*  See SF Bug # 771061 */
 118
 119 expand((vt . a^^(-1) . u+1)^^(-2));
 120 ((vt.a^^(-1).u)^^2+2*(vt.a^^(-1).u)+1)^^(-1)$
 121
 122 (a + 1)^^(-1) . (a + 1)^^(-1), expon:2;
 123 (a^^2+2*a+1)^^(-1)$
 124
 125 (a + b)^^2, expand;
 126 a^^2 + a.b + b . a + b^^2$
 127
 128 (a + b)^^-2, expand;
 129 (a^^2 + a . b + b . a + b^^2)^^-1$
 130
 131 (a - b)^^2, expand;
 132 a^^2 - a . b - b . a + b^^2$
 133
 134 (5 * a + 7 *b)^^2, expand;
 135 25 * a^^2 + 35 * a . b + 35 * b . a + 49 * b^^2$
 136
 137 (a . b)^^2, expand, dotexptsimp : false;
 138 ''(block([dotexptsimp : false], a . b . a . b))$
 139
 140 (a . b)^^-1, expand, dotexptsimp : false;
 141 (a . b)^^-1$
 142
 143 a.b.a.b, dotexptsimp : true;
 144 (a.b)^^2$
 145
 146 /*--- See SF Bug # 629716 ---*/
 147
 148 expand((x.y)^^2 . y - x.y.x.y^^2), dotexptsimp : false;
 149 0$
 150
 151 5 * a . b - a . (5 * b);
 152 0$
 153
 154 (5 * a) . b - 5 * (a . b);
 155 0$
 156
 157 a . (5 * b) - 5 *(a . b);
 158 0$
 159
 160 /*--- declared scalars---*/
 161
 162 (declare(alpha, scalar),0);
 163 0$
 164
 165 alpha . x, dotscrules : true;
 166 alpha * x$
 167
 168 (alpha . a) . b, dotscrules : true;
 169 alpha * a . b$
 170
 171 a . (alpha . b), dotscrules : true;
 172 alpha * a . b$
 173
 174 declare(".",antisymmetric);
 175 'done$
 176
 177 b . alpha, dotscrules : true;
 178 alpha * b$
 179
 180 a . (b . alpha), dotscrules : true;
 181 alpha * a . b$
 182
 183 a . (alpha * a), dotscrules : true;
 184 0;
 185
 186 remove(".",antisymmetric);
 187 'done$
 188
 189 /*--- dotassoc ---*/
 190
 191 is(equal(a . (b . c) , (a . b) . c)), dotassoc : true, prederror : false;
 192 true$
 193
 194 is(equal(a . (b . c) , (a . b) . c)), dotassoc : false, prederror : false;
 195 unknown$
 196
 197 (declare(".", commutative),0);
 198 0$
 199
 200 is(equal(a . b, b . a));
 201 true$
 202
 203 is(equal(a . (b . (c .d)), (a.b) . (c.d)));
 204 true$
 205
 206 (remove(".",commutative),0);
 207 0$
 208
 209 /* SF bug #3935: Noncommutative multiplication with string argument triggers "declare: argument must be a symbol" */
 210
 211 (kill(a, b),
 212  "a" . b);
 213 "a" . b;
 214
 215 string ("a" . b);
 216 "\"a\" . b";
 217
 218 a . "b";
 219 a . "b";
 220
 221 string (a . "b");
 222 "a . \"b\"";
 223
 224 "a" . "b";
 225 "a" . "b";
 226
 227 string ("a" . "b");
 228 "\"a\" . \"b\"";
 229
 230 (a: make_array (fixnum, 3, 2),
 231  b: make_array (fixnum, 2, 3),
 232  a . b);
 233 ''(a . b); /* this result just verifies there isn't an error */
 234
 235 /* try noncommutative exponent also */
 236 "a"^^-1;
 237 "a"^^-1;