Lines Matching refs:b
99 double x = 1.0, one = 1.0, b = 1.0;
102 while ( (x + b) == x ) b *= 2.0;
104 return ( (short) ((x + b) - x) );
162 double* diag, double* upper, double* b) in TriDiagGS() argument
264 b[i] -= lower[i] * b[i-1]; in TriDiagGS()
266 b[n-1] /= diag[n-1]; /* Rueckwaertselimination */ in TriDiagGS()
269 b[i] = ( b[i] - upper[i] * b[i+1] ) / diag[i]; in TriDiagGS()
294 double* upper, double* lowrow, double* ricol, double* b) in ZyklTriDiagGS() argument
415 b[0] /= diag[0]; /* Vorwaertselemination */ in ZyklTriDiagGS()
417 b[i] = ( b[i] - b[i-1] * lower[i] ) / diag[i]; in ZyklTriDiagGS()
420 temp -= lowrow[i] * b[i]; in ZyklTriDiagGS()
422 b[n-1] = ( b[n-1] + temp - lower[n-1] * b[n-2] ) / diag[n-1]; in ZyklTriDiagGS()
424 b[n-2] -= b[n-1] * upper[n-2]; /* Rueckwaertselimination */ in ZyklTriDiagGS()
427 b[i] -= upper[i] * b[i+1] + ricol[i] * b[n-1]; in ZyklTriDiagGS()
452 double* b, double* c, double* d) in NaturalSpline() argument
469 b[i]=h[i]; in NaturalSpline()
484 b[n-2]=b[n-2]-h[n-1]; in NaturalSpline()
507 error=TriDiagGS(sal_False,n-1,b,d,c,a); in NaturalSpline()
537 b[i]=(y[i+1]-y[i])/h[i]-h[i]*(c[i+1]+2.0*c[i])/3.0; in NaturalSpline()
559 double* b, double* c, double* d) in PeriodicSpline() argument
587 b[im1]=hl; in PeriodicSpline()
594 b[nm1]=hl; in PeriodicSpline()
599 Error=ZyklTriDiagGS(sal_False,n,b,d,c,lowrow,ricol,a); in PeriodicSpline()
612 b[i]=(y[i+1]-y[i])/hl; in PeriodicSpline()
613 b[i]=b[i]-hl*(c[i+1]+2.0*c[i])/3.0; in PeriodicSpline()