- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
-6ren">
我正在尝试用 C++ 实现 FFT。但是我收到此错误:
Undefined symbols for architecture x86_64:
"Complex::fft(std::vector<Complex*, std::allocator<Complex*> >*)", referenced from:
_main in ccMdeaRo.o
ld: symbol(s) not found for architecture x86_64
collect2: ld returned 1 exit status
这是我的代码:
// Complex.h -----------------------------------------------------------------
#include <cmath>
#include <string>
#include <sstream>
#include <vector>
using namespace std;
class Complex {
float _re,_im;
public:
static Complex ZERO;
static Complex ONE;
static const float TOLERANCE = 5E-6f;
Complex(float re, float im);
// Polar(float r, float theta);
// ~Complex();
float angle();
float magnitude();
bool equals(Complex* arg);
Complex* clone();
static vector<Complex*>* fft(vector<Complex*>* a);
// static Complex** fft(Complex** a, long n);
static Complex* omega(long n);
static Complex* polar(float r, float theta);
float re();
float im();
Complex* minus(Complex* that);
Complex* plus(Complex* that);
Complex* times(Complex* that);
string toString();
};
// Complex.cpp -----------------------------------------------------------------
#include "Complex.h"
Complex::Complex(float re, float im) {
this->_re = re;
this->_im = im;
}
float Complex::angle(){
return atan2(_im,_re);
}
float Complex::magnitude() {
return sqrt(this->_re * _re + this->_im * _im);
}
Complex* Complex::omega(long n) {
Complex* cplx = polar(1.0, 2*(M_PI)/n);
return cplx;
}
vector<Complex*>* fft(vector<Complex*>* a) {
long n = a->size();
if (n == 1) {
return a;
}
Complex* w_n = Complex::omega(n);
Complex* w = new Complex(1,0);
long half_n = n/2;
vector<Complex*>* aEvens = new vector<Complex*>(half_n);
vector<Complex*>* aOdds = new vector<Complex*>(half_n);
for (long i = 0; i < n; i++) {
if (i % 2 == 0) {
aEvens->push_back(a->at(i));
}
else {
aOdds->push_back(a->at(i));
}
}
vector<Complex*>* yEvens = fft(aEvens);
vector<Complex*>* yOdds = fft(aOdds);
vector<Complex*>* y = new vector<Complex*>(n);
y->resize(n, 0);
for (long k = 0; k < n/2; k++) {
Complex* prod = (yOdds->at(k))->times(w);
y->at(k) = yEvens->at(k)->plus(prod);
y->at(k+half_n) = yEvens->at(k + half_n)->minus(prod);
w = w->times(w_n);
}
delete w;
delete w_n;
return y;
}
Complex* Complex::polar (float r, float theta) {
Complex* cplx = new Complex(r*cos(theta),r*sin(theta));
return cplx;
}
bool Complex::equals(Complex* arg) {
double d_re = abs(this->_re - arg->_re);
double d_im = abs(this->_im - arg->_im);
return (d_re < TOLERANCE && d_im < TOLERANCE);
}
Complex* Complex::clone() {
Complex* theClone = new Complex(this->_re,this->_im);
return theClone;
}
float Complex::re() {
return this->_re;
}
float Complex::im() {
return _im;
}
Complex* Complex::minus(Complex* arg) {
float diff_re = this->_re - arg->_re;
float diff_im = this->_im - arg->_im;
Complex* diff = new Complex(diff_re, diff_im);
return diff;
}
Complex* Complex::plus(Complex* arg) {
float sum_re = this->_re + arg->_re;
float sum_im = this->_im + arg->_im;
Complex* sum = new Complex(sum_re, sum_im);
return sum;
}
Complex* Complex::times(Complex* arg) {
float magn = this->magnitude() * arg->magnitude();
float angl = this->angle() + arg->angle();
Complex* prod = polar(magn, angl);
return prod;
}
string Complex::toString() {
ostringstream s;
s << this->_re << " + " << this->_im << "i";
return s.str();
}
// main -----------------------------------------------------------------
#include "complex.h"
#include <iostream>
int main() {
Complex* a = new Complex(1,0);
Complex* b = new Complex(0,1);
Complex* c = new Complex(5,5);
Complex* d = new Complex(6,2);
Complex* e = a -> clone();
Complex* f = a->plus(b);
Complex* g = a->minus(b);
Complex* h = c->times(d);
Complex* j = Complex::omega(4);
Complex* k = Complex::polar(1, 2*(M_PI));
cout << "a:\t" << a->toString() << endl;
cout << "b:\t" << b->toString() << endl;
cout << "a->re():\t" << a->re() << endl;
cout << "a->im():\t" << a->im() << endl;
cout << "a plus b:\t" << f->toString() << endl;
cout << "a minus b:\t" << g->toString() << endl;
cout << "c times d:\t" << h->toString() << endl;
cout << "omega(4):\t" << j ->toString() << endl;
cout << "j->equals(b):\t" << j->equals(b) << endl;
cout << "j->equals(c):\t" << j->equals(c) << endl;
cout << "polar(1, 2PI):\t" << k->toString() << endl;
cout << "k->equals(a):\t" << k->equals(a) << endl;
cout << "k->equals(b):\t" << k->equals(b) << endl;
cout << "a->clone():\t" << e->toString() << endl;
Complex* one = new Complex(1,0);
Complex* two = new Complex(2,0);
Complex* three = new Complex(3,0);
Complex* four = new Complex(4,0);
vector<Complex*>* cplxVec = new vector<Complex *>;
cplxVec->push_back(one);
cplxVec->push_back(two);
cplxVec->push_back(three);
cplxVec->push_back(four);
cout << "\nPrinting out vector:\n" << endl;
for (unsigned int i = 0; i < cplxVec->size(); i++) {
cout << cplxVec->at(i) -> toString() << endl;
}
cout << endl;
Complex::fft(cplxVec);
return 0;
}
最佳答案
当您实现 fft
时,您忘记了 Complex::
:
vector<Complex*>* fft(vector<Complex*>* a) {
应该是:
vector<Complex*>* Complex::fft(vector<Complex*>* a) {
fft
有一个问题,它在调用自身时会在这里陷入无限循环:
vector<Complex*>* yEvens = fft(aEvens);
你减半一定不能让它们的大小减少到 1
。
关于C++ FFT 实现(体系结构 x86_64 的 undefined symbol ),我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/16826547/
问题故障解决记录 -- Java RMI Connection refused to host: x.x.x.x .... 在学习JavaRMI时,我遇到了以下情况 问题原因:可
我正在玩 Rank-N-type 并尝试输入 x x .但我发现这两个函数可以以相同的方式输入,这很不直观。 f :: (forall a b. a -> b) -> c f x = x x g ::
这个问题已经有答案了: How do you compare two version Strings in Java? (31 个回答) 已关闭 8 年前。 有谁知道如何在Java中比较两个版本字符串
这个问题已经有答案了: How do the post increment (i++) and pre increment (++i) operators work in Java? (14 个回答)
下面是带有 -n 和 -r 选项的 netstat 命令的输出,其中目标字段显示压缩地址 (127.1/16)。我想知道 netstat 命令是否有任何方法或选项可以显示整个目标 IP (127.1.
我知道要证明 : (¬ ∀ x, p x) → (∃ x, ¬ p x) 证明是: theorem : (¬ ∀ x, p x) → (∃ x, ¬ p x) := begin intro n
x * x 如何通过将其存储在“auto 变量”中来更改?我认为它应该仍然是相同的,并且我的测试表明类型、大小和值显然都是相同的。 但即使 x * x == (xx = x * x) 也是错误的。什么
假设,我们这样表达: someIQueryable.Where(x => x.SomeBoolProperty) someIQueryable.Where(x => !x.SomeBoolProper
我有一个字符串 1234X5678 我使用这个正则表达式来匹配模式 .X|..X|X. 我得到了 34X 问题是为什么我没有得到 4X 或 X5? 为什么正则表达式选择执行第二种模式? 最佳答案 这里
我的一个 friend 在面试时遇到了这个问题 找到使该函数返回真值的 x 值 function f(x) { return (x++ !== x) && (x++ === x); } 面试官
这个问题在这里已经有了答案: 10年前关闭。 Possible Duplicate: Isn't it easier to work with foo when it is represented b
我是 android 的新手,我一直在练习开发一个针对 2.2 版本的应用程序,我需要帮助了解如何将我的应用程序扩展到其他版本,即 1.x、2.3.x、3 .x 和 4.x.x,以及一些针对屏幕分辨率
为什么案例 1 给我们 :error: TypeError: x is undefined on line... //case 1 var x; x.push(x); console.log(x);
代码优先: # CASE 01 def test1(x): x += x print x l = [100] test1(l) print l CASE01 输出: [100, 100
我正在努力温习我的大计算。如果我有将所有项目移至 'i' 2 个空格右侧的函数,我有一个如下所示的公式: (n -1) + (n - 2) + (n - 3) ... (n - n) 第一次迭代我必须
给定 IP 字符串(如 x.x.x.x/x),我如何或将如何计算 IP 的范围最常见的情况可能是 198.162.1.1/24但可以是任何东西,因为法律允许的任何东西。 我要带198.162.1.1/
在我作为初学者努力编写干净的 Javascript 代码时,我最近阅读了 this article当我偶然发现这一段时,关于 JavaScript 中的命名空间: The code at the ve
我正在编写一个脚本,我希望避免污染 DOM 的其余部分,它将是一个用于收集一些基本访问者分析数据的第 3 方脚本。 我通常使用以下内容创建一个伪“命名空间”: var x = x || {}; 我正在
我尝试运行我的test_container_services.py套件,但遇到了以下问题: docker.errors.APIError:500服务器错误:内部服务器错误(“ b'{” message
是否存在这两个 if 语句会产生不同结果的情况? if(x as X != null) { // Do something } if(x is X) { // Do something } 编
我是一名优秀的程序员,十分优秀!