计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 14419 | 评论数: 15 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2017-7-10 21:54

正文摘要:

本帖最后由 my_aman 于 2017-7-10 22:03 编辑 请问各位老师大神: 我在用势函数+LBFGS,编写势函数势能量梯度代码,经过反复测试未得到预期的结构。 应该是梯度写的有问题,但是我就是没检查出来。现把代码贴出 ...

回复 Reply

my_aman 发表于 Post on 2017-7-11 10:29:44
呼叫 Sob老师,请求指点。
my_aman 发表于 Post on 2017-7-11 10:29:07
get-it 发表于 2017-7-11 10:11
讲道理,这种基础算法早就有库了,不用自己写的

做结构全局优化需要这个。心累
get-it 发表于 Post on 2017-7-11 10:11:07
my_aman 发表于 2017-7-11 09:51
代码里6C/... 有,不清楚到底哪里有问题?

讲道理,这种基础算法早就有库了,不用自己写的
my_aman 发表于 Post on 2017-7-11 09:53:37
卡开发发 发表于 2017-7-11 09:29
估计应该是这样的,应该也就是基本的矢量和微积分计算。

公式不带系数,代码里加上了。可是还是有问题。好崩溃呀
my_aman 发表于 Post on 2017-7-11 09:51:40

C:\Users\aman\Desktop

代码里6C/... 有,不清楚到底哪里有问题?
卡开发发 发表于 Post on 2017-7-11 09:29:39
本帖最后由 卡开发发 于 2017-7-11 09:31 编辑
my_aman 发表于 2017-7-11 08:01
请问*(rx)是什么意思呢,请详细解释下呗,代码应该怎么改呢?

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


QQ截图20170711083529.png (69.83 KB, 下载次数 Times of downloads: 56)

QQ截图20170711083529.png

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
get-it + 2 233333

查看全部评分 View all ratings

get-it 发表于 Post on 2017-7-11 09:25:59
本帖最后由 get-it 于 2017-7-11 09:28 编辑

d:\Pictures\剪贴板01.png

剪贴板01.png (25.02 KB, 下载次数 Times of downloads: 56)

剪贴板01.png

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
卡开发发 + 2 前后脚:-)

查看全部评分 View all ratings

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

请问*(rx)是什么意思呢,请详细解释下呗,代码应该怎么改呢?
my_aman 发表于 Post on 2017-7-11 07:58:02
zjxitcc 发表于 2017-7-10 22:36
讲真,这梯度公式我反正没看懂怎么推出来的··· 楼下解解惑,我都怀疑这梯度公式错了

r = sqrt(x^2 + y^2 + z^2),但是就是找不出哪里有问题
my_aman 发表于 Post on 2017-7-11 07:57:26
get-it 发表于 2017-7-10 23:40
G是一个矢量,r_ij = \sqrt(x^2 + y^2 + z^2)

是的,请问我写的代码哪里出问题了呢
my_aman 发表于 Post on 2017-7-11 07:56:59
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),请问,哪个地方有问题呢
coolrainbow 发表于 Post on 2017-7-11 05:43:10
-C/r^6    ===>   +6C/r^8 * (rx)
zjxitcc 发表于 Post on 2017-7-11 01:05:58
get-it 发表于 2017-7-10 23:40
G是一个矢量,r_ij = \sqrt(x^2 + y^2 + z^2)

我试过矢量,也没推出来长这样。。
get-it 发表于 Post on 2017-7-10 23:40:10
G是一个矢量,r_ij = \sqrt(x^2 + y^2 + z^2)

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 09:43 , Processed in 0.181371 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list