Question

急求解答！fortran运行时出现run-time error M6021：MATH - log:DOMAIN error,怎解答，

debug时出现：user breakpoint called from code at 0x7c92120e.求帮忙解答

Answer 1

这是很常见的算法错误，通常这种错误是不能再编译的时候发现的，因为者往往不是由语法错误引起的。M6201的错误信息可以到编译器的帮助文件中查找，大意就是自然对数的底超出范围，也就是说在计算log(N)的时候，括号里的N变成负数了，这时候运行的程序就会报错。这需要程序员对算法进行排查优化才能解决这种失误。所以，请自己检查并改进算法的代码吧。

追问

我在程序中，log（）加了abs（），不可能是负的

追答

你是怎么加的，把代码贴出来

追问

tao(j)=sigma(i)*(mu0+a*log(abs(yout(2)/v1))+psi(i))
就是这样加的

追答

还有一种情况可能导致log预算错误的就是，底数溢出。你的代码看起来所有变量应该都是单精度的，这样log后边括号里的数字就会有上限，你应该检查一下可能是v1变为0（或是一个非常小的数），也可能是yout(2)变为一个很大的数。试试看

追问

v1的数量级是很小e-9,那要怎么办

追答

使用双精度的变量试试，如果不行，就换下算法，把握住变量的变化程度。