Lines Matching refs:upper
162 double* diag, double* upper, double* b) in TriDiagGS() argument
257 diag[i] -= lower[i] * upper[i-1]; 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
383 lower[0] = upper[n-1] = 0.0; /* Matrix berechnen. */ in ZyklTriDiagGS()
388 upper[0] *= temp; /* MACH_EPS, so ex. keine */ in ZyklTriDiagGS()
392 { diag[i] -= lower[i] * upper[i-1]; in ZyklTriDiagGS()
395 upper[i] *= temp; in ZyklTriDiagGS()
399 diag[n-2] -= lower[n-2] * upper[n-3]; in ZyklTriDiagGS()
403 lowrow[i] = -lowrow[i-1] * upper[i-1]; in ZyklTriDiagGS()
405 lower[n-1] -= lowrow[n-3] * upper[n-3]; in ZyklTriDiagGS()
406 upper[n-2] = ( upper[n-2] - lower[n-2] * ricol[n-3] ) / diag[n-2]; in ZyklTriDiagGS()
410 diag[n-1] += temp - lower[n-1] * upper[n-2]; 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()