Merge branch 'master' into github-ci

Author: rtoy

ChangeLog

@@ -8,6 +8,8 @@ New items in core:
  * New command-line option --quit-on-error to make Maxima quit upon first
    runtime error
  * New command-line option --very-very-quiet to silence Maxima
+ * New command-line option --suppress-input-echo to suppress input echo
+   when running noninteractively
  * New global variable batch_answers_from_file
  * New Unicode display2d and global variable display2d_unicode
  * A new function random-name in plot.lisp that generates a random string
@@ -32,14 +34,19 @@ Changes in core:
  * In src/plot.lisp, clipping is no longer done when plot_format is xmaxima
  * The default y-range for plots is determined only by points within x-range
  * Function chinese renamed to solve_congruences
+ * Function 'system': capture output (e.g. with_stdout) from executed command
 
 Changes in share:
 -----------------
+ * openr_binary: work around bug in Allegro CL by omitting :ELEMENT-TYPE argument.
+ * readline: trim carriage return from end of string, if any.
 
 Changes in Xmaxima:
 -------------------
  * Line clipping in plot2d  and plotdf
  * Better integration algorithm for plotdf
+ * new options tstep and algorithm for plotdf
+ * plot linewidths and arrow lengths change when a plot window is resized
 
 Changes in the Windows installer:
 ---------------------------------
@@ -62,8 +69,10 @@ Bug fixes for numbered bugs:
  * \#2905: Assigning variable twice yields different results
  * \#3054: limit(exp(exp(2*log(x**5 + x)*log(log(x)))) / exp(exp(10*log(x)*log(log(x)))), x, inf)
  * \#3191: parsing problem: thru 3 for i in [a,b]
+ * \#3579: Overly eager simplification in trigreduce
  * \#3755: Insecure tmpdir usage
  * \#3803: Xmaxima should reopen the browser when needed
+ * \#3981: Trigreduce doesn't return principal value of atan(tan(...))
  * \#4117: Unexpected Behaviour of run_viewer and gnu_term dumb
  * \#4154: Some info nodes have bad headers (and no "Up"s)
  * \#4160: LISP error when integrating
@@ -74,6 +83,7 @@ Bug fixes for numbered bugs:
  * \#4191: limits of gamma_incomplete
  * \#4195: cmucl 21e intermittent lisp error in limit
  * \#4200: simplim%limit ignoring special variables
+ * \#4208: Error with complex numbers - example with trigreduce
  * \#4210: tlimit consults global assumptions for limit variable
  * \#4214: read_xpm not autoloading and not documented as requiring load("draw")
  * \#4215: limit(sqrt(-1+%i*x),x,0,minus)
@@ -103,6 +113,7 @@ Bug fixes for numbered bugs:
  * \#4388: tex() output under simp:false, incorrect placement of `+`
  * \#4408: letsimp doesn't simplify before returning result
  * \#4417: letsimp treats minus expressions as -1 times something
+ * \#4418: arctan and trigreduce error
  * \#4424: draw2d doesnt accept fill_density in order to create opaque shapes
  * \#4426: Remove repeated function resm1
  * \#4427: Turn number to cell in [email protected]
@@ -120,12 +131,16 @@ Bug fixes for numbered bugs:
 
 Unnumbered bugs fixed:
 ---------------------
- * Correct failing test in share/contrib/diffequations/tests
- * plotdf sliders hidden (reported on the mailing list on 2023/11/26)
  * mailing list 2023-11-25: Possibilities for basic image processing
+ * mailing list 2023-11-26: plotdf sliders hidden
  * mailing list 2023-09-25: Casting SPELs
- * recover option xfun of plotdf which was no longer working (2024-03-16)
  * mailing list 2024-12-24: Typographic Error (xmaxima)
+ * mailing list 2025-03-04: trigsolve requires to_poly_solve to work
+ * correct failing test in share/contrib/diffequations/tests
+ * recovers option xfun of plotdf which was no longer working
+ * fixes narrows option of plotdf (narrows should remain after window resizing)
+ * unreported bug: Clisp error when standard input is not writeable
+ * unreported bug: Function 'system': fix problems encountered when testing on Windows
 
 Documentation:
 --------------
@@ -137,3 +152,5 @@ Documentation:
 Build system:
 -------------
  * Install Maxima icons in an 'icons' subdirectory rather than 'pixmaps'
+ * maxima.asd: load defsystem.lisp in order to load share packages which use it.
+ * Various compatibility fixes for Allegro CL and LispWorks

doc/info/Numerical.texi

@@ -891,31 +891,26 @@ numerical method used is 4th order Runge-Kutta with variable time steps.
 
 @b{Plot window menu:}
 
-The menu bar of the plot window has the following seven icons:
+The menu bar of the plot window has the following five buttons:
 
-An X. Can be used to close the plot window.
+Close: can be used to close the plot window.
 
-A wrench and a screwdriver. Opens the configuration menu with several
+Config: opens the configuration menu with several
 fields that show the ODE(s) in use and various other settings. If a pair
 of coordinates are entered in the field @emph{Trajectory at} and the
 @key{enter} key is pressed, a new integral curve will be shown, in
 addition to the ones already shown.
 
-Two arrows following a circle. Replots the direction field with the
-new settings defined in the configuration menu and replots only the last
-integral curve that was previously plotted.
+Save: used to save a copy of the plot, in Postscript format, in the file
+specified in a field of the window that appears when that button is clicked.
 
-Hard disk drive with an arrow. Used to save a copy of the
-plot, in Postscript format, in the file specified in a field of the
-box that appears when that icon is clicked.
+Replot: replots the direction field with the new settings defined in the
+configuration menu and replots only the last integral curve that was
+previously plotted. If you just resized the plot window, the size and
+width of the arrows and curves will be adapted to the new size if you
+click on Replot.
 
