2017-09-26 26 views
-2

通过多次调用此函数将多次打印到文件时,我一直在运行内存泄漏。没有任何问题,除非我真的打印。通过多次调用该函数向文件打印一堆文件时发生内存泄露

void print_fields_weights_energies(MKL_Complex16 *stored_fields, 
            MKL_Complex16 *stored_energies, 
            MKL_Complex16 *weights, 
            int steps, 
            int_st ist, 
            cns_st cns) { 

    /*Print The Fields Out To File*/ 
    int walkers; 
    int fields; 
    int field_number; 
    double tau = steps*cns.dtau; 
    char message[150]; 
    FILE *pf; 

    /*Name and Make File*/ 
    sprintf(message, "fields-tau%f.dat", tau); 
    pf = fopen(message, "w+"); 

    for (walkers=0; walkers<cns.max_number_walkers; walkers++) { 
    for (fields=0; fields<ist.n_spin_orbitals_sq; fields++) { 
     field_number = walkers*ist.n_spin_orbitals_sq*3+fields*3; 
     fprintf(pf, "%f\t %d\t %d\t", tau, walkers, fields); 
     fprintf(pf, "%f\t %f\t %f\t %f+%fi\t\t", stored_energies[walkers*3].real, stored_energies[walkers*3+1].real, stored_energies[walkers*3+2].real, weights[walkers].real, weights[walkers].imag); 
     fprintf(pf, "%f\t %f+%fi\t %f+%fi\n", stored_fields[field_number].real, stored_fields[field_number+1].real, stored_fields[field_number+1].imag, stored_fields[field_number+2].real, stored_fields[field_number+2].imag); 
     } 
    } 
    fflush(pf); 
    return; 
} 
+4

你永远不会调用'fclose'? – AndyG

+0

抱歉错了标签。它写在C. – JY078

+0

是的。我打电话给fclose – JY078

回答

2

你需要一个FCLOSE()在你的函数的底部(你不需要fflush() - 当它被关闭该文件将得到刷新)。

在引用它之前,如果fopen()出现错误,您还应该检查(pf == NULL)。获取错误信息远比获取SEGV更有用。

对于未来的内存泄漏,我强烈推荐使用“valgrind”工具 - 它倾向于指出发生内存泄漏的确切位置。