The proposed algorithm combines the optimal convergence rate of the Bramble
-Pasciak-Xu (BPX) preconditioner and a natural treatment of local refinemen
t by the Hierarchical Basis Multigrid (HBM) method. The BPX-HBM combi natio
n blends the best features of both methods and is applied to 3D magnetostat
ic scalar potential problems with first order tetrahedral elements. Adaptiv
e local mesh refinement is incorporated in the overall computational proces
s.