0

我正在使用日晷数值求解器来求解 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 保存在文本文件(或其他有效格式)中?

4

0 回答 0