用QMOM求解出现负的Abscissas
-
用QMOM求解Pure aggregation中出现负的Abscissas，然后在公式（37）中就出现negative number cannot be raised to a fractional power 。有经验的大神怎么处理的呢。。
-
-
你是怎么发现的负值？无对流问题么？减小时间步长试试。
-
@李东岳 总是在一个时间点出现负的abscissas,提示错误“negative number cannot be raised to a fractional power ”。减小时间步长也不管用哟，，，。李老师这个是代码可以请您看看吗？
可能是算法的局限性？
-
@李东岳 没有对流。就是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对应的方程是什么？
-
@李东岳 恩，是源项。是图中的公式（37）和（38）
-
@李东岳 其中$\beta (L,\lambda) = C_0 = 1 $
-
@李东岳 现在只能算到1.4s左右，超过1.4s出现负特征尺寸了。
-
出现负值对应的m具体值是多少？
-
@李东岳
moments --->[0.574993, 0.670135, 0.805031, 1, 1.29021, 1.73242]
-
我试试，不过感觉你的算法没毛病，M3守恒。
你考虑过用Wheeler求节点么？比PD要强健。
