In this work, we propose an efficient matrix decomposition algorithm for the Method of Fundamental Solutions when applied to three–dimensional boundary value problems governed by elliptic systems of partial differential equations. In particular, we consider problems arising in linear elasticity in axisymmetric domains. The proposed algorithm exploits the block circulant structure of the coeffic...