Returns a finite difference representation of the Kuprov operator:
with the number of stencil points in the finite difference approximation specified by user. The resulting operator is a sparse matrix designed to act on the vectorisation of rho. The dimensions of rho are assumed to be ordered as [X Y Z]. For further information, see http://dx.doi.org/10.1039/C4CP03106G.
npoints - a three-element vector specifying the dimensions of the 3D cube of data that the operator will be acting on, in Angstroms. The dimensions are assu- med to be ordered as [X Y Z]. chi - the electron magnetic susceptibility tensor in cubic Angstroms, a symmetric 3x3 matrix. extents - a three-element vector specifying axis extents in Angstroms. The dimensions are assumed to be ordered as [X Y Z]. nstenc - number of finite-difference stencil points for the finite-difference approximation. Periodic boundary conditions are used.
K - a sparse matrix designed to act on the vectori- zation of the array. The dimensions are assumed to be ordered as [X Y Z].