Problem with gramschmidt and indexed matrix elements
- Subject: Problem with gramschmidt and indexed matrix elements
- From: Stavros Macrakis
- Date: Wed, 5 Jan 2011 10:45:10 -0500
This appears to be a bug in the general simplifier. Here's a simpler case:
%i85) display2d:false;
(%o85) false
(%i86) listarith:true;
(%o86) true
(%i87) a*[x,y];
(%o87) [a*x,a*y]
(%i88) a[1]*[x,y];
(%o88) a[1]*[x,y]
On Wed, Jan 5, 2011 at 10:32, Joshua Stults <joshua.stults at gmail.com> wrote:
> R : apply(matrix, makelist(makelist(concat(concat(r,i),j),i,1,3),j,1,3)) $
> rsubs :
> map("=",create_list(concat(concat(r,i),j),i,1,3,j,1,3),create_list(r[i,j],i,1,3,j,1,3))
> $
> R_indexed : subst(rsubs, R) $
> load("eigen") $
> gramschmidt(R);
> gramschmidt(R_indexed);
>