由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Programming版 - 数值计算
相关主题
C++ 做线性代数,方便使用的库?大家编程用得algebra lib是什么? (转载)
数组,结构,类在数值计算中哪个快请问:稀疏矩阵运算有什么好的C++库?
卫东大神来说说阿尔法狗横扫棋坛这事吧Matrix calculation in C++
Hi Thrust, I am a little worried about my privacyLAPACK/ARPACK under windows
我觉得王垠对Agile的看法很有道理啊Switch from Matlab to C(C++)?
各位轮子派?你们用到后面不觉得蛋疼?do you use blas/lapack?
有多少人觉得这一波AI是忽悠?Linux下运行lapack和blas的问题
求救:2个dense matrix的乘法How to use multi-core to speed Python program
相关话题的讨论汇总
话题: br话题: 计算话题: 数值话题: pde话题: 轨道
进入Programming版参与讨论
1 (共1页)
g****t
发帖数: 31659
1
数值计算开源选择其实还是非常少的。
一来要代码质量高。二来语法尽量描述性质的,
靠近Matlab,sql风格。
Numpy质量很高,这是python重要加分。
如果要快速开发。(我们假设你要研究如何用
ann算天气预报。)
那么Python/R/lua恐怕prototype
阶段还是仅有的选择了。
这三种不会有大的不可预测的问题。不会出现你开发了半年推倒重来的情况,或者憋死
一个月的
情况。
Rust,go,jvm上的各种线性代数,都没有大公司背书。
我觉得还是挺短视的。goog如果出钱弄好这块,其实我觉得
Go有潜力占住ML市场的。
有design pattern功力的多半不会自费武功去go.
做数值计算的多数是看不上design pattern等东西的。
思路很多其实还是fortran。
h**********c
发帖数: 4120
2
LAPACK,LINPACK, TNT这些美国高校大研究所已经用了多少年,
比较偏门DS tool, AUTO这些也都是开源,
你说那个ANN怎么做stability analysis,保证计算收敛?
可能就是用随机数来收敛吧。
估计你也就是在高校里,看看这各大公司热热闹闹,里面各种戏剧帮派。
w***g
发帖数: 5958
3
我以前要算一个啥规划, 只有NAG里有. 那时的数值算法真是数值算法.
现在的ANN, 就是乱棍打死老拳师.

【在 h**********c 的大作中提到】
: LAPACK,LINPACK, TNT这些美国高校大研究所已经用了多少年,
: 比较偏门DS tool, AUTO这些也都是开源,
: 你说那个ANN怎么做stability analysis,保证计算收敛?
: 可能就是用随机数来收敛吧。
: 估计你也就是在高校里,看看这各大公司热热闹闹,里面各种戏剧帮派。

h**********c
发帖数: 4120
4
现在没时间研究ANN,不是说NN这个东西完全没有用,没有严格的论证和可信的模型。
我老觉能一个编程语言静态分析系统,来避免常识性错误和reinvent wheels是有点用
处。不过觉得,你说那是knoledge base不知道区别在哪里,也许NN在路径规划和搜索
上有优势,没研究过。

【在 w***g 的大作中提到】
: 我以前要算一个啥规划, 只有NAG里有. 那时的数值算法真是数值算法.
: 现在的ANN, 就是乱棍打死老拳师.

g****t
发帖数: 31659
5
估计本版只有我一个知道你的id是Poincare书里定义的了,Lol
数值计算的基础部分,也是最早发展的部分是
Ode, pde求解。早先都是实验物理。
别的不和你多讲了,反正你也听不懂。


: LAPACK,LINPACK, TNT这些美国高校大研究所已经用了多少年,

: 比较偏门DS tool, AUTO这些也都是开源,

: 你说那个ANN怎么做stability analysis,保证计算收敛?

: 可能就是用随机数来收敛吧。

: 估计你也就是在高校里,看看这各大公司热热闹闹,里面各种戏剧帮派。



【在 h**********c 的大作中提到】
: 现在没时间研究ANN,不是说NN这个东西完全没有用,没有严格的论证和可信的模型。
: 我老觉能一个编程语言静态分析系统,来避免常识性错误和reinvent wheels是有点用
: 处。不过觉得,你说那是knoledge base不知道区别在哪里,也许NN在路径规划和搜索
: 上有优势,没研究过。

g****t
发帖数: 31659
6
现在是收获期。但是常见的ann技术已经出现降价的趋势了。
另外像我的工作,一个算法会装机上亿的量。
矩阵计算是不能出错的。massive industry products的软件以及算法的要求和别的行
业不同。
g****t
发帖数: 31659
7
大公司的优势在于scale out出去的时候成本低,以势压人。
但是难度就在于细节正确非常难。可以说extremely difficult。
产品走的量大了,什么毛病都会被市场打出来。质量不够,
就会影响share.例如手机前后摄像头的算法软件。还有audio写个小段子容易,量大了
就极难
d***a
发帖数: 13752
8
高性能数值计算暂时还是要靠FORTRAN和C,Java也凑合,再加上GPU编程,用MPI/PVM搞
并行。憋一个月不算什么,有些程序可以运行几个月才出结果,那么就值得多花点时间
,造跑得快一点的轮子。

【在 g****t 的大作中提到】
: 数值计算开源选择其实还是非常少的。
: 一来要代码质量高。二来语法尽量描述性质的,
: 靠近Matlab,sql风格。
: Numpy质量很高,这是python重要加分。
: 如果要快速开发。(我们假设你要研究如何用
: ann算天气预报。)
: 那么Python/R/lua恐怕prototype
: 阶段还是仅有的选择了。
: 这三种不会有大的不可预测的问题。不会出现你开发了半年推倒重来的情况,或者憋死
: 一个月的

z*********e
发帖数: 10149
9
数值计算最基础也最重要的是稳定性和收敛性
不管什么nn,无法证明稳定性的话都白搭

:估计本版只有我一个知道你的id是Poincare书里定义的了,Lol
h**********c
发帖数: 4120
10
我老最近几年比较相信事实比面子重要,
所以如果真的不懂被人揭短也不是太在意。有人肯给你辅导修改论文是你运气。
其实数学史上很重要一段在意大利,两个数学家PK就是解polynomial,解到3次5次在文
艺复兴时代(或更早)就是很难的了,而输家面临出局从此丢饭碗的下场。大家都很
serious,和陈真打擂台一样。所以为什么欧洲出身的数学家把galois的位置放的很高。
我也不是数学家,也不想和你打雷台。不过觉得你觉得你最近说话怪怪的,我的解释是
各为其主。
说东卫讲的乱棍打死老拳师很形像,不过老拳师混了那么多年没些弟子裙带,被打死活
该,不过接下来可能就是群欧。叔觉得押对了上一注,所以要求赌场cash out套现,享
受一下物质文明精神文明双丰收.
东卫所言的乱棍,没有明说,我觉的那个LAXIGILE,disruption mode,加上这个NN。你
可以说这是科学,证伪或证明都是远期。
实际上我们看到的是公交车上群欧司机,飞机上砸驾驶舱,所谓的学术技术会议霸住麦
克风不让别人说话。
这乱棍就是哲学上的虚无主义,没有真理,规律,谁霸住麦克风谁就随心所欲,甚至胡
说八道,指鹿为马。没有论证,设计,验收,没有标准,没有条例规范(这里是指工程
技术上的,政治上的已经太多)。但是事实还是人类长期总结的物理数学工程上的认识
,还是适用于人类对自然的再认识。

