势函数-LBFGS,程序代码问题

my_aman

请问各位老师大神： 我在用势函数+LBFGS，编写势函数势能量梯度代码，经过反复测试未得到预期的结构。
应该是梯度写的有问题，但是我就是没检查出来。现把代码贴出，请各位大神帮忙看看，应该怎么修改。
请各位一定要知无不言哦，先谢谢了。
上图是势函数，下图是梯度公式。其中代码在rm.c文件里，A,B,C,p为公式的参数。

zjxitcc

讲真，这梯度公式我反正没看懂怎么推出来的··· 楼下解解惑，我都怀疑这梯度公式错了

get-it

G是一个矢量，r_ij = \sqrt(x^2 + y^2 + z^2)

zjxitcc

get-it 发表于 2017-7-10 23:40
G是一个矢量，r_ij = \sqrt(x^2 + y^2 + z^2)

我试过矢量，也没推出来长这样。。

coolrainbow

-C/r^6    ===>   +6C/r^8 * (rx)

my_aman

get-it 发表于 2017-7-10 23:40
G是一个矢量，r_ij = \sqrt(x^2 + y^2 + z^2)

r=sqrt(x^2 + y^2 + z^2),请问，哪个地方有问题呢

my_aman

get-it 发表于 2017-7-10 23:40
G是一个矢量，r_ij = \sqrt(x^2 + y^2 + z^2)

是的，请问我写的代码哪里出问题了呢

my_aman

zjxitcc 发表于 2017-7-10 22:36
讲真，这梯度公式我反正没看懂怎么推出来的··· 楼下解解惑，我都怀疑这梯度公式错了

r = sqrt(x^2 + y^2 + z^2)，但是就是找不出哪里有问题

my_aman

coolrainbow 发表于 2017-7-11 05:43
-C/r^6    ===>   +6C/r^8 * (rx)

请问*（rx)是什么意思呢，请详细解释下呗，代码应该怎么改呢？

get-it

d:\Pictures\剪贴板01.png

卡开发发

my_aman 发表于 2017-7-11 08:01
请问*（rx)是什么意思呢，请详细解释下呗，代码应该怎么改呢？

估计应该是这样的，应该也就是基本的矢量和微积分计算。

my_aman

get-it 发表于 2017-7-11 09:25

C:\Users\aman\Desktop

代码里6C/... 有，不清楚到底哪里有问题？

my_aman

卡开发发 发表于 2017-7-11 09:29
估计应该是这样的，应该也就是基本的矢量和微积分计算。

公式不带系数，代码里加上了。可是还是有问题。好崩溃呀

get-it

my_aman 发表于 2017-7-11 09:51
代码里6C/... 有，不清楚到底哪里有问题？

讲道理，这种基础算法早就有库了，不用自己写的

my_aman

get-it 发表于 2017-7-11 10:11
讲道理，这种基础算法早就有库了，不用自己写的

做结构全局优化需要这个。心累