| 查看: 901 | 回复: 6 | |||
[求助]
求助：tinker程序包中，某段程序看不太懂，想知道其数学思想！
|
|
C compute the length of gradient and search direction g_norm = 0 s_norm = 0 do i = 1, nvar g_norm = g_norm+g(i) * g(i) ! here,g(i) is gradient s_norm = s_norm+s(i) * s(i) ! here,s(i) =-g(i) end do g_norm = sqrt( g_norm ) s_norm = sqrt( s_norm ) 这里，想请教高手，这里的数学思想是什么？ C store initial function, then normalize the search vector and find directional gradient sg_0 = 0 do i=1,nvar s(i) = s(i) / s_norm sg_0 = sg_0 + s(i) * g(i) end do C set the initail stepsize to the length of the passed search vector, or based on provious function decrease stpmax = 5.0 f_move = 0.5 * stpmax * g_norm step = 2.0 * abs ( f_move / sg_0) step = min ( step , s_norm ) 这两段，同样不太理解里面的数学思想，同时，不知道这是根据什么公式或者原理确定的步长 ！请求高手的帮助！ [ 来自小组 减肥吧 ] |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有21人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2013-12-03 22:43:11
3楼2013-12-04 08:38:02
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2013-12-04 08:55:14
5楼2013-12-04 09:23:13
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
6楼2013-12-04 09:57:57
7楼2013-12-04 10:23:31