python – 获取scipy的gmres迭代方法的迭代次数

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了python – 获取scipy的gmres迭代方法的迭代次数，小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含1028字，纯文字阅读大概需要2分钟。

内容图文

我想知道scipy.sparse.linalg.gmres需要多少次迭代才能收敛,但似乎没有任何争论.有一个maxiter参数,可以设置终止方法,但没有任何显示它正在采取的迭代次数.有人可以帮我弄这个吗？

解决方法:

为了说明@ cel的评论,你可以实现一个简单的计数器类,如下所示：

class gmres_counter(object):
    def __init__(self, disp=True):
        self._disp = disp
        self.niter = 0
    def __call__(self, rk=None):
        self.niter += 1
        if self._disp:
            print('iter %3i\trk = %s' % (self.niter, str(rk)))

你会像这样使用它：

import numpy as np
from scipy.sparse.linalg import gmres

A = np.random.randn(10, 10)
b = np.random.randn(10)

counter = gmres_counter()

x, info = gmres(A, b, callback=counter)
# iter   1        rk = 0.999558746968
# iter   2        rk = 0.960490282387
# iter   3        rk = 0.945887432101
# iter   4        rk = 0.931790454216
# iter   5        rk = 0.877655067142
# iter   6        rk = 0.739596963239
# iter   7        rk = 0.677886023198
# iter   8        rk = 0.52015135005
# iter   9        rk = 0.168298366785
# iter  10        rk = 9.22692033803e-16

print(counter.niter)
# 10

内容总结

以上是互联网集市为您收集整理的python – 获取scipy的gmres迭代方法的迭代次数全部内容，希望文章能够帮你解决python – 获取scipy的gmres迭代方法的迭代次数所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端