【在 g****t 的大作中提到】
: 估计本版只有我一个知道你的id是Poincare书里定义的了,Lol
: 数值计算的基础部分,也是最早发展的部分是
: Ode, pde求解。早先都是实验物理。
: 别的不和你多讲了,反正你也听不懂。
:
:
: LAPACK,LINPACK, TNT这些美国高校大研究所已经用了多少年,
:
: 比较偏门DS tool, AUTO这些也都是开源,
:
: 你说那个ANN怎么做stability analysis,保证计算收敛?
:
: 可能就是用随机数来收敛吧。
:
: 估计你也就是在高校里,看看这各大公司热热闹闹,里面各种戏剧帮派。

相关主题
各位轮子派?你们用到后面不觉得蛋疼?大家编程用得algebra lib是什么? (转载)
有多少人觉得这一波AI是忽悠?请问:稀疏矩阵运算有什么好的C++库?
求救:2个dense matrix的乘法Matrix calculation in C++
进入Programming版参与讨论
g****t
发帖数: 31659
11
我有PDE数值计算book chapter,论文,专利,NSF proposal。
毕业论文有几百页的Mathematica代码。
这些都是我十多年前干的事情。
当时方向比较偏。研究方向是数值计算的代数问题。
类似于微分代数。希望能对PDE找到高斯消去那样的计算scheme。
哦,本站那时候我还有不少与此有关的帖子。
你有这方面的高见,我仍然愿意听听.不然就别套近乎了。

【在 h**********c 的大作中提到】
: 我老最近几年比较相信事实比面子重要,
: 所以如果真的不懂被人揭短也不是太在意。有人肯给你辅导修改论文是你运气。
: 其实数学史上很重要一段在意大利,两个数学家PK就是解polynomial,解到3次5次在文
: 艺复兴时代(或更早)就是很难的了,而输家面临出局从此丢饭碗的下场。大家都很
: serious,和陈真打擂台一样。所以为什么欧洲出身的数学家把galois的位置放的很高。
: 我也不是数学家,也不想和你打雷台。不过觉得你觉得你最近说话怪怪的,我的解释是
: 各为其主。
: 说东卫讲的乱棍打死老拳师很形像,不过老拳师混了那么多年没些弟子裙带,被打死活
: 该,不过接下来可能就是群欧。叔觉得押对了上一注,所以要求赌场cash out套现,享
: 受一下物质文明精神文明双丰收.

g****t
发帖数: 31659
12
难点在coding。
不在于应用数学的那些皮毛.

【在 z*********e 的大作中提到】
: 数值计算最基础也最重要的是稳定性和收敛性
: 不管什么nn,无法证明稳定性的话都白搭
:
: :估计本版只有我一个知道你的id是Poincare书里定义的了,Lol
: :

m******r
发帖数: 1033
13
guest显然是搞数值计算的,曾经(也许现在)靠这个吃饭的。
heteroclinic似是数字信号处理的?
wdong是sde转人工智能外加培训和企业家的。
我呢,是首席民间科学家,没事喜欢瞎捉摸,但不靠这个吃饭, 所以旁观者门清。
n******g
发帖数: 2201
14
魏老师呢?早期国内中关村Delphi 程序员?
[在 magliner (magliner) 的大作中提到:]
:guest显然是搞数值计算的,曾经(也许现在)靠这个吃饭的。
:heteroclinic似是数字信号处理的?
:wdong是sde转人工智能外加培训和企业家的。
:我呢,是首席民间科学家,没事喜欢瞎捉摸,但不靠这个吃饭, 所以旁观者门清。
g****t
发帖数: 31659
15
我真的是neural network的master,然后F2来美国.
找不到奖学金才转pde/ode的。
前面不是有人问neural network的稳定性吗。
我多年前写过一个sequence2sequence的RNN estimator,用qualify/disqualify
好坏点的办法确保稳定性。但我只能做一维。
等5年后我告诉大家诀窍在哪。

【在 m******r 的大作中提到】
: guest显然是搞数值计算的,曾经(也许现在)靠这个吃饭的。
: heteroclinic似是数字信号处理的?
: wdong是sde转人工智能外加培训和企业家的。
: 我呢,是首席民间科学家,没事喜欢瞎捉摸,但不靠这个吃饭, 所以旁观者门清。

m******r
发帖数: 1033
16
teacherwei, i have no comment. 水平也许很高,赶走了goodbug, 我们看客也没得看
了。 好比看小说,两边正在捉对厮杀,大家都希望好人被打的鼻青脸肿,先写丢了性
命,可正义战胜了邪恶,最终还是赢了,至少和坏人同归于尽。 好比heteroclinic说
的<陈真>,最后一把刀插自己身上,和后面日本鬼子穿了葫芦串。
可惜到了最后好人却偏偏输了,还退出江湖,坏人还在这儿蹦跶呢。
连点悲剧美都没有, 令人不胜唏嘘啊。

【在 n******g 的大作中提到】
: 魏老师呢?早期国内中关村Delphi 程序员?
: [在 magliner (magliner) 的大作中提到:]
: :guest显然是搞数值计算的,曾经(也许现在)靠这个吃饭的。
: :heteroclinic似是数字信号处理的?
: :wdong是sde转人工智能外加培训和企业家的。
: :我呢,是首席民间科学家,没事喜欢瞎捉摸,但不靠这个吃饭, 所以旁观者门清。

