CompensatedBLAS: Higher precision BLAS with low precisions

Here I look at “Compensated” arithmetic and its application to a small side-project CompensatedBLAS . Compensated arithmetic allows you to implement higher-accuracy accumulation using expansions of lower-precision floating-point words. Two FP32 words carry at most roughly 48 significand bits before nonoverlap and implementation losses, so they do not in general reproduce the 53-bit significand or…