前提・実現したいこと
C言語でCSVへの書き出しを行いたいのですが,うまくいきません.
エラーメッセージとかは特になく,計算が終わったと思ったらファイルが作れていないという感じです.
わかる方がいらっしゃればご教示いただけないでしょうか.
該当のソースコード
#include < stdio.h >
#define _USE_MATH_DEFINES
#include < math.h >
#include < direct.h >
#define NX 100
#define NY 100
#define NZ 100
#define NT 100000
#define Tmax 3.0
void main(void)
{
double t[NX + 2][NY + 2][NZ + 2][2], f[NT + 2], dt, dx, r, a, time, pp, cp, lam, rho, tm, tinf, ss, phi, ep, sig, eta, dd, ff, p1, p2, rs, x, y, z, x1, y1, z1, xmax, ymax, zmax, t0, dy, dz;
int i, j, k, imax, jmax, kmax, n;
FILE* fp;transfer\thrust.csv", "w");
fopen_s(&fp, "C:\Users\Public\thrust.csv", "w");
dd = 1.63e-3;
xmax = 2.0e-3;
ymax = 2.0e-3;
zmax = 2.0e-3;
rs = dd / 2;
dt = Tmax / (double)NT;
dx = xmax / (double)NX;
dy = ymax / (double)NY;
dz = zmax / (double)NZ;
imax = NX;
jmax = NY;
kmax = NZ;
r = dt / dx / dx;
time = 0;
tm = 1850 + 273;
tinf = 1750 + 273;
t0 = 30 + 273;
phi = 0.7e-3;
ss = 0.25 * M_PI * pow(phi, 2);
ep = 0.32;
sig = 5.67e-8;
eta = 0.0464;
x = 0;
y = 0;
z = 0;
x1 = xmax/2;
y1 = ymax / 2;
z1 = zmax / 2;
fprintf(fp, "x,y,z,T\n");
//初期条件
//fprintf(fp, "0,");
for (i = 0; i <= imax; i++){
y = 0;
for (j = 0; j <= jmax; j++) {
z = 0;
for (k = 0; k <= kmax; k++) {
if (z >= z1 - sqrt(pow(rs, 2) - pow((y - y1), 2) - pow((x - x1), 2)) && z <= z1 + sqrt(pow(rs, 2) - pow((y - y1), 2) - pow((x - x1), 2))) t[i][j][k][0] = tinf;
else t[i][j][k][0]=t0; fprintf(fp, "%8.3e,", i * dx); fprintf(fp, "%8.3e,", j * dy); fprintf(fp, "%8.3e,", k * dz); fprintf(fp, "%8.3e\n", t[i][j][k][0]); z=z+dz; } y=y+dy; } x=x+dx; } fclose(fp); return 0;
}
試したこと
NX,NY,NZをそれぞれ10に設定したときはきちんとファイルが作れます.
また,fopen_sの部分でprintfによる書き出しを行ってみたのですが,書き出しがされなかったので,その前の部分(かなり冒頭部分)で計算が終わってしまっていると考えられます.
###環境
Visual studio 2019
回答3件
あなたの回答
tips
プレビュー