h**********c
发帖数: 4120
17
感觉您是我聊过最大牛的ID,
我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士我见过这人老
中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛。老板说过一
个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不过老板给过我
一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉得我不行,我
是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南大数学系的或
着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多嘴,水平可定
是有的,可能也在MMT混,有两三面之源,没有深交。
说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密,我还在雀斑
出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。
单说科学计算,我干过什么吧,解ODE,我做过两个模型,一个是CR3BP,通俗说是地球
月亮卫星三体,光这个模型的MATLAB code画图就的上千行MATLAB,python 加上C++/
OPENGL的code 就几万,行数不重要,重要的是理解什么东西干了什么。我觉得除了卖
力跑老板的工具以外,我认为我个人提出了通过能量守恒限定边界条件从使异构环轨道
的计算不是简单暴力/穷举。还做过一个lorenz manifold,帮助老板整理异构环轨道的
计算边界条件。两个模型都用不同的工具进行了大量的可视话。比如说过,康乃尔哪个
教授对lorenz模型的strange attractor的异构环轨道很感兴趣,不知道老板给他看了
没有。总之,这条路我没走通,混了个硕士,靠谢谢java, sql, node js混个饭吃。
从数学上讲,主要是理解gauss collocation的收敛特性,比起来,runge kuta,
implicit euler,一个不稳定,一个精度太低。用第三方工具计算floquet multiplier,
确定平衡点轨道的特性,卫星三体那个有三对,所以我问你明不明白还是个问题?过内
这个圈里,没见过什么paper.美过华人好象有一些名字。比较独门是Keller pseudo
arc length continuation,您要能说说这个东西干什么的,以后您在任何场合白话,i
don't give any sh*t. 条件是我不查资料评记忆,你也不能查资料。
叔在国内时候老板之一,酒店总经理兰方人讲,没有我不做的生意。和您套近乎奏是为
了生意,国内热钱很多很多。别人的博士论文我也看过百十来篇,SOSO。那真是一身功
夫,一樽还泪江月。
客观讲数的同事,力争吸引眼球和课读性,趣味性。

【在 g****t 的大作中提到】
: 我有PDE数值计算book chapter,论文,专利,NSF proposal。
: 毕业论文有几百页的Mathematica代码。
: 这些都是我十多年前干的事情。
: 当时方向比较偏。研究方向是数值计算的代数问题。
: 类似于微分代数。希望能对PDE找到高斯消去那样的计算scheme。
: 哦,本站那时候我还有不少与此有关的帖子。
: 你有这方面的高见,我仍然愿意听听.不然就别套近乎了。

g****t
发帖数: 31659
18
你看过证明洛伦兹系统是混沌的那个计算机程序证明吗?
拿过欧洲青年数学奖的。可能和你方向差不多。可能就是你说的康奈尔的教授。
你的三体程序在slow manifold怎么处理的?正常办法会一步只走一小点或者不走了。


: 感觉您是我聊过最大牛的ID,

: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士
我见过
这人老

: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛
。老板
说过一

: 个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不
过老板
给过我

: 一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉
得我不
行,我

: 是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南
大数学
系的或

: 着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多
嘴,水
平可定

: 是有的,可能也在MMT混,有两三面之源,没有深交。

: 说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密
,我还
在雀斑

: 出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。



【在 h**********c 的大作中提到】
: 感觉您是我聊过最大牛的ID,
: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士我见过这人老
: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛。老板说过一
: 个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不过老板给过我
: 一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉得我不行,我
: 是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南大数学系的或
: 着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多嘴,水平可定
: 是有的,可能也在MMT混,有两三面之源,没有深交。
: 说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密,我还在雀斑
: 出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。

w***g
发帖数: 5958
19
这个主页逼格很高!咱应该把贴主页的风气推起来。

【在 h**********c 的大作中提到】
: 感觉您是我聊过最大牛的ID,
: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士我见过这人老
: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛。老板说过一
: 个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不过老板给过我
: 一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉得我不行,我
: 是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南大数学系的或
: 着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多嘴,水平可定
: 是有的,可能也在MMT混,有两三面之源,没有深交。
: 说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密,我还在雀斑
: 出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。

l*********s
发帖数: 5409
20
同感,不过排版老式了点,不够好看。

【在 w***g 的大作中提到】
: 这个主页逼格很高!咱应该把贴主页的风气推起来。
相关主题
LAPACK/ARPACK under windowsLinux下运行lapack和blas的问题
Switch from Matlab to C(C++)?How to use multi-core to speed Python program
do you use blas/lapack?what kind of language this is and how to change it? (转载)
进入Programming版参与讨论
l*******m
发帖数: 1096
21
其实大牛可以搞搞Stochastic gradient descent on Riemannian manifolds,用到
deep learning上
还有hinton的capsules是不是可以和黎曼几何挂一挂?

【在 h**********c 的大作中提到】
: 感觉您是我聊过最大牛的ID,
: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士我见过这人老
: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛。老板说过一
: 个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不过老板给过我
: 一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉得我不行,我
: 是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南大数学系的或
: 着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多嘴,水平可定
: 是有的,可能也在MMT混,有两三面之源,没有深交。
: 说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密,我还在雀斑
: 出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。

g****t
发帖数: 31659
22
丘成桐已经进来做DL谈话了。前几天有个人转给我的。
你找找看。
有人说,最怕的就是外行大牛进来抢funding。
人家什么都不懂,照样糊的你一愣一愣。
还一个例子,还一个大牛除了让自己老婆开了个drive.ai偷资料,在百度屁都没干,一
年1200万美元。


: 其实大牛可以搞搞Stochastic gradient descent on Riemannian
manifolds,
用到

: deep learning上

: 还有hinton的capsules是不是可以和黎曼几何挂一挂?



【在 l*******m 的大作中提到】
: 其实大牛可以搞搞Stochastic gradient descent on Riemannian manifolds,用到
: deep learning上
: 还有hinton的capsules是不是可以和黎曼几何挂一挂?

m*****n
发帖数: 3575
23
你对PDE数值计算这么懂
请问抛物线型和椭圆型最方便的数值解法是什么?

【在 g****t 的大作中提到】
: 丘成桐已经进来做DL谈话了。前几天有个人转给我的。
: 你找找看。
: 有人说,最怕的就是外行大牛进来抢funding。
: 人家什么都不懂,照样糊的你一愣一愣。
: 还一个例子,还一个大牛除了让自己老婆开了个drive.ai偷资料,在百度屁都没干,一
: 年1200万美元。
:
:
: 其实大牛可以搞搞Stochastic gradient descent on Riemannian
: manifolds,
: 用到

h**********c
发帖数: 4120
24
我不是太理解这个slow manifold的概念,就说我论文里的计算,我导师自己给个概念
叫orbit continuation,在平衡点附近grow 一个周期轨道,不断算这个轨道的
continuation,因为连接两个不同平衡点的周期轨道的轨道也是异构环。我算的都是2
manifold.我又查了一下,你说的slow manifold不是homotopy,可能开会的时候听一个
人问过collision,是不是这个意思?
我这人挨溜号,搞计算的时候图形玩的多,当了马公,还老想着计算的东西。
感觉还是有点多动症,还好没有发展成自闭或弱智。

【在 g****t 的大作中提到】
: 你看过证明洛伦兹系统是混沌的那个计算机程序证明吗?
: 拿过欧洲青年数学奖的。可能和你方向差不多。可能就是你说的康奈尔的教授。
: 你的三体程序在slow manifold怎么处理的?正常办法会一步只走一小点或者不走了。
:
:
: 感觉您是我聊过最大牛的ID,
:
: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士
: 我见过
: 这人老
:
: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛
: 。老板

