Fast stray field computation on tensor grids
  A direct integration algorithm is described to compute the magnetostatic field and energy for given magnetization distributions on not necessarily uniform tensor grids. We use an analytically-based tensor approximation approach for function-related tensors, which reduces calculations to multilinear algebra operations. The algorithm scales with N^(4/3) for N computational cells used and with N^(2/3) (sublinear) when magnetization is given in canonical tensor format. In the final section we confirm our theoretical results concerning computing times and accuracy by means of numerical examples. 
Lukas Exl, Winfried Auzinger, Simon Bance, Markus Gusenbauer, Franz   Reichel, Thomas Schrefl
physics.comp-ph math.NA