Hybrid method for computing demagnetizing fields

A hybrid finite-element-boundary integral method is proposed for computing demagnetizing field. For a three-dimensional (3-D) mesh with N nodes in the magnetic region, the method requires storage O(N/sup 4/3/), and no nodes outside the magnetic region need to be considered. The method is suitable for nonlinear calculations, can be adapted to handle curved boundaries, and is especially convenient when the magnetic region consists of several parts whose relative positions are variable. >