Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F88854995
io_ascii.c
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Mon, Oct 21, 00:37
Size
1 KB
Mime Type
text/x-c
Expires
Wed, Oct 23, 00:37 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
21829386
Attached To
rSCINTROPARALLEL Poisson code for introduction to parallelism
io_ascii.c
View Options
#include "io_ascii.h"
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
/* -------------------------------------------------------------------------- */
int
write_to_file
(
int
n
,
float
**
x
,
unsigned
int
iter
,
float
minval
,
float
maxval
)
{
FILE
*
fp
;
char
filename
[
256
];
char
prefix
[]
=
"out_"
;
char
suffix
[]
=
".pgm"
;
char
number
[
5
];
int
i
,
j
,
lines
;
int
v
;
strcpy
(
filename
,
prefix
);
sprintf
(
number
,
"%u"
,
iter
);
while
(
iter
<
1000
)
{
strcat
(
filename
,
"0"
);
if
(
iter
==
0
)
iter
=
1
;
iter
*=
10
;
}
strcat
(
filename
,
number
);
strcat
(
filename
,
suffix
);
fp
=
fopen
(
filename
,
"w"
);
if
(
fp
==
NULL
)
{
fprintf
(
stderr
,
"Can't open output file %s!
\n
"
,
filename
);
exit
(
1
);
}
fprintf
(
fp
,
"%s
\n
%s
\n
"
,
"P2"
,
"# CREATOR: poisson program"
);
fprintf
(
fp
,
"%d %d
\n
"
,
n
,
n
);
fprintf
(
fp
,
"%d
\n
"
,
255
);
lines
=
4
;
for
(
i
=
0
;
i
<
n
;
i
++
)
{
for
(
j
=
0
;
j
<
n
;
j
++
)
{
v
=
(
int
)
roundf
(
255
*
(
x
[
i
][
j
]
-
minval
)
/
(
maxval
-
minval
));
v
=
v
<=
255
?
v
:
255
;
fprintf
(
fp
,
"%d
\n
"
,
(
int
)(
v
));
lines
++
;
}
}
fclose
(
fp
);
printf
(
"%d lines written to file %s
\n
"
,
lines
,
filename
);
return
0
;
}
/* -------------------------------------------------------------------------- */
int
colormap
(
int
n
)
{
FILE
*
fp
;
int
i
,
j
;
double
v
;
fp
=
fopen
(
"colormap.pgm"
,
"w"
);
fprintf
(
fp
,
"%s
\n
%s
\n
"
,
"P2"
,
"# CREATOR: poisson program"
);
fprintf
(
fp
,
"%d %d
\n
"
,
20
,
n
);
fprintf
(
fp
,
"%d
\n
"
,
255
);
for
(
j
=
0
;
j
<
n
;
j
++
)
{
v
=
(
int
)
roundf
(
255.0
*
j
/
(
n
-
1
));
for
(
i
=
0
;
i
<
20
;
i
++
)
{
fprintf
(
fp
,
"%d
\n
"
,
(
int
)(
v
));
}
}
fclose
(
fp
);
return
0
;
}
Event Timeline
Log In to Comment