v*****s
发帖数: 20290 | 1 你这个也是胡扯。啥叫最优，啥叫赢棋？狗肯定是每步都按胜率最高的来。你可以猜测
，在局面落后时，狗算出来的各个步数的胜率都差不多，所以乱下（虽然我觉得这个猜
测不太可能成立，再差不多也终会有细微的差别，我不相信狗会把胜率作roundoff）。
但这种情况下，你打算怎么定义最优？盘面目数最大，输目最小？这个改动可不小，他
得把每种情况下终盘的目数都存下来，然后做加权平均。 |
|
s******y
发帖数: 68 | 2 对呀, roundoff error 大了点. 呵呵 |
|
d****n
发帖数: 1637 | 3 /*working with larger numbers but still die */
/**C(80, 20) is okay, which was not ok in previous version**/
#include
#include
#include
float
gammln (float xx)
//Returns the value ln[Ãxx)] for xx > 0.
{
//Internal arithmetic will be done in double precision, a nicety that
//you can omit if five-figure
//accuracy is good enough.
double x, y, tmp, ser;
static double cof[6] = { 76.18009172947146, -86.50532032941677,
24.01409824083091, -1.231739572450155,
... 阅读全帖 |
|
d****n
发帖数: 1637 | 4 /*working with larger numbers but still die */
/**C(80, 20) is okay, which was not ok in previous version**/
#include
#include
#include
float
gammln (float xx)
//Returns the value ln[Ãxx)] for xx > 0.
{
//Internal arithmetic will be done in double precision, a nicety that
//you can omit if five-figure
//accuracy is good enough.
double x, y, tmp, ser;
static double cof[6] = { 76.18009172947146, -86.50532032941677,
24.01409824083091, -1.231739572450155,
... 阅读全帖 |
|
y**b
发帖数: 10166 | 5 我曾经用intel的编译器(好像也就它支持吧)来利用这个80bit，可是没看到什么效果。
好像是放内存里面的时候才80bit，取出来就64bit了，怎么回事有点记不清了。
但是用quad-precision就不同了，效果非常明显，控制并行计算的roundoff误差很有效。 |
|
c****n
发帖数: 2031 | 6 那就找某个reliable的算法，算出个很精确的解（比如取很密的mesh，保证error在
roundoff），让后用它作为精确解。 |
|
b****e
发帖数: 906 | 7 来自主题: Computation版 - 并行计算 考虑一下实际计算中的roundoff error,还有，完全的正定矩阵它的condition number
也可以很小，cg的收敛速度为1/condition number |
|
l*********2
发帖数: 95 | 8
很高兴终于有人回复了!
你说得有道理,我也想去问问两个系的弹性力学究竟有多大区别,我下周再去问一下我的
有限元老师吧.
我查了一下,下学期只有civil的弹性力学课程，航空的不开设。我想弹性力学应该算是
比较基础的固体力学课程吧，应该尽早修读吧？
至于数学课，今日一位学长简单向我推荐MA 51400 - Numerical Analysis,说对有限元
什么有帮助.你觉得呢?
我刚看了看Numerical Analysis和MA527的课程介绍，感觉MA527讲的topic以前都学过
，但是或许会比较深吧。而且下学期不开527，只有MA528.
我打算下周把这些课程资料整理好，找我有限元老师问一下吧。当然还是希望能听听板
上大侠们的意见！
Numerical Analysis
Iterative methods for solving nonlinear, linear difference equations,
applications to solution of polynomial equations; differentiation and
integration form... 阅读全帖 |
|
l*********2
发帖数: 95 | 9
很高兴终于有人回复了!
你说得有道理,我也想去问问两个系的弹性力学究竟有多大区别,我下周再去问一下我的
有限元老师吧.
我查了一下,下学期只有civil的弹性力学课程，航空的不开设。我想弹性力学应该算是
比较基础的固体力学课程吧，应该尽早修读吧？
至于数学课，今日一位学长简单向我推荐MA 51400 - Numerical Analysis,说对有限元
什么有帮助.你觉得呢?
我刚看了看Numerical Analysis和MA527的课程介绍，感觉MA527讲的topic以前都学过
，但是或许会比较深吧。而且下学期不开527，只有MA528.
我打算下周把这些课程资料整理好，找我有限元老师问一下吧。当然还是希望能听听板
上大侠们的意见！
Numerical Analysis
Iterative methods for solving nonlinear, linear difference equations,
applications to solution of polynomial equations; differentiation and
integration form... 阅读全帖 |
|
S*********g
发帖数: 5298 | 10 >> help rem
REM Remainder after division.
REM(x,y) is x - n.*y where n = fix(x./y) if y ~= 0. If y is not an
integer and the quotient x./y is within roundoff error of an integer,
then n is that integer. The inputs x and y must be real arrays of the
same size, or real scalars.
By convention:
REM(x,0) is NaN.
REM(x,x), for x~=0, is 0.
REM(x,y), for x~=y and y~=0, has the same sign as x.
Note: MOD(x,y), for x~=y and y~=0, has the same sign as y.
|
|
A*******s
发帖数: 3942 | 11 大侠这个macro我用过，非常不错。不过有个小feedback，就是里面的对score和weight
的roundoff会带来一些误差，我觉得也没有必要。 |
|
