假设你有一个std::complex<double> array[N];
您如何在 for 循环中将 all-inf
和values 设置为 0?inf
std::isinf
对我不起作用编译器早于 C++11
假设你有一个std::complex<double> array[N];
您如何在 for 循环中将 all-inf
和values 设置为 0?inf
std::isinf
对我不起作用编译器早于 C++11
试试这个:std::numeric_limits<double>::infinity()
将提供一个值,告诉您一个数字是否为无穷大。
这是一个值的示例vector
代码double
;您可以更改repl_inf
为使用complex<double>
。
#include <iostream>
#include <algorithm>
#include <vector>
#include <limits>
using namespace std;
int repl_inf (double v) { return (v== std::numeric_limits<double>::infinity() || -v== std::numeric_limits<double>::infinity()) ? 0.0 : v; }
int main() {
vector<double> v;
v.push_back(1);
v.push_back(1.0/0.0);
v.push_back(2);
v.push_back(-1.0/0.0);
transform (v.begin(), v.end(), v.begin(), repl_inf);
for (int i = 0 ; i != v.size() ; i++) {
cout << v[i] << endl;
}
return 0;
}
在 ideone:链接。
对于complex<double>
:
bool isInf(double v) {
return v== std::numeric_limits<double>::infinity()
|| -v== std::numeric_limits<double>::infinity();
}
bool isInf(complex<double> c) {
return isInf(c.real) || isInf(c.imag);
}
std::replace( array, array + N, std::complex( std::numeric_limits<double>::infinity(), 0.0), std::complex(0.0) );
std::replace( array, array + N, std::complex( -std::numeric_limits<double>::infinity(), 0.0), std::complex(0.0) );
或者
bool isinf( std::complex<double> d )
{
return std::abs( d.real() ) == std::numeric_limits<double>::infinity()
}
std::replace_if( array, array + N, isinf, 0.0 );