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

[Maxima-bugs] [maxima:bugs] #4647 Maxima error in --preload file causes Lisp error

[Raymond T. <rt...@us...>]

- **status**: open --> closed
- **assigned_to**: Raymond Toy
- **Comment**:

Fixed in [ce3f0a]



---

**[bugs:#4647] Maxima error in --preload file causes Lisp error**

**Status:** closed
**Group:** None
**Labels:** command line batch error handling 
**Created:** Wed Dec 10, 2025 05:48 PM UTC by Robert Dodier
**Last Updated:** Wed Dec 10, 2025 10:57 PM UTC
**Owner:** Raymond Toy


Working with a recent build (1c6a095). A Maxima error in a file loaded via `--preload` causes a Lisp error, and puts the user in the Lisp debugger. 

E.g. let `foo.mac` contain just `error("Foo");`. Then:

```
sh maxima-local -l sbcl -p foo.mac
read and interpret /home/dodier/sourceforge/maxima-code/foo.mac
Foo
 -- an error. To debug this try: debugmode(true);
Maxima encountered a Lisp error:

 attempt to THROW to a tag that does not exist: MACSYMA-QUIT

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

debugger invoked on a SB-INT:SIMPLE-CONTROL-ERROR @52A00940 in thread
#<THREAD "main thread" RUNNING {10028F4103}>:
  attempt to THROW to a tag that does not exist: RETURN-FROM-DEBUGGER

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

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

(THROW)
0]
```
I think whatever mechanism is loading `--preload` files needs to catch MACSYMA-QUIT and just go on to the next command line argument.


---

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: #4647 Maxima error in --preload file causes Lisp error

[Raymond T. <rt...@us...>]

In a bit.  But this fix also does something reasonable for lisp files.  With your foo.mac, we now get:
```
$ ./maxima-local --preload foo.mac
read and interpret /home/toy/src/sourceforge/maxima/foo.mac
Foo
 -- an error. To debug this try: debugmode(true);
```
With foo.lisp containing just `(error "Foo")`, we get:
```
$ ./maxima-local --preload foo.lisp
loadfile: loading foo.lisp.
loadfile: failed to load foo.lisp
 -- an error. To debug this try: debugmode(true);
```


---

**[bugs:#4647] Maxima error in --preload file causes Lisp error**

**Status:** open
**Group:** None
**Labels:** command line batch error handling 
**Created:** Wed Dec 10, 2025 05:48 PM UTC by Robert Dodier
**Last Updated:** Wed Dec 10, 2025 10:05 PM UTC
**Owner:** nobody


Working with a recent build (1c6a095). A Maxima error in a file loaded via `--preload` causes a Lisp error, and puts the user in the Lisp debugger. 

E.g. let `foo.mac` contain just `error("Foo");`. Then:

```
sh maxima-local -l sbcl -p foo.mac
read and interpret /home/dodier/sourceforge/maxima-code/foo.mac
Foo
 -- an error. To debug this try: debugmode(true);
Maxima encountered a Lisp error:

 attempt to THROW to a tag that does not exist: MACSYMA-QUIT

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

debugger invoked on a SB-INT:SIMPLE-CONTROL-ERROR @52A00940 in thread
#<THREAD "main thread" RUNNING {10028F4103}>:
  attempt to THROW to a tag that does not exist: RETURN-FROM-DEBUGGER

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

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

(THROW)
0]
```
I think whatever mechanism is loading `--preload` files needs to catch MACSYMA-QUIT and just go on to the next command line argument.


---

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: #4647 Maxima error in --preload file causes Lisp error

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

Looks good, can you commit that?


---

**[bugs:#4647] Maxima error in --preload file causes Lisp error**

**Status:** open
**Group:** None
**Labels:** command line batch error handling 
**Created:** Wed Dec 10, 2025 05:48 PM UTC by Robert Dodier
**Last Updated:** Wed Dec 10, 2025 07:06 PM UTC
**Owner:** nobody


Working with a recent build (1c6a095). A Maxima error in a file loaded via `--preload` causes a Lisp error, and puts the user in the Lisp debugger. 

E.g. let `foo.mac` contain just `error("Foo");`. Then:

```
sh maxima-local -l sbcl -p foo.mac
read and interpret /home/dodier/sourceforge/maxima-code/foo.mac
Foo
 -- an error. To debug this try: debugmode(true);
Maxima encountered a Lisp error:

 attempt to THROW to a tag that does not exist: MACSYMA-QUIT

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

debugger invoked on a SB-INT:SIMPLE-CONTROL-ERROR @52A00940 in thread
#<THREAD "main thread" RUNNING {10028F4103}>:
  attempt to THROW to a tag that does not exist: RETURN-FROM-DEBUGGER

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

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

(THROW)
0]
```
I think whatever mechanism is loading `--preload` files needs to catch MACSYMA-QUIT and just go on to the next command line argument.


---

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] #4648 autoload problem with hstep

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

---

**[bugs:#4648] autoload problem with hstep**

**Status:** open
**Group:** None
**Labels:** autoload hstep 
**Created:** Wed Dec 10, 2025 07:50 PM UTC by Barton Willis
**Last Updated:** Wed Dec 10, 2025 07:50 PM UTC
**Owner:** nobody


Wrong: should be 1, not `hstep(5)`
~~~
(%i1)	'hstep(5);
(%o1)	hstep(5)
~~~
But after entering `hstep(5)`, the `hstep` package  autoload  and after that `'hstep(5)` simplifies correctly
~~~
(%i2)	hstep(5);
(%o2)	1

(%i3)	'hstep(5);
(%o3)	1
~~~


---

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] #4647 Maxima error in --preload file causes Lisp error

[Raymond T. <rt...@us...>]

I modified the action for `--preload` to
```
(let (($loadprint (not *maxima-quiet*)))
  (catch 'macsyma-quit
    ($load file)))
```
This works.


---

**[bugs:#4647] Maxima error in --preload file causes Lisp error**

**Status:** open
**Group:** None
**Labels:** command line batch error handling 
**Created:** Wed Dec 10, 2025 05:48 PM UTC by Robert Dodier
**Last Updated:** Wed Dec 10, 2025 05:48 PM UTC
**Owner:** nobody


Working with a recent build (1c6a095). A Maxima error in a file loaded via `--preload` causes a Lisp error, and puts the user in the Lisp debugger. 

E.g. let `foo.mac` contain just `error("Foo");`. Then:

```
sh maxima-local -l sbcl -p foo.mac
read and interpret /home/dodier/sourceforge/maxima-code/foo.mac
Foo
 -- an error. To debug this try: debugmode(true);
Maxima encountered a Lisp error:

 attempt to THROW to a tag that does not exist: MACSYMA-QUIT

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

debugger invoked on a SB-INT:SIMPLE-CONTROL-ERROR @52A00940 in thread
#<THREAD "main thread" RUNNING {10028F4103}>:
  attempt to THROW to a tag that does not exist: RETURN-FROM-DEBUGGER

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

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

(THROW)
0]
```
I think whatever mechanism is loading `--preload` files needs to catch MACSYMA-QUIT and just go on to the next command line argument.


---

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] #4647 Maxima error in --preload file causes Lisp error

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

---

**[bugs:#4647] Maxima error in --preload file causes Lisp error**

**Status:** open
**Group:** None
**Labels:** command line batch error handling 
**Created:** Wed Dec 10, 2025 05:48 PM UTC by Robert Dodier
**Last Updated:** Wed Dec 10, 2025 05:48 PM UTC
**Owner:** nobody


Working with a recent build (1c6a095). A Maxima error in a file loaded via `--preload` causes a Lisp error, and puts the user in the Lisp debugger. 

E.g. let `foo.mac` contain just `error("Foo");`. Then:

```
sh maxima-local -l sbcl -p foo.mac
read and interpret /home/dodier/sourceforge/maxima-code/foo.mac
Foo
 -- an error. To debug this try: debugmode(true);
Maxima encountered a Lisp error:

 attempt to THROW to a tag that does not exist: MACSYMA-QUIT

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

debugger invoked on a SB-INT:SIMPLE-CONTROL-ERROR @52A00940 in thread
#<THREAD "main thread" RUNNING {10028F4103}>:
  attempt to THROW to a tag that does not exist: RETURN-FROM-DEBUGGER

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

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

(THROW)
0]
```
I think whatever mechanism is loading `--preload` files needs to catch MACSYMA-QUIT and just go on to the next command line argument.


---

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: #4644 bfloat((x+1)/x includes unnecessary 1.0b0 factor

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

This error occurs after modifying `timesbigfloat` to not introduce the `1.0b0` factor.


---

**[bugs:#4644] bfloat((x+1)/x includes unnecessary 1.0b0 factor**

**Status:** open
**Group:** None
**Labels:** bfloat nfloat 
**Created:** Wed Dec 03, 2025 07:15 PM UTC by Stavros Macrakis
**Last Updated:** Mon Dec 08, 2025 05:58 PM UTC
**Owner:** nobody


~~~
bfloat((x+2)/x) => (1.0b0*(x+2.0b0))/x
~~~
What is the ``1.0b0`` factor doing there?

Set as priority 1 (which I hope means lowest, not highest) because this probably isn't doing any harm.



---

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] #4644 bfloat((x+1)/x includes unnecessary 1.0b0 factor

[Raymond T. <rt...@us...>]

I always thought priority 1 was the highest, but I don't actually know.


---

**[bugs:#4644] bfloat((x+1)/x includes unnecessary 1.0b0 factor**

**Status:** open
**Group:** None
**Labels:** bfloat nfloat 
**Created:** Wed Dec 03, 2025 07:15 PM UTC by Stavros Macrakis
**Last Updated:** Mon Dec 08, 2025 05:20 PM UTC
**Owner:** nobody


~~~
bfloat((x+2)/x) => (1.0b0*(x+2.0b0))/x
~~~
What is the ``1.0b0`` factor doing there?

Set as priority 1 (which I hope means lowest, not highest) because this probably isn't doing any harm.



---

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: #4644 bfloat((x+1)/x includes unnecessary 1.0b0 factor

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

I can't seem to reproduce that error in Maxima 5.48.1 SBCL 2.5.7.

~~~
map('nfloat,
    [1,1.0,1.0b0,2/3,x,1.0*x,1.0b0*x,2/3*x,
     x/y,2*x/y,2/3*x/y,1.0*x/y,1.0b0*x/y,
     %i*x,%i*1.0*x,%i*1.0b0*x,%i*2/3*x,
     %i*x/y,%i*2*x/y,%i*2/3*x/y,%i*1.0*x/y,%i*1.0b0*x/y,
     1+%i*x,1+%i*1.0*x,1+%i*1.0b0*x,1+%i*2/3*x,1+%i*x/y,
     1+%i*2*x/y,1+%i*2/3*x/y,1+%i*1.0*x/y,1+%i*1.0b0*x/y])
       => no errors
~~~
Maybe a bug has been introduced into `nfloat` recently?

Of course, it's hard to know what's a bug and what isn't because `nfloat` isn't documented!



---

**[bugs:#4644] bfloat((x+1)/x includes unnecessary 1.0b0 factor**

**Status:** open
**Group:** None
**Labels:** bfloat nfloat 
**Created:** Wed Dec 03, 2025 07:15 PM UTC by Stavros Macrakis
**Last Updated:** Mon Dec 08, 2025 05:17 PM UTC
**Owner:** nobody


~~~
bfloat((x+2)/x) => (1.0b0*(x+2.0b0))/x
~~~
What is the ``1.0b0`` factor doing there?

Set as priority 1 (which I hope means lowest, not highest) because this probably isn't doing any harm.



---

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] #4644 bfloat((x+1)/x includes unnecessary 1.0b0 factor

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

- **labels**: bfloat --> bfloat, nfloat



---

**[bugs:#4644] bfloat((x+1)/x includes unnecessary 1.0b0 factor**

**Status:** open
**Group:** None
**Labels:** bfloat nfloat 
**Created:** Wed Dec 03, 2025 07:15 PM UTC by Stavros Macrakis
**Last Updated:** Sun Dec 07, 2025 05:21 PM UTC
**Owner:** nobody


~~~
bfloat((x+2)/x) => (1.0b0*(x+2.0b0))/x
~~~
What is the ``1.0b0`` factor doing there?

Set as priority 1 (which I hope means lowest, not highest) because this probably isn't doing any harm.



---

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] #4646 minor hstep problems and missing features

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

---

**[bugs:#4646] minor hstep problems and missing features**

**Status:** open
**Group:** None
**Created:** Mon Dec 08, 2025 02:50 PM UTC by Barton Willis
**Last Updated:** Mon Dec 08, 2025 02:50 PM UTC
**Owner:** Barton Willis


This is OK
~~~
(%i1) hstep(x^2+1);
(%o1)                                  1
~~~
But this is not OK--it should be a `hstep` nounform:
~~~
(%i2) declare(x,complex)$
(%i3) hstep(x^2+1);
(%o3)                                  1
~~~

Not wrong, but `pz` is better:
~~~
(%i1) sign(hstep(x));
(%o1)                                 pnz
~~~

And we should give `hstep` an antiderivative property:
~~~
(%i3) integrate(hstep(x),x);

(%o3) 'integrate(hstep(x),x)
~~~

I have a putative fix for these deficiencies. 



---

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] #4644 bfloat((x+1)/x includes unnecessary 1.0b0 factor

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

Removing the `1.0b0` factor breaks the `nfloat` function:

```
nfloat(1.0b0*%i);
BIGFLOAT: unable to convert %i 1 to a CL or BIGFLOAT number.
 -- an error. To debug this try: debugmode(true);
```


---

**[bugs:#4644] bfloat((x+1)/x includes unnecessary 1.0b0 factor**

**Status:** open
**Group:** None
**Labels:** bfloat 
**Created:** Wed Dec 03, 2025 07:15 PM UTC by Stavros Macrakis
**Last Updated:** Sun Dec 07, 2025 04:39 PM UTC
**Owner:** nobody


~~~
bfloat((x+2)/x) => (1.0b0*(x+2.0b0))/x
~~~
What is the ``1.0b0`` factor doing there?

Set as priority 1 (which I hope means lowest, not highest) because this probably isn't doing any harm.



---

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] #4644 bfloat((x+1)/x includes unnecessary 1.0b0 factor

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

Actually, it's any product. `bfloat(x*y)` also introduces the `1.0b0` factor. It's caused by `timesbigfloat`.


---

**[bugs:#4644] bfloat((x+1)/x includes unnecessary 1.0b0 factor**

**Status:** open
**Group:** None
**Labels:** bfloat 
**Created:** Wed Dec 03, 2025 07:15 PM UTC by Stavros Macrakis
**Last Updated:** Sat Dec 06, 2025 03:33 AM UTC
**Owner:** nobody


~~~
bfloat((x+2)/x) => (1.0b0*(x+2.0b0))/x
~~~
What is the ``1.0b0`` factor doing there?

Set as priority 1 (which I hope means lowest, not highest) because this probably isn't doing any harm.



---

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] #4644 bfloat((x+1)/x includes unnecessary 1.0b0 factor

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

The same thing happens with `bfloat(x/y)`.


---

**[bugs:#4644] bfloat((x+1)/x includes unnecessary 1.0b0 factor**

**Status:** open
**Group:** None
**Labels:** bfloat 
**Created:** Wed Dec 03, 2025 07:15 PM UTC by Stavros Macrakis
**Last Updated:** Wed Dec 03, 2025 07:15 PM UTC
**Owner:** nobody


~~~
bfloat((x+2)/x) => (1.0b0*(x+2.0b0))/x
~~~
What is the ``1.0b0`` factor doing there?

Set as priority 1 (which I hope means lowest, not highest) because this probably isn't doing any harm.



---

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: #4645 rectform/polarform complex division overflow

[Raymond T. <rt...@us...>]

I'm not aware of any Lisp that does anything fancy with complex division other than taking care not to overflow or underflow unnessarily.  My limited understanding is for the worst cases, you need to use extra precision., as your bfloat result shows.

A Google search returns this [paper on HW complex division that is correctly rounded](http://arith.cs.ucla.edu/publications/CmplexDiv-ASAP03.pdf).  I very briefly skimmed it. It's pretty complicated and I can't say I understand it.


---

**[bugs:#4645] rectform/polarform complex division overflow**

**Status:** open
**Group:** None
**Labels:** rectform 
**Created:** Wed Dec 03, 2025 07:23 PM UTC by Stavros Macrakis
**Last Updated:** Thu Dec 04, 2025 06:19 PM UTC
**Owner:** nobody


~~~
fpprintprec:5$

q: (1.13e-100*%i+5.43e-10)/(5.7e-312*%i+1.2e-311)
rectform(q) => OVERFLOW
float(rectform(bfloat(q))) => 3.692e+301-1.7537e+301*%i

polarform(q) => OVERFLOW
polarform(bfloat(q)) => 4.0873b301*%e^-(4.4345b-1*%i)
~~~



---

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] #4645 rectform/polarform complex division overflow

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

I tried using SBCL native complexes (following Barton's idea) and I don't get overflow, but the answer is off by a hundred ULPs, which is the sort of thing to expect with "naive" complex float division.
~~~
(defun $cldiv (a b)
  (let* ((recta (risplit a))
	 (ac (complex (car recta) (cdr recta))) ; will err out if not numeric
	 (rectb (risplit b))
	 (bc (complex (car rectb) (cdr rectb)))
	 (q (/ ac bc)))
    (add (realpart q) (mul '$%i (imagpart q)))))

bfdiv(a,b):=block([fpprec:32],float(rectform(bfloat(a)/bfloat(b))));
bfcabs(a):=block([fpprex:100],float(cabs(bfloat(a))));

clq: cldiv(1.13e-100*%i+5.43e-10,5.7e-312*%i+1.2e-311);
   => 3.691993880674702e301 - 1.7536970933199895e301 %i

bfq: bfdiv(1.13e-100*%i+5.43e-10,5.7e-312*%i+1.2e-311);
   =>3.6919938806746145e301 - 1.7536970933199478e301 %i
   
bfdiv(clq,bfq)
   => 1.0000000000000238 - 1.0515920681324065e-17 %i

qr: bfcabs(%) => 1.0000000000000238
(qr-1)/unit_in_last_place(qr) => 107
~~~




---

**[bugs:#4645] rectform/polarform complex division overflow**

**Status:** open
**Group:** None
**Labels:** rectform 
**Created:** Wed Dec 03, 2025 07:23 PM UTC by Stavros Macrakis
**Last Updated:** Wed Dec 03, 2025 07:23 PM UTC
**Owner:** nobody


~~~
fpprintprec:5$

q: (1.13e-100*%i+5.43e-10)/(5.7e-312*%i+1.2e-311)
rectform(q) => OVERFLOW
float(rectform(bfloat(q))) => 3.692e+301-1.7537e+301*%i

polarform(q) => OVERFLOW
polarform(bfloat(q)) => 4.0873b301*%e^-(4.4345b-1*%i)
~~~



---

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: #4641 is does not work

[Raymond T. <rt...@us...>]

Thanks.  I rarely use matrices for anything, so I never remember all the associated functions.  I did now matrices were basically a list of lists.


---

**[bugs:#4641] is does not work**

**Status:** not-a-bug
**Group:** None
**Created:** Mon Dec 01, 2025 04:03 PM UTC by dan hayes
**Last Updated:** Wed Dec 03, 2025 10:43 PM UTC
**Owner:** nobody


build_info()  or bug_report()
"branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks")

(t:matrix([1],[2]),t,t[2],is( t[2]>0));   

and it comes back with unknown



---

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: #4641 is does not work

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

> I do not know why t[2] doesn't return a matrix.

Maxima defines a matrix as being composed of rows, each of which is a vector/list.
t[2] is the second row, and is a list, not a matrix.
If you want the second row as a matrix, you can use row(M,2), which gives a nx1 matrix, not a list. col(M,2) gives an 1xn matrix. But, beware! You can modify the parent matrix via the rows, but not via the columns.

I'd like to say that this is all beautifully designed, but frankly it was probably just an expedient way of implementing basic matrix functionality.




---

**[bugs:#4641] is does not work**

**Status:** not-a-bug
**Group:** None
**Created:** Mon Dec 01, 2025 04:03 PM UTC by dan hayes
**Last Updated:** Wed Dec 03, 2025 10:33 PM UTC
**Owner:** nobody


build_info()  or bug_report()
"branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks")

(t:matrix([1],[2]),t,t[2],is( t[2]>0));   

and it comes back with unknown



---

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] #4641 is does not work

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

You don't seem to be paying attention to what Raymond and I are telling you.
The list [0.1] is not the same as the number 0.1. 
Maxima does allow arithmetic on lists, so that 0.1*[2] yields [0.2]. But that is not equal to 0.2.



---

**[bugs:#4641] is does not work**

**Status:** not-a-bug
**Group:** None
**Created:** Mon Dec 01, 2025 04:03 PM UTC by dan hayes
**Last Updated:** Wed Dec 03, 2025 10:33 PM UTC
**Owner:** nobody


build_info()  or bug_report()
"branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks")

(t:matrix([1],[2]),t,t[2],is( t[2]>0));   

and it comes back with unknown



---

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: #4641 is does not work

[Raymond T. <rt...@us...>]

I think that's caused by two things.  `t[2]` is a list, not matrix.  And `listarith` idefaults to `true`, so list arithmetic is done.  If you set `listarith` to `false`,  then I get:
```
(%i17) listarith:false;
(%o17)                               false
(%i18) .1/(t[2]-1);
                                      0.1
(%o18)                              ───────
                                    [2] - 1
```
I do not know why `t[2]` doesn't return a matrix.


---

**[bugs:#4641] is does not work**

**Status:** not-a-bug
**Group:** None
**Created:** Mon Dec 01, 2025 04:03 PM UTC by dan hayes
**Last Updated:** Wed Dec 03, 2025 03:37 PM UTC
**Owner:** nobody


build_info()  or bug_report()
"branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks")

(t:matrix([1],[2]),t,t[2],is( t[2]>0));   

and it comes back with unknown



---

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: #4641 is does not work

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

 The reason i think it should work because(t:matrix([1],[2]),t,t[2],.1/( t[2]-1));does come back with the correct answer using it as i should think without any t[2][1] or first(t[2]) etc.so it should do so in(t:matrix([1],[2]),t,t[2],is( t[2]=2));rather than giving false and (t:matrix([1],[2]),is(abs(t[2])>0));should come back true rather than unknown. To be consistent if it workedcorrectly in the numeric case it should do the same in these other cases. It is not consistent if it didn't.I would have written this in the forum but i could not find it in bugs. I could not even get a list of bugs so i wrote it here.
    On Wednesday, December 3, 2025 at 09:37:42 AM CST, Stavros Macrakis via Maxima-bugs <max...@li...> wrote:  
 
 
Perhaps you should look at the value of flatten(t)[2]:
t:matrix([1],[2])$
flatten(t)[2] => [2]

In some systems (e.g. R), a singleton value is the same as a vector of length 1. That is not true in Maxima or most other systems, as you can verify:
is([2]=2) => false

"Of coourse it should work like it is supposed to without either any of these extra things."
Comments like this are inappropriate in a bug report, for multiple reasons:
* It is obnoxiously phrased. "Of course" is presumptuous.
* It presumes that your (mis)understanding of the system is what the system is "supposed" to do.
* It suggests that using the correct functions is "extra things".

Please try to use more neutral language and recognize that we are here to help you. Attacking us as incompetent is not helpful.

[bugs:#4641] is does not work

Status: not-a-bug
Group: None
Created: Mon Dec 01, 2025 04:03 PM UTC by dan hayes
Last Updated: Wed Dec 03, 2025 02:51 PM UTC

Owner: nobody
build_info() or bug_report()
"branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks")

(t:matrix([1],[2]),t,t[2],is( t[2]>0)); 

and it comes back with unknown

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 mailing list
Max...@li...
https://lists.sourceforge.net/lists/listinfo/maxima-bugs
  


---

**[bugs:#4641] is does not work**

**Status:** not-a-bug
**Group:** None
**Created:** Mon Dec 01, 2025 04:03 PM UTC by dan hayes
**Last Updated:** Wed Dec 03, 2025 03:37 PM UTC
**Owner:** nobody


build_info()  or bug_report()
"branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks")

(t:matrix([1],[2]),t,t[2],is( t[2]>0));   

and it comes back with unknown



---

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: #4641 is does not work

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

 The reason i think it should work because(t:matrix([1],[2]),t,t[2],.1/( t[2]-1));does come back with the correct answer using it as i should think without any t[2][1] or first(t[2]) etc.so it should do so in(t:matrix([1],[2]),t,t[2],is( t[2]=2));rather than giving false should be true and (t:matrix([1],[2]),is(abs(t[2])>0));should come back true rather than unknown. To be consistent if it workedcorrectly in the numeric case it should do the same in these other cases. It is not consistent if it didn't.I would have written this in the forum but i could not find it in bugs. I could not even get a list of bugs so i wrote it here.

    On Monday, December 1, 2025 at 10:10:50 AM CST, Raymond Toy <rt...@us...> wrote:  
 
 
With your example t[2] is [2], a 1x1 matrix. (Well, a list of one element). I think unknown is a good answer. Perhaps you really wanted is(t[2][1]>0). This returns true.

[bugs:#4641] is does not work

Status: open
Group: None
Created: Mon Dec 01, 2025 04:03 PM UTC by dan hayes
Last Updated: Mon Dec 01, 2025 04:03 PM UTC
Owner: nobody

build_info() or bug_report()
"branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks")

(t:matrix([1],[2]),t,t[2],is( t[2]>0)); 

and it comes back with unknown

Sent from sourceforge.net because you indicated interest in https://sourceforge.net/p/maxima/bugs/4641/

To unsubscribe from further messages, please visit https://sourceforge.net/auth/subscriptions/
     


---

**[bugs:#4641] is does not work**

**Status:** not-a-bug
**Group:** None
**Created:** Mon Dec 01, 2025 04:03 PM UTC by dan hayes
**Last Updated:** Wed Dec 03, 2025 03:37 PM UTC
**Owner:** nobody


build_info()  or bug_report()
"branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks")

(t:matrix([1],[2]),t,t[2],is( t[2]>0));   

and it comes back with unknown



---

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] #4645 rectform/polarform complex division overflow

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

---

**[bugs:#4645] rectform/polarform complex division overflow**

**Status:** open
**Group:** None
**Labels:** rectform 
**Created:** Wed Dec 03, 2025 07:23 PM UTC by Stavros Macrakis
**Last Updated:** Wed Dec 03, 2025 07:23 PM UTC
**Owner:** nobody


~~~
fpprintprec:5$

q: (1.13e-100*%i+5.43e-10)/(5.7e-312*%i+1.2e-311)
rectform(q) => OVERFLOW
float(rectform(bfloat(q))) => 3.692e+301-1.7537e+301*%i

polarform(q) => OVERFLOW
polarform(bfloat(q)) => 4.0873b301*%e^-(4.4345b-1*%i)
~~~



---

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] #4644 bfloat((x+1)/x includes unnecessary 1.0b0 factor

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

---

**[bugs:#4644] bfloat((x+1)/x includes unnecessary 1.0b0 factor**

**Status:** open
**Group:** None
**Labels:** bfloat 
**Created:** Wed Dec 03, 2025 07:15 PM UTC by Stavros Macrakis
**Last Updated:** Wed Dec 03, 2025 07:15 PM UTC
**Owner:** nobody


~~~
bfloat((x+2)/x) => (1.0b0*(x+2.0b0))/x
~~~
What is the ``1.0b0`` factor doing there?

Set as priority 1 (which I hope means lowest, not highest) because this probably isn't doing any harm.



---

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: #4641 is does not work

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

Perhaps you should look at the value of flatten(t)[2]:
~~~
t:matrix([1],[2])$
flatten(t)[2] => [2]
~~~
In some systems (e.g. R), a singleton value is the same as a vector of length 1. That is not true in Maxima or most other systems, as you can verify:
~~~
is([2]=2) => false
~~~

"Of coourse it should work like it is supposed to without either any of these extra things."
Comments like this are inappropriate in a bug report, for multiple reasons:
* It is obnoxiously phrased. "Of course" is presumptuous.
* It presumes that your (mis)understanding of the system is what the system is "supposed" to do.
* It suggests that using the correct functions is "extra things".

Please try to use more neutral language and recognize that we are here to help you. Attacking us as incompetent is not helpful.


---

**[bugs:#4641] is does not work**

**Status:** not-a-bug
**Group:** None
**Created:** Mon Dec 01, 2025 04:03 PM UTC by dan hayes
**Last Updated:** Wed Dec 03, 2025 02:51 PM UTC
**Owner:** nobody


build_info()  or bug_report()
"branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks")

(t:matrix([1],[2]),t,t[2],is( t[2]>0));   

and it comes back with unknown



---

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.