我該如何引用std :: sin（const valarray <double>＆）？

Question

我有一些valarray函數指針代碼的麻煩：

double (*fp)(double) = sin; 
valarray<double> (*fp)(const valarray<double> &) = sin; 

第一個彙總，第二給出：

error: no matches converting function 'sin' to type 'class std::valarray<double> (*)(const class std::valarray<double>&)' 

Answer 1

這將編譯。看起來GCC的valarray使用表達式模板來延遲竇的計算。但是這會使sin模板的返回類型不完全是valarray<T>，而是一些奇怪的複雜類型。

#include <valarray> 

template<typename T> struct id { typedef T type; }; 
int main() { 
    using std::valarray; 
    using std::sin; 

    id<__typeof__(sin(valarray<double>()))>::type (*fp)(const valarray<double> &) = sin; 
} 

編輯：見AProgrammer的標準報價爲什麼GCC是罰款這樣做。

編輯：符合標準的解決方法

在嚴格符合標準的方式這樣做沒有__typeof__是有點棘手。您將需要獲得返回類型sin。您可以像這樣使用條件運算符，如Eric Niebler has shown。它的作用是讓sin函數實際上不被調用，但只能進行類型檢查。通過嘗試另一個分支條件運算符來同類型的（這實際上是評估的一個）轉換，我們可以生成一個虛擬參數只是爲了能夠推斷出函數指針的類型：

#include <valarray> 

using std::valarray; 

template<typename T> struct id { 
    typedef T type; 
}; 

struct ded_ty { 
    template<typename T> 
    operator id<T>() { return id<T>(); } 
}; 

template<typename E, typename T> 
id<T(*)(valarray<E> const&)> genFTy(T t) { 
    return id<T(*)(valarray<E> const&)>(); 
} 

template<typename T> 
void work(T fp, id<T>) { 
    // T is the function pointer type, fp points 
    // to the math function. 
} 

int main() { 
    work(std::sin, 1 ? ded_ty() : genFTy<double>(std::sin(valarray<double>()))); 
} 

如果您想立即獲取地址，則可以編寫work，以便再次返回fp。現在

template<typename T> 
T addy(T fp, id<T>) { return fp; } 

，你終於可以編寫一個宏來封裝有條件的經營者弄虛作假，當你想要得到任何這樣的數學函數的地址，使用它。

#define DEDUCE(FN,Y) (1 ? ded_ty() : genFTy<Y>(FN(std::valarray<Y>()))) 

要獲取地址，並將它傳遞給一些通用的功能，下面的作品，然後

std::transform(v1.begin(), v1.end(), v1.begin(), 
    addy(std::sin, DEDUCE(std::sin, double))); 
std::transform(v2.begin(), v2.end(), v2.begin(), 
    addy(std::cos, DEDUCE(std::cos, double))); 

Answer 2

你說的稱號約std::sin，但隨後分配::sin。

valarray<double> (*fp)(const valarray<double> &) = std::sin; 

這應該有效。請注意，您應該限定sin的所有用途，儘管大多數實現將名稱注入到全局名稱空間，即使包含<cmath>（這是非標準行爲）。

編輯：不幸的是，你運氣不好。該標準說約sin(valarray<T> const &)以下（26.3.3.3）。

該函數將返回一個值，該值是類型T的或它可以是明確地 轉換爲類型T.與gcc執行

優化由標準理所當然的。上面的代碼不能保證工作。

Answer 3

26 3.1/3

任何函數返回的valarray被允許返回另一個類型的對象， 提供的所有const成員的valarray的功能也適用於這種類型。

的目的是允許使用模板表達式來優化結果（即循環整個陣列做各次的計算在一個時間，直接分配到所得的valarray <>而不是建立一個臨時的）。

z = sin(x+y); 

可以使用__typeof__ GCC擴展進行優化，以

for (i = 0; i < N; ++i) 
    z[i] = sin(x[i] + y[i]);