I am trying to add a spherical harmonics solver model (Keaveny and Maxey) for magnetic force calculation between paramagnetic particles in a uniform magnetic field. Currently, I have the solver written in python using many scipy libraries like linear system solver, combination formula, and associated Legendre functions. Is there a simple way to write these in c++ for LIGGGHTS, or can I include some other libraries like Eigen which may help me?
Eric E. Keaveny, Martin R. Maxey,
Modeling the magnetic interactions between paramagnetic beads in magnetorheological fluids,
Journal of Computational Physics,
Volume 227, Issue 22,