Арифметические операции с числами четверной точности (binary 128)

Требуется реализовать возможность работы с числами повышенной точности без использования готовых библиотек, например, quadmath (__float128). Может быть, существуют какие-то дельные описания алгоритмов как производить сложение и умножение таких экзотических чисел? Несколько дней гугления пока не привели к положительным результатам :( Был бы признателен за любые идеи или наброски решений в виде кода.