share/amatrix/amatrix.dem

   1 if ?get ('amatrix, 'present) = false then load (amatrix) else 'already_loaded;
   2
   3 "... Hilbert matrix ..."$
   4
   5 foo : make_matrix (7, 5)$
   6 for i thru foo@nr do for j thru foo@nc do foo [i, j] : 1/(i + j - 1);
   7 foo;
   8
   9 "... Integer subscripts ..."$
  10
  11 foo [1, 1];
  12
  13 foo [1, 1] : 42;
  14
  15 foo;
  16
  17 "... 'all subscripts ..."$
  18
  19 foo [all, 1];
  20 foo [1, all];
  21
  22 foo [all, all];
  23
  24 bar : foo [all, 1];
  25
  26 mumble : foo [1, all];
  27
  28 "... Collective assignment ..."$
  29
  30 foo [all, 2] : %pi;
  31
  32 foo;
  33
  34 foo [all, 5] : foo [all, 1];
  35
  36 foo;
  37
  38 "... Implicit subscript=1 when #rows=1 or #columns=1 ..."$
  39
  40 bar;
  41
  42 bar [1];
  43
  44 bar [all];
  45
  46 "... Copy-on-write policy ..."$
  47
  48 bar [1];
  49
  50 bar [1] : 1729;
  51
  52 bar;
  53
  54 mumble;
  55
  56 "... List of integers subscripts ..."$
  57
  58 foo [[1, 5, 3], 1];
  59
  60 foo [all, [5, 4, 3]];
  61
  62 foo [[5, 6, 7], [2, 5]];
  63
  64 "... amatrix of Boolean values subscripts ..."$
  65
  66 baz : amatrixmap (lambda ([x], is (x > 3/10)), bar);
  67
  68 quux : amatrixmap (lambda ([x], is (x < 1/4)), bar);
  69
  70 blurf : amatrixmap (lambda ([x], is (x < 1/2)), mumble);
  71
  72 snort : amatrixmap (lambda ([x], is (x >= 1/3)), mumble);
  73
  74 foo [baz, 4];
  75
  76 foo [6, blurf];
  77
  78 foo [baz, all];
  79
  80 foo [all, blurf];
  81
  82 foo [baz, [1, 5]];
  83
  84 foo [[2, 3, 6], blurf];
  85
  86 foo [baz, blurf];
  87
  88 "... Relational expressions subscripts ..."$
  89
  90 bar;
  91
  92 bar [bar > 2/3];
  93
  94 foo [bar >= 1/4, 4];
  95
  96 foo [6, 3/4 >= mumble];
  97
  98 foo [bar > 1/2, all];
  99
 100 foo [all, mumble < 1/2];
 101
 102 foo [bar # 1/2, [5, 2]];
 103
 104 foo [[4, 3, 2], mumble < 3/7];
 105
 106 foo [bar <= 7/8, blurf];
 107
 108 foo [baz, mumble > 1/8];
 109
 110 foo [bar > 2/7, 8/9 >= mumble];
 111
 112 "... Boolean expressions subscripts ..."$
 113
 114 foo [not baz, 3];
 115
 116 foo [7, blurf and snort];
 117
 118 foo [baz and not quux, all];
 119
 120 foo [all, not blurf or snort];
 121
 122 foo [baz or quux, [1, 2, 4]];
 123
 124 foo [[3, 2, 5], not blurf and snort];
 125
 126 foo [quux and baz, snort];
 127
 128 foo [quux, not blurf];
 129
 130 foo [not quux and not baz, mumble < 0.95];
 131
 132 foo [bar > 0.15, not blurf or not snort];
 133
 134 foo [quux or baz, blurf or not snort];
 135
 136