2

我有一些代码可以在单位球面上找到一个点。回想一下，对于一个单位球体：

1 = sqrt( x^2 + y^2 + z^2 )

该算法在零和一之间选择两个随机点（x 和 y 坐标）。如果它们的大小小于 1，我们就有空间通过求解上述 z 方程来定义第三个坐标。

void pointOnSphere(double *point){
double x, y;

do {
    x = 2*randf() - 1;
    y = 2*randf() - 1;
} while (x*x + y*y > 1);

double mag = sqrt(fabs(1 - x*x - y*y));

point[0] = 2*(x*mag);
point[1] = 2*(y*mag);
point[2] = 1 - 2*(mag*mag);
}

从技术上讲，我继承了这段代码。前任所有者使用“无视严格的标准合规性”的 -Ofast 编译。TL;DR 这意味着您的代码不需要遵循严格的 IEEE 标准。因此，当我尝试在没有优化的情况下进行编译时，我遇到了错误。

 undefined reference to `sqrt'

什么是IEEE标准？好吧，因为计算机不能以无限精度存储浮点数，如果你不小心，在某些计算过程中会出现舍入错误。

经过一番谷歌搜索后，我遇到了这个问题，这让我走上了正确使用 IEEE 材料的正确轨道。我什至读过这篇关于浮点数的文章（我推荐）。不幸的是，它没有回答我的问题。

我想在我的函数中使用 sqrt() 而不是像Newton Iteration这样的东西。我知道我的算法中的问题可能来自我可能（即使不是真的）将负数传递给 sqrt() 函数的事实。我只是不太确定如何解决这个问题。感谢所有的帮助！

哦，如果相关的话，我正在使用Mersenne Twister数字生成器。

澄清一下，我将 libm 与 -lm 链接起来！我还确认它指向正确的库。

4

3 回答 3

1

至于未定义的引用，sqrt您需要链接 with libm，通常是 with-lm或类似选项。

另请注意

如果它们的大小小于 1，我们就有空间通过求解上述 z 方程来定义第三个坐标。

是错的。x 和 y 必须满足 x * x + y * y <= 1 才能有 z 的解。

于 2017-11-01T21:16:02.707 回答
0

为了确保点满足条件，测试条件本身作为while循环的一部分，而不是条件的推导。

// functions like `sqrt(), hypot()` benefit with declaration before use
//   and without it may generate "undefined reference to `sqrt'"
// Some functions like `sqrt()` are understood and optimized out by a smart compiler.
// Still, best to always declare them.
#include <math.h>

void pointOnSphere(double *point){
  double x, y, z;
  do {
    x = 2*randf() - 1;
    y = 2*randf() - 1;
    double zz = 1.0 - hypot(x,y); 
    if (zz < 0.) continue;  // On rare negative values due to imprecision
    z = sqrt(zz);
    if (rand()%2) z = -z;  // Flip z half the time
  } while (x*x + y*y + z*z > 1);  // Must meet this condition

  point[0] = x;
  point[1] = y;
  point[2] = z;
}
于 2017-11-01T21:27:20.717 回答
0

我会使用球坐标

theta = randf()*M_PI;
phi = randf()*2*M_PI;
r = 1.0;
x = r*sin(theta)*cos(phi);
y = r*sin(theta)*sin(phi);
z = r*cos(theta);
于 2017-11-01T21:30:46.230 回答