A fast algorithm to predict elastic fields due to arbitrarily shaped eigenstrains in an elastic, isotropic half-space is advanced in this paper. The inclusion domain is partitioned in a set of cuboids of uniform eigenstrains, ...