Lines Matching refs:csum
2515 double x, scale, oldcsum, csum = 1.0, frac1 = 0.0, frac2 = 0.0, frac3 = 0.0; in vector_norm() local
2548 assert(fabs(csum) >= fabs(x)); in vector_norm()
2549 oldcsum = csum; in vector_norm()
2550 csum += x; in vector_norm()
2551 frac1 += (oldcsum - csum) + x; in vector_norm()
2554 assert(fabs(csum) >= fabs(x)); in vector_norm()
2555 oldcsum = csum; in vector_norm()
2556 csum += x; in vector_norm()
2557 frac2 += (oldcsum - csum) + x; in vector_norm()
2559 assert(csum + lo * lo == csum); in vector_norm()
2562 h = sqrt(csum - 1.0 + (frac1 + frac2 + frac3)); in vector_norm()
2571 assert(fabs(csum) >= fabs(x)); in vector_norm()
2572 oldcsum = csum; in vector_norm()
2573 csum += x; in vector_norm()
2574 frac1 += (oldcsum - csum) + x; in vector_norm()
2577 assert(fabs(csum) >= fabs(x)); in vector_norm()
2578 oldcsum = csum; in vector_norm()
2579 csum += x; in vector_norm()
2580 frac2 += (oldcsum - csum) + x; in vector_norm()
2583 assert(fabs(csum) >= fabs(x)); in vector_norm()
2584 oldcsum = csum; in vector_norm()
2585 csum += x; in vector_norm()
2586 frac3 += (oldcsum - csum) + x; in vector_norm()
2588 x = csum - 1.0 + (frac1 + frac2 + frac3); in vector_norm()
2600 assert(fabs(csum) >= fabs(x)); in vector_norm()
2601 oldcsum = csum; in vector_norm()
2602 csum += x; in vector_norm()
2603 frac1 += (oldcsum - csum) + x; in vector_norm()
2605 return max * sqrt(csum - 1.0 + frac1); in vector_norm()