2018年2月10日
stiff chemistry模型出现NaN错误
摘要： 通过定位可以看到，是usr_rates.f中出现了奇异值，因为我的代码中有这样一句话： 而同时我的ConH2在声明后没有赋初值，因此，当X_g(IJK,H2) < c_Limiter后，ConH2就是一个随机值，到后面参与计算就会出现奇异值的问题，所以只需要在前面声明之后赋初值 ConH2 = ZE 阅读全文
posted @ 2018-02-10 21:41 JayYin 阅读(307) 评论(0) 推荐(0) 编辑