h**********c
发帖数: 4120
25
我找过一阵子手机pad friendly的framework,后来太忙,不了了之

【在 w***g 的大作中提到】
: 这个主页逼格很高!咱应该把贴主页的风气推起来。
h**********c
发帖数: 4120
26
随机的东西,好象正好是homotopy/continuation这些人相反的一面,不太清楚。

【在 l*******m 的大作中提到】
: 其实大牛可以搞搞Stochastic gradient descent on Riemannian manifolds,用到
: deep learning上
: 还有hinton的capsules是不是可以和黎曼几何挂一挂?

g****t
发帖数: 31659
27
我不跟人social


: 你对PDE数值计算这么懂

: 请问抛物线型和椭圆型最方便的数值解法是什么?



【在 m*****n 的大作中提到】
: 你对PDE数值计算这么懂
: 请问抛物线型和椭圆型最方便的数值解法是什么?

g****t
发帖数: 31659
28
我的意思很简单,就是空间里有一块地方,simulation很多步,然而结果走的极慢。但
是它又不是invariant manifold,所以会有数值困难。


: 我不是太理解这个slow manifold的概念,就说我论文里的计算,我导师自己给
个概念

: 叫orbit continuation,在平衡点附近grow 一个周期轨道,不断算这个轨道的

: continuation,因为连接两个不同平衡点的周期轨道的轨道也是异构环。我算的
都是2

: manifold.我又查了一下,你说的slow manifold不是homotopy,可能开会的时候
听一个

: 人问过collision,是不是这个意思?

: 我这人挨溜号,搞计算的时候图形玩的多,当了马公,还老想着计算的东西。

: 感觉还是有点多动症,还好没有发展成自闭或弱智。



【在 h**********c 的大作中提到】
: 随机的东西,好象正好是homotopy/continuation这些人相反的一面,不太清楚。
h**********c
发帖数: 4120
29
你有具体的例题吗?
我个人的理解,你在一个有洞的空间里,使用依赖于连续特性的数学法则。好象是一个
悖论。也许这就是为什么我看到的实际的计算模型维持在低纬度以保证有限参数的扰动
保证流形的不变和连续。

【在 g****t 的大作中提到】
: 我的意思很简单,就是空间里有一块地方,simulation很多步,然而结果走的极慢。但
: 是它又不是invariant manifold,所以会有数值困难。
:
:
: 我不是太理解这个slow manifold的概念,就说我论文里的计算,我导师自己给
: 个概念
:
: 叫orbit continuation,在平衡点附近grow 一个周期轨道,不断算这个轨道的
:
: continuation,因为连接两个不同平衡点的周期轨道的轨道也是异构环。我算的
: 都是2
:
: manifold.我又查了一下,你说的slow manifold不是homotopy,可能开会的时候
: 听一个

g****t
发帖数: 31659
30
最简单的就是一个信号发生器。
输出是3个频率a,b,c的指数函数叠加一些噪声。如果一个a比b,c两个small 1000倍。
如何根据输出估计a,b,c
你可以试验下
y = A*exp(-iat)+B*exp(-ibt)+C*exp(-ict) + noise
Given y, find best estimation of a

【在 h**********c 的大作中提到】
: 你有具体的例题吗?
: 我个人的理解,你在一个有洞的空间里,使用依赖于连续特性的数学法则。好象是一个
: 悖论。也许这就是为什么我看到的实际的计算模型维持在低纬度以保证有限参数的扰动
: 保证流形的不变和连续。

相关主题
问一个C++下的Bug(Linux下)数组,结构,类在数值计算中哪个快
C++里面 class的数据成员的顺序是什么样的?卫东大神来说说阿尔法狗横扫棋坛这事吧
C++ 做线性代数,方便使用的库?Hi Thrust, I am a little worried about my privacy
进入Programming版参与讨论
h**********c
发帖数: 4120
31
我老当年学传热学时候公式比这复杂多,不过老师也不明白,或着她认为我们不需要明
白,知道算什么查那本书,不过传热的计算多数可以在工程许用范围内验证,不是蒙人
的随机数。工程上科一半要做实验,验证理论,每个人都能重复。
你不能建立你的模型象auto 手册,搞成invariant Ode模型,search auto07p.pdf 或
auto2000.pdf
比如象bruta problem,Sandstede 的模型.
我说的仅供参考。

【在 g****t 的大作中提到】
: 最简单的就是一个信号发生器。
: 输出是3个频率a,b,c的指数函数叠加一些噪声。如果一个a比b,c两个small 1000倍。
: 如何根据输出估计a,b,c
: 你可以试验下
: y = A*exp(-iat)+B*exp(-ibt)+C*exp(-ict) + noise
: Given y, find best estimation of a

N*****r
发帖数: 94
32

我对这类的数值稳定性的问题也很感兴趣
原来比较关心的是天体的长时间轨道模拟的数值稳定性,

【在 g****t 的大作中提到】
: 最简单的就是一个信号发生器。
: 输出是3个频率a,b,c的指数函数叠加一些噪声。如果一个a比b,c两个small 1000倍。
: 如何根据输出估计a,b,c
: 你可以试验下
: y = A*exp(-iat)+B*exp(-ibt)+C*exp(-ict) + noise
: Given y, find best estimation of a

N*****r
发帖数: 94
33

不过你这个如果只是1000倍,假定 noise 量级小的话, 取样点够多, fft应该是还是
能给出a的

【在 g****t 的大作中提到】
: 最简单的就是一个信号发生器。
: 输出是3个频率a,b,c的指数函数叠加一些噪声。如果一个a比b,c两个small 1000倍。
: 如何根据输出估计a,b,c
: 你可以试验下
: y = A*exp(-iat)+B*exp(-ibt)+C*exp(-ict) + noise
: Given y, find best estimation of a

g****t
发帖数: 31659
34
数值计算开源选择其实还是非常少的。
一来要代码质量高。二来语法尽量描述性质的,
靠近Matlab,sql风格。
Numpy质量很高,这是python重要加分。
如果要快速开发。(我们假设你要研究如何用
ann算天气预报。)
那么Python/R/lua恐怕prototype
阶段还是仅有的选择了。
这三种不会有大的不可预测的问题。不会出现你开发了半年推倒重来的情况,或者憋死
一个月的
情况。
Rust,go,jvm上的各种线性代数,都没有大公司背书。
我觉得还是挺短视的。goog如果出钱弄好这块,其实我觉得
Go有潜力占住ML市场的。
有design pattern功力的多半不会自费武功去go.
做数值计算的多数是看不上design pattern等东西的。
思路很多其实还是fortran。
h**********c
发帖数: 4120
35
LAPACK,LINPACK, TNT这些美国高校大研究所已经用了多少年,
比较偏门DS tool, AUTO这些也都是开源,
你说那个ANN怎么做stability analysis,保证计算收敛?
可能就是用随机数来收敛吧。
估计你也就是在高校里,看看这各大公司热热闹闹,里面各种戏剧帮派。
w***g
发帖数: 5958
36
我以前要算一个啥规划, 只有NAG里有. 那时的数值算法真是数值算法.
现在的ANN, 就是乱棍打死老拳师.

