通过定位可以看到,是usr_rates.f中出现了奇异值,因为我的代码中有这样一句话:
而同时我的ConH2在声明后没有赋初值,因此,当X_g(IJK,H2) < c_Limiter后,ConH2就是一个随机值,到后面参与计算就会出现奇异值的问题,所以只需要在前面声明之后赋初值
ConH2 = ZERO即可
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步