我是一個C++的新人,我想做一個模板函數。我遇到了模板返回值的問題。的代碼是這樣的模板函數返回問題
#include <iostream>
#include <cmath>
#include <gsl/gsl_rng.h>
#define pi 3.1415926535
using namespace std;
template <class T>
T cpradius(T a,T b,T p, int n)
{
const gsl_rng_type *R;
gsl_rng *r;
gsl_rng_env_setup();
R = gsl_rng_default;
r = gsl_rng_alloc(R);
T p1[n],p3[n],p2 ;
T radius[n] ;
for (int i = 0; i<n; i++)
{
p1[i] = gsl_rng_uniform(r);
p2 += p1[i];
}
for (int j = 0; j<n; j++)
{
p3[j] = p1[j]/p2;
radius[j] = sqrt(p3[j]*a*b*p/pi);
//cout << radius[j] << endl;
}
return radius[n];
}
int main(){
double r[30] = {0};
r[30] = cpradius(30.0,30.0,0.6,30);
for (int i = 0;i<30;i++){
cout << r[i] <<endl;
}
return 0;
}
然後,編譯:
g++ -Wall -I/usr/local/include/ tst3.cpp -lgsl -lgslcblas
結果被1點式柱,30零:
0
0
0
...
0
似乎最初的陣列沒有更新,是否有人幫我?謝謝!
有效數組索引爲'[0,n-1]',因此'radius [n]'和'r [30]'調用未定義的行爲。 – 2013-02-21 14:41:45