// Initializes a convolution to evaluate the demag field for the given mesh geometry.
func NewMFM(mesh *data.Mesh, lift, tipsize float64) *MFMConvolution {
k := mag.MFMKernel(mesh, lift, tipsize)
size := mesh.Size()
c := new(MFMConvolution)
c.size = size
c.kern = k
c.kernSize = k[X].Size()
c.init()
c.mesh = mesh
return c
}
func (c *MFMConvolution) Reinit(lift, tipsize float64) {
c.kern = mag.MFMKernel(c.mesh, lift, tipsize)
c.initFFTKern3D()
}