模板函数返回问题

时间:2013-02-21 14:36:51

标签: c++ templates return gsl

我是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

似乎初始数组没有更新,有人帮助我吗?谢谢!

1 个答案:

答案 0 :(得分:2)

这与模板无关。具有30个元素的数组在indizes 0..29处具有元素。您在r[30]中存储了一些不受限制的内容。

此外,您正在创建非对象长度(p1,p3和半径)的数组,这是无效的C ++。从你的返回我假设你想要返回整个数组,包含30个值。你要做的是返回第31个(!)值,或者至少你尝试访问result[n]。这是未定义的行为,这意味着任何事情都可能发生:没什么,崩溃,你可能期望的任何东西,你的地址的在线披萨订单......

看起来您仍然需要对C ++数组以及如何处理它们有一些基本的了解。

<强> PS: 作为第一个提示,std::vector可能是您正在寻找的,因为它是处理可变长度数组的C ++方法。

<强> PS2: 不要将#define用于常量。改为使用固定类型的常量变量。它永远不会受到伤害,但有时会为您节省一些痛苦的调试。