Looks like a pari bug to me since their docs suggest that flag=1 and flag=4 have the same output, just a slightly different algorithm, and flag =1 does not include the columns of 0s.
John On 23 January 2012 19:01, Volker Braun <vbraun.n...@gmail.com> wrote: > For the record, its clear that pari4 does output the unwanted zero rows > (columns in pari), but is this a bug or documented somewhere in pari? > > > [vbraun@volker-laptop-two keys]$ sage -gp > > GP/PARI CALCULATOR Version 2.5.0 (development svn-13228) > > amd64 running linux (x86-64/GMP-5.0.1 kernel) 64-bit version > > compiled: Jan 18 2012, gcc-4.6.2 20111027 (Red Hat 4.6.2-1) (GCC) > > (readline v6.2 enabled, extended help enabled) > > > Copyright (C) 2000-2011 The PARI Group > > PARI/GP is free software, covered by the GNU General Public License, and > comes WITHOUT ANY WARRANTY WHATSOEVER. > > Type ? for help, \q to quit. > Type ?12 for how to get moral (and possibly technical) support. > > parisize = 8000000, primelimit = 500509 > <1, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 1, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 1, 0, 0, > 0, 0; 0, 0, 0, 0, 0, 0, 1, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0, 1, 0, 0; 0, 0, 0, > 0, 0, 0, 0, 0, 1, 0; 0, 0, 0, 0, 0, 0, 0, 0, 0, 1] > %1 = > [1 1 -1 -5 0 -1 0 1 1 0] > [-1 1 0 3 12 2 0 0 2 -2] > [-1 0 -2 -4 -1 0 -1 6 1 3] > [-19 0 -1 -3 -2 -7 -2 0 1 1] > [1 1 1 1 -12 -14 -3 3 -3 -1] > [1 0 0 0 0 0 0 0 0 0] > [0 1 0 0 0 0 0 0 0 0] > [0 0 1 0 0 0 0 0 0 0] > [0 0 0 1 0 0 0 0 0 0] > [0 0 0 0 1 0 0 0 0 0] > [0 0 0 0 0 1 0 0 0 0] > [0 0 0 0 0 0 1 0 0 0] > [0 0 0 0 0 0 0 1 0 0] > [0 0 0 0 0 0 0 0 1 0] > [0 0 0 0 0 0 0 0 0 1] > > ? mathnf(m~,flag=4)[1] > %2 = > [0 0 0 0 0 1 0 0 0 0 0 0 0 0 0] > [0 0 0 0 0 0 1 0 0 0 0 0 0 0 0] > [0 0 0 0 0 0 0 1 0 0 0 0 0 0 0] > [0 0 0 0 0 0 0 0 1 0 0 0 0 0 0] > [0 0 0 0 0 0 0 0 0 1 0 0 0 0 0] > [0 0 0 0 0 0 0 0 0 0 1 0 0 0 0] > [0 0 0 0 0 0 0 0 0 0 0 1 0 0 0] > [0 0 0 0 0 0 0 0 0 0 0 0 1 0 0] > [0 0 0 0 0 0 0 0 0 0 0 0 0 1 0] > [0 0 0 0 0 0 0 0 0 0 0 0 0 0 1] > > -- > To post to this group, send an email to sage-devel@googlegroups.com > To unsubscribe from this group, send an email to > sage-devel+unsubscr...@googlegroups.com > For more options, visit this group at > http://groups.google.com/group/sage-devel > URL: http://www.sagemath.org -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
