Dear ELK developers and users, I am calculating the spontaneous electric polarization of NiI₂, a spin-spiral multiferroic, using ELK version 10.9.5. My setup uses spinsprl .true. with a spin-spiral q-vector, DFT+U on Ni (FLL, U=4.2 eV, J=0.8 eV), GGA, spin-orbit coupling, and a 1×1×1 primitive cell. After reading the full task list and grepping the source code, I found that ELK does not appear to have a standalone task for computing spontaneous Berry phase polarization. The subroutine polar() in...