Scattering amplitudes computed at a fixed loop order, along with any other object computed in perturbative QFT, can be expressed as a linear combination of a finite basis of loop integrals. To compute loop amplitudes in practise, such a basis of integrals must be determined. In this talk I introduce a new algorithm for finding bases of loop integrals and discuss its implementation in the publically available package Azurite.