From: abpetrov <abp...@uf...> - 2013-08-12 16:28:32
|
Hi, I just have tried make linear transformation for fermi operators. But after such transformation result has strange property. When I set switch factor on some sum's members dissapears. When I set switch factor off they appears again. I suppose this is a bad thing because in many situations it can lead to illegal interpretation calculation results. So, questions is 1) Why dissapears some sum's members 2) How can I avoid that situations? I can add, that when I made new fermi operators a2 and a2c (see below) as commutative variables they don't disappears. The program with non-commutative variables below. In program I have fermi operators a,ac with rules and a2 and a2c with rules. After that I defined operator H, in which I made linear transformation. Last two outputs of operator H2 is different because of switch factor on. off lower$ on intstr$ load_package noncom2; operator a,ac; noncom a,a; noncom ac,ac; noncom a,ac; for all i let a(i)*ac(i) = 1 - ac(i)*a(i); for all i,j such that i neq j let a(i)*ac(j) = -ac(j)*a(i); for all i,j such that ordp(i,j) let a(i)*a(j) = -a(j)*a(i); for all i,j such that ordp(i,j) let ac(i)*ac(j) = -ac(j)*ac(i); operator a2,a2c; noncom a2,a2; noncom a2c,a2c; noncom a2,a2c; for all i let a2(i)*a2c(i) = 1 - a2c(i)*a2(i); for all i,j such that i neq j let a2(i)*a2c(j) = -a2c(j)*a2(i); for all i,j such that ordp(i,j) let a2(i)*a2(j) = -a2(j)*a2(i); for all i,j such that ordp(i,j) let a2c(i)*a2c(j) = -a2c(j)*a2c(i); operator H; H := J1*ac(i)*a(i+1) + J1*ac(i+1)*a(i); operator U,UT,V,VT; sub_fermi := { a(~i) => U(i,j1)*a2(j1) + V(i,j1)*a2c(j1), ac(~i) => VT(i,j2)*a2(j2) + UT(i,j2)*a2c(j2) }; H2 := (H where sub_fermi); on factor; H2; Best regards, Petrov Alexander |
From: abpetrov <abp...@uf...> - 2013-08-13 15:57:08
|
> Hi, > I just have tried make linear transformation for fermi operators. > But after such transformation result has strange property. > When I set switch factor on some sum's members dissapears. > When I set switch factor off they appears again. I suppose this is a bad > thing because in many situations it can lead to illegal interpretation > calculation results. > So, questions is > 1) Why dissapears some sum's members > 2) How can I avoid that situations? > I can add, that when I made new fermi operators a2 and a2c (see below) > as commutative variables they don't disappears. > The program with non-commutative variables below. > In program I have fermi operators a,ac with rules and > a2 and a2c with rules. After that I defined operator H, > in which I made linear transformation. > Last two outputs of operator H2 is different because of switch > factor on. > > > off lower$ > on intstr$ > > load_package noncom2; > > operator a,ac; > noncom a,a; > noncom ac,ac; > noncom a,ac; > > for all i let a(i)*ac(i) = 1 - ac(i)*a(i); > for all i,j such that i neq j let a(i)*ac(j) = -ac(j)*a(i); > for all i,j such that ordp(i,j) let a(i)*a(j) = -a(j)*a(i); > for all i,j such that ordp(i,j) let ac(i)*ac(j) = -ac(j)*ac(i); > > operator a2,a2c; > noncom a2,a2; > noncom a2c,a2c; > noncom a2,a2c; > > for all i let a2(i)*a2c(i) = 1 - a2c(i)*a2(i); > for all i,j such that i neq j let a2(i)*a2c(j) = -a2c(j)*a2(i); > for all i,j such that ordp(i,j) let a2(i)*a2(j) = -a2(j)*a2(i); > for all i,j such that ordp(i,j) let a2c(i)*a2c(j) = -a2c(j)*a2c(i); > > > operator H; > H := J1*ac(i)*a(i+1) + J1*ac(i+1)*a(i); > > operator U,UT,V,VT; > > sub_fermi := { a(~i) => U(i,j1)*a2(j1) + V(i,j1)*a2c(j1), > ac(~i) => VT(i,j2)*a2(j2) + UT(i,j2)*a2c(j2) }; > H2 := (H where sub_fermi); > on factor; > H2; > > > Best regards, Petrov Alexander > > ------------------------------------------------------------------------------ > Get 100% visibility into Java/.NET code with AppDynamics Lite! > It's a free troubleshooting tool designed for production. > Get down to code-level detail for bottlenecks, with <2% overhead. > Download for free and get started troubleshooting in minutes. > http://pubads.g.doubleclick.net/gampad/clk?id=48897031&iu=/4140/ostg.clktrk > _______________________________________________ > Reduce-algebra-developers mailing list > Red...@li... > https://lists.sourceforge.net/lists/listinfo/reduce-algebra-developers > I just have rewritten result for H2 as input expression and switched on factor. Result same as in prior program. So, linear transformaition hasn't relation to that behaviour of Reduce. Possible, something wrong with non-commutative operators or with on factor working with non-commutative operators. The new program, which show different view H2 depending from factor switch is: off lower$ load_package noncom2; operator a2,a2c; noncom a2,a2; noncom a2c,a2c; noncom a2,a2c; operator U,UT,V,VT; H2 := J1*(U(i + 1,j1)*UT(i,j2)*a2c(j2)*a2(j1) + U(i + 1,j1)*VT(i,j2)*a2(j2)*a2(j1) + U (i,j1)*UT(i + 1,j2)*a2c(j2)*a2(j1) + U(i,j1)*VT(i + 1,j2)*a2(j2)*a2(j1) + UT(i + 1,j2)*V(i,j1)*a2c(j2)*a2c(j1) + UT(i,j2)*V(i + 1,j1)*a2c(j2)*a2c(j1) - V(i + 1, j1)*VT(i,j2)*a2c(j1)*a2(j2) - V(i,j1)*VT(i + 1,j2)*a2c(j1)*a2(j2)); on factor; H2; Where is 4 first summand? Best regards, Petrov Alexander |
From: Rainer S. <rai...@gm...> - 2013-08-14 18:28:13
|
Hello, I had a quick look on your problem. First of all, it does not occur when you remove the line off lower$ It seems to me that the error comes from the noncom2 package. This package was written 20 years ago for Reduce 3.3, and redefines some of the inner parts of Reduce. In particular, it changes the internal ordering of expressions, which doesn't work well with the current version. I'll look at this in the next few days. Rainer |
From: abpetrov <abp...@uf...> - 2013-08-15 07:50:36
|
thank you 14.08.2013 18:28, Rainer Schöpf пишет: > Hello, > > I had a quick look on your problem. > > First of all, it does not occur when you remove the line > > off lower$ > > It seems to me that the error comes from the noncom2 package. This package was > written 20 years ago for Reduce 3.3, and redefines some of the inner parts of > Reduce. In particular, it changes the internal ordering of expressions, which > doesn't work well with the current version. > > I'll look at this in the next few days. > > > Rainer > |
From: Rainer S. <rai...@gm...> - 2013-08-16 20:56:22
|
Hello, I believe I have corrected the error in the subversion repository. Please update your working copy and rebuild. Rainer |
From: Rainer S. <rai...@gm...> - 2013-08-17 09:13:18
|
On Sat, 17 Aug 2013 at 13:09 -0000, abpetrov wrote: > Hello, > Nothing changed. But I am not sure that, may be I got another repository? > I saw directory trunk. As I can see, old file scripts/findos.sh is in > repository without Gentoo, but it must be here according to our > correspondence. > I got repository according documentation with command > > svn co https://reduce-algebra.svn.sourceforge.net/svnroot/reduce-algebra reduce-algebra This is an old repository. sourceforge.net changed repositories a while ago, but didn't remove the old one. Where did you find this URL? The correct one is shown on the REDUCE download page http://reduce-algebra.sourceforge.net/downloading.html It is svn checkout http://svn.code.sf.net/p/reduce-algebra/code/trunk reduce-algebra Rainer |
From: abpetrov <abp...@uf...> - 2013-08-18 15:29:25
|
Now it works. Dissapearing sum's members is absents. Link for svn I got from the page http://reduce-algebra.com/downloading.htm Thank you. Best regards, Petrov Alexander. 17.08.2013 09:13, Rainer Schöpf пишет: > On Sat, 17 Aug 2013 at 13:09 -0000, abpetrov wrote: > > > Hello, > > Nothing changed. But I am not sure that, may be I got another repository? > > I saw directory trunk. As I can see, old file scripts/findos.sh is in > > repository without Gentoo, but it must be here according to our > > correspondence. > > I got repository according documentation with command > > > > svn co https://reduce-algebra.svn.sourceforge.net/svnroot/reduce-algebra reduce-algebra > > This is an old repository. sourceforge.net changed repositories a while ago, > but didn't remove the old one. > > Where did you find this URL? > > The correct one is shown on the REDUCE download page > > http://reduce-algebra.sourceforge.net/downloading.html > > It is > > svn checkout http://svn.code.sf.net/p/reduce-algebra/code/trunk reduce-algebra > > Rainer > |