add 2025 BLAS sources

git-svn-id: https://svn.r-project.org/R/trunk@87891 00db46b3-68df-0310-9c12-caf00c1e9a41

Author: ripley

src/include/Defn.h

@@ -1714,7 +1714,7 @@ int Rf_FixupDigits(SEXP, warn_type);
 int Rf_FixupWidth (SEXP, warn_type);
 SEXP Rf_installDDVAL(int i);
 SEXP Rf_installS3Signature(const char *, const char *);
-Rboolean Rf_isFree(SEXP);
+bool Rf_isFree(SEXP);
 Rboolean Rf_isUnmodifiedSpecSym(SEXP sym, SEXP env);
 SEXP Rf_matchE(SEXP, SEXP, int, SEXP);
 // void Rf_setSVector(SEXP*, int, SEXP);
@@ -1914,8 +1914,8 @@ int	R_ShowFiles(int, const char **, const char **, const char *,
 int     R_EditFiles(int, const char **, const char **, const char *);
 int	R_ChooseFile(int, char *, int);
 char	*R_HomeDir(void);
-Rboolean R_FileExists(const char *);
-Rboolean R_HiddenFile(const char *);
+bool    R_FileExists(const char *);
+bool    R_HiddenFile(const char *);
 double	R_FileMtime(const char *);
 int	R_GetFDLimit(void);
 int	R_EnsureFDLimit(int);
@@ -2284,7 +2284,7 @@ int utf8clen(char c);
 int Rf_AdobeSymbol2ucs2(int n);
 double R_strtod5(const char *str, char **endptr, char dec,
 		 Rboolean NA, int exact);
-SEXP R_listCompact(SEXP s, Rboolean keep_initial);
+SEXP R_listCompact(SEXP s, bool keep_initial);
 
 typedef unsigned short R_ucs2_t;
 size_t mbcsToUcs2(const char *in, R_ucs2_t *out, int nout, int enc);
@@ -2321,7 +2321,7 @@ int Rasprintf_malloc(char **str, const char *fmt, ...)
 
 SEXP fixup_NaRm(SEXP args); /* summary.c */
 void invalidate_cached_recodings(void);  /* from sysutils.c */
-void resetICUcollator(Rboolean disable); /* from util.c */
+void resetICUcollator(bool disable); /* from util.c */
 void dt_invalidate_locale(void); /* from Rstrptime.h */
 extern int R_OutputCon; /* from connections.c */
 

src/main/complex.c

@@ -1,6 +1,6 @@
 /*
  *  R : A Computer Language for Statistical Data Analysis
- *  Copyright (C) 2000-2022  The R Core Team
+ *  Copyright (C) 2000-2025  The R Core Team
  *  Copyright (C) 2005       The R Foundation
  *  Copyright (C) 1995-1997  Robert Gentleman and Ross Ihaka
  *
@@ -585,18 +585,18 @@ static double complex z_atanh(double complex z)
     return -I * z_atan(z * I);
 }
 
-static Rboolean cmath1(double complex (*f)(double complex),
+static bool cmath1(double complex (*f)(double complex),
 		       const Rcomplex *x, Rcomplex *y, R_xlen_t n)
 {
     R_xlen_t i;
-    Rboolean naflag = FALSE;
+    bool naflag = false;
     for (i = 0 ; i < n ; i++) {
 	if (ISNA(x[i].r) || ISNA(x[i].i)) {
 	    y[i].r = NA_REAL; y[i].i = NA_REAL;
 	} else {
 	    SET_C99_COMPLEX(y, i, f(toC99(x + i)));
 	    if ( (ISNAN(y[i].r) || ISNAN(y[i].i)) &&
-		!(ISNAN(x[i].r) || ISNAN(x[i].i)) ) naflag = TRUE;
+		!(ISNAN(x[i].r) || ISNAN(x[i].i)) ) naflag = true;
 	}
     }
     return naflag;
@@ -606,7 +606,7 @@ attribute_hidden SEXP complex_math1(SEXP call, SEXP op, SEXP args, SEXP env)
 {
     SEXP x, y;
     R_xlen_t n;
-    Rboolean naflag = FALSE;
+    bool naflag = false;
 
     PROTECT(x = CAR(args));
     n = XLENGTH(x);
@@ -696,7 +696,7 @@ attribute_hidden SEXP complex_math2(SEXP call, SEXP op, SEXP args, SEXP env)
     Rcomplex ai, bi, *y;
     const Rcomplex *a, *b;
     SEXP sa, sb, sy;
-    Rboolean naflag = FALSE;
+    bool naflag = false;
     cm2_fun f;
 
     switch (PRIMVAL(op)) {
@@ -734,7 +734,7 @@ attribute_hidden SEXP complex_math2(SEXP call, SEXP op, SEXP args, SEXP env)
 	    f(&y[i], &ai, &bi);
 	    if ( (ISNAN(y[i].r) || ISNAN(y[i].i)) &&
 		 !(ISNAN(ai.r) || ISNAN(ai.i) || ISNAN(bi.r) || ISNAN(bi.i)) )
-		naflag = TRUE;
+		naflag = true;
 	}
     });
     if (naflag)
@@ -787,12 +787,12 @@ attribute_hidden SEXP do_complex(SEXP call, SEXP op, SEXP args, SEXP rho)
 }
 
 static void R_cpolyroot(double *opr, double *opi, int *degree,
-			double *zeror, double *zeroi, Rboolean *fail);
+			double *zeror, double *zeroi, bool *fail);
 
 attribute_hidden SEXP do_polyroot(SEXP call, SEXP op, SEXP args, SEXP rho)
 {
     SEXP z, zr, zi, r, rr, ri;
-    Rboolean fail;
+    bool fail;
     int degree, i, n;
 
     checkArity(op, args);
@@ -904,10 +904,10 @@ attribute_hidden SEXP do_polyroot(SEXP call, SEXP op, SEXP args, SEXP rho)
 
 #include <Rmath.h> /* for R_pow_di */
 
-static void calct(Rboolean *);
-static Rboolean fxshft(int, double *, double *);
-static Rboolean vrshft(int, double *, double *);
-static void nexth(Rboolean);
+static void calct(bool *);
+static bool fxshft(int, double *, double *);
+static bool vrshft(int, double *, double *);
+static void nexth(bool);
 static void noshft(int);
 
 static void polyev(int, double, double,
@@ -931,30 +931,30 @@ static const double mre = 2. * M_SQRT2 * /* eta, i.e. */DBL_EPSILON;
 static const double infin = DBL_MAX;
 
 static void R_cpolyroot(double *opr, double *opi, int *degree,
-			double *zeror, double *zeroi, Rboolean *fail)
+			double *zeror, double *zeroi, bool *fail)
 {
     static const double smalno = DBL_MIN;
     static const double base = (double)FLT_RADIX;
     static int d_n, i, i1, i2;
     static double zi, zr, xx, yy;
     static double bnd, xxx;
-    Rboolean conv;
+    bool conv;
     int d1;
     double *tmp;
     static const double cosr =/* cos 94 */ -0.06975647374412529990;
     static const double sinr =/* sin 94 */  0.99756405025982424767;
     xx = M_SQRT1_2;/* 1/sqrt(2) = 0.707.... */
 
     yy = -xx;
-    *fail = FALSE;
+    *fail = false;
 
     nn = *degree;
     d1 = nn - 1;
 
     /* algorithm fails if the leading coefficient is zero. */
 
     if (opr[0] == 0. && opi[0] == 0.) {
-	*fail = TRUE;
+	*fail = true;
 	return;
     }
 
@@ -1037,7 +1037,7 @@ static void R_cpolyroot(double *opr, double *opi, int *degree,
 	/* the zerofinder has failed on two major passes */
 	/* return empty handed */
 
-	*fail = TRUE;
+	*fail = true;
 	vmaxset(vmax);
 	return;
 
@@ -1113,7 +1113,7 @@ static void noshft(int l1)
  *  initiates a variable-shift iteration and returns with the
  *  approximate zero if successful.
  */
-static Rboolean fxshft(int l2, double *zr, double *zi)
+static bool fxshft(int l2, double *zr, double *zi)
 {
 /*  l2	  - limit of fixed shift steps
  *  zr,zi - approximate zero if convergence (result TRUE)
@@ -1123,7 +1123,7 @@ static Rboolean fxshft(int l2, double *zr, double *zi)
  * Uses global (sr,si), nn, pr[], pi[], .. (all args of polyev() !)
 */
 
-    Rboolean pasd, h_s_0, test;
+    bool pasd, h_s_0, test;
     static double svsi, svsr;
     static int i, j, n;
     static double oti, otr;
@@ -1134,8 +1134,8 @@ static Rboolean fxshft(int l2, double *zr, double *zi)
 
     polyev(nn, sr, si, pr, pi, qpr, qpi, &pvr, &pvi);
 
-    test = TRUE;
-    pasd = FALSE;
+    test = true;
+    pasd = false;
 
     /* calculate first t = -p(s)/h(s). */
 
@@ -1160,10 +1160,10 @@ static Rboolean fxshft(int l2, double *zr, double *zi)
 
 	if (!h_s_0 && test && j != l2) {
 	    if (hypot(tr - otr, ti - oti) >= hypot(*zr, *zi) * 0.5) {
-		pasd = FALSE;
+		pasd = false;
 	    }
 	    else if (! pasd) {
-		pasd = TRUE;
+		pasd = true;
 	    }
 	    else {
 
@@ -1179,14 +1179,14 @@ static Rboolean fxshft(int l2, double *zr, double *zi)
 		svsr = sr;
 		svsi = si;
 		if (vrshft(10, zr, zi)) {
-		    return TRUE;
+		    return true;
 		}
 
 		/* the iteration failed to converge. */
 		/* turn off testing and restore */
 		/* h, s, pv and t. */
 
-		test = FALSE;
+		test = false;
 		for (i=1 ; i<=n ; i++) {
 		    hr[i-1] = shr[i-1];
 		    hi[i-1] = shi[i-1];
@@ -1208,7 +1208,7 @@ static Rboolean fxshft(int l2, double *zr, double *zi)
 
 /* carries out the third stage iteration.
  */
-static Rboolean vrshft(int l3, double *zr, double *zi)
+static bool vrshft(int l3, double *zr, double *zi)
 {
 /*  l3	    - limit of steps in stage 3.
  *  zr,zi   - on entry contains the initial iterate;
@@ -1218,12 +1218,12 @@ static Rboolean vrshft(int l3, double *zr, double *zi)
  *
  * Assign and uses  GLOBAL sr, si
 */
-    Rboolean h_s_0, b;
+    bool h_s_0, b;
     static int i, j;
     static double r1, r2, mp, ms, tp, relstp;
     static double omp;
 
-    b = FALSE;
+    b = false;
     sr = *zr;
     si = *zi;
 
@@ -1254,7 +1254,7 @@ static Rboolean vrshft(int l3, double *zr, double *zi)
 		/* one zero to dominate. */
 
 		tp = relstp;
-		b = TRUE;
+		b = true;
 		if (relstp < eta)
 		    tp = eta;
 		r1 = sqrt(tp);
@@ -1275,7 +1275,7 @@ static Rboolean vrshft(int l3, double *zr, double *zi)
 		/* increases significantly. */
 
 		if (mp * .1 > omp)
-		    return FALSE;
+		    return false;
 	    }
 	}
 	omp = mp;
@@ -1292,15 +1292,15 @@ static Rboolean vrshft(int l3, double *zr, double *zi)
 	    si += ti;
 	}
     }
-    return FALSE;
+    return false;
 
 L_conv:
     *zr = sr;
     *zi = si;
-    return TRUE;
+    return true;
 }
 
-static void calct(Rboolean *h_s_0)
+static void calct(bool *h_s_0)
 {
     /* computes	 t = -p(s)/h(s).
      * h_s_0   - logical, set true if h(s) is essentially zero.	*/
@@ -1322,7 +1322,7 @@ static void calct(Rboolean *h_s_0)
     }
 }
 
-static void nexth(Rboolean h_s_0)
+static void nexth(bool h_s_0)
 {
     /* calculates the next shifted h polynomial.
      * h_s_0 :	if TRUE  h(s) is essentially zero

src/main/identical.c

@@ -103,7 +103,7 @@ attribute_hidden SEXP do_identical(SEXP call, SEXP op, SEXP args, SEXP env)
 #define EXTPTR_AS_REF   (flags & IDENT_EXTPTR_AS_REF)
 
 /* do the two objects compute as identical?
-   Also used in unique.c */
+   Also used in unique.c, eval.c, relop.c, dotcode.c, and exported */
 Rboolean
 R_compute_identical(SEXP x, SEXP y, int flags)
 {

src/main/objects.c

@@ -1028,7 +1028,7 @@ static SEXP inherits3(SEXP x, SEXP what, SEXP which)
 
     if( !isLogical(which) || (LENGTH(which) != 1) )
 	error(_("'which' must be a length 1 logical vector"));
-    Rboolean isvec = asRbool(which, R_NilValue);
+    bool isvec = asRbool(which, R_NilValue);
 
     if(isvec)
 	PROTECT(rval = allocVector(INTSXP, nwhat));
@@ -1231,7 +1231,7 @@ static SEXP R_isMethodsDispatchOn(SEXP onOff)
     R_stdGen_ptr_t old = R_get_standardGeneric_ptr();
     int ival =  !NOT_METHODS_DISPATCH_PTR(old);
     if(length(onOff) > 0) {
-	Rboolean onOffValue = asRbool(onOff, R_NilValue);
+	bool onOffValue = asRbool(onOff, R_NilValue);
 	if(onOffValue == NA_INTEGER)
 	    error(_("'onOff' must be TRUE or FALSE"));
 	else if(onOffValue == FALSE)
@@ -1782,7 +1782,7 @@ SEXP R_do_new_object(SEXP class_def)
     }
     PROTECT(e = R_do_slot(class_def, s_className));
     PROTECT(value = duplicate(R_do_slot(class_def, s_prototype)));
-    Rboolean xDataType = TYPEOF(value) == ENVSXP || TYPEOF(value) == SYMSXP ||
+    bool xDataType = TYPEOF(value) == ENVSXP || TYPEOF(value) == SYMSXP ||
 	TYPEOF(value) == EXTPTRSXP;
     if((TYPEOF(value) == OBJSXP || getAttrib(e, R_PackageSymbol) != R_NilValue) &&
        !xDataType)

src/main/plot.c

@@ -86,7 +86,7 @@ SEXP CreateAtVector(double axp[], const double usr[], int nint, Rboolean logflag
 	}
     }
     else { /* ------ log axis ----- */
-	Rboolean reversed = FALSE;
+	bool reversed = false;
 	double
 	    umin = usr[0],
 	    umax = usr[1];