【在 h**********c 的大作中提到】
: LAPACK,LINPACK, TNT这些美国高校大研究所已经用了多少年,
: 比较偏门DS tool, AUTO这些也都是开源,
: 你说那个ANN怎么做stability analysis,保证计算收敛?
: 可能就是用随机数来收敛吧。
: 估计你也就是在高校里,看看这各大公司热热闹闹,里面各种戏剧帮派。

h**********c
发帖数: 4120
37
现在没时间研究ANN,不是说NN这个东西完全没有用,没有严格的论证和可信的模型。
我老觉能一个编程语言静态分析系统,来避免常识性错误和reinvent wheels是有点用
处。不过觉得,你说那是knoledge base不知道区别在哪里,也许NN在路径规划和搜索
上有优势,没研究过。

【在 w***g 的大作中提到】
: 我以前要算一个啥规划, 只有NAG里有. 那时的数值算法真是数值算法.
: 现在的ANN, 就是乱棍打死老拳师.

g****t
发帖数: 31659
38
估计本版只有我一个知道你的id是Poincare书里定义的了,Lol
数值计算的基础部分,也是最早发展的部分是
Ode, pde求解。早先都是实验物理。
别的不和你多讲了,反正你也听不懂。


: LAPACK,LINPACK, TNT这些美国高校大研究所已经用了多少年,

: 比较偏门DS tool, AUTO这些也都是开源,

: 你说那个ANN怎么做stability analysis,保证计算收敛?

: 可能就是用随机数来收敛吧。

: 估计你也就是在高校里,看看这各大公司热热闹闹,里面各种戏剧帮派。



【在 h**********c 的大作中提到】
: 现在没时间研究ANN,不是说NN这个东西完全没有用,没有严格的论证和可信的模型。
: 我老觉能一个编程语言静态分析系统,来避免常识性错误和reinvent wheels是有点用
: 处。不过觉得,你说那是knoledge base不知道区别在哪里,也许NN在路径规划和搜索
: 上有优势,没研究过。

g****t
发帖数: 31659
39
现在是收获期。但是常见的ann技术已经出现降价的趋势了。
另外像我的工作,一个算法会装机上亿的量。
矩阵计算是不能出错的。massive industry products的软件以及算法的要求和别的行
业不同。
g****t
发帖数: 31659
40
大公司的优势在于scale out出去的时候成本低,以势压人。
但是难度就在于细节正确非常难。可以说extremely difficult。
产品走的量大了,什么毛病都会被市场打出来。质量不够,
就会影响share.例如手机前后摄像头的算法软件。还有audio写个小段子容易,量大了
就极难
相关主题
Hi Thrust, I am a little worried about my privacy有多少人觉得这一波AI是忽悠?
我觉得王垠对Agile的看法很有道理啊求救:2个dense matrix的乘法
各位轮子派?你们用到后面不觉得蛋疼?大家编程用得algebra lib是什么? (转载)
进入Programming版参与讨论
d***a
发帖数: 13752
41
高性能数值计算暂时还是要靠FORTRAN和C,Java也凑合,再加上GPU编程,用MPI/PVM搞
并行。憋一个月不算什么,有些程序可以运行几个月才出结果,那么就值得多花点时间
,造跑得快一点的轮子。

【在 g****t 的大作中提到】
: 数值计算开源选择其实还是非常少的。
: 一来要代码质量高。二来语法尽量描述性质的,
: 靠近Matlab,sql风格。
: Numpy质量很高,这是python重要加分。
: 如果要快速开发。(我们假设你要研究如何用
: ann算天气预报。)
: 那么Python/R/lua恐怕prototype
: 阶段还是仅有的选择了。
: 这三种不会有大的不可预测的问题。不会出现你开发了半年推倒重来的情况,或者憋死
: 一个月的

z*********e
发帖数: 10149
42
数值计算最基础也最重要的是稳定性和收敛性
不管什么nn,无法证明稳定性的话都白搭

:估计本版只有我一个知道你的id是Poincare书里定义的了,Lol
h**********c
发帖数: 4120
43
我老最近几年比较相信事实比面子重要,
所以如果真的不懂被人揭短也不是太在意。有人肯给你辅导修改论文是你运气。
其实数学史上很重要一段在意大利,两个数学家PK就是解polynomial,解到3次5次在文
艺复兴时代(或更早)就是很难的了,而输家面临出局从此丢饭碗的下场。大家都很
serious,和陈真打擂台一样。所以为什么欧洲出身的数学家把galois的位置放的很高。
我也不是数学家,也不想和你打雷台。不过觉得你觉得你最近说话怪怪的,我的解释是
各为其主。
说东卫讲的乱棍打死老拳师很形像,不过老拳师混了那么多年没些弟子裙带,被打死活
该,不过接下来可能就是群欧。叔觉得押对了上一注,所以要求赌场cash out套现,享
受一下物质文明精神文明双丰收.
东卫所言的乱棍,没有明说,我觉的那个LAXIGILE,disruption mode,加上这个NN。你
可以说这是科学,证伪或证明都是远期。
实际上我们看到的是公交车上群欧司机,飞机上砸驾驶舱,所谓的学术技术会议霸住麦
克风不让别人说话。
这乱棍就是哲学上的虚无主义,没有真理,规律,谁霸住麦克风谁就随心所欲,甚至胡
说八道,指鹿为马。没有论证,设计,验收,没有标准,没有条例规范(这里是指工程
技术上的,政治上的已经太多)。但是事实还是人类长期总结的物理数学工程上的认识
,还是适用于人类对自然的再认识。

【在 g****t 的大作中提到】
: 估计本版只有我一个知道你的id是Poincare书里定义的了,Lol
: 数值计算的基础部分,也是最早发展的部分是
: Ode, pde求解。早先都是实验物理。
: 别的不和你多讲了,反正你也听不懂。
:
:
: LAPACK,LINPACK, TNT这些美国高校大研究所已经用了多少年,
:
: 比较偏门DS tool, AUTO这些也都是开源,
:
: 你说那个ANN怎么做stability analysis,保证计算收敛?
:
: 可能就是用随机数来收敛吧。
:
: 估计你也就是在高校里,看看这各大公司热热闹闹,里面各种戏剧帮派。

