k*********g 发帖数: 791 | |
j****x 发帖数: 943 | |
t*****z 发帖数: 812 | 3 当然不是万能丹。总有个适用范围的
【在 k*********g 的大作中提到】 : 不解释;
|
k*********g 发帖数: 791 | 4 即便在有限元适用的范围、自觉得很牛p的fields,也有一个数值方法比它快10~100倍
左右;
【在 t*****z 的大作中提到】 : 当然不是万能丹。总有个适用范围的
|
t*****z 发帖数: 812 | 5 这个有点意思
能不能举个例子。最好是实际的工业中的例子。简单的例子解析解最快了:)
【在 k*********g 的大作中提到】 : 即便在有限元适用的范围、自觉得很牛p的fields,也有一个数值方法比它快10~100倍 : 左右;
|
k*********g 发帖数: 791 | 6 需要举例?
等到你能看到例子的时候;
世界上就有很多人知道这些了;
真正做过数值研究的、同时也会思考的;
就会体会到;
至少这么一点:
一个产生的矩阵有几百个非零元素的数值方法;
比一个bandwidth只有几个元素的方法;
要慢10到100倍;
(如果其他情况一样的话); |
k****a 发帖数: 4162 | 7 你这说了等于没说嘛
【在 k*********g 的大作中提到】 : 需要举例? : 等到你能看到例子的时候; : 世界上就有很多人知道这些了; : 真正做过数值研究的、同时也会思考的; : 就会体会到; : 至少这么一点: : 一个产生的矩阵有几百个非零元素的数值方法; : 比一个bandwidth只有几个元素的方法; : 要慢10到100倍; : (如果其他情况一样的话);
|
v****r 发帖数: 11 | 8 给你source加个delta function,你弄出个bandwidth解解看 |
k*********g 发帖数: 791 | 9 关于数值方法之间的性能差别,你们这些没有亲手做过数值实验,进行过比较的人,哪
来的那么多自信?
数值计算不是物理学那样well documented,数值计算的好东西像典型工程那样装在实
践者、研究者的头脑里
教材大多数是由2、3流研究者写的,paper多数是垃圾,光看这些是没有用的。 |
O*******d 发帖数: 69 | |
|
|
c*******g 发帖数: 475 | 11 只能说明你傻, 不解释.
【在 k*********g 的大作中提到】 : 即便在有限元适用的范围、自觉得很牛p的fields,也有一个数值方法比它快10~100倍 : 左右;
|
b*******h 发帖数: 2585 | |
k*********g 发帖数: 791 | 13 谱元法是好方法。
采用正交完备函数,比拉格朗日多项式要快10~100倍。 |
k*********g 发帖数: 791 | 14 教材大多数是由2、3流研究者写的,paper多数是垃圾,光看这些是没有用的。 |
k*********g 发帖数: 791 | 15 dirac delta 可以用正交完备函数(即本征函数)expand,所以它对数值计算不造成任
何困难。
这些是最基本的 mathematical methods in physics的内容。
而且,source项往往是Ax=b的右边相,跟矩阵性态有神马关系?
这些东东都不知道,也研发个p数值计算?
现在的孩子们数值计算水平很烂,主要是他们的老师水平太烂。
我算是中国、美国的数值牛人见识了不少的。一开始非常敬畏他们,因为我自己啥也做
不出来,他们至少张牙舞爪看上去挺newbee的。后来博士期间推了1千页数学、写了20
万训练程序后,发觉原来他们啥也不懂,根本没有吃透数值计算的精髓。
【在 v****r 的大作中提到】 : 给你source加个delta function,你弄出个bandwidth解解看
|
o******l 发帖数: 828 | 16 开馆收徒吧。我也是看他们写的不知所云。
20
【在 k*********g 的大作中提到】 : dirac delta 可以用正交完备函数(即本征函数)expand,所以它对数值计算不造成任 : 何困难。 : 这些是最基本的 mathematical methods in physics的内容。 : 而且,source项往往是Ax=b的右边相,跟矩阵性态有神马关系? : 这些东东都不知道,也研发个p数值计算? : 现在的孩子们数值计算水平很烂,主要是他们的老师水平太烂。 : 我算是中国、美国的数值牛人见识了不少的。一开始非常敬畏他们,因为我自己啥也做 : 不出来,他们至少张牙舞爪看上去挺newbee的。后来博士期间推了1千页数学、写了20 : 万训练程序后,发觉原来他们啥也不懂,根本没有吃透数值计算的精髓。
|
i**g 发帖数: 134 | 17 请问 delta 函数用正交完备基展开的时候,
展开项的个数是有限的还是无限的?
20
【在 k*********g 的大作中提到】 : dirac delta 可以用正交完备函数(即本征函数)expand,所以它对数值计算不造成任 : 何困难。 : 这些是最基本的 mathematical methods in physics的内容。 : 而且,source项往往是Ax=b的右边相,跟矩阵性态有神马关系? : 这些东东都不知道,也研发个p数值计算? : 现在的孩子们数值计算水平很烂,主要是他们的老师水平太烂。 : 我算是中国、美国的数值牛人见识了不少的。一开始非常敬畏他们,因为我自己啥也做 : 不出来,他们至少张牙舞爪看上去挺newbee的。后来博士期间推了1千页数学、写了20 : 万训练程序后,发觉原来他们啥也不懂,根本没有吃透数值计算的精髓。
|
k*********g 发帖数: 791 | 18 无限的;
不记得那个0到1的跳跃可以用无限傅立叶级数展开吗?
【在 i**g 的大作中提到】 : 请问 delta 函数用正交完备基展开的时候, : 展开项的个数是有限的还是无限的? : : 20
|
k*********g 发帖数: 791 | 19 green函数可以用正交完备函数展开;
什么是green?
不就是一个奇点吗 |
i**g 发帖数: 134 | 20 需要无限项的话,那就不是对数值算法不造成任何影响了
【在 k*********g 的大作中提到】 : 无限的; : 不记得那个0到1的跳跃可以用无限傅立叶级数展开吗?
|
|
|
k*********g 发帖数: 791 | 21 ???
【在 i**g 的大作中提到】 : 需要无限项的话,那就不是对数值算法不造成任何影响了
|
k*********g 发帖数: 791 | 22 我看你还是赶紧做我的门徒吧,你的老师啥也没有教会你。
【在 i**g 的大作中提到】 : 需要无限项的话,那就不是对数值算法不造成任何影响了
|
k*********g 发帖数: 791 | 23 并行计算
并行计算在科学、工程计算中是不可或缺的。通常人们认为,并行计算的瓶颈在于如何
实现不同处理器之间数据通信的时间远少于在当地处理器上的计算时间。本人的数值实
验和经验表明,这个根本不是并行计算的最大难点,这个其实这个很容易实现。假设我
们解决一个2维问题,N为一个方向上的未知数,把2维区域分成m^2个模块,每一个模块
由一个处理器管理。那么当地的计算量和N/m的平方成正比,而数据交换只跟N/m成正比
。对于典型的并行计算问题来说,N远比m要大,故通常的担忧、观点是偏颇的,没有抓
住问题要害。
数值实验表明,已经解决了人们关心的那个所谓瓶颈问题的现有软件的并行效果没有一
个是理想的。本人的研究表明,并行计算的瓶颈是,并行过程改变了原先的单处理器计
算的迭代结构,这种改变通常是往不利的方向改变的,并且跟处理器之间的数据依赖程
度有关,数据依赖程度越高,迭代的恶化程度越厉害。在有限元中,外来处理器在当地
处理器的渗透程度很高,而在谱元法中,外来处理器只允许影响本地处理器最表面的一
层。本人曾经的实验室同事的数值实验表明,谱元法能达到前所未有的并行效果。 |
k*********g 发帖数: 791 | 24 有限元是不堪一击的,我对她很merciful,而且因为::
搞有限元的对数值计算的推广是有很大贡献的,这一点要表扬。
数值计算communicity在科技、工业中的marketshare太小太小,造成很多人没工作没饭
吃。
搞数值计算应该多到处游说、鼓吹,把这块饼做大。 |
C**********B 发帖数: 24 | 25 个人意见,不是来打架的哈。我觉得每个方法有自己的特点和缺点,有自己的适用范围
和问题。就算是对某一方法的适用问题,还要看用的人是是否是结合自己的项目用好了
,否真的懂了,是否真的会用。而且每个方法的弱点都有很多人在做研究试图去弥补。
楼主仅仅凭借自己的经验和看法,一概而论一个被广泛使用得到了很多人实践认可的方
法,个人觉得是不太合适的,这里多数人的回帖也说明了这一点。而且这还是在楼主经
验丰富的前提下说出来的,个人有点难以理解。所以个人觉得楼主有点发泄不满的嫌疑
,呵呵。
我觉得有限元最大的弱点就是上手慢。算法,分析和写程序相对复杂,当然还有其他的
弱点。付出了代价自然要有额外收获才行。很多人用有限元,自然是有他们的需求。如
果一个问题不需要有限元所拥有的那些优势点,使用有限元恐怕就事倍功半了。其实据
我所知,研究和用有限元的人远少于差分法有限体,也有这些原因在里面。
没有听说什么大事,我一般不来mitbbs,这次是给侄女想读的健康管理专业找点信息,
才来的。突然看到还有个科学计算版,就进来晃一晃,觉得这个题目挺有意思,所以发
表一点小小看法。以后有空常来和大家交流交流。。。。。如果有知道健康管理专业方
面信息的朋友,也请不吝赐教。 |
c***r 发帖数: 1570 | 26 看到ls的,我笑了,
还没看出来lz不过是写过几个课程小project的跳梁小丑,
然后鼓吹自己是 '10万行公式20万行代码'
‘而且这还是在楼主经验丰富的前提下说出来的,‘ -- 哪里看出来经验丰富了,我倒
是觉得经验很少,你让楼主写个最简单的newton solver出来看看?
眼睛擦亮点,别把小丑当大师。
之前还有人说要拜师,这种人还是不要做科学计算了,没潜质,连简单的bool类型true
和false都分辨不了,怎么写if-statement。 |
C**********B 发帖数: 24 | 27 ls的: 不管lz是跳梁小丑还是经验丰富,不管我有没有看出来他是在干什么,也不管想拜他为师的人怎么想,在这个自由发表意见的地方,你要想得到被人的尊重,不想被别人攻击,在表达自己观点的时候就要尊重别人的观点,不要去攻击他人。否则讨论就变成了吵架,而你自己也是得不到尊重的。当然,在论坛上这种事情经常发生,怕你也习以为常了。我对你这种行径的反感是我不大上论坛的主要原因之一。如我前面所说,要不是需要给侄女找专业信息,我也不会偶然看到这个版面,来跟你胡搅蛮缠。发第四个帖子就是跟人针锋相对,感觉很不爽。
你如果对他人观点不屑一顾,可以一笑了之。但如果你要发表言论,就算是反对的言论
,也应该要以尊重他人为基础。lz的确有不少漏洞,我也不敢苟同他这个帖子和其他几个帖子里的很多观点。但他自称经验丰富,我无法也无时间去确认,也不愿去争这个口舌,当时更是没有时间去浏览其他帖子了解别人和lz在其他帖子里的争论。所以为了避免吵架也就假设他的确经验丰富,然后来表达我的不同观点。难道我不攻击人就是把他当大师了?笑话!如果你也知道一点要尊重他人,你有什么好笑的?在你叫别人擦亮眼睛的时候,你有没有意识到你的眼睛已经被自己的自大蒙上了一层猪油?你和lz在其他地方的争论我是刚刚才看到的,恐怕你也有发泄不满的嫌疑。就算lz如你所说是只写过一点程序,有半壶水响叮当的嫌疑,或者是有意想宣传谱方法,你要客观指出来可以,但公开搞藐视嘲笑攻击,甚至想踩着他人显摆自己,也只能得到别人的鄙视。正如你现在从我这里得到的一样!不尊重他人的人永远也得不到别人的尊重。当然,现在我已经在做我不希望做的事情,攻击你不尊重你了,所以没有指望你的尊重。这也是为什么我感觉很不爽的原因。。。。。。其实我不是个谦虚的人,毛病不少,脾气有点急,话也多(从这个帖子怕也看的出来^-^),以前比较狂比较傲比较爱吹牛。但我也知道在论坛上不要随便攻击人,因为我希望我的观点不管是否被大家接受,至少能够被大家公正的对待,而不是跟大家吵架。当然你的观点可能不一样。那我们就只能各持己见了。
你不要把自己粉得那么经验丰富,一眼就能看出别人是什么人。网络上的事情说不清楚,同一个人发不同帖子时候的不同出发点导致的结果可能让你对同一个人做出两个不同的判断,所以我不觉得你有这么大本事。我不知道你到底多有经验,所以我不敢说我在做研究项目写程序做分析上比你有经验,但你也不一定敢这么说。而且就算你经验比这里所有人都丰富,你就有资格攻击人了么?如果你真是经验丰富,是这个行当的老人,你这就叫为老不尊!如果不是,那就叫狂妄自大!
每个人看问题都有不同的出发点。在他现在看来对的事情,在他以后,或者在其他人看来可能就是个小儿科的错误。也许他的观点得不到大家的同意,但也不是活该被大家攻击的。我不知道他在其他地方有没有攻击他人,至少在这里他虽然口气傲慢自大,还在贬低我也搞了十来年的有限元,但他没有像你这么攻击人,所以自然我也不愿攻击他。你觉得你比他厉害,他是跳梁小丑,还有比你更厉害人也觉得你也是个跳梁小丑而已。而且你自己的行为也不过把自己描述成了一个跳梁小丑。将心比心,你自己怕也不希望被别人这么贬低吧。今天翻帖子时候看到你不少回帖,看得出你是这里的常客。虽然现在我在和你针锋相对,但也得承认你比较专业,很多回帖也没有攻击人。但是,你自己在这里这么随意攻击他人,就算你比lz更专业,能找出lz的错误,也只能得到别人的鄙视,无他。 |
i**g 发帖数: 134 | 28 ft
既然delta函数需要用无限项的正交基展开,
你用有限阶数的Chebyshev 或者 Legendre
谱元法,居然敢说对数值计算没有任何影响?
而且第一次见到做谱元的敢诋毁有限元,
谱元只不过是使用了不同的基函数,实现
了质量矩阵对角化,但是所基于的数学基础
完全是同有限元来的,披了个马甲就真当
是本质区别了?
【在 k*********g 的大作中提到】 : 我看你还是赶紧做我的门徒吧,你的老师啥也没有教会你。
|
k*********g 发帖数: 791 | 29 my gosh,你这话说得。
后面的项越来越小,即便把后面的无穷多项项加起来,也很小,也可以ignored。可以
这样的。
delta函数用当然可以用各种 正交函数 表示,这个随便抓一本 mathematical methods
in physics 就说了,
the following is excellent text, most popular one among american graduate
programs in physics & in engineering:
mathematical methods for physicists
sixth edition, by Arfken & Weber
【在 i**g 的大作中提到】 : ft : 既然delta函数需要用无限项的正交基展开, : 你用有限阶数的Chebyshev 或者 Legendre : 谱元法,居然敢说对数值计算没有任何影响? : 而且第一次见到做谱元的敢诋毁有限元, : 谱元只不过是使用了不同的基函数,实现 : 了质量矩阵对角化,但是所基于的数学基础 : 完全是同有限元来的,披了个马甲就真当 : 是本质区别了?
|
s******g 发帖数: 38 | |
|
|
i**g 发帖数: 134 | 31 请不要随意假设别人的观点,我重来没说过delta函数不可以用完备正交函数展开。
我只是对你说的delta函数对数值计算,尤其你推崇的SEM,不造成任何影响,这个结论
不同意。
delta在各个频谱上的幅值是相同的,你用什么正交函数展开可以实现后面的项越来越
小,以至后面无穷多项加起来,也可以忽略?
methods
【在 k*********g 的大作中提到】 : my gosh,你这话说得。 : 后面的项越来越小,即便把后面的无穷多项项加起来,也很小,也可以ignored。可以 : 这样的。 : delta函数用当然可以用各种 正交函数 表示,这个随便抓一本 mathematical methods : in physics 就说了, : the following is excellent text, most popular one among american graduate : programs in physics & in engineering: : mathematical methods for physicists : sixth edition, by Arfken & Weber
|
t*****n 发帖数: 4908 | 32 坑兄20万程序化了多久时间写完的?有没有成型产品?
20
【在 k*********g 的大作中提到】 : dirac delta 可以用正交完备函数(即本征函数)expand,所以它对数值计算不造成任 : 何困难。 : 这些是最基本的 mathematical methods in physics的内容。 : 而且,source项往往是Ax=b的右边相,跟矩阵性态有神马关系? : 这些东东都不知道,也研发个p数值计算? : 现在的孩子们数值计算水平很烂,主要是他们的老师水平太烂。 : 我算是中国、美国的数值牛人见识了不少的。一开始非常敬畏他们,因为我自己啥也做 : 不出来,他们至少张牙舞爪看上去挺newbee的。后来博士期间推了1千页数学、写了20 : 万训练程序后,发觉原来他们啥也不懂,根本没有吃透数值计算的精髓。
|
k*********g 发帖数: 791 | 33 20万程序 was due to research, most are eventually erased.
【在 t*****n 的大作中提到】 : 坑兄20万程序化了多久时间写完的?有没有成型产品? : : 20
|
k*********g 发帖数: 791 | 34 the MAIN problem with spectral polynomials handling delta function is:
Gibbs phenomena
by the way, the Gibbs paradox kills the possibility of concurrent
statistical mechanics to be the most fundamental theory for the universe.
【在 i**g 的大作中提到】 : 请不要随意假设别人的观点,我重来没说过delta函数不可以用完备正交函数展开。 : 我只是对你说的delta函数对数值计算,尤其你推崇的SEM,不造成任何影响,这个结论 : 不同意。 : delta在各个频谱上的幅值是相同的,你用什么正交函数展开可以实现后面的项越来越 : 小,以至后面无穷多项加起来,也可以忽略? : : methods
|
t*****n 发帖数: 4908 | 35 坑兄是应用数学的吗?如果是research代码,能否公开?我们也可以拜读拜读。
【在 k*********g 的大作中提到】 : 20万程序 was due to research, most are eventually erased.
|
t*****n 发帖数: 4908 | 36 放狗搜之原来坑兄是搞物理的。佩服佩服。怪不得理论一套一套的。让我们学工的好崇
拜。
【在 t*****n 的大作中提到】 : 坑兄是应用数学的吗?如果是research代码,能否公开?我们也可以拜读拜读。
|
x*******1 发帖数: 28835 | |
b****e 发帖数: 906 | 38 Gibbs paradox 不是考虑全同粒子就解决了吗?
【在 k*********g 的大作中提到】 : the MAIN problem with spectral polynomials handling delta function is: : Gibbs phenomena : by the way, the Gibbs paradox kills the possibility of concurrent : statistical mechanics to be the most fundamental theory for the universe.
|