-Magnifying glass with a plus sign. Zooms in the plot.
-
-Magnifying glass with a minus sign. Zooms out the plot. The plot can be
-displaced by holding down the right mouse button while the mouse is
-moved.
-
-Icon of a plot. Opens another window with a plot of the two variables
+Time plot: creates two new window showing the plots of the two variables
 in terms of time, for the last integral curve that was plotted.
 
 @b{Plot options:}
@@ -977,6 +972,12 @@ interactively using slider buttons, and the range of variation of those
 parameters. The names and ranges of the parameters must be given in a
 string with a comma-separated sequence of elements @code{name=min:max}
 
+@item
+@dfn{tstep} sets the value of the time intervals used in the integration
+algorithm. It must be a floating point number; you might have to adjust
+its value to get good results for the integral curves. If not given, a
+default value will be chosen according to the region to be plotted.
+
 @item
 @dfn{xfun} defines a string with semi-colon-separated sequence
 of functions of @var{x} to be displayed, on top of the direction field.

doc/info/commandline-options.texi

@@ -138,10 +138,6 @@ The following command line options are available for Maxima:
         Do not verify on startup that the set of html topics is consistent with
         text topics.
 
-@need 150
-@item     --suppress-input-echo
-        Do not print input expressions when processing noninteractively.
-
 
 @end table
 

interfaces/xmaxima/Tkmaxima/Plotdf.tcl

@@ -4,7 +4,7 @@
 # For distribution under GNU public License.  See COPYING. #
 #                                                          #
 #     Modified by Jaime E. Villate                         #
-#     Time-stamp: "2025-02-01 16:35:14 villate"            #
+#     Time-stamp: "2025-05-27 10:54:55 villate"            #
 ############################################################
 
 global plotdfOptions
@@ -107,15 +107,15 @@ proc doIntegrate { win x0 y0 } {
     set linewidth [expr {$linewidth*[vectorlength $width $height]/1000.}]
     set arrowshape [scalarTimesVector $linewidth {3 5 2}]
 
-    # method can be rungeKutta, rungeKuttaA or adamsMoulton
-    set method {adamsMoulton}
+    # integrator can be rungeKutta, rungeKuttaA or adamsMoulton
+    set integrator {adamsMoulton}
     oset $win trajectory_at [format "%.10g  %.10g" $x0 $y0]
     lappend trajectoryStarts [list $x0 $y0]
     set didLast {}
     # puts "doing at $trajectory_at"
-    # A reasonabel value of tstep has already been set up in drawDF by
-    # using the maximum length of the field vectors. This is just in case.
-    if {$tstep eq {}} {set tstep 0.1}
+    # Default value for tstep equal to the plot box's diagonal divided by 400
+    set steps $nsteps
+    set h $tstep
 
     set todo {1}
     switch -- $direction {
@@ -150,7 +150,7 @@ proc doIntegrate { win x0 y0 } {
                         set arrow {last}
                         set coords {}}}}
             set h [expr {$sgn*$tstep}]
-            set form [list $method xff yff $tinitial $x0 $y0 $h $nsteps]
+            set form [list $integrator xff yff $tinitial $x0 $y0 $h $nsteps]
 
             # puts "doing: $form"
             # pts will be a list with values of t, x and y, at the initial
@@ -263,7 +263,7 @@ proc drawArrowScreen { c atx aty dfx dfy color } {
 proc drawDF { win tinitial } {
     global axisGray
     makeLocal  $win xmin xmax xcenter ycenter c ymin ymax transform vectors \
-        xaxislabel yaxislabel nobox axes width height narrows tstep
+        xaxislabel yaxislabel nobox axes width height narrows
     set rtosx rtosx$win
     set rtosy rtosy$win
     set storx storx$win
@@ -295,7 +295,6 @@ proc drawDF { win tinitial } {
 		append all " $len $dfx $dfy "
 		if { $min > $len } {set min $len}
                         if { $max < $len } {set max $len}}}
-        if {$tstep eq {}} {oset $win tstep [expr {$stepsize/(10.0*$max)}]}
 	set arrowmin [expr {0.25*$stepsize}]
 	set arrowrange [expr {0.85*$stepsize - $arrowmin}]
 	set s1 [expr {($arrowrange*$min+$arrowmin*$min-$arrowmin*$max)/($min-$max)}]
@@ -399,7 +398,10 @@ proc plotdf { args } {
     getOptions $plotdfOptions $args -usearray [oarray $win]
     oset $win didLast {}
     # Makes extra vertical space for sliders
-    linkLocal $win sliders height
+    linkLocal $win sliders height tstep xradius yradius
+    if { "$tstep" == "" } {
+	set tstep [expr {[vectorlength $xradius $yradius] / 200.0}]
+    }
     if {[string length $sliders] > 0} {
         oset $win height [expr {$height + 40*[llength [split $sliders ,]]}]}
 
@@ -471,7 +473,7 @@ proc doConfigdf { win } {
     pack $frdydx.dxdt  $frdydx.dydt -side bottom  -fill x -expand 1
     pack $frdydx.dydxbut $frdydx.dydtbut -side left -fill x -expand 1
 
-    foreach w {narrows parameters xfun linewidth xradius yradius xcenter ycenter tinitial versus_t nsteps direction curves vectors fieldlines } {
+    foreach w {narrows parameters xfun linewidth xradius yradius xcenter ycenter tinitial versus_t tstep nsteps direction curves vectors fieldlines } {
 	mkentry $wb1.$w [oloc $win $w] $w $buttonFont
 	pack $wb1.$w -side bottom -expand 1 -fill x
     }