我目前正在为径向对称系统实施经典密度泛函理论。用数学术语来说，我正在寻找一个函数来最小化一个函数。的最简单部分如下所示：，但还有其他涉及 . 我正在尝试使用nlopt的数值实现，它在平面坐标中适用于相同的功能。$f(r)$$\Omega[f]$$\Omega$$\int\limits_0^R r^2 f(r)(\ln(f(r))-1)\mathrm{d}r$$f$$\Omega$

$f(r)$表示为向量在均匀间隔点处的样本。这些点从Δr开始，一直到R。f(0)无关紧要。有一个常数边界条件f(r>R)=f_∞。$f_i$$f$$r_i$$Δr$$R$$f(0)$$f(r>R)=f_∞$

但是，在径向坐标中，由于因子r^2（雅可比行列式） ，小$r$的收敛性比大$r$如果我只关心\Omega的值，这不会有问题，但在这种情况下，\Omega只是获得f的工具。$r^2$$\Omega$$\Omega$$f$

是否有更适合此类问题的库？还是有其他方法可以解决问题并继续使用 nlopt？