Fast matrix methods for calculating the relaxation about defects in crystals