compute aux with less roundoff discrepancies