mathematica的bug
mathematica吧
全部回复
仅看楼主
level 4
LatticeData里的An的矩阵好像给成了Bn
E6 E7没有数据
是不是bug
2022年09月03日 06点09分 1
level 4
c[x_] = Sum[x^2^i, {i, 0, Infinity}]
N[ResidueSum[{c[x], Abs[x] < 0.1}, x]]
\:6B63\:5728\:8BA1\:7B97In[5]:= ResidueSum::ncsol: 无法求出在 \!\(\*UnderoverscriptBox[\(\[Sum]\), \(i = 0\), \(\[Infinity]\)]\*SuperscriptBox[\(x\), SuperscriptBox[\(2\), \(i\)]]\) 在解集 Abs[x]<0.1 中的的残差和.
\:6B63\:5728\:8BA1\:7B97In[5]:= NSum::nsnum: 被加数(或者它的导数)x^2^i 在点 i = 46661 不是数值的.
\:6B63\:5728\:8BA1\:7B97In[5]:= NSum::nsnum: 被加数(或者它的导数)x^2^i 在点 i = 46661 不是数值的.
\:6B63\:5728\:8BA1\:7B97In[5]:= NSum::nsnum: 被加数(或者它的导数)x^2^i 在点 i = 46661 不是数值的.
\:6B63\:5728\:8BA1\:7B97In[5]:= General::stop: 在本次计算中,NSum::nsnum 的进一步输出将被抑制.
\:6B63\:5728\:8BA1\:7B97In[5]:= ResidueSum::ncsol: 无法求出在 NSum[x^2^i,{i,0,\[Infinity]}] 在解集 Abs[x]<0.1 中的的残差和.
\:6B63\:5728\:8BA1\:7B97In[5]:= ResidueSum::ncsol: 无法求出在 NSum[x^2^i,{i,0,\[Infinity]}] 在解集 Abs[x]<0.1 中的的残差和.
\:6B63\:5728\:8BA1\:7B97In[5]:= General::stop: 在本次计算中,ResidueSum::ncsol 的进一步输出将被抑制.
Out[5]= ResidueSum[{NSum[x^2^i, {i, 0, \[Infinity]}], Abs[x] < 0.1},
x]
2022年09月03日 07点09分 2
level 4
In[7]:= FunctionPoles[{c[x], Abs[x] < 0.1}, x]
\:6B63\:5728\:8BA1\:7B97In[7]:= FunctionPoles::ncsol: 无法在 Abs[x]<0.1 的解集内求出 \!\(\*UnderoverscriptBox[\(\[Sum]\), \(i = 0\), \(\[Infinity]\)]\*SuperscriptBox[\(x\), SuperscriptBox[\(2\), \(i\)]]\) 的极点.
Out[7]= FunctionPoles[{\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(i = 0\), \(\[Infinity]\)]
\*SuperscriptBox[\(x\),
SuperscriptBox[\(2\), \(i\)]]\), Abs[x] < 0.1}, x]
2022年09月03日 07点09分 3
level 4
这个没有极点 但就是求不出来
2022年09月03日 07点09分 4
吧务
level 15
1楼的问题太模糊,而且我也不懂晶体学,就不发表意见了。
关于2楼和3楼:算不出来一般不认为是bug,就我所知,目前Mathematica里应该没有哪个代数变换函数可以对
未计算的合式进行符号分析。(至少文档里没有相关例子。)——当然了,实践中确实发现过个别孤例,比如SE帖子《Backslide of Limit》(mathematica.stackexchange.com/q/64633/1871)下面的这个例子(注意,在较新的版本里这个现象已经无法重现):
但这个例子没有推广价值。
2022年09月03日 08点09分 5
报错太奇怪了
2022年09月03日 08点09分
@贴吧用户_aJAJeEJ 首先,我想ResidueSum::ncsol这条警告还是很明确的。(“无法求出”。)至于后面那条NSum的警告,这是因为当N[…]内部存在未计算的Sum时,软件会自动调用NSum。这一点在NSum自带帮助的“更多信息与选项”部分提及了。
2022年09月03日 08点09分
level 4
LatticeData["A3", "GeneratorMatrix"] // Normal
输出的是B3的
和英文维基百科写的不一样(Root System)
2022年09月03日 08点09分 6
可能那个An的是在n+1维中 然后被截断了
2022年09月03日 08点09分
有办法在n维中表示 但是有无理数
2022年09月03日 08点09分
level 4
还有个问题 给2f那个函数在单位圆内绘图(把求和上限改成24~36)时内存占用极大
2022年09月03日 08点09分 7
虽然你又没给具体代码,但我好像猜出来是怎么回事了。参9楼。
2022年09月03日 09点09分
吧务
level 15
还有,不相关的问题请尽量别挤一个帖里,你这样搞,今后别人想搜相关内容可能会出现麻烦。
2022年09月03日 08点09分 8
吧务
level 15
虽然7楼说得不清楚,但我好像猜出来LZ指的是啥了。这个问题或许该算是对极端情形的优化不足。总之,下面这样是没问题的:
但是,如果把c的定义中的 2. 改成 2 (注意这里的区别!)就会出现内存狂飙的问题,这似乎和Plot系函数的符号预处理机制有关——使用Echo包裹x,会发现ComplexPlot直接使用了边界值-1-I(注意,是准确值,不是机器精度的-1. - 1. I)去试探c[z]的取值,而c[z]在-1-I处的
准确取值嘛……
上面所说的应该不是唯一的原因,因为把ComplexPlot的第二参数改成{z, -1. - I, 1. + I}虽然也可以回避内存问题,但代码需要执行12秒左右。原因一时搞不清。
2022年09月03日 09点09分 9
c(-1-i)不存在
2022年09月03日 10点09分
2022年09月03日 10点09分
@xzcyr 无穷的不存在
2022年09月03日 10点09分
level 4
Sum[x^f[i],{i,0,Infinity}] 其中f:N->N 为严格递增函数
在单位圆内必定解析
mma还是不够强 不能判断级数的收敛半径
2022年09月03日 10点09分 10
level 4
一个点就能改变巨大的内存 离谱
2022年09月03日 10点09分 11
level 4
为什么时间变长了?(12700H)
2022年09月03日 10点09分 12
13.1
2022年09月03日 10点09分
……嗬,居然有这么严重的性能退化。我测试用的是12.3。
2022年09月03日 11点09分
@xzcyr mma什么破毛病 12.0在以前的破电脑上占用内存较少 时间也短(最开始那个 不优化)
2022年09月03日 12点09分
@xzcyr 电脑配置?
2022年09月03日 12点09分
吧务
level 15
看样子13.1版以降(也许是13.0开始的?)ComplexPlot出现了严重的性能退化。这似乎依旧是试探导致的问题。(可是明明都设了RegionFunction,ComplexPlot还在域外疯狂试探是想干嘛呢……)总之,还能抢救一下:
Clear[c];
c[x_] := Total@Table[x^2.^i, {i, 0, 36}]
ComplexPlot[If[Abs[z] < 1, c[z], 7.2 10^284], {z, -1 - I, 1 + I},
RegionFunction -> Function[{z}, Abs[z] <= 1], Exclusions -> None,
PlotLegends -> Automatic] // AbsoluteTiming
2022年09月03日 12点09分 13
此bug已在版本13.2修复。
2024年04月06日 10点04分
吧务
level 15
刚检查了一下,这个直接试探左边界的奇怪机制在版本9时代不存在。具体哪版引入的不明。
可以使用以下代码进行检查:
Plot[Sin[Echo@x], {x, 123/457, 7}, MaxRecursion -> 0, PlotPoints -> 2]
2022年09月03日 12点09分 14
Echo是10.3版引入的,早期版本自己写一个也容易(当然,未实现全部功能):If[$VersionNumber<10,System`Echo=(Print@#;#)&]
2022年09月03日 12点09分
1