g****t
发帖数: 31659
44
我有PDE数值计算book chapter,论文,专利,NSF proposal。
毕业论文有几百页的Mathematica代码。
这些都是我十多年前干的事情。
当时方向比较偏。研究方向是数值计算的代数问题。
类似于微分代数。希望能对PDE找到高斯消去那样的计算scheme。
哦,本站那时候我还有不少与此有关的帖子。
你有这方面的高见,我仍然愿意听听.不然就别套近乎了。

【在 h**********c 的大作中提到】
: 我老最近几年比较相信事实比面子重要,
: 所以如果真的不懂被人揭短也不是太在意。有人肯给你辅导修改论文是你运气。
: 其实数学史上很重要一段在意大利,两个数学家PK就是解polynomial,解到3次5次在文
: 艺复兴时代(或更早)就是很难的了,而输家面临出局从此丢饭碗的下场。大家都很
: serious,和陈真打擂台一样。所以为什么欧洲出身的数学家把galois的位置放的很高。
: 我也不是数学家,也不想和你打雷台。不过觉得你觉得你最近说话怪怪的,我的解释是
: 各为其主。
: 说东卫讲的乱棍打死老拳师很形像,不过老拳师混了那么多年没些弟子裙带,被打死活
: 该,不过接下来可能就是群欧。叔觉得押对了上一注,所以要求赌场cash out套现,享
: 受一下物质文明精神文明双丰收.

g****t
发帖数: 31659
45
难点在coding。
不在于应用数学的那些皮毛.

【在 z*********e 的大作中提到】
: 数值计算最基础也最重要的是稳定性和收敛性
: 不管什么nn,无法证明稳定性的话都白搭
:
: :估计本版只有我一个知道你的id是Poincare书里定义的了,Lol
: :

m******r
发帖数: 1033
46
guest显然是搞数值计算的,曾经(也许现在)靠这个吃饭的。
heteroclinic似是数字信号处理的?
wdong是sde转人工智能外加培训和企业家的。
我呢,是首席民间科学家,没事喜欢瞎捉摸,但不靠这个吃饭, 所以旁观者门清。
n******g
发帖数: 2201
47
魏老师呢?早期国内中关村Delphi 程序员?
[在 magliner (magliner) 的大作中提到:]
:guest显然是搞数值计算的,曾经(也许现在)靠这个吃饭的。
:heteroclinic似是数字信号处理的?
:wdong是sde转人工智能外加培训和企业家的。
:我呢,是首席民间科学家,没事喜欢瞎捉摸,但不靠这个吃饭, 所以旁观者门清。
g****t
发帖数: 31659
48
我真的是neural network的master,然后F2来美国.
找不到奖学金才转pde/ode的。
前面不是有人问neural network的稳定性吗。
我多年前写过一个sequence2sequence的RNN estimator,用qualify/disqualify
好坏点的办法确保稳定性。但我只能做一维。
等5年后我告诉大家诀窍在哪。

【在 m******r 的大作中提到】
: guest显然是搞数值计算的,曾经(也许现在)靠这个吃饭的。
: heteroclinic似是数字信号处理的?
: wdong是sde转人工智能外加培训和企业家的。
: 我呢,是首席民间科学家,没事喜欢瞎捉摸,但不靠这个吃饭, 所以旁观者门清。

m******r
发帖数: 1033
49
teacherwei, i have no comment. 水平也许很高,赶走了goodbug, 我们看客也没得看
了。 好比看小说,两边正在捉对厮杀,大家都希望好人被打的鼻青脸肿,先写丢了性
命,可正义战胜了邪恶,最终还是赢了,至少和坏人同归于尽。 好比heteroclinic说
的<陈真>,最后一把刀插自己身上,和后面日本鬼子穿了葫芦串。
可惜到了最后好人却偏偏输了,还退出江湖,坏人还在这儿蹦跶呢。
连点悲剧美都没有, 令人不胜唏嘘啊。

【在 n******g 的大作中提到】
: 魏老师呢?早期国内中关村Delphi 程序员?
: [在 magliner (magliner) 的大作中提到:]
: :guest显然是搞数值计算的,曾经(也许现在)靠这个吃饭的。
: :heteroclinic似是数字信号处理的?
: :wdong是sde转人工智能外加培训和企业家的。
: :我呢,是首席民间科学家,没事喜欢瞎捉摸,但不靠这个吃饭, 所以旁观者门清。

h**********c
发帖数: 4120
50
感觉您是我聊过最大牛的ID,
我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士我见过这人老
中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛。老板说过一
个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不过老板给过我
一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉得我不行,我
是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南大数学系的或
着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多嘴,水平可定
是有的,可能也在MMT混,有两三面之源,没有深交。
说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密,我还在雀斑
出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。
单说科学计算,我干过什么吧,解ODE,我做过两个模型,一个是CR3BP,通俗说是地球
月亮卫星三体,光这个模型的MATLAB code画图就的上千行MATLAB,python 加上C++/
OPENGL的code 就几万,行数不重要,重要的是理解什么东西干了什么。我觉得除了卖
力跑老板的工具以外,我认为我个人提出了通过能量守恒限定边界条件从使异构环轨道
的计算不是简单暴力/穷举。还做过一个lorenz manifold,帮助老板整理异构环轨道的
计算边界条件。两个模型都用不同的工具进行了大量的可视话。比如说过,康乃尔哪个
教授对lorenz模型的strange attractor的异构环轨道很感兴趣,不知道老板给他看了
没有。总之,这条路我没走通,混了个硕士,靠谢谢java, sql, node js混个饭吃。
从数学上讲,主要是理解gauss collocation的收敛特性,比起来,runge kuta,
implicit euler,一个不稳定,一个精度太低。用第三方工具计算floquet multiplier,
确定平衡点轨道的特性,卫星三体那个有三对,所以我问你明不明白还是个问题?过内
这个圈里,没见过什么paper.美过华人好象有一些名字。比较独门是Keller pseudo
arc length continuation,您要能说说这个东西干什么的,以后您在任何场合白话,i
don't give any sh*t. 条件是我不查资料评记忆,你也不能查资料。
叔在国内时候老板之一,酒店总经理兰方人讲,没有我不做的生意。和您套近乎奏是为
了生意,国内热钱很多很多。别人的博士论文我也看过百十来篇,SOSO。那真是一身功
夫,一樽还泪江月。
客观讲数的同事,力争吸引眼球和课读性,趣味性。

【在 g****t 的大作中提到】
: 我有PDE数值计算book chapter,论文,专利,NSF proposal。
: 毕业论文有几百页的Mathematica代码。
: 这些都是我十多年前干的事情。
: 当时方向比较偏。研究方向是数值计算的代数问题。
: 类似于微分代数。希望能对PDE找到高斯消去那样的计算scheme。
: 哦,本站那时候我还有不少与此有关的帖子。
: 你有这方面的高见,我仍然愿意听听.不然就别套近乎了。

