• 本期要文(可全文下载) •

### 高精度曲面建模的一种快速算法

1. 中国科学院地理科学与资源研究所 资源与环境信息系统国家重点实验室,北京 100101
• 收稿日期:2011-11-07 修回日期:2012-04-18 出版日期:2012-06-25 发布日期:2012-06-25
• 通讯作者: 岳天祥(1963-),男,甘肃庆阳人,博士生导师,研究员。主要研究方向:资源环境模型与系统模拟。 E-mail:yue@lreis.ac.cn E-mail:yue@lreis.ac.cn
• 作者简介:赵 娜(1986-),女,山东莱芜人,博士研究生。主要研究方向:生态模型与系统模拟、气候变化的模拟与研究 。E-mail: zhaon@lreis.ac.cn
• 基金资助:

国家杰出青年科学基金项目(40825003);全球变化研究国家重大科学研究计划项目课题(2010CB950904);国家自然科学基金重点项目(41023010)。

### Fast Methods for High Accuracy Surface Moldeling

ZHAO Na, YUE Tianxiang

1. State Key Laboratory of Resources and Environment Information System,Institute of Geographic Sciences and Natural Resources Research, CAS, Beijing 100101, China
• Received:2011-11-07 Revised:2012-04-18 Online:2012-06-25 Published:2012-06-25

Abstract:

As a novel surface modeling method, the whole computing process of HASM can be divided into three parts: deriving finite difference approximations to differential equations, establishing the sampling point equations and solving the algebra equations. The well known conjugate gradient method has been used extensively to solve this symmetric positive definite system instead of Gaussian's elimination based direct methods, especially for very large systems when parallel solution environment is preferred. However, a difficulty associated with the method of conjugate gradients is that it works well on matrices that are either well conditioned or have just a few distinct eigenvalues and the coefficient matrix of the algebra equation in HASM is ill-conditioned. In this paper, we show how to preprocess a linear system so that the matrix of coefficients assumes one of these nice forms. We give two other preconditioners, i.e. incomplete Cholesky decomposition conjugate gradient method (ICCG) and symmetric successive over relaxation-preconditioned conjugate gradient method (SSORCG), so as to improve the convergence rate of HASM. Furthermore, we give adequate consideration in storage scheme of the large sparse matrix and optimize the performance of sparse matrix-vector multiplication. The cost of the computation is also considered in each iteration. We implement and test the proposed method on a Dell OptiPlex 990MT machine. Numerical tests show that ICCG has the fastest convergence rate of HASM. We also find that both ICCG and SSORCG have much faster convergence rates than others.