Efficient implementation of the domain-integrated field relations method for quasi-static magnetic fields