在表中预先计算一次cos()和sin()

最简单的解决方案是编写一个单独的程序,创建 .cc 带有矢量定义的文件。

例如:

#include <iostream>
#include <cmath>

int main()
{
    std::ofstream out("values.cc");

    out << "#include \"static_values.h\"\n"; 
    out << "#include <vector>\n";

    out << "std::vector<float> pi_values = {\n";
    out << std::precision(10);

    // We only need to compute the range from 0 to PI/2, and use trigonometric
    // transformations for values outside this range.
    double range = 3.141529 / 2;
    unsigned num_results = 250000;

    for (unsigned i = 0; i < num_results; i++) {
        double value = (range / num_results) * i;
        double res = std::sin(value);

        out << "    " << res << ",\n";
    }
    out << "};\n"
    out.close();
}

请注意,这不太可能提高性能,因为这种大小的表可能不适合您的二级缓存。这意味着很大一部分三角计算将需要访问RAM;每个这样的访问花费大约几百个CPU周期。

顺便说一句,你看过近似SSE SIMD三角函数库吗。这看起来是一个很好的用例。

您可以使用预计算,而不是将它们存储在可执行文件中:

double precomputed_sin[65536];

struct table_filler {
    table_filler() {
        for (int i=0; i<65536; i++) {
            precomputed_sin[i] = sin(i*2*3.141592654/65536);
        }
    }
} table_filler_instance;

这样,在程序启动时只计算一次表,并且它仍然位于固定的内存地址。之后 tsin 和 tcos 可以内联实现为

inline double tsin(int x) { return precomputed_sin[x & 65535]; }
inline double tcos(int x) { return precomputed_sin[(x + 16384) & 65535]; }

对于这类问题,通常的答案是写一个 生成C++源文件的程序,该文件中的值 一个表,并将其编译到您的DLL中。如果你正在考虑 具有128000个条目的表(128000个双是1MB), 编译器中可能会遇到一些内部限制。 在这种情况下,您可能会考虑将值写入 作为内存转储的文件,以及 mmap 加载时正在处理此文件 DLL。(在窗户下面,我想你甚至可以把这一秒放进去 文件转换为DLL文件的第二个流,这样您就不会 以分发第二个文件。)