计算化学公社

标题: 势函数-LBFGS,程序代码问题 [打印本页]

作者
Author: my_aman 时间: 2017-7-10 21:54
标题: 势函数-LBFGS,程序代码问题
本帖最后由 my_aman 于 2017-7-10 22:03 编辑

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

作者
Author: zjxitcc 时间: 2017-7-10 22:36
讲真，这梯度公式我反正没看懂怎么推出来的···

楼下解解惑，我都怀疑这梯度公式错了

作者
Author: get-it 时间: 2017-7-10 23:40
G是一个矢量，r_ij = \sqrt(x^2 + y^2 + z^2)

作者
Author: zjxitcc 时间: 2017-7-11 01:05

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

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

作者
Author: coolrainbow 时间: 2017-7-11 05:43
-C/r^6 ===> +6C/r^8 * (rx)

作者
Author: my_aman 时间: 2017-7-11 07:56

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),请问，哪个地方有问题呢

作者
Author: my_aman 时间: 2017-7-11 07:57

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

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

作者
Author: my_aman 时间: 2017-7-11 07:58

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

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

作者
Author: my_aman 时间: 2017-7-11 08:01

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

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

作者
Author: get-it 时间: 2017-7-11 09:25
本帖最后由 get-it 于 2017-7-11 09:28 编辑

d:\Pictures\剪贴板01.png

作者
Author: 卡开发发 时间: 2017-7-11 09:29
本帖最后由卡开发发于 2017-7-11 09:31 编辑

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

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

作者
Author: my_aman 时间: 2017-7-11 09:51

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

C:\Users\aman\Desktop

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

作者
Author: my_aman 时间: 2017-7-11 09:53

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

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

作者
Author: get-it 时间: 2017-7-11 10:11

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

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

作者
Author: my_aman 时间: 2017-7-11 10:29

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

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

作者
Author: my_aman 时间: 2017-7-11 10:29
呼叫 Sob老师，请求指点。

欢迎光临计算化学公社 (http://ccc.keinsci.com/)