fn sub_limbs(a: &[u64; 4], b: &[u64; 4]) -> ([u64; 4], bool)
Subtracts two 4-limb numbers, returning (result, borrow).