stiff chemistry模型出现NaN错误
Posted JayYin
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了stiff chemistry模型出现NaN错误相关的知识，希望对你有一定的参考价值。
通过定位可以看到，是usr_rates.f中出现了奇异值，因为我的代码中有这样一句话：
而同时我的ConH2在声明后没有赋初值，因此，当X_g(IJK,H2) < c_Limiter后，ConH2就是一个随机值，到后面参与计算就会出现奇异值的问题，所以只需要在前面声明之后赋初值
ConH2 = ZERO即可
以上是关于stiff chemistry模型出现NaN错误的主要内容，如果未能解决你的问题，请参考以下文章
线性模型函数 lm() 错误：外部函数调用中的 NA/NaN/Inf (arg 1)