我是一名有用的助手,可以翻译文本。
我遇到了使用 malloc 的问题。
我有一个名为 jacobi_gpu
的函数需要被多次调用:
int main(int argc, char* argv[]){
/* ... */
int totalrot=0;
while(nrot>0){
iter++;
nrot=jacobi_gpu(a,q, tol, dimmat);
totalrot+=nrot;
printf("iter =%3d nrot=%3d\n",iter, nrot);
}
/* ... */
}
参数 a、q、tol 和 dimmat 已正确初始化。A 和 Q 是两个方阵,dimmat 是它们的维度。
这是我的代码:
int jacobi_gpu(double A[], double Q[], double tol, long int dim){
int nrot, p, q, k, tid;
double c, s;
double *mc, *vc;
printf("jacobi begins \n");
mc = (double *)malloc(2 * dim * sizeof(double));
vc = (double *)malloc(2 * dim * sizeof(double));
if( mc == NULL || vc == NULL){
fprintf(stderr, "pb allocation matricre\n");
exit(1);
}
nrot = 0;
for(k = 0; k < dim - 1; k++){
eye(mc, dim);
eye(vc, dim);
for(tid = 0; tid < floor(dim /2); tid++){
p = (tid + k)%(dim - 1);
if(tid != 0)
q = (dim - tid + k - 1)%(dim - 1);
else
q = dim - 1;
//printf("p = %d | q = %d\n", p, q);
if(fabs(A[p + q*dim]) > tol){
nrot++;
symschur2(A, dim, p, q, &c, &s);
mc[2*tid] = p; vc[2 * tid] = c;
mc[2*tid + 1] = q; vc[2*tid + 1] = -s;
mc[2*tid + 2*(dim - 2*tid) - 2] = p;
vc[2*tid + 2*(dim - 2*tid) - 2 ] = s;
mc[2*tid + 2*(dim - 2*tid) - 1] = q;
vc[2 * tid + 2*(dim - 2*tid) - 1 ] = c;
}
}
affiche(mc,dim,2,"Matrice creuse");
affiche(vc,dim,2,"Valeur creuse");
}
printf("end\n");
free(mc);
free(vc);
return nrot;
}
我的问题出现在对mc变量进行malloc调用时:
*** glibc detected *** ./jacobi_gpu: double free or corruption (!prev): 0x00000000022944a0 ***
*** glibc detected *** ./jacobi_gpu: malloc(): memory corruption: 0x0000000002294580 ***
有什么建议吗?
[编辑]
- 函数eye初始化一个单位矩阵。
- 函数affiche用行和列显示矩阵。第一个参数是矩阵,第二个参数是行数,第三个参数是列数。
更多解释
矩阵mc的目的是存储变量p和q。这些变量包含列索引。 矩阵vc的目的是存储这些列中包含的值。 例如,如果矩阵mc的第一行是0和5(p = 0,q = 5),那么意味着矩阵vc中的值将在列0和5中。 如果矩阵的第五行是2 3(p = 2,q = 3),那么意味着矩阵vc中第五行的值将在列2和3中。
希望这次我表述更清楚了。
谢谢你的帮助。
eye()
函数是做什么的?它在什么情况下会失败?即,它是在第一次调用malloc()
后失败还是在多次调用后失败? - e.Jamesaffiche()
函数是做什么用的? - e.Jamesaffiche()
被翻译为海报,因此我猜测它是一个输出包装器。 - fbstj