期刊文献+
共找到1,714篇文章
< 1 2 86 >
每页显示 20 50 100
一种基于敏捷集群计算系统的并行GMRES方法
1
作者 何康馨 席国江 陈颖 《无线电通信技术》 北大核心 2024年第1期162-167,共6页
随着通信系统和人工智能的飞速发展,以智慧城市、智慧工厂和智能制造等为代表的多种新型应用场景不断涌现,使得通信、感知和计算等系统的一体化成为技术发展的新趋势。人工智能新型应用场景对大规模高效敏捷计算提出了新的要求,基于敏... 随着通信系统和人工智能的飞速发展,以智慧城市、智慧工厂和智能制造等为代表的多种新型应用场景不断涌现,使得通信、感知和计算等系统的一体化成为技术发展的新趋势。人工智能新型应用场景对大规模高效敏捷计算提出了新的要求,基于敏捷集群计算系统,提出了一种并行广义最小残差(Generalized Minimal Residual, GMRES)方法,主要通过并行矩阵向量乘法和并行高瘦矩阵QR(Tall and Skinny QR,TSQR)分解实现Krylov子空间的高效并行构造,充分利用集群计算系统的计算和通信性能,实现大规模线性方程组Ax=b的快速求解,其中A为一个n×n的矩阵,在工程实践中,n可达数十万甚至百万规模。通过求解二维泊松方程的有限元离散得到的刚度方程,验证了算法的有效性。 展开更多
关键词 敏捷集群计算 并行广义最小残差方法 KRYLOV子空间 大规模线性方程组
下载PDF
Some Results on the Range-Restricted GMRES Method
2
作者 Yiqin Lin 《Journal of Applied Mathematics and Physics》 2023年第12期3902-3908,共7页
In this paper we reconsider the range-restricted GMRES (RRGMRES) method for solving nonsymmetric linear systems. We first review an important result for the usual GMRES method. Then we give an example to show that the... In this paper we reconsider the range-restricted GMRES (RRGMRES) method for solving nonsymmetric linear systems. We first review an important result for the usual GMRES method. Then we give an example to show that the range-restricted GMRES method does not admit such a result. Finally, we give a modified result for the range-restricted GMRES method. We point out that the modified version can be used to show that the range-restricted GMRES method is also a regularization method for solving linear ill-posed problems. 展开更多
关键词 Nonsymmetric Linear System Krylov Subspace method Arnoldi Process gmres RRgmres
下载PDF
A Novel Method for Linear Systems of Fractional Ordinary Differential Equations with Applications to Time-Fractional PDEs
3
作者 Sergiy Reutskiy Yuhui Zhang +1 位作者 Jun Lu Ciren Pubu 《Computer Modeling in Engineering & Sciences》 SCIE EI 2024年第5期1583-1612,共30页
This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering a... This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering and science.An approximate solution of the system is sought in the formof the finite series over the Müntz polynomials.By using the collocation procedure in the time interval,one gets the linear algebraic system for the coefficient of the expansion which can be easily solved numerically by a standard procedure.This technique also serves as the basis for solving the time-fractional partial differential equations(PDEs).The modified radial basis functions are used for spatial approximation of the solution.The collocation in the solution domain transforms the equation into a system of fractional ordinary differential equations similar to the one mentioned above.Several examples have verified the performance of the proposed novel technique with high accuracy and efficiency. 展开更多
关键词 System of FODEs numerical solution müntz polynomial basis time fractional PDE BSm collocation method
下载PDF
求解具有多个右端项线性方程组的HBFGl-GMRES方法
4
作者 温荣丽 张建华 《江西科学》 2023年第2期217-224,共8页
重启灵活预处理总体广义极小残差(FGl-GMRES)方法是求解具有多个右端项线性方程组的经典迭代方法之一。然而重启会丢失旧循环Krylov子空间产生的信息,因此重启Gl-GMRES方法和FGl-GMRES方法在求解一些具有多个右端项线性方程组时常常出... 重启灵活预处理总体广义极小残差(FGl-GMRES)方法是求解具有多个右端项线性方程组的经典迭代方法之一。然而重启会丢失旧循环Krylov子空间产生的信息,因此重启Gl-GMRES方法和FGl-GMRES方法在求解一些具有多个右端项线性方程组时常常出现收敛速度慢或不收敛的情形。为了加快重启Gl-GMRES方法和FGl-GMRES方法的收敛速度,使用Heavy ball技术,提出了求解具有多个右端项线性方程组的Heavy ball Gl-GMRES方法(HBGl-GMRES)和Heavy ball FGl-GMRES方法(HBFGl-GMRES)。最后,数值实验结果验证了新方法的有效性。 展开更多
关键词 Gl-gmres方法 预处理 灵活Gl-gmres方法 Heavy ball技术
下载PDF
Fully coupled flow-induced vibration of structures under small deformation with GMRES method 被引量:2
5
作者 张立翔 郭亚昆 张洪明 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 2010年第1期87-96,共10页
Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompress... Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompressible viscous fluid flow. The spatial discretization of the formulations is based on the multi-linear interpolating functions by using the finite element method for both the fluid and solid structures. The generalized trapezoidal rule is used to obtain apparently non-symmetric linear equations in an incremental form for the variables of the flow and vibration. The nonlinear convective term and time factors are contained in the non-symmetric coefficient matrix of the equations. The generalized minimum residual (GMRES) method is used to solve the incremental equations. A new stable algorithm of GMRES-Hughes-Newmark is developed to deal with the flow-induced vibration with dynamical fluid-structure interaction in complex geometries. Good agreement between the simulations and laboratory measurements of the pressure and blade vibration accelerations in a hydro turbine passage was obtained, indicating that the GiViRES-Hughes-Newmark algorithm presented in this paper is suitable for dealing with the flow-induced vibration of structures under small deformation. 展开更多
关键词 flow-induced vibration fluid-structure interaction generalized variationalprinciple numerical methods generalized minimum residual gmres method
下载PDF
A Note on the Adaptive Simpler Block GMRES Method
6
作者 Qiaohua Liu Lei Yao Aijing Liu 《Communications on Applied Mathematics and Computation》 2019年第3期435-447,共13页
The adaptive simpler block GMRES method was investigated by Zhong et al.(J Comput Appl Math 282:139-156, 2015) where the condition number of the adaptively chosen basis for the Krylov subspace was evaluated. In this p... The adaptive simpler block GMRES method was investigated by Zhong et al.(J Comput Appl Math 282:139-156, 2015) where the condition number of the adaptively chosen basis for the Krylov subspace was evaluated. In this paper, the new upper bound for the condition number is investigated. Numerical tests show that the new upper bound is tighter. 展开更多
关键词 Large-scale NONSYmmETRIC linear system Krylov SUBSPACE method ADAPTIVE simpler BLOCK gmres Condition number
下载PDF
以对称反对称分裂预条件处理GMRES(m)的不精确牛顿法潮流计算 被引量:11
7
作者 刘凯 陈红坤 +1 位作者 向铁元 高志新 《电网技术》 EI CSCD 北大核心 2009年第19期123-126,共4页
针对大规模电力系统修正方程式高度稀疏的特点,研究了一种基于对称反对称预处理的不精确牛顿法。利用矩阵的对称反对称分裂,提出一种新的预处理子,并将其与GMRES(m)算法相结合,改进潮流计算的收敛性和收敛速度。IEEE300节点系统的计算... 针对大规模电力系统修正方程式高度稀疏的特点,研究了一种基于对称反对称预处理的不精确牛顿法。利用矩阵的对称反对称分裂,提出一种新的预处理子,并将其与GMRES(m)算法相结合,改进潮流计算的收敛性和收敛速度。IEEE300节点系统的计算结果验证了所提算法的有效性。 展开更多
关键词 潮流计算 对称反对称分裂 广义极小残余法(gmres(m)) 预条件处理
下载PDF
基于改进Jacobian-Free Newton-GMRES(m)的电力系统分布式潮流计算 被引量:77
8
作者 陈颖 沈沉 +1 位作者 梅生伟 卢强 《电力系统自动化》 EI CSCD 北大核心 2006年第9期5-8,37,共5页
为了提高互联电网分布式潮流算法的实用性,应尽可能减少协调计算中的交换信息。文中采用隐函数方式表示分布式潮流的边界协调方程,并采用具有自适应预处理能力的NewtonGMRES(m)方法构建分解协调算法。该方法具有Jacobian-Free特性,协调... 为了提高互联电网分布式潮流算法的实用性,应尽可能减少协调计算中的交换信息。文中采用隐函数方式表示分布式潮流的边界协调方程,并采用具有自适应预处理能力的NewtonGMRES(m)方法构建分解协调算法。该方法具有Jacobian-Free特性,协调计算中只需要交换边界节点状态信息。以IEEE标准系统和大规模实际系统为算例的测试结果表明,新的分布式潮流算法具有较高的收敛性,数据交换接口简单,实用性强,适合于求解连续变化的分布式潮流问题。 展开更多
关键词 分布式潮流计算 Newton-gmres(m) 预处理
下载PDF
大区域地下水模拟的预优并行GMRES(m)算法研究 被引量:3
9
作者 黄林显 邵景力 +2 位作者 崔亚莉 程汤培 李玲 《现代地质》 CAS CSCD 北大核心 2012年第3期614-619,共6页
大区域研究区由于涉及范围大、水文地质参数复杂多变,一直是进行地下水数值模拟的热点和难点。针对大区域地下水模拟的特点,在MPI环境中对Krylov子空间GMRES(m)算法的并行性进行分析,提出基于区域分解法的并行实现策略,并对不同的预条... 大区域研究区由于涉及范围大、水文地质参数复杂多变,一直是进行地下水数值模拟的热点和难点。针对大区域地下水模拟的特点,在MPI环境中对Krylov子空间GMRES(m)算法的并行性进行分析,提出基于区域分解法的并行实现策略,并对不同的预条件子的加速效果进行比较。数值实验结果表明:并行GMRES(m)算法在求解大区域三维地下水模型时可以显著的加快求解速度,且具有较好的可扩展性。另外,Jacobi预条件子与GMRES算法的组合具有更优的加速比和执行效率,是一种求解大型化、复杂化地下水水流问题的可行方案。 展开更多
关键词 大区域地下水模拟 区域分解法 并行gmres(m)算法 预条件技术
下载PDF
GMRES(m)算法在离散不适定问题中的应用 被引量:3
10
作者 张海燕 闵涛 刘相国 《科技导报》 CAS CSCD 2007年第13期54-59,共6页
基于投影方法的规划算法——Krylov子空间技术,研究了离散不适定正则化和Krylov子空间广义极小残余算法(GMRES(m))的基本理论,特别是残余向量与Krylov子空间的关系。利用离散不适定正则化方法,将不适定问题转化为适定问题,利用广义极小... 基于投影方法的规划算法——Krylov子空间技术,研究了离散不适定正则化和Krylov子空间广义极小残余算法(GMRES(m))的基本理论,特别是残余向量与Krylov子空间的关系。利用离散不适定正则化方法,将不适定问题转化为适定问题,利用广义极小残余算法对此适定问题进行数值求解。数值结果表明该算法是可靠和有效的。 展开更多
关键词 gmres(m)算法 不适定 ARNOLDI 正则化
下载PDF
集群系统中基于MPI的并行GMRES(m)计算通信的研究及应用 被引量:2
11
作者 杨爱民 刘韧 +1 位作者 赵广华 崔玉环 《微电子学与计算机》 CSCD 北大核心 2009年第9期129-131,135,共4页
针对求解大型稠密线性方程组的GMRES(m)算法的内在并行性,应用可移植消息传递标准MPI的集群通信机制在分布式存储并行系统上,设计了一种粗粒度、低通信开销的并行算法,并且应用于边界元求解的大型弹性问题的计算中.通过与串行算法进行比... 针对求解大型稠密线性方程组的GMRES(m)算法的内在并行性,应用可移植消息传递标准MPI的集群通信机制在分布式存储并行系统上,设计了一种粗粒度、低通信开销的并行算法,并且应用于边界元求解的大型弹性问题的计算中.通过与串行算法进行比较,设计的并行算法具有较高的计算精度和计算效率. 展开更多
关键词 边界元 gmres(m)算法 并行算法 通信开销
下载PDF
基于Jacobian-Free Newton-GMRES(m)方法的电力系统分布式暂态仿真算法 被引量:18
12
作者 陈颖 沈沉 +1 位作者 梅生伟 卢强 《电力系统自动化》 EI CSCD 北大核心 2006年第10期12-16,共5页
分布式暂态仿真是实现市场环境下互联电力系统在线一体化仿真分析的有效途径。文中研究了电力系统分布式暂态仿真计算模型,提出了基于Jacobian-Free Newton-GMRES(m)方法的协调求解算法。该算法只需要交换边界母线状态信息,接口简单,实... 分布式暂态仿真是实现市场环境下互联电力系统在线一体化仿真分析的有效途径。文中研究了电力系统分布式暂态仿真计算模型,提出了基于Jacobian-Free Newton-GMRES(m)方法的协调求解算法。该算法只需要交换边界母线状态信息,接口简单,实用性强。为了提高算法收敛速度,减少协调求解所需通信次数,提出了自适应预处理和连续修正预处理矩阵、预估边界条件初值及多时步同时协调等改进方法并将其应用于新算法中。测试结果表明,新的暂态仿真分解协调算法收敛快,通信次数少,非常适合在基于广域网络的分布式环境中实现。 展开更多
关键词 分布式计算 暂态仿真 电力网格 Jacobian—Free Newton-gmres(m)
下载PDF
变电站关键设备工频电场计算的预条件处理GMRES(m)边界元法 被引量:5
13
作者 张占龙 邓军 +4 位作者 许焱 何为 毛玉星 肖冬萍 韦军 《重庆大学学报(自然科学版)》 EI CAS CSCD 北大核心 2010年第1期78-82,93,共6页
在计算大尺度变电站关键设备工频电场时,传统方法效率低、性能差,计算困难。针对常规方法在大尺度工频电场计算中的瓶颈问题,提出了一种提高变电站关键设备三维电场分布计算效率的预条件GMRES(m)边界元法。阐述了预条件GMRES(m)迭代边... 在计算大尺度变电站关键设备工频电场时,传统方法效率低、性能差,计算困难。针对常规方法在大尺度工频电场计算中的瓶颈问题,提出了一种提高变电站关键设备三维电场分布计算效率的预条件GMRES(m)边界元法。阐述了预条件GMRES(m)迭代边界元法的基本原理及实现方法,并针对500kV变电站中部分关键设备周围电场分布进行了计算与比较分析。结果表明,预条件GMRES(m)边界元法经过预条件处理电位系数矩阵后,收敛速度快、残值收敛速度快、迭代次数少;在不降低计算精度的前提下,计算时间明显优越于直接迭代法;在满足工程误差和提高计算效率的同时,预条件GMRES(m)边界元法更适合于计算大尺度变电站关键设备的工频电场。 展开更多
关键词 边界元法 变电站 工频电场 广义极小残值算法(gmres) 预条件
下载PDF
基于GMRES(m)法的双连通区域数值保角变换的计算法 被引量:5
14
作者 吕毅斌 赖富明 +1 位作者 王樱子 武德安 《数学杂志》 CSCD 北大核心 2016年第5期1028-1034,共7页
本文研究了基于模拟电荷法的双连通区域的数值保角变换问题.利用限制Krylov子空间最大维数的算法–GMRES(m)算法,求解基于模拟电荷法的双连通区域数值保角变换中的约束方程,获得了模拟电荷和变换半径,构造了近似保角变换函数.数值实验... 本文研究了基于模拟电荷法的双连通区域的数值保角变换问题.利用限制Krylov子空间最大维数的算法–GMRES(m)算法,求解基于模拟电荷法的双连通区域数值保角变换中的约束方程,获得了模拟电荷和变换半径,构造了近似保角变换函数.数值实验表明了本文算法的有效性. 展开更多
关键词 模拟电荷法 双连通区域 KRYLOV子空间 gmres(m)法
下载PDF
基于预条件处理GMRES的不精确牛顿法潮流计算 被引量:16
15
作者 胡博 周家启 +1 位作者 刘洋 陈炜骏 《电工技术学报》 EI CSCD 北大核心 2007年第2期98-104,共7页
结合大规模电力系统修正方程组高维超稀疏性以及短向量的特点,提出以Krylov子空间方法研究电力系统方程计算问题。针对牛顿法潮流计算,采用预条件处理的GMRES方法求解高维稀疏的修正方程组,提出一种完整的基于预条件处理GMRES的不精确... 结合大规模电力系统修正方程组高维超稀疏性以及短向量的特点,提出以Krylov子空间方法研究电力系统方程计算问题。针对牛顿法潮流计算,采用预条件处理的GMRES方法求解高维稀疏的修正方程组,提出一种完整的基于预条件处理GMRES的不精确牛顿潮流算法,设计实现不同的预条件子,并以此为基础详细比较各类预条件子的预处理效果。通过对IEEE30、IEEE118和多个合成的大规模电力系统进行潮流计算,结果表明ILU预条件子比其他预条件子需要更少的迭代次数和浮点运算次数,当系统规模达到3000节点左右时,基于ILU预条件子的不精确牛顿法与传统的LU直接分解法相比,浮点运算次数减少了50%,内存使用量减少了将近10%,并且随着系统规模的增大,浮点运算次数基本上保持在LU直接法的50%左右,对大规模电力系统的潮流计算极为有利。 展开更多
关键词 潮流计算 KRYLOV子空间方法 不精确牛顿法 gmres 预条件处理
下载PDF
基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用 被引量:5
16
作者 龚小权 贾洪印 +2 位作者 陈江涛 赵辉 周桂宇 《空气动力学学报》 CSCD 北大核心 2019年第1期121-132,共12页
为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进... 为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进一步提高GMRES的计算效率,发展了方程组右端项残值雅可比精确计算方法,针对无黏通量Roe格式和黏性通量BR2(Bassi Rebay 2)黏性计算方法,分别解析给出其对守恒变量多项式自由度的雅可比矩阵。基于建立的方法首先采用NACA0012翼型研究了GMRES的重启次数及收敛参数对方法收敛性影响,然后采用无黏及黏性算例对比研究了基于雅可比矩阵不同计算方法的GMRES计算效率,同时对比研究了雅可比矩阵完全近似求解下GMRES和LU-SGS(Lower Upper-Symmetric Gauss-Seidel)的计算效率。结果表明,建立的基于右端项残值雅可比矩阵精确求解的GMRES方法能够大幅提高不同精度DG方法的CFL(CourantFriedrichs-Lewy)数,相比前面提到的其它方法具有更高的计算效率,其收敛速度实现量级以上的提高。 展开更多
关键词 隐式迭代 间断Galerkin有限元方法 gmres 精确雅可比矩阵 计算效率
下载PDF
基于Beowulf集群的大规模电力系统牛顿法潮流求解的并行GMRES方法 被引量:12
17
作者 胡博 谢开贵 曹侃 《电工技术学报》 EI CSCD 北大核心 2011年第4期145-152,共8页
大规模电力系统牛顿法潮流计算中,修正方程组的系数矩阵具有高维、稀疏、非对称的特点,结合该特点,提出基于预条件GMRES的并行牛顿法潮流计算方法。其中,对块Jacobi预条件子矩阵而言,根据处理器数确定其分块数,依此设计出高效的准对角... 大规模电力系统牛顿法潮流计算中,修正方程组的系数矩阵具有高维、稀疏、非对称的特点,结合该特点,提出基于预条件GMRES的并行牛顿法潮流计算方法。其中,对块Jacobi预条件子矩阵而言,根据处理器数确定其分块数,依此设计出高效的准对角并行预条件子矩阵;通过对Jacobi矩阵更新过程的矢量化处理,结合并行稀疏矩阵向量运算技术,提出Jacobi矩阵更新的并行化计算方法。对7 680节点、12 000节点等多个大规模电力系统进行潮流计算,结果表明:随着系统规模的增大(达到3 000节点及以上时),本文提出的并行潮流计算方法比传统并行LU分解法在并行加速比、并行效率等方面有明显优势。 展开更多
关键词 潮流计算 准对角预条件子矩阵 Jacobi矩阵更新 并行计算 gmres方法
下载PDF
误差向量与Krylov子空间对GMRES(m)算法收敛速度的影响 被引量:3
18
作者 吴果林 王晟 《广西科学》 CAS 2011年第3期214-217,221,共5页
从广义极小残量法GMRES(m)的结构出发,分析其误差向量与Krylov子空间对该算法收敛速度的影响,推导出误差向量与Krylov子空间第1个向量和第m+1个向量的方向余弦关系,并用数值算例验证其合理性.当误差向量rk+1在Krylov子空间向量v1的投影... 从广义极小残量法GMRES(m)的结构出发,分析其误差向量与Krylov子空间对该算法收敛速度的影响,推导出误差向量与Krylov子空间第1个向量和第m+1个向量的方向余弦关系,并用数值算例验证其合理性.当误差向量rk+1在Krylov子空间向量v1的投影较大而在向量vm+1的投影较小时,GMRES(m)算法收敛速度较慢,反之亦反.算例结果与理论结果相符. 展开更多
关键词 线性方程 迭代方法 广义极小残量法 KRYLOV子空间
下载PDF
基于辛Gauss方法及预处理GMRES方法的暂态稳定性并行计算 被引量:3
19
作者 温柏坚 胡佳怡 +2 位作者 郭文鑫 汪芳宗 李钦 《电力系统保护与控制》 EI CSCD 北大核心 2012年第22期19-24,共6页
将s级2s阶的辛Gauss方法用于电力系统暂态稳定性计算,提出了一种新的并行计算方法。该算法首先将微分—代数方程组经多级差分后转化为大规模非线性方程组,并利用牛顿法对其进行求解。在此基础上,利用矩阵分解方法将整体计算任务分解为... 将s级2s阶的辛Gauss方法用于电力系统暂态稳定性计算,提出了一种新的并行计算方法。该算法首先将微分—代数方程组经多级差分后转化为大规模非线性方程组,并利用牛顿法对其进行求解。在此基础上,利用矩阵分解方法将整体计算任务分解为两部分:一部分计算任务可按相应的级数或在不同的时间点上进行'解耦',因而具有完全的时间并行性;对剩下的一部分计算任务,采用预处理GMRES方法对其进行空间并行求解,并为此提出了一种新的预处理方法。利用三个不同规模的算例系统,对所提算法的收敛性进行了测试,并在GPU上对算法进行了实际测试。测试结果表明,该算法可以获得很高的加速比,可以用于大规模电网暂态稳定性的实时分析计算。 展开更多
关键词 暂态稳定性 辛Gauss算法 并行计算 gmres方法 W-变换 预处理 GPU
下载PDF
预条件GMRES(m)法迭代求解大规模边界元弹性问题 被引量:5
20
作者 陈泽军 肖宏 《固体力学学报》 CAS CSCD 北大核心 2006年第S1期50-55,共6页
提出了一种基于单元节点的块雅可比预条件方法,扩大了边界元法的计算规模,使之可用于大规模工程问题的求解.数值实验说明了这种预条件技术的有效性,表明预条件GMRES(m)算法具有较好的收敛特性,适合于求解大规模问题边界元弹性问题所形... 提出了一种基于单元节点的块雅可比预条件方法,扩大了边界元法的计算规模,使之可用于大规模工程问题的求解.数值实验说明了这种预条件技术的有效性,表明预条件GMRES(m)算法具有较好的收敛特性,适合于求解大规模问题边界元弹性问题所形成的稠密非对称线性方程组. 展开更多
关键词 广义极小残值算法(gmres) 预条件 边界元法 弹性问题 迭代法
下载PDF
上一页 1 2 86 下一页 到第
使用帮助 返回顶部