Rapid computation of magnetic anomalies with demagnetization included, for arbitrarily shaped magnetic bodies