用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要强健。