maxima-bugs — Reports from the SF bug-tracking system

[Maxima-bugs] [maxima:bugs] #4636 signum(ind) is an error

[Barton W. <wil...@us...>]

One more reason: I wanted to make `signum(2.0b0) = 1.0b0, signum(2/3) = 1, ...`.  the silly code was an easy way to achieve that


---

**[bugs:#4636] signum(ind) is an error**

**Status:** open
**Group:** None
**Labels:** extended real signum limit 
**Created:** Thu Nov 27, 2025 12:41 PM UTC by Barton Willis
**Last Updated:** Fri Nov 28, 2025 01:21 PM UTC
**Owner:** nobody


Either a  `signum` nounform or `ind` is a better result than an error:
~~~
(%i4)	signum(ind);
sign: sign of ind is undefined.
~~~
Should the general simplifier, or the one‑argument limit function, handle` F(extended-real)`? We’ve discussed this—do we have a consensus?



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4636 signum(ind) is an error

[Barton W. <wil...@us...>]

I inserted the code  `($expand (div x ... ) 1 0)`. I didn't like it too much, but it I did it to avoid a check if the imaginary part was zero and make sure that the complex number case returns a number in rectangular form.  

I'll revise this bit of code.


---

**[bugs:#4636] signum(ind) is an error**

**Status:** open
**Group:** None
**Labels:** extended real signum limit 
**Created:** Thu Nov 27, 2025 12:41 PM UTC by Barton Willis
**Last Updated:** Fri Nov 28, 2025 01:14 PM UTC
**Owner:** nobody


Either a  `signum` nounform or `ind` is a better result than an error:
~~~
(%i4)	signum(ind);
sign: sign of ind is undefined.
~~~
Should the general simplifier, or the one‑argument limit function, handle` F(extended-real)`? We’ve discussed this—do we have a consensus?



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] Re: #4636 signum(ind) is an error

[Stavros M. <mac...@us...>]

Not sure when this code got in here:
   ($expand (div x (ftake 'mexpt (add (mul re re) (mul im im)) 1//2)) 1 0))
but that seems like a silly way to handle the real case (and makes the rest of the code unnecessary for real numerical inputs). Presumably this should not be used when im=0.


---

**[bugs:#4636] signum(ind) is an error**

**Status:** open
**Group:** None
**Labels:** extended real signum limit 
**Created:** Thu Nov 27, 2025 12:41 PM UTC by Barton Willis
**Last Updated:** Fri Nov 28, 2025 01:01 PM UTC
**Owner:** nobody


Either a  `signum` nounform or `ind` is a better result than an error:
~~~
(%i4)	signum(ind);
sign: sign of ind is undefined.
~~~
Should the general simplifier, or the one‑argument limit function, handle` F(extended-real)`? We’ve discussed this—do we have a consensus?



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4636 signum(ind) is an error

[Barton W. <wil...@us...>]

I simple fix is to errcatch sign.  I addition to a few minor edits, I used the fact that `complex-number-p` returns multiple values to make the complex number case a bit more efficient.

Let me know what you all think of this revision.

~~~
;; This is missing signum(x^n) = signsum(x)^n, where n is an integer--possibly should be restricted to a positive integer.
(def-simplifier signum (x)
  (cond
    ;; Case: zero--this makes signum(0.0) -> 0.0 and signum(0.0b0) -> 0.0b0
    ((zerop1 x) x)
    (t
     (multiple-value-bind (is-complex-number re im)
         (complex-number-p x #'mnump)
       (cond
         ;; Case: complex number
         (is-complex-number
          ($expand (div x (ftake 'mexpt (add (mul re re) (mul im im)) 1//2)) 1 0))
         ;; Idempotent: signum(signum(z)) = signum(z)
         ((and (consp x) (eq '%signum (caar x))) x)       
         (t
          ;; We use the errcatch to avoid signum(ind) throwing an error
          (let ((sgn (car (errcatch ($csign x)))))
            (cond
              ;; signum(negative) -> -1
              ((eq sgn '$neg) -1)
              ;; signum(zero) -> 0
              ((eq sgn '$zero) 0)
              ;; signum(positive) -> 1
              ((eq sgn '$pos) 1)
              ;; Multiplicative: signum(ab) = signum(a) * signum(b)
              ((mtimesp x)
               (fapply 'mtimes (mapcar #'(lambda (s) (ftake '%signum s)) (cdr x))))
              ;; Reflection rule: signum(-x) -> -signum(x)
              ((great (neg x) x)
               (neg (ftake '%signum (neg x))))
              ;; Give up
              (t (give-up))))))))))
~~~


---

**[bugs:#4636] signum(ind) is an error**

**Status:** open
**Group:** None
**Labels:** extended real signum limit 
**Created:** Thu Nov 27, 2025 12:41 PM UTC by Barton Willis
**Last Updated:** Thu Nov 27, 2025 12:41 PM UTC
**Owner:** nobody


Either a  `signum` nounform or `ind` is a better result than an error:
~~~
(%i4)	signum(ind);
sign: sign of ind is undefined.
~~~
Should the general simplifier, or the one‑argument limit function, handle` F(extended-real)`? We’ve discussed this—do we have a consensus?



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4637 error partition_set

[dan h. <zm...@us...>]

mistake it should read jave read partition_set(3,length(list)=2)  but i now see that there is a better option than this and should have used integer_partitions(3,2) instead. So cancel this bug report.


---

**[bugs:#4637] error partition_set**

**Status:** not-a-bug
**Group:** None
**Created:** Thu Nov 27, 2025 03:49 PM UTC by dan hayes
**Last Updated:** Thu Nov 27, 2025 04:17 PM UTC
**Owner:** nobody


partition_set(3,length(lis)=2)

instead of giving the correct answer of [[0,3], [ 1,2]) 
it gives error  like length argument cannot be a symbol.....
ok so it doesn't print out the 0' so i do
partition_set(3,length(lis)<3)
still gives same error


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4637 error partition_set

[Barton W. <wil...@us...>]

- **status**: open --> not-a-bug
- **Comment**:

The first argument to  must be a set, and the second must be a predicate (a boolean valued function).

For the example you gave, the first entry is an integer, not a set, and the second is not even a function. 

The user documentation for `partition_set`gives some examples that should clarify how to use this function.

Possibly, you want to use the function 'integer_partitions`. 

Closing this ticket as `not-a-bug`. 




---

**[bugs:#4637] error partition_set**

**Status:** not-a-bug
**Group:** None
**Created:** Thu Nov 27, 2025 03:49 PM UTC by dan hayes
**Last Updated:** Thu Nov 27, 2025 03:49 PM UTC
**Owner:** nobody


partition_set(3,length(lis)=2)

instead of giving the correct answer of [[0,3], [ 1,2]) 
it gives error  like length argument cannot be a symbol.....
ok so it doesn't print out the 0' so i do
partition_set(3,length(lis)<3)
still gives same error


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4637 error partition_set

[dan h. <zm...@us...>]

---

**[bugs:#4637] error partition_set**

**Status:** open
**Group:** None
**Created:** Thu Nov 27, 2025 03:49 PM UTC by dan hayes
**Last Updated:** Thu Nov 27, 2025 03:49 PM UTC
**Owner:** nobody


partition_set(3,length(lis)=2)

instead of giving the correct answer of [[0,3], [ 1,2]) 
it gives error  like length argument cannot be a symbol.....
ok so it doesn't print out the 0' so i do
partition_set(3,length(lis)<3)
still gives same error


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4636 signum(ind) is an error

[Barton W. <wil...@us...>]

---

**[bugs:#4636] signum(ind) is an error**

**Status:** open
**Group:** None
**Labels:** extended real signum limit 
**Created:** Thu Nov 27, 2025 12:41 PM UTC by Barton Willis
**Last Updated:** Thu Nov 27, 2025 12:41 PM UTC
**Owner:** nobody


Either a  `signum` nounform or `ind` is a better result than an error:
~~~
(%i4)	signum(ind);
sign: sign of ind is undefined.
~~~
Should the general simplifier, or the one‑argument limit function, handle` F(extended-real)`? We’ve discussed this—do we have a consensus?



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4629 box expression missing simp flags & expand applied to boxed objects

[Robert D. <rob...@us...>]

- **labels**:  --> box, expand, simplification



---

**[bugs:#4629] box expression missing simp flags & expand applied to boxed objects**

**Status:** open
**Group:** None
**Labels:** box expand simplification 
**Created:** Tue Nov 04, 2025 12:45 PM UTC by Barton Willis
**Last Updated:** Tue Nov 04, 2025 12:45 PM UTC
**Owner:** nobody


The user documentation doesn't say what `expand` is supposed to do to a boxed expression, but I'm pretty sure that we don't want `expand` to strip a boxed expression of some of its `simp` flags:

~~~
(%i25) xxx : box(x/2)$

(%i26) ?print([xxx,expand(xxx)]);

((MLIST SIMP) ((MBOX) ((MTIMES SIMP) ((RAT SIMP) 1 2) $X))
                                        ((MBOX) ((MTIMES) ((RAT) 1 2) $X)))
~~~
In the expanded expression, `rat` lost its  `simp` flag. Additionally,  the `mbox` operators are missing a `simp` flag.

For a labeled box, all is well:
~~~
(%i32) xxx : box(x/2,0)$

(%i33) ?print([xxx,expand(xxx)]);

((MLIST SIMP) ((MLABOX SIMP) ((MTIMES SIMP) ((RAT SIMP) 1 2) $X) 0)
                                        ((MLABOX SIMP) ((MTIMES SIMP) ((RAT SIMP) 1 2) $X) 0))
(%o33) [box(x/2,0),box(x/2,0)]
~~~

This was done with a recent Maxima build. With Maxima 5.48.1, I think all is well:

~~~
(%i1)	xxx : box(x/2);
(xxx)	box(x/2)
(%i2)	?print(%);
((MBOX SIMP) ((MTIMES SIMP) ((RAT SIMP) 1 2) $X)) 

(%o2)	box(x/2)
(%i3)	expand(xxx);
(%o3)	box(x/2)
(%i4)	?print(%);
((MBOX SIMP) ((MTIMES SIMP) ((RAT SIMP) 1 2) $X)) 

(%o4)	box(x/2)
~~~




---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4634 Divergent integral simplifies to -1/4 with simplify_sum

[Robert D. <rob...@us...>]

- **labels**:  --> simplify_sum, limit, sum



---

**[bugs:#4634] Divergent integral simplifies to -1/4 with simplify_sum**

**Status:** open
**Group:** None
**Labels:** simplify_sum limit sum 
**Created:** Tue Nov 18, 2025 12:58 PM UTC by kcrisman
**Last Updated:** Tue Nov 18, 2025 01:45 PM UTC
**Owner:** nobody


In the [Sage support list](https://groups.google.com/g/sage-support/c/pCYoD6x2Y8M) we have the following message about a divergent integral giving a solution with the optional package `simplify_sum`.

"In maxima 5.47.0 I get:

```
(%i1) load(simplify_sum)$
(%i2) simplify_sum( sum(n*(-1)^n,n,1,inf));
-1/4
```

This happens both on sage's ECL based one and on the system maxima (running on SBCL). "

That this is a reversion is seen from the following:

```
Maxima 5.42.2 http://maxima.sourceforge.net
using Lisp ECL 16.1.2
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) load(simplify_sum)$

(%i2) simplify_sum( sum(n*(-1)^n,n,1,inf));
(%o2)                                 ind
```



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4628 %solve(abs(sin(x))=1, x) misses solutions

[David S. <tom...@us...>]

Interesting! One would think that solving `abs(f(x)) = c` is done by solving `f(x) = c` and `-f(x) = c` and then returning the union of the solutions. Apparently, this is not what's happening, because solving `sin(x) = 1` and `sin(x) = -1` yields only the real solutions.


---

**[bugs:#4628] %solve(abs(sin(x))=1,x) misses solutions**

**Status:** open
**Group:** None
**Labels:** %solve 
**Created:** Sat Nov 01, 2025 08:50 PM UTC by Stavros Macrakis
**Last Updated:** Sat Nov 01, 2025 08:50 PM UTC
**Owner:** Barton Willis


<pre>
%solve(abs(sin(x))=1,x) =>
   %union([x = 2*%pi*%z3405-%i*log(sqrt(2)-1)],
                     [x = 2*%pi*%z3454-%i*log(sqrt(2)+1)])
</pre>
This is missing all the real solutions (namely %pi*(%z+1/2)) as well as many other complex solutions.


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4635 one argument limit of linear comb of inf and minf

[Barton W. <wil...@us...>]

---

**[bugs:#4635] one argument limit of linear comb of inf and minf**

**Status:** open
**Group:** None
**Created:** Thu Nov 20, 2025 07:21 PM UTC by Barton Willis
**Last Updated:** Thu Nov 20, 2025 07:21 PM UTC
**Owner:** nobody


This is wrong--Maxima needs to ask about the sign of `a-b`, not the sign of `a` and `b` separately:
~~~
(%i1) limit(a*inf + b*minf);
Is a positive, negative or zero?
p;
Is b positive, negative or zero?
n;
(%o1)                                 inf
~~~



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4633 integrate(exp(- t) log(t), t, 0, 1) --> integral nounform

[Barton W. <wil...@us...>]

- **status**: open --> closed
- **Comment**:

Fixed by Commit [73bf03] master.  Closing this ticket. 



---

**[bugs:#4633] integrate(exp(- t) log(t), t, 0, 1) --> integral nounform**

**Status:** closed
**Group:** None
**Labels:** limit gamma_incomplete defint 
**Created:** Fri Nov 14, 2025 01:06 AM UTC by Barton Willis
**Last Updated:** Fri Nov 14, 2025 01:06 AM UTC
**Owner:** nobody


We get a noun form because Maxima can't find a limit of an incomplete gamma function.  Repair this limit limitation, and I think this definite integral will give a correct value:
~~~
(%i3) integrate(exp(-t) * log(t), t, 0, 1);

(%o3) 'integrate(%e^-t*log(t),t,0,1)
~~~
This is a test in `rtest_log` (problem 120 , line 433) 


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #3041 limit(inf*(zeroa+inf)) => und, should be inf

[Barton W. <wil...@us...>]

- **status**: open --> closed
- **Comment**:

Fixed by Commit [31a229] master.  The test suggested in this ticket was included in `rtest_limit_extra` some time ago.  Closing this ticket.



---

**[bugs:#3041] limit(inf*(zeroa+inf)) => und, should be inf**

**Status:** closed
**Group:** None
**Labels:** limit 
**Created:** Sun Oct 25, 2015 07:30 PM UTC by Stavros Macrakis
**Last Updated:** Tue Nov 11, 2025 01:04 PM UTC
**Owner:** nobody


limit(inf(zeroa+inf)) => und, should be inf

single-arg limit should be true regardless of how the various parts reach the limits, i.e., in this case it is equivalent to the *multiple* limit:

      limit (a->inf, b->0+, c->inf) a*(b+c)
      
it looks like instead it is expanding to 

     inf * zeroa + infinf => und + inf => und


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4634 Divergent integral simplifies to -1/4 with simplify_sum

[Barton W. <wil...@us...>]

Specifically, this is a bug in `simplimplus`:

~~~
(%i66) limit(((-1)^N*(2*N+1))/4-1/4,N,inf);

1 Call   simplimplus [((-1)^N*N)/2+(-1)^N/4-1/4]
1 Return simplimplus ind
(%o66) ind
~~~


---

**[bugs:#4634] Divergent integral simplifies to -1/4 with simplify_sum**

**Status:** open
**Group:** None
**Created:** Tue Nov 18, 2025 12:58 PM UTC by kcrisman
**Last Updated:** Tue Nov 18, 2025 01:18 PM UTC
**Owner:** nobody


In the [Sage support list](https://groups.google.com/g/sage-support/c/pCYoD6x2Y8M) we have the following message about a divergent integral giving a solution with the optional package `simplify_sum`.

"In maxima 5.47.0 I get:

```
(%i1) load(simplify_sum)$
(%i2) simplify_sum( sum(n*(-1)^n,n,1,inf));
-1/4
```

This happens both on sage's ECL based one and on the system maxima (running on SBCL). "

That this is a reversion is seen from the following:

```
Maxima 5.42.2 http://maxima.sourceforge.net
using Lisp ECL 16.1.2
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) load(simplify_sum)$

(%i2) simplify_sum( sum(n*(-1)^n,n,1,inf));
(%o2)                                 ind
```



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4634 Divergent integral simplifies to -1/4 with simplify_sum

[Barton W. <wil...@us...>]

With Maxima 5.48.1, we're back to the* incorrect* value of `ind`

~~~
(%i11)	simplify_sum(sum(n*(-1)^n,n,1,inf));
(%o11)	ind
~~~

This result is wrong and it is due to a bug in the limit 

~~~
(%i12)	simplify_sum(sum(n*(-1)^n,n,1,N));
(%o12)	((-1)^N*(2*N+1))/4-1/4

(%i13)	limit(%,N,inf);
(%o13)	ind
~~~
The output `%o12` is correct, but the  output `%o13` should be `und`, not `ind`.

Thanks for this reporting this bug.




---

**[bugs:#4634] Divergent integral simplifies to -1/4 with simplify_sum**

**Status:** open
**Group:** None
**Created:** Tue Nov 18, 2025 12:58 PM UTC by kcrisman
**Last Updated:** Tue Nov 18, 2025 12:58 PM UTC
**Owner:** nobody


In the [Sage support list](https://groups.google.com/g/sage-support/c/pCYoD6x2Y8M) we have the following message about a divergent integral giving a solution with the optional package `simplify_sum`.

"In maxima 5.47.0 I get:

```
(%i1) load(simplify_sum)$
(%i2) simplify_sum( sum(n*(-1)^n,n,1,inf));
-1/4
```

This happens both on sage's ECL based one and on the system maxima (running on SBCL). "

That this is a reversion is seen from the following:

```
Maxima 5.42.2 http://maxima.sourceforge.net
using Lisp ECL 16.1.2
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) load(simplify_sum)$

(%i2) simplify_sum( sum(n*(-1)^n,n,1,inf));
(%o2)                                 ind
```



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4634 Divergent integral simplifies to -1/4 with simplify_sum

[kcrisman <kcr...@us...>]

---

**[bugs:#4634] Divergent integral simplifies to -1/4 with simplify_sum**

**Status:** open
**Group:** None
**Created:** Tue Nov 18, 2025 12:58 PM UTC by kcrisman
**Last Updated:** Tue Nov 18, 2025 12:58 PM UTC
**Owner:** nobody


In the [Sage support list](https://groups.google.com/g/sage-support/c/pCYoD6x2Y8M) we have the following message about a divergent integral giving a solution with the optional package `simplify_sum`.

"In maxima 5.47.0 I get:

```
(%i1) load(simplify_sum)$
(%i2) simplify_sum( sum(n*(-1)^n,n,1,inf));
-1/4
```

This happens both on sage's ECL based one and on the system maxima (running on SBCL). "

That this is a reversion is seen from the following:

```
Maxima 5.42.2 http://maxima.sourceforge.net
using Lisp ECL 16.1.2
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) load(simplify_sum)$

(%i2) simplify_sum( sum(n*(-1)^n,n,1,inf));
(%o2)                                 ind
```



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #3041 limit(inf*(zeroa+inf)) => und, should be inf

[Barton W. <wil...@us...>]

The function  `toplevel-$limit` calls `($expand exp 1 0)`  before sending the limit to the main limit code.  Changing the the call to `($expand exp 0 0)` (for pure simplification) results in several testsuite failures, some of them serious.  Previously, I've attempted, but  failed to fix these bugs. 

As a workaround,  I'm experimenting with calling `($expand exp 0 0)` for expressions that are free of the limit variable and  otherwise call `($expand exp 1 0)`.

There is a source code comment about all this




---

**[bugs:#3041] limit(inf*(zeroa+inf)) => und, should be inf**

**Status:** open
**Group:** None
**Labels:** limit 
**Created:** Sun Oct 25, 2015 07:30 PM UTC by Stavros Macrakis
**Last Updated:** Sun Oct 25, 2015 07:31 PM UTC
**Owner:** nobody


limit(inf(zeroa+inf)) => und, should be inf

single-arg limit should be true regardless of how the various parts reach the limits, i.e., in this case it is equivalent to the *multiple* limit:

      limit (a->inf, b->0+, c->inf) a*(b+c)
      
it looks like instead it is expanding to 

     inf * zeroa + infinf => und + inf => und


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4632 gentran file management functions mishandle streams

[Robert D. <rob...@us...>]

- **status**: open --> closed
- **Comment**:

Fixed by commit [ d9dcda4 ].



---

**[bugs:#4632] gentran file management functions mishandle streams**

**Status:** closed
**Group:** None
**Labels:** gentran share packages 
**Created:** Sat Nov 08, 2025 06:40 AM UTC by Robert Dodier
**Last Updated:** Sat Nov 08, 2025 06:40 AM UTC
**Owner:** nobody


`gentran` has functions to handle a list of output files; `gentranin`, `gentranout`, etc. Looks like those functions assume  their arguments are always file names, and not streams. `rtest_gentran.mac` has a helper function which tries apply one or the other of those functions to a string output stream, so `batch("rtest_gentran.mac")` provokes a Lisp error.


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4632 gentran file management functions mishandle streams

[Robert D. <rob...@us...>]

---

**[bugs:#4632] gentran file management functions mishandle streams**

**Status:** open
**Group:** None
**Labels:** gentran share packages 
**Created:** Sat Nov 08, 2025 06:40 AM UTC by Robert Dodier
**Last Updated:** Sat Nov 08, 2025 06:40 AM UTC
**Owner:** nobody


`gentran` has functions to handle a list of output files; `gentranin`, `gentranout`, etc. Looks like those functions assume  their arguments are always file names, and not streams. `rtest_gentran.mac` has a helper function which tries apply one or the other of those functions to a string output stream, so `batch("rtest_gentran.mac")` provokes a Lisp error.


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4631 Lisp error in gentran with gentranlang = c (attempting symbol operation on non-symbol)

[Robert D. <rob...@us...>]

- **status**: open --> closed
- **Comment**:

Fixed by commit [ 6b7c453 ].



---

**[bugs:#4631] Lisp error in gentran with gentranlang = c (attempting symbol operation on non-symbol)**

**Status:** closed
**Group:** None
**Labels:** share packages gentran 
**Created:** Thu Nov 06, 2025 06:04 AM UTC by Robert Dodier
**Last Updated:** Thu Nov 06, 2025 06:05 AM UTC
**Owner:** nobody


Followup to bug #4625. Trying the example originally shown there provokes a Lisp error -- SYMBOL-PLIST is being called on a non-symbol.  The bit about "t0 undeclared variable" is just an advisory; it's not an error.
```
(%i1) load (gentran) $
(%i2) gentranlang:c $
(%i3) var:-((sqrt(m_piM^4+W^8*cost^2-W^8+4*P0^2*W^6)+((cost^2-1)*m_nuc^2+W^2*cost^2-W^2)*m_piM^4+sqrt(m_piM^2+P0^2)*(4*W^2*m_nuc*m_piM^2-4*W^2*m_nuc^3+4*W^4*m_nuc)+((2-2*cost^2)*m_nuc^4-4*W^2*cost^2*m_nuc^2-2*W^4*cost^2-2*W^4)*m_piM^2+(cost^2-1)*m_nuc^6+(W^2-W^2*cost^2)*m_nuc^4+(W^4-W^4*cost^2)*m_nuc^2+W^6*cost^2-W^6)/((cost^2-1)*m_piM^4+((2-2*cost^2)*m_nuc^2-2*W^2*cost^2-2*W^2)*m_piM^2+(cost^2-1)*m_nuc^4+(2*W^2-2*W^2*cost^2)*m_nuc^2+W^4*cost^2-W^4)) $

(%i4) gentran (eval (var));
gentran: t0 undeclared variable

debugger invoked on a TYPE-ERROR @B800C26BD7 in thread
#<THREAD tid=32336 "main thread" RUNNING {1001760003}>:
  The value
    (SETQ T0 (MINUS (PLUS # # # # # # # # #)))
  is not of type
    SYMBOL
  when binding SYMBOL

Type HELP for debugger help, or (SB-EXT:EXIT) to exit from SBCL.

restarts (invokable by number or by possibly-abbreviated name):
  0: [MACSYMA-QUIT] Maxima top-level
  1: [CONTINUE    ] Ignore runtime option --eval "(cl-user::run)".
  2: [ABORT       ] Skip rest of --eval and --load options.
  3:                Skip to toplevel READ/EVAL/PRINT loop.
  4: [EXIT        ] Exit SBCL (calling #'EXIT, killing the process).

(SYMBOL-PLIST (SETQ T0 (MINUS (PLUS (EXPT (PLUS # # # #) (QUOTIENT 1.0d0 2.0d0)) (TIMES (PLUS # # #) (EXPT |m_piM| 4.0d0)) (TIMES (EXPT # #) (PLUS # # #)) (TIMES (PLUS # # # #) (EXPT |m_piM| 2.0d0)) (TIMES (PLUS # -1) (EXPT M_NUC 6.0d0)) (TIMES (PLUS # #) (EXPT M_NUC 4.0d0)) (TIMES (PLUS # #) (EXPT M_NUC 2.0d0)) (TIMES (EXPT |w| 6.0d0) (EXPT COST 2.0d0)) (MINUS (EXPT |w| 6.0d0)))))) [external]
0] 
```
I think the problem is that gentran has created a new assignment statement (to assign a partial result to the made-up variable `t0`) and appended it to the list of expressions to be returned, and there is one level of nested too many. E.g. putting `(list foo)` onto the return list instead of just `foo`, something like that. Not 100% sure.


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4603 Control stack regression with abs_integrate / 5.48.0

[Michael O. <mjo...@us...>]

Thank you both!


---

**[bugs:#4603] Control stack regression with abs_integrate / 5.48.0**

**Status:** closed
**Group:** None
**Labels:** abs_integrate defint limit 
**Created:** Thu Aug 28, 2025 01:46 PM UTC by Michael Orlitzky
**Last Updated:** Thu Nov 06, 2025 05:25 PM UTC
**Owner:** nobody


https://sourceforge.net/p/maxima/bugs/3144/

This works in 5.47.0:

```
Maxima 5.47.0 https://maxima.sourceforge.io
using Lisp ECL 24.5.10
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) integrate(sin(k*x)/x*erf(x^2),x,0,inf);
                            inf
                           /                  2
                           [    sin(k x) erf(x )
(%o1)                      I    ---------------- dx
                           ]           x
                           /
                            0
(%i2)  load(abs_integrate);  
(%o2) /usr/share/maxima/5.47.0/share/contrib/integration/abs_integrate.mac
(%i3) integrate(sin(k*x)/x*erf(x^2),x,0,inf);
                            inf
                           /                  2
                           [    sin(k x) erf(x )
(%o3)                      I    ---------------- dx
                           ]           x
                           /
                            0
```

But has regressed in 5.48.1:

```
Maxima 5.48.1 https://maxima.sourceforge.io
using Lisp ECL 24.5.10
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) integrate(sin(k*x)/x*erf(x^2),x,0,inf);
                            inf
                           ⌠                  2
                           ⎮    sin(k x) erf(x )
(%o1)                      ⎮    ──────────────── dx
                           ⎮           x
                           ⌡
                            0
(%i2) load(abs_integrate);
(%o2) /usr/share/maxima/5.48.1/share/contrib/integration/abs_integrate.mac
(%i3) integrate(sin(k*x)/x*erf(x^2),x,0,inf);

Maxima encountered a Lisp error:

 BINDING-STACK overflow at size 10240. Stack can probably be resized.
Proceed with caution.

Automatically continuing.
To enable the Lisp debugger set *debugger-hook* to nil.
```


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4603 Control stack regression with abs_integrate / 5.48.0

[Barton W. <wil...@us...>]

- **status**: open --> closed
- **Comment**:

Fixed by Commit [0c11d9] master.  The result is an integrate nounform, but it nolonger gives a stack overflow. 

Append test to `rtest_abs_integrate`.  Closing this ticket.  Thanks to Robert who identified the reason for this bug.



---

**[bugs:#4603] Control stack regression with abs_integrate / 5.48.0**

**Status:** closed
**Group:** None
**Labels:** abs_integrate defint limit 
**Created:** Thu Aug 28, 2025 01:46 PM UTC by Michael Orlitzky
**Last Updated:** Thu Oct 09, 2025 06:23 PM UTC
**Owner:** nobody


https://sourceforge.net/p/maxima/bugs/3144/

This works in 5.47.0:

```
Maxima 5.47.0 https://maxima.sourceforge.io
using Lisp ECL 24.5.10
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) integrate(sin(k*x)/x*erf(x^2),x,0,inf);
                            inf
                           /                  2
                           [    sin(k x) erf(x )
(%o1)                      I    ---------------- dx
                           ]           x
                           /
                            0
(%i2)  load(abs_integrate);  
(%o2) /usr/share/maxima/5.47.0/share/contrib/integration/abs_integrate.mac
(%i3) integrate(sin(k*x)/x*erf(x^2),x,0,inf);
                            inf
                           /                  2
                           [    sin(k x) erf(x )
(%o3)                      I    ---------------- dx
                           ]           x
                           /
                            0
```

But has regressed in 5.48.1:

```
Maxima 5.48.1 https://maxima.sourceforge.io
using Lisp ECL 24.5.10
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) integrate(sin(k*x)/x*erf(x^2),x,0,inf);
                            inf
                           ⌠                  2
                           ⎮    sin(k x) erf(x )
(%o1)                      ⎮    ──────────────── dx
                           ⎮           x
                           ⌡
                            0
(%i2) load(abs_integrate);
(%o2) /usr/share/maxima/5.48.1/share/contrib/integration/abs_integrate.mac
(%i3) integrate(sin(k*x)/x*erf(x^2),x,0,inf);

Maxima encountered a Lisp error:

 BINDING-STACK overflow at size 10240. Stack can probably be resized.
Proceed with caution.

Automatically continuing.
To enable the Lisp debugger set *debugger-hook* to nil.
```


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4619 limit(inf = inf) causes stack overflow

[Barton W. <wil...@us...>]

- **status**: open --> closed
- **Comment**:

Fixed by Commit [0c11d9] master. Closing this ticket.



---

**[bugs:#4619] limit(inf = inf) causes stack overflow**

**Status:** closed
**Group:** None
**Labels:** limit 
**Created:** Thu Oct 09, 2025 06:22 PM UTC by Robert Dodier
**Last Updated:** Thu Nov 06, 2025 03:52 PM UTC
**Owner:** nobody


Looking at bug report #4603, looks like it is triggered by a bug in `limit`. Here is the simplest case I found:
```
(%i1) limit(inf=inf);
INFO: Binding stack guard page unprotected
Binding stack guard page temporarily disabled: proceed with caution

debugger invoked on a SB-KERNEL::BINDING-STACK-EXHAUSTED in thread
#<THREAD tid=4380 "main thread" RUNNING {1001760003}>:
  Binding stack exhausted.
```
This other example also causes a stack overflow: `limit(minf = minf)`.

Bug #4603 is triggered by `limit(foo - inf = -inf)`, where `foo` is a gensym (and the bug is also triggered when `foo` is just `foo`).

Looks like the stack overflow originates in INFSIMP2. I didn't investigate beyond that.

When this bug gets fixed, we should go back to #4603 and verify that bug is also fixed.


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4619 limit(inf = inf) causes stack overflow

[Barton W. <wil...@us...>]

The function `toplevel-$limit` does a check  (using the predicate limunknown) for a limit expression that should return a limit nounform.  I think this check should be moved to down to the function `limit`.  That's because the limit code sometimes calls `limit` without going through `toplevel-$limit`--when that happens limits such as `limit(inf=inf)` confuse the code.

Additionally, I think I'll alter `limunknow` to pass the limit variable instead of using the special. 

Just this change, not a modification of `infsimp` is all that is needed for a fix to this bug. It also fixes bug #4603.

The suggestion to change `(or (not (free...)))` with a call to among is good--maybe I'll do that too, but not as a part of this bug fix.




---

**[bugs:#4619] limit(inf = inf) causes stack overflow**

**Status:** open
**Group:** None
**Labels:** limit 
**Created:** Thu Oct 09, 2025 06:22 PM UTC by Robert Dodier
**Last Updated:** Sun Oct 12, 2025 02:30 PM UTC
**Owner:** nobody


Looking at bug report #4603, looks like it is triggered by a bug in `limit`. Here is the simplest case I found:
```
(%i1) limit(inf=inf);
INFO: Binding stack guard page unprotected
Binding stack guard page temporarily disabled: proceed with caution

debugger invoked on a SB-KERNEL::BINDING-STACK-EXHAUSTED in thread
#<THREAD tid=4380 "main thread" RUNNING {1001760003}>:
  Binding stack exhausted.
```
This other example also causes a stack overflow: `limit(minf = minf)`.

Bug #4603 is triggered by `limit(foo - inf = -inf)`, where `foo` is a gensym (and the bug is also triggered when `foo` is just `foo`).

Looks like the stack overflow originates in INFSIMP2. I didn't investigate beyond that.

When this bug gets fixed, we should go back to #4603 and verify that bug is also fixed.


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.