tests/rtest_scalarp.mac

   1 (kill (all), 0);
   2 0;
   3
   4 (remvalue(a,b,c,d,x,y),0);
   5 0$
   6
   7 (declare([a,b],scalar),0);
   8 0$
   9
  10 (sc : [-5 ,0, 8, 7.8, 1.23b0, %pi, %phi, %i, sqrt(1+%pi), a, a*b, a+b, a/b,a^3, b^a, a^sqrt(2)],0);
  11 0$
  12
  13 (buggy : [],0);
  14 0$
  15
  16 (for si in sc do if not(scalarp(si)) then buggy : cons(si,buggy),0);
  17 0$
  18
  19 buggy;
  20 []$
  21
  22 (nsc : [x,x+y,x*y,x^3, x = y, x < y, x^^5, cosh(x), sqrt(x), [5,x], matrix([1,2],[3,4])],0);
  23 0$
  24
  25 (buggy : [],0);
  26 0$
  27
  28 (for si in nsc do if scalarp(si) then buggy : cons(si,buggy),0);
  29 0$
  30
  31 buggy;
  32 []$
  33
  34 determinant(matrix([1,2],[3,4]));
  35 -2$
  36
  37 determinant(matrix([a,b],[c,d]));
  38 a*d - b * c$
  39
  40 determinant(a);
  41 a$
  42
  43 is(verbify(op(determinant(x))) = 'determinant);
  44 true$
  45
  46 determinant(sqrt(5));
  47 sqrt(5)$
  48
  49 determinant(7*%i + 5);
  50 7*%i + 5$
  51
  52 block([matrix_element_add:'f, matrix_element_mult:'g],
  53   determinant( matrix([a,b],[c,d]) ) );
  54 f(f(0,g(1,a,d)),g(-1,b,c))$
  55
  56 (remove([a,b],scalar),0);
  57 0$
  58
  59 (remvalue(sc, nsc,buggy),0);
  60 0$
  61
  62
  63
  64
  65
  66
  67
  68
  69
  70