小编Tim*_*Tim的帖子

我正在寻找一些帮助或提示来加速我的代码。

我已经实现了一个例程，用于根据一组球谐系数 C_{n,m} 和 S_{n,m} 计算空间中某个点 (r,phi,lambda) 的引力势。等式显示在下面的链接中：

并且包括与纬度 (phi) 相关的关联勒让德多项式 P_{n,m} 的递归计算，从前两个值 P{0,0} 和 P_{1,1} 开始。

起初，我将其作为 MATLAB C-MEX 代码实现，只有我代码的核心部分是用 C 语言编写的。我想做一个纯C程序，但发现代码运行速度慢了3-5倍，这让我想知道为什么。这可能是我定义结构并使用指向中央代码中的指针的方式吗？看起来它是需要额外时间的核心计算部分，但该部分没有改变，尽管在我将指针直接传递给变量之前，现在我在结构内部使用指针。任何帮助表示赞赏！

在下文中，我将尝试解释我的代码并展示一些摘录：

在程序开始时，我定义了三个结构。一个保存球谐系数，C_{n,m} 和 S_{n,m}, (ggm_struct)，一个保存计算坐标 (comp_struct)，一个保存结果 (func_struct)：

// Define constant variables
const double deg2rad = M_PI/180.0;                 // degrees to radians conversion factor
const double sfac = 1.0000E-280;                   // scaling factor
const double sqrt2 = 1.414213562373095;            // sqrt(2)
const double sqrt3 = 1.732050807568877;            // sqrt(3)

// Define structure to hold geopotential model
struct ggm_struct {
    char   product_type[100], modelname[100], …