辛几何算法在射线追踪中的应用

被引：12

作者：

陈景波

秦孟兆

机构：

[1] 中国科学院数学与系统科学研究院计算数学与科学工程计算研究所

来源：

数值计算与计算机应用 | 2000年 / 04期

关键词：

Ray tracing; Symplectic algorithm; Maslov theory;

D O I：

暂无

中图分类号：

O186 [微分几何、积分几何];

学科分类号：

摘要：

Ray tracing is a basic aspect in tomography. To solve the caustic problem in inhomogeneous media using Maslov asymptotic theory, we need to calculate the position and slowness vector at every point. Therefore, ray tracing must rely on the ray equations in Hamiltonian form. In this paper, fourth order symplectic scheme and nonsymplectic Runge-Kutta scheme are compared in ray tracing for sinusoidal velocity model. The result indicates that ray paths obtained by two schemes are almost the same. But on keeping Hamilton quantities, the symplectic scheme is far better than the Runge-Kutta scheme. On computing travel time for Htamiltonian system with T parameter, we use trapezoid formula for numerical integration. The result coincides with that obtained using Hamiltonian system with t parameter.

引用

页码：255 / 265

页数：11