Re: [fricas-devel] better integration result -- avoid complicated roots

2023-11-19 Thread Waldek Hebisch
On Mon, Nov 20, 2023 at 08:48:30AM +0800, Qian Yun wrote: > > > On 11/20/23 02:19, Waldek Hebisch wrote: > > > > "atan doubling" makes thing harder (increases degrees) and after > > doubling Rioboo was failing. Doubling may happen both before > > irexpand and after. > > > > What is "atan

Re: [fricas-devel] better integration result -- avoid complicated roots

2023-11-19 Thread Waldek Hebisch
On Sun, Nov 19, 2023 at 07:17:14PM +0800, Qian Yun wrote: > > The core idea is that for polynomials with real coefficient, > its roots come in pairs (symmetric over the real axis). > So in the following computation of the rootSum of ilogs, > many terms can cancel with each other. Now, some

Re: [fricas-devel] better integration result -- avoid complicated roots

2023-11-19 Thread Qian Yun
On 11/20/23 02:19, Waldek Hebisch wrote: "atan doubling" makes thing harder (increases degrees) and after doubling Rioboo was failing. Doubling may happen both before irexpand and after. What is "atan doubling"? There are several changes in 'mapleok.input'. Did you look that they all

Re: [fricas-devel] real? sqrt (-sqrt 2)

2023-11-19 Thread Qian Yun
Does RealClosure work with Expression? (e.g doing integration) Currently we can't do so in FriCAS. Is that possible in theory? - Qian On 11/19/23 23:12, Waldek Hebisch wrote: On Sun, Nov 19, 2023 at 10:50:17PM +0800, oldk1331 wrote: Then is there a way to specify the algebraic number which

Re: [fricas-devel] better integration result -- avoid complicated roots

2023-11-19 Thread Waldek Hebisch
On Sun, Nov 19, 2023 at 08:47:26PM +0800, Qian Yun wrote: > > > On 11/19/23 20:28, Waldek Hebisch wrote: > I was talking about section 2.8 "Riboo's Algorithm for Real Rational > Functions", page 65 to 70 for second edition of "Symbolic > Integration I". > > > AFAICS most significant reason for

Re: [fricas-devel] real? sqrt (-sqrt 2)

2023-11-19 Thread Waldek Hebisch
On Sun, Nov 19, 2023 at 10:50:17PM +0800, oldk1331 wrote: > Then is there a way to specify the algebraic number > which square is 2 and is near +1.41421? (1) -> rC := RealClosure(Fraction(Integer)) (1) RealClosure(Fraction(Integer))

Re: [fricas-devel] real? sqrt (-sqrt 2)

2023-11-19 Thread oldk1331
Then is there a way to specify the algebraic number which square is 2 and is near +1.41421? - Qian -- You received this message because you are subscribed to the Google Groups "FriCAS - computer algebra system" group. To unsubscribe from this group and stop receiving emails from it, send an

Re: [fricas-devel] real? sqrt (-sqrt 2)

2023-11-19 Thread Waldek Hebisch
On Sun, Nov 19, 2023 at 09:49:10PM +0800, Qian Yun wrote: > > > On 11/19/23 21:36, Waldek Hebisch wrote: > > On Sun, Nov 19, 2023 at 07:58:50PM +0800, Qian Yun wrote: > > > "real? sqrt (-sqrt 2)" returns "true" by default. > > > > You mean should return false? We get this when working > > with

Re: [fricas-devel] real? sqrt (-sqrt 2)

2023-11-19 Thread Qian Yun
On 11/19/23 21:36, Waldek Hebisch wrote: On Sun, Nov 19, 2023 at 07:58:50PM +0800, Qian Yun wrote: "real? sqrt (-sqrt 2)" returns "true" by default. You mean should return false? We get this when working with expressions: (2) -> real? sqrt (-sqrt 2::EXPR(INT)) (2) false Because

