我正在使用日晷数值求解器来求解 DAE 方程。我能够创建一个工作代码,但我不确定有效保存解决方案向量的最佳方法是什么。下面是求解方程并使用for
循环将解保存在文本文件中的代码片段。
FILE* file;
int i;
tout = dt;
file = fopen("solution.txt", "w");
/*
for (iout = 1; iout < NOUT; iout++) {
tout = iout * dt; // dt : Time step
retval = IDASolve(mem, tout, &tret, uu, up, IDA_NORMAL); // Solve equation
if (retval < 0) break;
/* Save the solution vector in a text file*/
// MGRID : length of vector, uu
for (i = 0; i < MGRID; i++) {
fprintf(file, "%lf\t", NV_Ith_S(uu, i)); // access vector uu element using macro NV_Ith_S
}
fprintf(file, "\n");
}
fclose(file);
有没有其他方法可以在没有内部“for”循环的情况下将向量 uu 保存在文本文件(或其他有效格式)中?