share/vector/vector_rebuild.mac

   1 /*
   2   Some utilities for working with vectors.
   3
   4   Copyright (C)  Nov. 2008  Volker van Nek
   5
   6   For examples see vector_rebuild.usg .
   7
   8   This program is free software; you can redistribute it and/or modify
   9   it under the terms of the GNU General Public License as published by
  10   the Free Software Foundation; either version 2 of the License, or
  11   (at your option) any later version.
  12
  13   This program is distributed in the hope that it will be useful,
  14   but WITHOUT ANY WARRANTY; without even the implied warranty of
  15   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16   GNU General Public License for more details.
  17
  18   You should have received a copy of the GNU General Public License
  19   along with this program; if not, write to the Free Software
  20   Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
  21   MA 02110-1301, USA.
  22 */
  23
  24 eval_when([translate,batch,demo,load], load("vector_rebuild.lisp"))$
  25
  26 eval_when([batch,demo,load], ttyoff:true)$
  27
  28 listarith: false$
  29 doallmxops: false$
  30 domxmxops: false$
  31
  32 dotexptsimp: false$
  33 /* dotexptsimp: false avoids the following (same for matrices):
  34 (%i1) (listarith: false, x: [1,2].[1,2], listarith: true, ev(x));
  35 (%o1)                                [ 5 ]
  36 */
  37
  38 domxplus: false$ /* default */
  39 scalarmatrixp: true$ /* default */
  40 mx0simp: true$ /* default, if false then 0*matrix([1,2]); --> 0 */
  41
  42 /* 1. define mtimesq (vector_rebuild.lisp)
  43    2. set listarith & friends
  44    3. then we are ready for tellsimpafter
  45 */
  46 ttyoff: true$
  47 lsmx_or_abs1_p(x):=
  48   is( listp(x) or matrixp(x) or is(equal(abs(x),1))=true )$
  49 matchdeclare(
  50   nlsmx, lambda([x], not lsmx_or_abs1_p(x)),
  51   tmsq, ?mtimesqp,
  52   lsmx, lambda([x],listp(x) or matrixp(x)) )$
  53 simp:false$
  54 tellsimpafter(
  55   nlsmx*tmsq,
  56     buildq([nlm:nlsmx, tq1:args(tmsq)[1], tq2:args(tmsq)[2]],
  57       nlm*tq1 ?mtimesq tq2) )$
  58 tellsimpafter(
  59   nlsmx*lsmx, buildq([nlm:nlsmx, lm:lsmx], nlm ?mtimesq lm) )$
  60 tellsimpafter(
  61   (-1)*nlsmx*lsmx, buildq([nlm:nlsmx, lm:lsmx], - nlm ?mtimesq lm) )$
  62 tellsimpafter(
  63   (-1)*lsmx, buildq([lm:lsmx], - lm) )$
  64 simp:true$
  65 ttyoff: false$
  66
  67 /* 2 lines copied from share/matrix/eigen.mac: */
  68 columnvector(x):=transpose(matrix(x))$
  69 covect(x):=columnvector(x)$
  70
  71 eval_when([batch,demo,load], ttyoff:false)$