[Maxima-commits] [git] Maxima CAS branch master updated. branch-5_48-base-185-g3a9a37f76

[peterpall <pet...@us...>]

This is an automated email from the git hooks/post-receive script. It was
generated because a ref change was pushed to the repository containing
the project "Maxima CAS".

The branch, master has been updated
       via  3a9a37f764757dddbc0028f6abde43631f80be9e (commit)
      from  0f29b7f5a0ae8eda3bbf81b882e1dbad793aa6e2 (commit)

Those revisions listed above that are new to this repository have
not appeared on any other notification email; so we list those
revisions in full, below.

- Log -----------------------------------------------------------------
commit 3a9a37f764757dddbc0028f6abde43631f80be9e
Author: Gunter Königsmann <gu...@pe...>
Date:   Thu Sep 11 09:29:47 2025 +0200

    wrstcse: taught the rss methods to deal with asymmetric limits
    
    If a diff() returns 'ind we now use the limit with the bigger abs()
    value instead of generating an error message.

diff --git a/share/contrib/wrstcse.mac b/share/contrib/wrstcse.mac
index b95edae5c..7b9d65cb5 100644
--- a/share/contrib/wrstcse.mac
+++ b/share/contrib/wrstcse.mac
@@ -215,7 +215,20 @@ wc_getrssparams(x):=block([tols:[],derivates_linearized,derivates,sigmaval,meanv
     meanval:wc_typicalvalues(x),
     /* Create a list of derivatives, linearized over the point that all tol[x] are 0.*/
     derivates:makelist(block([derivate:limit(%wc_typicalvalues(diff(x,i),i),i,0)],
-        if derivate='ind then error("wc_getrssparams: Cannot get meaningful normal distribution params for ",i),
+        if derivate='ind then block([wc_poslimit,wc_neglimit],
+           /* We have two derivates. As we are performing a worst case
+              calculation we make a worst-case assumption and use the bigger
+              one ignoring its sign */
+           wc_poslimit:limit(wc_diffofwctol,wc_tol,0,'plus),
+           wc_neglimit:limit(wc_diffofwctol,wc_tol,0,'minus),
+           derivate:max(abs(wc_poslimit),abs(wc_neglimit)),
+           /* If we can determine which of the limits is bigger we can
+              do without the abs() */
+           if abs(wc_poslimit) > abs(wc_neglimit) then
+              derivate:wc_poslimit,
+           if abs(wc_poslimit) < abs(wc_neglimit) then
+              derivate:wc_neglimit
+           ),
         derivate),
         i,tols),
     derivates_linearized:makelist(taytorat(taylor(i)),i,derivates),

-----------------------------------------------------------------------

Summary of changes:
 share/contrib/wrstcse.mac | 15 ++++++++++++++-
 1 file changed, 14 insertions(+), 1 deletion(-)


hooks/post-receive
-- 
Maxima CAS