gmres MATLAB中的relres

计算科学 matlab 迭代法 格瑞斯
2021-12-18 04:42:55

我认为 MATLAB 中的 relres 是 relres = norm(M(bA x))/norm(M\b) 的形式,当它小于 tol 时停止迭代。我想知道如何将 relres 更改为 norm((bA x))/norm(b)。或者使用其他代码而不是内置代码

1个回答

该定义的优点是计算 relres 来自 GMRES 迭代“免费”。您可以在没有预处理器的情况下切换到其他定义,但是您必须使用更多操作来计算它。

我不认为你能比在每一步从头开始据我在文档中看到的,Matlab没有回调来引入自定义停止标准,因此您必须自己修改其源代码(可用 --- 类型来阅读它)。bAxgmresedit gmres