Re: [fricas-devel] real? sqrt (-sqrt 2)

2023-11-19 Thread Waldek Hebisch
On Sun, Nov 19, 2023 at 07:58:50PM +0800, Qian Yun wrote: > "real? sqrt (-sqrt 2)" returns "true" by default. You mean should return false? We get this when working with expressions: (2) -> real? sqrt (-sqrt 2::EXPR(INT)) (2) false > Because interpreter prefers CTRIGMNP over TRIGMNIP > as

Re: [fricas-devel] better integration result -- avoid complicated roots

2023-11-19 Thread Qian Yun
On 11/19/23 20:28, Waldek Hebisch wrote: Until recently I found that Rioboo's algorithm described in Manuel Bronstein's book, at the end of section 2.8, covers this situation. Hmm, I am not sure what you want to say here. We use Rioboo's algorithm, but in some cases it does not give us

Re: [fricas-devel] better integration result -- avoid complicated roots

2023-11-19 Thread Waldek Hebisch
On Sun, Nov 19, 2023 at 07:17:14PM +0800, Qian Yun wrote: > Some of the integration results are pages long -- > because IntegrationResult often has a 'logpart', aka > a rootSum form. > > 'lg2func' in IR2F uses 'zeroOf' which can introduce > large radicals into the result. > > For example > >

[fricas-devel] real? sqrt (-sqrt 2)

2023-11-19 Thread Qian Yun
"real? sqrt (-sqrt 2)" returns "true" by default. Because interpreter prefers CTRIGMNP over TRIGMNIP as shown by ")set mes bot on": Function Selection for real? Arguments: AN [1] signature: EXPR(COMPLEX(INT)) -> BOOLEAN implemented: slot (Boolean)(Expression (Complex

Re: [fricas-devel] [PATCH] fix build with CLISP (no ffi)

2023-11-19 Thread Qian Yun
On 11/19/23 18:27, Waldek Hebisch wrote: On Sun, Nov 19, 2023 at 10:05:10AM +0800, Qian Yun wrote: Currently, the only FFI function used during build (without-x11) is |remove_directory| in |erase_lib|, and there has a Lisp fallback. So the autotools variable

[fricas-devel] better integration result -- avoid complicated roots

2023-11-19 Thread Qian Yun
Some of the integration results are pages long -- because IntegrationResult often has a 'logpart', aka a rootSum form. 'lg2func' in IR2F uses 'zeroOf' which can introduce large radicals into the result. For example internalIntegrate(1/(z^4+1),z) internalIntegrate((1+z^2)^4/ (1+z^4)^4,z)

Re: [fricas-devel] [PATCH] fix build with CLISP (no ffi)

2023-11-19 Thread Waldek Hebisch
On Sun, Nov 19, 2023 at 06:35:14PM +0800, Qian Yun wrote: > In your commit 055f3a40 (or svn r490), 2009-01-18, said > "Allow build even if Clisp lacks FFI.". Old mail says it > was needed for sage interface (back when they were using > CLISP instead of ECL, it seems). Yes, but for many years sage

Re: [fricas-devel] [PATCH] fix build with CLISP (no ffi)

2023-11-19 Thread Qian Yun
In your commit 055f3a40 (or svn r490), 2009-01-18, said "Allow build even if Clisp lacks FFI.". Old mail says it was needed for sage interface (back when they were using CLISP instead of ECL, it seems). Poplog doesn't support FFI, I've not tried it so I don't know if it will build. I plan to

Re: [fricas-devel] [PATCH] fix build with CLISP (no ffi)

2023-11-19 Thread Waldek Hebisch
On Sun, Nov 19, 2023 at 10:05:10AM +0800, Qian Yun wrote: > Currently, the only FFI function used during build (without-x11) > is |remove_directory| in |erase_lib|, and there has a Lisp fallback. > > So the autotools variable "fricas_has_directory_operations" > means we have those C functions,