用QMOM求解出现负的Abscissas

[Astremon]

用QMOM求解Pure aggregation中出现负的Abscissas，然后在公式（37）中就出现negative number cannot be raised to a fractional power 。有经验的大神怎么处理的呢。。

[Astremon]

@astremon 文献传送门

[东岳]

你是怎么发现的负值？无对流问题么？减小时间步长试试。

[Astremon]

@李东岳 总是在一个时间点出现负的abscissas,提示错误“negative number cannot be raised to a fractional power ”。减小时间步长也不管用哟，，，。李老师这个是代码可以请您看看吗？
可能是算法的局限性？

[Astremon]

@李东岳 没有对流。就是homogeneous 常数聚并核函数。

[东岳]

73, 73行是你的源项么？

temp_sum1[i2,i3] = w_i[i2]*w_i[i3]*(r_i[i2]**3 + r_i[i3]**3)**(m/3.)  ##
temp_sum2[i2,i3] = w_i[i2]*w_i[i3]*r_i[i2]**m

对应的方程是什么？

[Astremon]

@李东岳 恩，是源项。是图中的公式（37）和（38）

[Astremon]

@李东岳 其中$\beta (L,\lambda) = C_0 = 1 $

[Astremon]

@李东岳 现在只能算到1.4s左右，超过1.4s出现负特征尺寸了。

[东岳]

出现负值对应的m具体值是多少？

[Astremon]

@李东岳

moments  --->[0.574993, 0.670135, 0.805031, 1, 1.29021, 1.73242]

[东岳]

我试试，不过感觉你的算法没毛病，M3守恒。
你考虑过用Wheeler求节点么？比PD要强健。