相关主题
请问:稀疏矩阵运算有什么好的C++库?Switch from Matlab to C(C++)?
Matrix calculation in C++do you use blas/lapack?
LAPACK/ARPACK under windowsLinux下运行lapack和blas的问题
进入Programming版参与讨论
g****t
发帖数: 31659
51
你看过证明洛伦兹系统是混沌的那个计算机程序证明吗?
拿过欧洲青年数学奖的。可能和你方向差不多。可能就是你说的康奈尔的教授。
你的三体程序在slow manifold怎么处理的?正常办法会一步只走一小点或者不走了。


: 感觉您是我聊过最大牛的ID,

: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士
我见过
这人老

: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛
。老板
说过一

: 个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不
过老板
给过我

: 一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉
得我不
行,我

: 是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南
大数学
系的或

: 着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多
嘴,水
平可定

: 是有的,可能也在MMT混,有两三面之源,没有深交。

: 说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密
,我还
在雀斑

: 出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。



【在 h**********c 的大作中提到】
: 感觉您是我聊过最大牛的ID,
: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士我见过这人老
: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛。老板说过一
: 个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不过老板给过我
: 一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉得我不行,我
: 是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南大数学系的或
: 着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多嘴,水平可定
: 是有的,可能也在MMT混,有两三面之源,没有深交。
: 说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密,我还在雀斑
: 出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。

w***g
发帖数: 5958
52
这个主页逼格很高!咱应该把贴主页的风气推起来。

【在 h**********c 的大作中提到】
: 感觉您是我聊过最大牛的ID,
: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士我见过这人老
: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛。老板说过一
: 个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不过老板给过我
: 一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉得我不行,我
: 是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南大数学系的或
: 着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多嘴,水平可定
: 是有的,可能也在MMT混,有两三面之源,没有深交。
: 说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密,我还在雀斑
: 出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。

l*********s
发帖数: 5409
53
同感,不过排版老式了点,不够好看。

【在 w***g 的大作中提到】
: 这个主页逼格很高!咱应该把贴主页的风气推起来。
l*******m
发帖数: 1096
54
其实大牛可以搞搞Stochastic gradient descent on Riemannian manifolds,用到
deep learning上
还有hinton的capsules是不是可以和黎曼几何挂一挂?

【在 h**********c 的大作中提到】
: 感觉您是我聊过最大牛的ID,
: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士我见过这人老
: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛。老板说过一
: 个人在做navier stoke,可能打错了。我这么专的人,当然很感兴趣。不过老板给过我
: 一个小模型,让我解解看,数值解,ODE,我两个晚上没弄出来,可能觉得我不行,我
: 是纯CS,偏重计算机图形,本加硕都是加拿大的。老板那个硕士据说是南大数学系的或
: 着计算机双修的,奥数金牌等等,没考证过,门清者自门清,不需要我多嘴,水平可定
: 是有的,可能也在MMT混,有两三面之源,没有深交。
: 说到我做过什么基本都在我个人主页我以前公开过,从来也不是什么秘密,我还在雀斑
: 出没,所以你如果www.heteroclinic.net有什么有兴趣的东西,我寻找VC。

g****t
发帖数: 31659
55
丘成桐已经进来做DL谈话了。前几天有个人转给我的。
你找找看。
有人说,最怕的就是外行大牛进来抢funding。
人家什么都不懂,照样糊的你一愣一愣。
还一个例子,还一个大牛除了让自己老婆开了个drive.ai偷资料,在百度屁都没干,一
年1200万美元。


: 其实大牛可以搞搞Stochastic gradient descent on Riemannian
manifolds,
用到

: deep learning上

: 还有hinton的capsules是不是可以和黎曼几何挂一挂?



【在 l*******m 的大作中提到】
: 其实大牛可以搞搞Stochastic gradient descent on Riemannian manifolds,用到
: deep learning上
: 还有hinton的capsules是不是可以和黎曼几何挂一挂?

m*****n
发帖数: 3575
56
你对PDE数值计算这么懂
请问抛物线型和椭圆型最方便的数值解法是什么?

【在 g****t 的大作中提到】
: 丘成桐已经进来做DL谈话了。前几天有个人转给我的。
: 你找找看。
: 有人说,最怕的就是外行大牛进来抢funding。
: 人家什么都不懂,照样糊的你一愣一愣。
: 还一个例子,还一个大牛除了让自己老婆开了个drive.ai偷资料,在百度屁都没干,一
: 年1200万美元。
:
:
: 其实大牛可以搞搞Stochastic gradient descent on Riemannian
: manifolds,
: 用到

h**********c
发帖数: 4120
57
我不是太理解这个slow manifold的概念,就说我论文里的计算,我导师自己给个概念
叫orbit continuation,在平衡点附近grow 一个周期轨道,不断算这个轨道的
continuation,因为连接两个不同平衡点的周期轨道的轨道也是异构环。我算的都是2
manifold.我又查了一下,你说的slow manifold不是homotopy,可能开会的时候听一个
人问过collision,是不是这个意思?
我这人挨溜号,搞计算的时候图形玩的多,当了马公,还老想着计算的东西。
感觉还是有点多动症,还好没有发展成自闭或弱智。

【在 g****t 的大作中提到】
: 你看过证明洛伦兹系统是混沌的那个计算机程序证明吗?
: 拿过欧洲青年数学奖的。可能和你方向差不多。可能就是你说的康奈尔的教授。
: 你的三体程序在slow manifold怎么处理的?正常办法会一步只走一小点或者不走了。
:
:
: 感觉您是我聊过最大牛的ID,
:
: 我这点东西也七八年不碰了就是晒晒吧,PDE我没做过,老板有一个硕士
: 我见过
: 这人老
:
: 中,给他做过PDE网格,三角化,我理解是通过网格的精细化来保证收敛
: 。老板

h**********c
发帖数: 4120
58
我找过一阵子手机pad friendly的framework,后来太忙,不了了之

【在 w***g 的大作中提到】
: 这个主页逼格很高!咱应该把贴主页的风气推起来。
h**********c
发帖数: 4120
59
随机的东西,好象正好是homotopy/continuation这些人相反的一面,不太清楚。

【在 l*******m 的大作中提到】
: 其实大牛可以搞搞Stochastic gradient descent on Riemannian manifolds,用到
: deep learning上
: 还有hinton的capsules是不是可以和黎曼几何挂一挂?

g****t
发帖数: 31659
60
我不跟人social


: 你对PDE数值计算这么懂

: 请问抛物线型和椭圆型最方便的数值解法是什么?



【在 m*****n 的大作中提到】
: 你对PDE数值计算这么懂
: 请问抛物线型和椭圆型最方便的数值解法是什么?

相关主题
How to use multi-core to speed Python programC++里面 class的数据成员的顺序是什么样的?
what kind of language this is and how to change it? (转载)C++ 做线性代数,方便使用的库?
问一个C++下的Bug(Linux下)数组,结构,类在数值计算中哪个快
进入Programming版参与讨论
g****t
发帖数: 31659
61
我的意思很简单,就是空间里有一块地方,simulation很多步,然而结果走的极慢。但
是它又不是invariant manifold,所以会有数值困难。


