高手指路啊～～～～～～～～

[allendragon]

下面是一个模版函数
template <class T, int N>
static inline bool ldltdc(T **A, T *rdiag)
{
T v[N-1];
int xx = N; //here is the problem, why N = N+1????
for (int i = 0; i < N; i++) {
for (int k = 0; k < i; k++)
v[k] = A[i][k] * rdiag[k];
for (int j = i; j < N; j++) {
T sum = A[i][j];
for (int k = 0; k < i; k++)
sum -= v[k] * A[j][k];
if (i == j) {
if (unlikely(sum <= T(0)))
return false;
rdiag[i] = T(1) / sum;
} else {
A[j][i] = sum;
}
}
}

return true;
}

在mfc中调用出现下面的情况：
ldltdc<float,3>(w, diag)；
调试的时候模版函数中的N＝4而不是N＝3。
同样的调用在console程序中没有问题，N＝3。

这是为什么啊？

[allendragon]

又发现下面这个函数好使
// LU decomposition
template <class T, int N>
static inline bool ludcmp(T **a, int *indx, T *d = NULL)
{
int i, j, k;
T vv[N];
if (d)
*d = 1;
for (i = 0; i < N; i++)
。。。。。。
}
撞鬼了啊～～～～～

[huaren801007]

在学习mfc对com的支持的时候,记忆当中是mfc只要使用类厂对象创建一个com组件的实例,就会把它的引用计数设为1,然后客户端引用的时候,这时引用计数已经变为2了

以上两个函数的不同在于第一个函数用到了N作为右值,问题应该是出在这里

说的不对的话,见谅!

[allendragon]

int xx = N; //here is the problem, why N = N+1????
你说的是上面这句吗？这句是为了调试加上的，本来的函数中是没有的，即使没有这句仍然有不对啊。