You can subscribe to this list here.
2009 |
Jan
(2) |
Feb
(5) |
Mar
|
Apr
|
May
(2) |
Jun
(8) |
Jul
(4) |
Aug
|
Sep
|
Oct
(2) |
Nov
(6) |
Dec
|
---|---|---|---|---|---|---|---|---|---|---|---|---|
2010 |
Jan
(1) |
Feb
(1) |
Mar
(3) |
Apr
(2) |
May
(2) |
Jun
(2) |
Jul
(18) |
Aug
(13) |
Sep
(7) |
Oct
|
Nov
|
Dec
(2) |
2011 |
Jan
|
Feb
(11) |
Mar
|
Apr
(4) |
May
|
Jun
(1) |
Jul
(18) |
Aug
(16) |
Sep
(12) |
Oct
(12) |
Nov
(19) |
Dec
(42) |
2012 |
Jan
(16) |
Feb
(3) |
Mar
(8) |
Apr
(14) |
May
(30) |
Jun
(5) |
Jul
(7) |
Aug
(3) |
Sep
(10) |
Oct
(4) |
Nov
(10) |
Dec
(1) |
2013 |
Jan
(14) |
Feb
(8) |
Mar
(5) |
Apr
(3) |
May
(9) |
Jun
(19) |
Jul
|
Aug
(27) |
Sep
(5) |
Oct
(18) |
Nov
(12) |
Dec
(8) |
2014 |
Jan
(5) |
Feb
(8) |
Mar
(20) |
Apr
(22) |
May
(28) |
Jun
(9) |
Jul
(6) |
Aug
(46) |
Sep
(40) |
Oct
(15) |
Nov
(8) |
Dec
(34) |
2015 |
Jan
(20) |
Feb
(15) |
Mar
(18) |
Apr
(20) |
May
(3) |
Jun
(13) |
Jul
(10) |
Aug
(19) |
Sep
(8) |
Oct
(31) |
Nov
(26) |
Dec
(13) |
2016 |
Jan
(13) |
Feb
(4) |
Mar
(14) |
Apr
(28) |
May
(19) |
Jun
(7) |
Jul
(1) |
Aug
|
Sep
(19) |
Oct
(5) |
Nov
(4) |
Dec
(9) |
2017 |
Jan
(4) |
Feb
(30) |
Mar
|
Apr
(5) |
May
(1) |
Jun
(1) |
Jul
(3) |
Aug
(2) |
Sep
(11) |
Oct
(3) |
Nov
(1) |
Dec
(6) |
2018 |
Jan
(5) |
Feb
(12) |
Mar
(5) |
Apr
(12) |
May
(22) |
Jun
(86) |
Jul
(7) |
Aug
(5) |
Sep
(6) |
Oct
(17) |
Nov
(1) |
Dec
(3) |
2019 |
Jan
(17) |
Feb
(4) |
Mar
(2) |
Apr
(7) |
May
(1) |
Jun
(2) |
Jul
(7) |
Aug
(9) |
Sep
|
Oct
(11) |
Nov
(20) |
Dec
(24) |
2020 |
Jan
(13) |
Feb
(1) |
Mar
(9) |
Apr
(2) |
May
(6) |
Jun
(6) |
Jul
(4) |
Aug
(2) |
Sep
(4) |
Oct
(1) |
Nov
(2) |
Dec
(6) |
2021 |
Jan
(10) |
Feb
(49) |
Mar
(26) |
Apr
(2) |
May
(1) |
Jun
|
Jul
(4) |
Aug
(6) |
Sep
|
Oct
(8) |
Nov
(5) |
Dec
(11) |
2022 |
Jan
|
Feb
|
Mar
(14) |
Apr
(19) |
May
(14) |
Jun
(4) |
Jul
|
Aug
|
Sep
(6) |
Oct
(4) |
Nov
|
Dec
(1) |
2023 |
Jan
|
Feb
(4) |
Mar
(6) |
Apr
|
May
|
Jun
(6) |
Jul
|
Aug
|
Sep
(13) |
Oct
(1) |
Nov
|
Dec
(16) |
2024 |
Jan
(66) |
Feb
(13) |
Mar
(5) |
Apr
(4) |
May
|
Jun
|
Jul
|
Aug
|
Sep
|
Oct
|
Nov
|
Dec
|
From: Rainer S. <rai...@gm...> - 2015-04-22 18:58:13
|
On Wed, 22 Apr 2015 at 13:53 -0400, Kostas Oikonomou wrote: > Thanks. I don't pretend to understand what is going on, but I also > discovered this: > > 1: z := exp(2*x)/x^2; > > 2*x > e > z := ------ > 2 > x > > 2: limit(z,x,infinity); > > 0 Oh dear, the limit code needs a thorough code review. Also this particular problem can be traced to this one: 1: y := x^2*e^(1/x^2); 2 1/x 2 y := e *x 2: sub(x=0,y); 0 Rainer |
From: Kostas O. <ko...@re...> - 2015-04-22 17:53:52
|
Thanks. I don't pretend to understand what is going on, but I also discovered this: 1: z := exp(2*x)/x^2; 2*x e z := ------ 2 x 2: limit(z,x,infinity); 0 On 04/22/2015 12:45, Rainer Schöpf wrote: > On Tue, 21 Apr 2015 at 15:26 -0400, Kostas Oikonomou wrote: > > > I just found the following problem with integrals: > > > > Reduce (Free CSL version), 06-Apr-15 ... > > > > 1: int(e^(-x^2),x,-infinity,infinity); > > > > sqrt(pi) > > > > > > 2: int(e^(-(x-mu)^2/(2*sigma^2)),x,-infinity,infinity); > > > > 2 > > ***** sigma invalid as number > > > It's more a problem with the limits package. Try this: > > y := e^(-(x-mu)^2/(2*sigma^2)); > > z := df(y,x)/x; > > limit(z,x,infinity); > > on complex; > > limit(z,x,infinity); > > Rainer |
From: Hearn, A. <he...@ra...> - 2015-04-22 17:17:39
|
-----Original Message----- From: scarpaflorent [mailto:sca...@ip...] Sent: Tuesday, April 07, 2015 10:40 AM To: re...@ra... Subject: A question about factorization in Reduce To whom it may concern, My name is Florent Scarpa, and I am using Reduce to realize calculations in physics. It seems very efficient, but I have a problem that really bothers me. It is about factorizing my expressions ; indeed, all functions doing factorizations, substitutions, limits, Taylor developments... only accept kernels as argument to be used as the common factor to find in expressions. My problem is that I would need to be able to factorize with a factor which is a product, or quotient, of kernels. Unfortunately, I did not find any solution in the user's manual, and I can not find a way to replace my product of kernels by a simple one; the consequence is that I am always confronted to the error message : "... invalid as kernel". Would you have any solution to propose for this problem ? Thank you in advance. Yours faithfully, Florent Scarpa __________________________________________________________________________ This email message is for the sole use of the intended recipient(s) and may contain confidential information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not the intended recipient, please contact the sender by reply email and destroy all copies of the original message. |
From: Rainer S. <rai...@gm...> - 2015-04-22 16:46:09
|
On Tue, 21 Apr 2015 at 15:26 -0400, Kostas Oikonomou wrote: > I just found the following problem with integrals: > > Reduce (Free CSL version), 06-Apr-15 ... > > 1: int(e^(-x^2),x,-infinity,infinity); > > sqrt(pi) > > > 2: int(e^(-(x-mu)^2/(2*sigma^2)),x,-infinity,infinity); > > 2 > ***** sigma invalid as number It's more a problem with the limits package. Try this: y := e^(-(x-mu)^2/(2*sigma^2)); z := df(y,x)/x; limit(z,x,infinity); on complex; limit(z,x,infinity); Rainer |
From: Kostas O. <ko...@re...> - 2015-04-21 19:27:34
|
I just found the following problem with integrals: Reduce (Free CSL version), 06-Apr-15 ... 1: int(e^(-x^2),x,-infinity,infinity); sqrt(pi) 2: int(e^(-(x-mu)^2/(2*sigma^2)),x,-infinity,infinity); 2 ***** sigma invalid as number 2 ***** sigma invalid as number 2 ***** sigma invalid as number Kostas |
From: Rainer S. <rai...@gm...> - 2015-04-17 10:39:26
|
On Thu, 9 Apr 2015 at 08:31 -0000, Andrea Magnoni wrote: > Somebody are able to compile correctly this package written for a old version of Reduce ? > package: http://www.kfki.hu/~kadlec/sw/riccir/riccir.tar.gz > instructions: http://wwwold.kfki.hu/cnc/szhkpub/riccir/riccir.html The package build process has changed in the years since the riccir package was written. You might try this: 1. In subsmsq.red, move the powlis1!* variable from the global to the fluid declaration. 2. In ini, comment out the first in statement: %in "riccir.red"$ Regards, Rainer |
From: Andrea M. <and...@ya...> - 2015-04-09 12:15:28
|
He doesn't know if the package is still operating with the newer version of Reduce "That's a very old package for Reduce and I have no idea whether it still works or not. You can download it from http://www.kfki.hu/~kadlec/sw/riccir/riccir.tar.gz Best regards, Jozsef -- E-mail : kad...@wi... PGP key: http://www.kfki.hu/~kadlec/pgp_public_key.txt Address: Wigner Research Centre for Physics, Hungarian Academy of Sciences H-1525 Budapest 114, POB. 49, Hungary" Il Giovedì 9 Aprile 2015 10:31, Andrea Magnoni <and...@ya...> ha scritto: Somebody are able to compile correctly this package written for a old version of Reduce ? package: http://www.kfki.hu/~kadlec/sw/riccir/riccir.tar.gz instructions: http://wwwold.kfki.hu/cnc/szhkpub/riccir/riccir.html Regards Andrea Magnoni |
From: Eberhard S. <esc...@ca...> - 2015-04-09 10:18:01
|
It is probably the best to contact the author of RICCIR. To my knowledge Jozsef can be reached at the addresses above. Eberhard On 04/09/2015 10:31 AM, Andrea Magnoni wrote: > Somebody are able to compile correctly this package written for a old > version of Reduce ? > > package: http://www.kfki.hu/~kadlec/sw/riccir/riccir.tar.gz > <http://www.kfki.hu/%7Ekadlec/sw/riccir/riccir.tar.gz> > > instructions: http://wwwold.kfki.hu/cnc/szhkpub/riccir/riccir.html > > > Regards > > Andrea Magnoni > > > ------------------------------------------------------------------------------ > BPM Camp - Free Virtual Workshop May 6th at 10am PDT/1PM EDT > Develop your own process in accordance with the BPMN 2 standard > Learn Process modeling best practices with Bonita BPM through live exercises > http://www.bonitasoft.com/be-part-of-it/events/bpm-camp-virtual- event?utm_ > source=Sourceforge_BPM_Camp_5_6_15&utm_medium=email&utm_campaign=VA_SF > > > _______________________________________________ > Reduce-algebra-developers mailing list > Red...@li... > https://lists.sourceforge.net/lists/listinfo/reduce-algebra-developers |
From: Andrea M. <and...@ya...> - 2015-04-09 08:32:24
|
Somebody are able to compile correctly this package written for a old version of Reduce ? package: http://www.kfki.hu/~kadlec/sw/riccir/riccir.tar.gz instructions: http://wwwold.kfki.hu/cnc/szhkpub/riccir/riccir.html Regards Andrea Magnoni |
From: Raymond R. <ray...@gm...> - 2015-04-07 14:38:39
|
I asked on the TeXmacs site but they thought that I should ask here. First I would like to say that the the texmacs/reduce version works very well. There is one small thing with the TeXmacs display: with "on zb_trace" some of the long outputs terminate prematurely. Copying the same input into the reduce monitor app continues the output to an answer. In fact the premature termination in texmacs is erratic and usually after three or four tries sometimes displays the whole output. Attached is my particular program; if you like I will try to shrink it but I don't think the problem is with the input; but some random thing in the display module. You can switch the problem on and off by switching "on zb_tracs" and "off zb_trace" at the first line. The problem is with the output of the last line 719 of input. Thanks -- Two views on life: life is an art not to be learned by observation. George Santayana:Interpretations of Poetry and Religion It's kinda nice to participate in your life Raymond Rogers |
From: Raffaele V. <raf...@un...> - 2015-04-07 13:55:53
|
On 07/04/15 13:44, Alan Barnes wrote: > Recently I have experienced a few problems with the handling of large > amounts of output in the GUI version of redcsl. > > If output is directed to a file there is no problem. However, problems > occur with direct output. When scrolling back the first part of the > output is corrupt: wrong values and/or raw TeX displayed. This often > leads to a subsequent crash. > > The problems occur on Ubuntu with a fairly recent version of redcsl > (mid-March) -- I haven't tried redcsl on Windows to see if the same > problem occurs. No problems with rfcsl or with PSL-based versions as > far as I can see. > > One way to reproduce this problem fairly predictably (but not to crash > the system immediately) is to load the file 2015-04-07-invfn3.tst from > packages/regressions and then scroll back to the beginning of the output. > > Has anyone experienced similar problems? I tried your file but I do not have the same problem (at least for the moment). However, the CSL GUI has some known bugs, eg: if you try to display a 'big' matrix it will refuse with an error; cut-and-paste menu items are broken; etc. I am trying to help Arthur Norman to spot the bugs, I suggest you to ask him directly. Best regards, Raffaele Vitolo. > Alan Barnes > > > ------------------------------------------------------------------------------ > BPM Camp - Free Virtual Workshop May 6th at 10am PDT/1PM EDT > Develop your own process in accordance with the BPMN 2 standard > Learn Process modeling best practices with Bonita BPM through live exercises > http://www.bonitasoft.com/be-part-of-it/events/bpm-camp-virtual- event?utm_ > source=Sourceforge_BPM_Camp_5_6_15&utm_medium=email&utm_campaign=VA_SF > _______________________________________________ > Reduce-algebra-developers mailing list > Red...@li... > https://lists.sourceforge.net/lists/listinfo/reduce-algebra-developers > |
From: Alan B. <ala...@gm...> - 2015-04-07 11:44:52
|
Recently I have experienced a few problems with the handling of large amounts of output in the GUI version of redcsl. If output is directed to a file there is no problem. However, problems occur with direct output. When scrolling back the first part of the output is corrupt: wrong values and/or raw TeX displayed. This often leads to a subsequent crash. The problems occur on Ubuntu with a fairly recent version of redcsl (mid-March) -- I haven't tried redcsl on Windows to see if the same problem occurs. No problems with rfcsl or with PSL-based versions as far as I can see. One way to reproduce this problem fairly predictably (but not to crash the system immediately) is to load the file 2015-04-07-invfn3.tst from packages/regressions and then scroll back to the beginning of the output. Has anyone experienced similar problems? Alan Barnes |
From: Kostas O. <ko...@re...> - 2015-03-31 18:44:28
|
Great, thanks very much. Kostas On 03/31/2015 03:27, Rainer Schöpf wrote: > > > > 6: rank matrix_augment(A,B); > > > > +++ Error undefined function (2 args): matrix_augment > > I believe I corrected this problem as well. Changes are committed. > > Rainer |
From: Rainer S. <rai...@gm...> - 2015-03-31 07:27:29
|
Hi Kostas, > Rainer, thank you. This is what happens: > > Reduce (Free CSL version), 23-Mar-15 ... > > 1: load_package linalg; > > 2: put('matrix_augment,'rtypefn,'quotematrix); > > quotematrix > > 3: A := mat((a),(b))$ B := mat((c),(d))$ > > 5: tp matrix_augment(A,B); > > [a b] > [c d] > > 6: rank matrix_augment(A,B); > > +++ Error undefined function (2 args): matrix_augment I believe I corrected this problem as well. Changes are committed. Rainer |
From: Kostas O. <ko...@re...> - 2015-03-30 22:57:58
|
Rainer, thank you. This is what happens: Reduce (Free CSL version), 23-Mar-15 ... 1: load_package linalg; 2: put('matrix_augment,'rtypefn,'quotematrix); quotematrix 3: A := mat((a),(b))$ B := mat((c),(d))$ 5: tp matrix_augment(A,B); [a b] [c d] 6: rank matrix_augment(A,B); +++ Error undefined function (2 args): matrix_augment 7: Kostas On 03/30/2015 04:59, Rainer Schöpf wrote: > On Sun, 29 Mar 2015 at 19:08 -0400, Kostas Oikonomou wrote: > > > I would expect that I could write > > > > 39: x := mat((a),(b))$ > > > > 40: y := mat((c),(d))$ > > > > and then ask for the transpose: > > > > 41: tp matrix_augment(x,y); > > > > ***** tp(matrix_augment(x,y)) invalid as scalar > > Ah, well, that's one of the reasons why linalg is listed as contributed > package... > > Seriously, please try this > > put('matrix_augment,'rtypefn,'quotematrix); > > (and likewise for matrix_stack). > > If this works to your satisfaction, I'll commit the change. > > Rainer |
From: Tony R. <ant...@ad...> - 2015-03-30 19:16:22
|
-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi Rainer, many thanks, it appears to be working now for my programs, Tony On 30/03/2015 5:54 pm, Rainer Schöpf wrote: > Hi Tony, > > > have been using a 2012 version of Reduce, but just upgraded and now > > all my let rules invoking expressions stored in arrays break. The > > simplest example is the following code > > > array u(9); > > operator v; > > let v(~p)=>u(p); > > end; > > > which generates the output > > > 1: in "temp.red"; > > array u(9); > > operator v; > > let v(~p)=>u(p); > > > ***** ~p invalid as array index > > Cont? (Y or N) > > > Is this a bug to be corrected? or do I have to change all my code? > > I think it is a consequence of a change made in August 2014, to properly handle lists as array elements. > > I'm just committed a correction to allow free rule variables in array indices on the right side of rules. Your example does no longer produce the error; please test. > > Rainer > > > > > > I cannot find anything in the manual saying that the construct is > > invalid---the manual says that arrays on the left-hand side are > > evaluated, but my arrays are on the right-hand side. > > > Tony > > > Rainer Schöpf - -- - --------------------------------------------------------------------- Professor A.J. Roberts School of Mathematical Sciences phone: +61 8 8313 3035 University of Adelaide fax: +61 8 8313 3696 South Australia 5005. mailto:ant...@ad... http://www.maths.adelaide.edu.au/anthony.roberts/ ==.0000001000000100000110001000011010001111110010111011101000010000== -----BEGIN PGP SIGNATURE----- Version: GnuPG v1 Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/ iEYEARECAAYFAlUZoPcACgkQ7TX8dTbro1vs6QCgjXYOba4qisFzlpFjfK7T1bSK ruAAnj6njoc7o1n8v4gV6sgZFC/9knKw =Ab91 -----END PGP SIGNATURE----- |
From: Rainer S. <rai...@gm...> - 2015-03-30 08:59:27
|
On Sun, 29 Mar 2015 at 19:08 -0400, Kostas Oikonomou wrote: > I would expect that I could write > > 39: x := mat((a),(b))$ > > 40: y := mat((c),(d))$ > > and then ask for the transpose: > > 41: tp matrix_augment(x,y); > > ***** tp(matrix_augment(x,y)) invalid as scalar Ah, well, that's one of the reasons why linalg is listed as contributed package... Seriously, please try this put('matrix_augment,'rtypefn,'quotematrix); (and likewise for matrix_stack). If this works to your satisfaction, I'll commit the change. Rainer |
From: Rainer S. <rai...@gm...> - 2015-03-30 07:24:12
|
Hi Tony, > have been using a 2012 version of Reduce, but just upgraded and now > all my let rules invoking expressions stored in arrays break. The > simplest example is the following code > > array u(9); > operator v; > let v(~p)=>u(p); > end; > > which generates the output > > 1: in "temp.red"; > array u(9); > operator v; > let v(~p)=>u(p); > > ***** ~p invalid as array index > Cont? (Y or N) > > Is this a bug to be corrected? or do I have to change all my code? I think it is a consequence of a change made in August 2014, to properly handle lists as array elements. I'm just committed a correction to allow free rule variables in array indices on the right side of rules. Your example does no longer produce the error; please test. Rainer > > I cannot find anything in the manual saying that the construct is > invalid---the manual says that arrays on the left-hand side are > evaluated, but my arrays are on the right-hand side. > > Tony > > -- > --------------------------------------------------------------------- > Professor A.J. Roberts > School of Mathematical Sciences phone: +61 8 8313 3035 > University of Adelaide fax: +61 8 8313 3696 > South Australia 5005. mailto:ant...@ad... > http://www.maths.adelaide.edu.au/anthony.roberts/ > ==.0000001000000100000110001000011010001111110010111011101000010000== > > [ PGP Signature check FAILED - Mo 30. Mär 08:58:04 CEST 2015 ] > Rainer Schöpf |
From: Tony R. <ant...@ad...> - 2015-03-30 06:00:50
|
-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi all, have been using a 2012 version of Reduce, but just upgraded and now all my let rules invoking expressions stored in arrays break. The simplest example is the following code array u(9); operator v; let v(~p)=>u(p); end; which generates the output 1: in "temp.red"; array u(9); operator v; let v(~p)=>u(p); ***** ~p invalid as array index Cont? (Y or N) Is this a bug to be corrected? or do I have to change all my code? I cannot find anything in the manual saying that the construct is invalid---the manual says that arrays on the left-hand side are evaluated, but my arrays are on the right-hand side. Tony - -- - --------------------------------------------------------------------- Professor A.J. Roberts School of Mathematical Sciences phone: +61 8 8313 3035 University of Adelaide fax: +61 8 8313 3696 South Australia 5005. mailto:ant...@ad... http://www.maths.adelaide.edu.au/anthony.roberts/ ==.0000001000000100000110001000011010001111110010111011101000010000== -----BEGIN PGP SIGNATURE----- Version: GnuPG v1 Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/ iEYEARECAAYFAlUY5oIACgkQ7TX8dTbro1v3dACfcRdgOc69ATIWQssmSr5RNwZJ mQcAn2d93zNROf46N6n12Snv2xGjcozy =5Nys -----END PGP SIGNATURE----- |
From: Kostas O. <ko...@re...> - 2015-03-29 23:08:23
|
I would expect that I could write 39: x := mat((a),(b))$ 40: y := mat((c),(d))$ and then ask for the transpose: 41: tp matrix_augment(x,y); ***** tp(matrix_augment(x,y)) invalid as scalar 42: But this doesn't work. However, this works 48: tp (M := matrix_augment(A,B)); [a b] [c d] but I don't want to have to create an extra matrix. Apparently matrix_augment doesn't quite return a matrix? Same problem occurs with matrix_stack(). Kostas Kostas |
From: Raymond R. <ray...@gm...> - 2015-03-11 20:52:24
|
Hello, I joined this list because I wanted to do some integrations over multi-variable hypergeometric functions. To make a long story short; none of the programs on sourceforge would work right. The precompiled ones would bring up a window/terminal but wouldn't allow me to edit directory names or file names. The compilation process also failed to work. But! the svn svn co http://svn.code.sf.net/p/reduce-algebra/code/trunk reduce-algebra ./configure --with-csl make worked flawlessly. Rather than submit a bug report or some such I thought that you would like to just update the sourceforge files? Of course if you want I can send log files, straces, and such; but since a solution is in hand I thought it would be a waste of time. Ubuntu 14.04 (trusty) Gnome 3.8.3 Kernel 3.13.0-46-generic AMD dual proccessor etc..... BTW: the latex display is great! And it looks like the Gosper/Zeilberger packages will make my work easy (thank you). Ray -- Two views on life: life is an art not to be learned by observation. George Santayana:Interpretations of Poetry and Religion It's kinda nice to participate in your life Raymond Rogers |
From: Rainer S. <rai...@gm...> - 2015-03-05 18:35:07
|
Here's my proposal for the trigonometric rules, with sin as example for all others: Replace sin( (~~w + ~~k*pi)/~~d) => (if evenp fix(k/d) then 1 else -1) * sin((w + remainder(k,d)*pi)/d) when w freeof pi and ratnump(k/d) and fixp k and abs(k/d) >= 1, by sin( (~~w + ~~k*pi)/~~d) => (if evenp fix(k/d) then 1 else -1) * sin(w/d + ((k/d)-fix(k/d))*pi) when w freeof pi and ratnump(k/d) and abs(k/d) >= 1, and sin((~x + ~~k*pi)/~d) => sign(k/d)*cos(x/d) when x freeof pi and abs(k/d) = 1/2, by sin((~x + ~~k*pi)/~~d) => sign(k/d)*cos(x/d) when x freeof pi and abs(k/d) = 1/2, This makes sure that the rules match with both OFF and ON RATIONAL. Rainer |
From: Rainer S. <rai...@gm...> - 2015-03-04 18:18:19
|
On Wed, 4 Mar 2015 at 12:17 -0000, Alan Barnes wrote: > The simplification of expressions such as sin(x+5*pi/2) to cos x fails > after ON RATIONAL. Expressions with any rational multiple of PI are not > folded towards zero, but the simplification of expressions involving > integer multiples of PI is not affected. OFF MCD has much the same effect. > > It seems that the pattern matching in opmtch is messed up. It's not really messed up. The problem is this rule: sin( (~~w + ~~k*pi)/~~d) => (if evenp fix(k/d) then 1 else -1) * sin((w + remainder(k,d)*pi)/d) when w freeof pi and ratnump(k/d) and fixp k and abs(k/d) >= 1, First of all, ratnump doesn't work correctly with ON RATIONAL, it returns nil when given a rational number. I'll commit a correction later. Secondly, with ON RATIONAL, the above rule would match with d=1 and k a rational number, which violates the condition (fixp k being false). A second rule needs to be defined for this case: sin ( ~~w + ~~r*pi) => (if evenp fix(r) then 1 else -1) * sin((w + (r-fix(r))*sign(r)*pi)) when w freeof pi and ratnump r and abs(r) >=1; Note the use of (r-fix(r)) instead of remainder. It is certainly possible to fold the two rules into one. Rainer |
From: Alan B. <ala...@gm...> - 2015-03-04 12:17:45
|
The simplification of expressions such as sin(x+5*pi/2) to cos x fails after ON RATIONAL. Expressions with any rational multiple of PI are not folded towards zero, but the simplification of expressions involving integer multiples of PI is not affected. OFF MCD has much the same effect. It seems that the pattern matching in opmtch is messed up. Is this a known feature? I couldn't immediately see an easy way of generalising the relevant substitution rule to cater for the case when RATIONAL is on, but my knowledge of the pattern matcher is currently rather limited. Whilst investigating the above feature, I noticed that in CSL-based versions that I couldn't arrange to simultaneously trace calls to simpiden and opmtch. In PSL-based versions the tracing worked fine. Alan Barnes |
From: Rainer S. <rai...@gm...> - 2015-03-03 06:50:30
|
On Mon, 2 Mar 2015 at 15:44 -0000, Alan Barnes wrote: > Currently Reduce algebraic mode does not (as far as I can see) support > trig functions with arguments in degrees, although symbolic mode does! > > However, there is a lot of code in arith/rdelem.red to support such > functionality which is currently commented out. This code seems to work > well although presumably there is some (perhaps subtle) reason why it > was commented out (back in Reduce 3.7 days at least). Looking at the historical code I find that it was commented out when the bigfloat package was converted from base 10 to base 2 arithmetic. I cannot remember why; I suspect this was done since sind and friends are not defined in algebraic mode. In my opinion, defining them is worthwhile. Rainer |