Implementation aspects of 3d lattice-bgk: Boundaries, accuracy and a new fast relaxation method