3333
3434#define TRUE 1
3535#define FALSE 0
36- #define M 64
37- #define N 64
36+ #define M 16
37+ #define N 16
3838#define M_LEN (M + 1)
3939#define N_LEN (N + 1)
4040#define SIZE ((M_LEN)*(N_LEN))
41- #define ITMAX 4000
41+ #define ITMAX 1
4242#define L_OUT TRUE
4343#define VAL_OUT TRUE
4444
4545extern double wtime ();
4646extern void dswap (double * * a , double * * b );
4747extern void write_to_file (double * array , int tM , int tN , char * filename );
48+ extern void print_to_file (double * array , int tM , int tN , char * filename );
4849
4950//! Benchmark weather prediction program for comparing the
5051//! preformance of current supercomputers. The model is
@@ -363,9 +364,12 @@ int main(int argc, char **argv) {
363364
364365 // Output p, u, v fields and run times.
365366 if (VAL_OUT ){
366- write_to_file (pnew , M_LEN , N_LEN , "p.txt" );
367- write_to_file (unew , M_LEN , N_LEN , "u.txt" );
368- write_to_file (vnew , M_LEN , N_LEN , "v.txt" );
367+ write_to_file (pnew , M_LEN , N_LEN , "p.bin" );
368+ write_to_file (unew , M_LEN , N_LEN , "u.bin" );
369+ write_to_file (vnew , M_LEN , N_LEN , "v.bin" );
370+ print_to_file (pnew , M_LEN , N_LEN , "p.txt" );
371+ print_to_file (unew , M_LEN , N_LEN , "u.txt" );
372+ print_to_file (vnew , M_LEN , N_LEN , "v.txt" );
369373 }
370374
371375 if (L_OUT ) {
@@ -431,7 +435,7 @@ void dswap(double **pA, double **pB)
431435 * pB = pTemp ;
432436}
433437
434- void write_to_file (double * array , int tM , int tN , char * filename ) {
438+ void print_to_file (double * array , int tM , int tN , char * filename ) {
435439 FILE * file = fopen (filename , "w" );
436440 if (file == NULL ) {
437441 printf ("Error opening file %s\n" , filename );
@@ -445,3 +449,18 @@ void write_to_file(double *array, int tM, int tN, char *filename) {
445449 }
446450 fclose (file );
447451}
452+
453+ void write_to_file (double * array , int tM , int tN , char * filename ) {
454+ FILE * file = fopen (filename , "wb" );
455+ if (file == NULL ) {
456+ printf ("Error opening file %s\n" , filename );
457+ return ;
458+ }
459+ for (int i = 0 ; i < tM ; i ++ ) {
460+ for (int j = 0 ; j < tN ; j ++ ) {
461+ fwrite (& array [(i * (tN ))+ j ], sizeof (double ), 1 , file );
462+ }
463+ fprintf (file , "\n" );
464+ }
465+ fclose (file );
466+ }
0 commit comments