This paper presents derivation of the three-dimensional (3-D) fast multipol
e method for use in calculating magnetostatic fields. Many other derivation
s use complex, inefficient, and awkward function normalizations. In this de
rivation, tenser functions and an operator are defined that greatly simplif
y the theory and implementation of the algorithm. This paper provides proof
of O(N) scaling, implementation, and parallelization tips, performance res
ults, and an error analysis.