components/scientific/gdl/patches/08-brent.cpp.patch

   1 --- gdl-1.0.2/src/brent.cpp.orig        2023-01-06 18:39:47.000000000 +0100
   2 +++ gdl-1.0.2/src/brent.cpp     2023-02-05 15:14:16.437488573 +0100
   3 @@ -41,7 +41,7 @@
   4  double praxis ( double t0, double h0, int n, int prin, double x[],
   5    double f ( double x[], int n ) );
   6  void print2 ( int n, double x[], int prin, double fx, int nf, int nl );
   7 -void quad ( int n, double f ( double x[], int n ), double x[], double t,
   8 +void quad2 ( int n, double f ( double x[], int n ), double x[], double t,
   9    double h, double v[], double q0[], double q1[], int &nl, int &nf, double dmin,
  10    double ldt, double &fx, double &qf1, double &qa, double &qb, double &qc,
  11    double &qd0, double &qd1 );
  12 @@ -1413,7 +1413,7 @@
  13  //
  14  //  Try quadratic extrapolation in case we are in a curved valley.
  15  //
  16 -    quad ( n, f, x, t, h, v, q0, q1, nl, nf, dmin, ldt, fx, qf1,
  17 +    quad2 ( n, f, x, t, h, v, q0, q1, nl, nf, dmin, ldt, fx, qf1,
  18        qa, qb, qc, qd0, qd1 );
  19
  20      for ( j = 0; j < n; j++ )
  21 @@ -1653,7 +1653,7 @@
  22  }
  23  //****************************************************************************80
  24
  25 -void quad ( int n, double f ( double x[], int n ), double x[], double t,
  26 +void quad2 ( int n, double f ( double x[], int n ), double x[], double t,
  27    double h, double v[], double q0[], double q1[], int &nl, int &nf, double dmin,
  28    double ldt, double &fx, double &qf1, double &qa, double &qb, double &qc,
  29    double &qd0, double &qd1 )