: 我不是太理解这个slow manifold的概念,就说我论文里的计算,我导师自己给
个概念

: 叫orbit continuation,在平衡点附近grow 一个周期轨道,不断算这个轨道的

: continuation,因为连接两个不同平衡点的周期轨道的轨道也是异构环。我算的
都是2

: manifold.我又查了一下,你说的slow manifold不是homotopy,可能开会的时候
听一个

: 人问过collision,是不是这个意思?

: 我这人挨溜号,搞计算的时候图形玩的多,当了马公,还老想着计算的东西。

: 感觉还是有点多动症,还好没有发展成自闭或弱智。



【在 h**********c 的大作中提到】
: 随机的东西,好象正好是homotopy/continuation这些人相反的一面,不太清楚。
h**********c
发帖数: 4120
62
你有具体的例题吗?
我个人的理解,你在一个有洞的空间里,使用依赖于连续特性的数学法则。好象是一个
悖论。也许这就是为什么我看到的实际的计算模型维持在低纬度以保证有限参数的扰动
保证流形的不变和连续。

【在 g****t 的大作中提到】
: 我的意思很简单,就是空间里有一块地方,simulation很多步,然而结果走的极慢。但
: 是它又不是invariant manifold,所以会有数值困难。
:
:
: 我不是太理解这个slow manifold的概念,就说我论文里的计算,我导师自己给
: 个概念
:
: 叫orbit continuation,在平衡点附近grow 一个周期轨道,不断算这个轨道的
:
: continuation,因为连接两个不同平衡点的周期轨道的轨道也是异构环。我算的
: 都是2
:
: manifold.我又查了一下,你说的slow manifold不是homotopy,可能开会的时候
: 听一个

g****t
发帖数: 31659
63
最简单的就是一个信号发生器。
输出是3个频率a,b,c的指数函数叠加一些噪声。如果一个a比b,c两个small 1000倍。
如何根据输出估计a,b,c
你可以试验下
y = A*exp(-iat)+B*exp(-ibt)+C*exp(-ict) + noise
Given y, find best estimation of a

【在 h**********c 的大作中提到】
: 你有具体的例题吗?
: 我个人的理解,你在一个有洞的空间里,使用依赖于连续特性的数学法则。好象是一个
: 悖论。也许这就是为什么我看到的实际的计算模型维持在低纬度以保证有限参数的扰动
: 保证流形的不变和连续。

h**********c
发帖数: 4120
64
我老当年学传热学时候公式比这复杂多,不过老师也不明白,或着她认为我们不需要明
白,知道算什么查那本书,不过传热的计算多数可以在工程许用范围内验证,不是蒙人
的随机数。工程上科一半要做实验,验证理论,每个人都能重复。
你不能建立你的模型象auto 手册,搞成invariant Ode模型,search auto07p.pdf 或
auto2000.pdf
比如象bruta problem,Sandstede 的模型.
我说的仅供参考。

【在 g****t 的大作中提到】
: 最简单的就是一个信号发生器。
: 输出是3个频率a,b,c的指数函数叠加一些噪声。如果一个a比b,c两个small 1000倍。
: 如何根据输出估计a,b,c
: 你可以试验下
: y = A*exp(-iat)+B*exp(-ibt)+C*exp(-ict) + noise
: Given y, find best estimation of a

N*****r
发帖数: 94
65

我对这类的数值稳定性的问题也很感兴趣
原来比较关心的是天体的长时间轨道模拟的数值稳定性,

【在 g****t 的大作中提到】
: 最简单的就是一个信号发生器。
: 输出是3个频率a,b,c的指数函数叠加一些噪声。如果一个a比b,c两个small 1000倍。
: 如何根据输出估计a,b,c
: 你可以试验下
: y = A*exp(-iat)+B*exp(-ibt)+C*exp(-ict) + noise
: Given y, find best estimation of a

N*****r
发帖数: 94
66

不过你这个如果只是1000倍,假定 noise 量级小的话, 取样点够多, fft应该是还是
能给出a的

【在 g****t 的大作中提到】
: 最简单的就是一个信号发生器。
: 输出是3个频率a,b,c的指数函数叠加一些噪声。如果一个a比b,c两个small 1000倍。
: 如何根据输出估计a,b,c
: 你可以试验下
: y = A*exp(-iat)+B*exp(-ibt)+C*exp(-ict) + noise
: Given y, find best estimation of a

h**********c
发帖数: 4120
67
这是数学问题,法语大学里基本是研究生6年纪有足够的基础搞清楚。
具体是operator theory, function analysis这些课。我也没学过,半拉卡几的听了一
些。具体的技术就是spectral radius什么的。
这几年混饭吃人太多,老老实实学这些的人都上不了岗。

【在 N*****r 的大作中提到】
:
: 不过你这个如果只是1000倍,假定 noise 量级小的话, 取样点够多, fft应该是还是
: 能给出a的

h**********c
发帖数: 4120
68
忘了说是法语大学数学专业,嘴上不能示弱,
叔稍微一眨眼,就会成阿比那夫食谱上的一道recipe.
于是我们每天有固定的几个小时在机房的角落里roar each other.尽管机房里有噪音很
大的服务器,不是服务器噪音都很大,有噪音很大的服务器,其它人都会意的点头,
that is what they are supposed to do.
finally, both beasts flet felt tired and retreated under the cover of the
dearkness of the jungle's nightfall.

【在 h**********c 的大作中提到】
: 这是数学问题,法语大学里基本是研究生6年纪有足够的基础搞清楚。
: 具体是operator theory, function analysis这些课。我也没学过,半拉卡几的听了一
: 些。具体的技术就是spectral radius什么的。
: 这几年混饭吃人太多,老老实实学这些的人都上不了岗。

1 (共1页)
进入Programming版参与讨论
相关主题
How to use multi-core to speed Python program我觉得王垠对Agile的看法很有道理啊
what kind of language this is and how to change it? (转载)各位轮子派?你们用到后面不觉得蛋疼?
问一个C++下的Bug(Linux下)有多少人觉得这一波AI是忽悠?
C++里面 class的数据成员的顺序是什么样的?求救:2个dense matrix的乘法
C++ 做线性代数,方便使用的库?大家编程用得algebra lib是什么? (转载)
数组,结构,类在数值计算中哪个快请问:稀疏矩阵运算有什么好的C++库?
卫东大神来说说阿尔法狗横扫棋坛这事吧Matrix calculation in C++
Hi Thrust, I am a little worried about my privacyLAPACK/ARPACK under windows
相关话题的讨论汇总
话题: br话题: 计算话题: 数值话题: pde话题: 轨道