通过定位可以看到,是usr_rates.f中出现了奇异值,因为我的代码中有这样一句话:
而同时我的ConH2在声明后没有赋初值,因此,当X_g(IJK,H2) < c_Limiter后,ConH2就是一个随机值,到后面参与计算就会出现奇异值的问题,所以只需要在前面声明之后赋初值
ConH2 = ZERO即可
通过定位可以看到,是usr_rates.f中出现了奇异值,因为我的代码中有这样一句话:
而同时我的ConH2在声明后没有赋初值,因此,当X_g(IJK,H2) < c_Limiter后,ConH2就是一个随机值,到后面参与计算就会出现奇异值的问题,所以只需要在前面声明之后赋初值
ConH2 = ZERO即可
转载于:https://www.cnblogs.com/Jay-CFD/p/8439824.html