Field smoothing and adaptive mesh generation

A cheap smoothing process is described which both improves results from a standard first order solver and forms the basis of an effective adaptive solution system. The resulting system is demonstrated on several devices and is shown to develop accurate results in the presence of mixed media and with a very coarse starting discretization. The adaptive solver can be made to converge on a specified output result; either global, e.g., inductance or force, or local, e.g., the flux density in a particular region. While the criterion for measuring error seems somewhat simplistic, it can generate good results, particularly when combined with a convergence criterion based on the desired results of the analysis