| 1 | greg | 1.1 | #ifndef lint | 
| 2 | greg | 2.58 | static const char       RCSid[] = "$Id: rv2.c,v 2.57 2008/08/29 00:30:51 greg Exp $"; | 
| 3 | greg | 1.1 | #endif | 
| 4 |  |  | /* | 
| 5 |  |  | *  rv2.c - command routines used in tracing a view. | 
| 6 |  |  | * | 
| 7 | greg | 2.38 | *  External symbols declared in rpaint.h | 
| 8 |  |  | */ | 
| 9 |  |  |  | 
| 10 | schorsch | 2.41 | #include "copyright.h" | 
| 11 |  |  |  | 
| 12 | schorsch | 2.40 | #include  <ctype.h> | 
| 13 | schorsch | 2.41 | #include  <string.h> | 
| 14 | greg | 1.1 |  | 
| 15 | schorsch | 2.40 | #include  "platform.h" | 
| 16 | greg | 1.1 | #include  "ray.h" | 
| 17 | greg | 2.46 | #include  "source.h" | 
| 18 | schorsch | 2.43 | #include  "ambient.h" | 
| 19 | greg | 1.1 | #include  "otypes.h" | 
| 20 |  |  | #include  "rpaint.h" | 
| 21 |  |  |  | 
| 22 | greg | 2.38 | extern int  psample;                    /* pixel sample size */ | 
| 23 |  |  | extern double  maxdiff;                 /* max. sample difference */ | 
| 24 | greg | 1.1 |  | 
| 25 | greg | 2.4 | #define  CTRL(c)        ((c)-'@') | 
| 26 | greg | 1.1 |  | 
| 27 | greg | 2.6 | #ifdef  SMLFLT | 
| 28 |  |  | #define  sscanvec(s,v)  (sscanf(s,"%f %f %f",v,v+1,v+2)==3) | 
| 29 |  |  | #else | 
| 30 |  |  | #define  sscanvec(s,v)  (sscanf(s,"%lf %lf %lf",v,v+1,v+2)==3) | 
| 31 |  |  | #endif | 
| 32 |  |  |  | 
| 33 | greg | 2.20 | extern char  rifname[128];              /* rad input file name */ | 
| 34 | greg | 2.18 |  | 
| 35 | greg | 1.1 | extern char  *progname; | 
| 36 | greg | 1.29 | extern char  *octname; | 
| 37 | greg | 1.1 |  | 
| 38 |  |  |  | 
| 39 | greg | 2.55 | void | 
| 40 | schorsch | 2.43 | getframe(                               /* get a new frame */ | 
| 41 |  |  | char  *s | 
| 42 |  |  | ) | 
| 43 | greg | 1.1 | { | 
| 44 | greg | 1.12 | if (getrect(s, &pframe) < 0) | 
| 45 | greg | 1.1 | return; | 
| 46 |  |  | pdepth = 0; | 
| 47 |  |  | } | 
| 48 |  |  |  | 
| 49 |  |  |  | 
| 50 | greg | 2.55 | void | 
| 51 | schorsch | 2.43 | getrepaint(                             /* get area and repaint */ | 
| 52 |  |  | char  *s | 
| 53 |  |  | ) | 
| 54 | greg | 1.12 | { | 
| 55 |  |  | RECT  box; | 
| 56 |  |  |  | 
| 57 |  |  | if (getrect(s, &box) < 0) | 
| 58 |  |  | return; | 
| 59 | greg | 2.55 | paintrect(&ptrunk, &box); | 
| 60 | greg | 1.12 | } | 
| 61 |  |  |  | 
| 62 |  |  |  | 
| 63 | greg | 2.55 | void | 
| 64 | schorsch | 2.43 | getview(                                /* get/show view parameters */ | 
| 65 |  |  | char  *s | 
| 66 |  |  | ) | 
| 67 | greg | 1.1 | { | 
| 68 |  |  | FILE  *fp; | 
| 69 |  |  | char  buf[128]; | 
| 70 | greg | 1.10 | char  *fname; | 
| 71 | greg | 1.3 | int  change = 0; | 
| 72 | greg | 2.44 | VIEW  nv = ourview; | 
| 73 | greg | 1.1 |  | 
| 74 | greg | 2.9 | while (isspace(*s)) | 
| 75 |  |  | s++; | 
| 76 |  |  | if (*s == '-') {                        /* command line parameters */ | 
| 77 |  |  | if (sscanview(&nv, s)) | 
| 78 |  |  | newview(&nv); | 
| 79 |  |  | else | 
| 80 |  |  | error(COMMAND, "bad view option(s)"); | 
| 81 |  |  | return; | 
| 82 |  |  | } | 
| 83 | greg | 1.1 | if (sscanf(s, "%s", buf) == 1) {        /* write parameters to a file */ | 
| 84 | greg | 1.10 | if ((fname = getpath(buf, NULL, 0)) == NULL || | 
| 85 | greg | 1.1 | (fp = fopen(fname, "a")) == NULL) { | 
| 86 |  |  | sprintf(errmsg, "cannot open \"%s\"", buf); | 
| 87 |  |  | error(COMMAND, errmsg); | 
| 88 |  |  | return; | 
| 89 |  |  | } | 
| 90 | greg | 1.2 | fputs(progname, fp); | 
| 91 | greg | 1.1 | fprintview(&ourview, fp); | 
| 92 | greg | 1.19 | fputs(sskip(s), fp); | 
| 93 | greg | 2.18 | putc('\n', fp); | 
| 94 | greg | 1.1 | fclose(fp); | 
| 95 |  |  | return; | 
| 96 |  |  | } | 
| 97 |  |  | sprintf(buf, "view type (%c): ", ourview.type); | 
| 98 |  |  | (*dev->comout)(buf); | 
| 99 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 100 | greg | 2.4 | if (buf[0] == CTRL('C')) return; | 
| 101 | greg | 1.3 | if (buf[0] && buf[0] != ourview.type) { | 
| 102 | greg | 1.1 | nv.type = buf[0]; | 
| 103 | greg | 1.3 | change++; | 
| 104 | greg | 2.44 | } | 
| 105 | greg | 1.1 | sprintf(buf, "view point (%.6g %.6g %.6g): ", | 
| 106 |  |  | ourview.vp[0], ourview.vp[1], ourview.vp[2]); | 
| 107 |  |  | (*dev->comout)(buf); | 
| 108 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 109 | greg | 2.4 | if (buf[0] == CTRL('C')) return; | 
| 110 | greg | 2.6 | if (sscanvec(buf, nv.vp)) | 
| 111 | greg | 1.3 | change++; | 
| 112 | greg | 1.1 | sprintf(buf, "view direction (%.6g %.6g %.6g): ", | 
| 113 | greg | 2.45 | ourview.vdir[0]*ourview.vdist, | 
| 114 |  |  | ourview.vdir[1]*ourview.vdist, | 
| 115 |  |  | ourview.vdir[2]*ourview.vdist); | 
| 116 | greg | 1.1 | (*dev->comout)(buf); | 
| 117 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 118 | greg | 2.4 | if (buf[0] == CTRL('C')) return; | 
| 119 | greg | 2.45 | if (sscanvec(buf, nv.vdir)) { | 
| 120 |  |  | nv.vdist = 1.; | 
| 121 | greg | 1.3 | change++; | 
| 122 | greg | 2.45 | } | 
| 123 | greg | 1.1 | sprintf(buf, "view up (%.6g %.6g %.6g): ", | 
| 124 |  |  | ourview.vup[0], ourview.vup[1], ourview.vup[2]); | 
| 125 |  |  | (*dev->comout)(buf); | 
| 126 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 127 | greg | 2.4 | if (buf[0] == CTRL('C')) return; | 
| 128 | greg | 2.6 | if (sscanvec(buf, nv.vup)) | 
| 129 | greg | 1.3 | change++; | 
| 130 | greg | 1.1 | sprintf(buf, "view horiz and vert size (%.6g %.6g): ", | 
| 131 |  |  | ourview.horiz, ourview.vert); | 
| 132 |  |  | (*dev->comout)(buf); | 
| 133 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 134 | greg | 2.4 | if (buf[0] == CTRL('C')) return; | 
| 135 | greg | 1.3 | if (sscanf(buf, "%lf %lf", &nv.horiz, &nv.vert) == 2) | 
| 136 |  |  | change++; | 
| 137 | greg | 2.26 | sprintf(buf, "fore and aft clipping plane (%.6g %.6g): ", | 
| 138 |  |  | ourview.vfore, ourview.vaft); | 
| 139 |  |  | (*dev->comout)(buf); | 
| 140 |  |  | (*dev->comin)(buf, NULL); | 
| 141 |  |  | if (buf[0] == CTRL('C')) return; | 
| 142 |  |  | if (sscanf(buf, "%lf %lf", &nv.vfore, &nv.vaft) == 2) | 
| 143 |  |  | change++; | 
| 144 | greg | 1.14 | sprintf(buf, "view shift and lift (%.6g %.6g): ", | 
| 145 |  |  | ourview.hoff, ourview.voff); | 
| 146 | greg | 1.1 | (*dev->comout)(buf); | 
| 147 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 148 | greg | 2.4 | if (buf[0] == CTRL('C')) return; | 
| 149 | greg | 1.14 | if (sscanf(buf, "%lf %lf", &nv.hoff, &nv.voff) == 2) | 
| 150 | greg | 1.3 | change++; | 
| 151 |  |  | if (change) | 
| 152 |  |  | newview(&nv); | 
| 153 | greg | 1.1 | } | 
| 154 |  |  |  | 
| 155 |  |  |  | 
| 156 | greg | 2.55 | void | 
| 157 | schorsch | 2.43 | lastview(                               /* return to a previous view */ | 
| 158 |  |  | char  *s | 
| 159 |  |  | ) | 
| 160 | greg | 1.1 | { | 
| 161 |  |  | char  buf[128]; | 
| 162 | greg | 1.10 | char  *fname; | 
| 163 | greg | 1.1 | int  success; | 
| 164 |  |  | VIEW  nv; | 
| 165 |  |  |  | 
| 166 |  |  | if (sscanf(s, "%s", buf) == 1) {        /* get parameters from a file */ | 
| 167 | schorsch | 2.42 | nv = stdview; | 
| 168 | greg | 2.19 | if ((fname = getpath(buf, "", R_OK)) == NULL || | 
| 169 | greg | 1.33 | (success = viewfile(fname, &nv, NULL)) == -1) { | 
| 170 | greg | 1.1 | sprintf(errmsg, "cannot open \"%s\"", buf); | 
| 171 |  |  | error(COMMAND, errmsg); | 
| 172 |  |  | return; | 
| 173 |  |  | } | 
| 174 |  |  | if (!success) | 
| 175 |  |  | error(COMMAND, "wrong file format"); | 
| 176 |  |  | else | 
| 177 |  |  | newview(&nv); | 
| 178 |  |  | return; | 
| 179 |  |  | } | 
| 180 | greg | 1.20 | if (oldview.type == 0) {        /* no old view! */ | 
| 181 | greg | 1.1 | error(COMMAND, "no previous view"); | 
| 182 |  |  | return; | 
| 183 |  |  | } | 
| 184 | schorsch | 2.42 | nv = ourview; | 
| 185 |  |  | ourview = oldview; | 
| 186 |  |  | oldview = nv; | 
| 187 | greg | 2.55 | newimage(NULL); | 
| 188 | greg | 2.18 | } | 
| 189 |  |  |  | 
| 190 |  |  |  | 
| 191 | greg | 2.55 | void | 
| 192 | schorsch | 2.43 | saveview(                               /* save view to rad file */ | 
| 193 |  |  | char  *s | 
| 194 |  |  | ) | 
| 195 | greg | 2.18 | { | 
| 196 |  |  | char  view[64]; | 
| 197 |  |  | char  *fname; | 
| 198 |  |  | FILE  *fp; | 
| 199 |  |  |  | 
| 200 | greg | 2.19 | if (*atos(view, sizeof(view), s)) { | 
| 201 |  |  | if (isint(view)) { | 
| 202 |  |  | error(COMMAND, "cannot write view by number"); | 
| 203 |  |  | return; | 
| 204 |  |  | } | 
| 205 | greg | 2.18 | s = sskip(s); | 
| 206 | greg | 2.19 | } | 
| 207 | greg | 2.20 | while (isspace(*s)) | 
| 208 |  |  | s++; | 
| 209 |  |  | if (*s) | 
| 210 |  |  | atos(rifname, sizeof(rifname), s); | 
| 211 |  |  | else if (rifname[0] == '\0') { | 
| 212 | greg | 2.18 | error(COMMAND, "no previous rad file"); | 
| 213 |  |  | return; | 
| 214 |  |  | } | 
| 215 |  |  | if ((fname = getpath(rifname, NULL, 0)) == NULL || | 
| 216 |  |  | (fp = fopen(fname, "a")) == NULL) { | 
| 217 |  |  | sprintf(errmsg, "cannot open \"%s\"", rifname); | 
| 218 |  |  | error(COMMAND, errmsg); | 
| 219 |  |  | return; | 
| 220 |  |  | } | 
| 221 |  |  | fputs("view= ", fp); | 
| 222 |  |  | fputs(view, fp); | 
| 223 |  |  | fprintview(&ourview, fp); | 
| 224 |  |  | putc('\n', fp); | 
| 225 |  |  | fclose(fp); | 
| 226 |  |  | } | 
| 227 |  |  |  | 
| 228 |  |  |  | 
| 229 | greg | 2.55 | void | 
| 230 | schorsch | 2.43 | loadview(                               /* load view from rad file */ | 
| 231 |  |  | char  *s | 
| 232 |  |  | ) | 
| 233 | greg | 2.18 | { | 
| 234 |  |  | char  buf[512]; | 
| 235 |  |  | char  *fname; | 
| 236 |  |  | FILE  *fp; | 
| 237 |  |  | VIEW  nv; | 
| 238 |  |  |  | 
| 239 |  |  | strcpy(buf, "rad -n -s -V -v "); | 
| 240 |  |  | if (sscanf(s, "%s", buf+strlen(buf)) == 1) | 
| 241 |  |  | s = sskip(s); | 
| 242 |  |  | else | 
| 243 |  |  | strcat(buf, "1"); | 
| 244 | greg | 2.20 | if (*s) | 
| 245 |  |  | atos(rifname, sizeof(rifname), s); | 
| 246 |  |  | else if (rifname[0] == '\0') { | 
| 247 | greg | 2.18 | error(COMMAND, "no previous rad file"); | 
| 248 |  |  | return; | 
| 249 |  |  | } | 
| 250 | greg | 2.19 | if ((fname = getpath(rifname, "", R_OK)) == NULL) { | 
| 251 | greg | 2.18 | sprintf(errmsg, "cannot access \"%s\"", rifname); | 
| 252 |  |  | error(COMMAND, errmsg); | 
| 253 |  |  | return; | 
| 254 |  |  | } | 
| 255 |  |  | sprintf(buf+strlen(buf), " %s", fname); | 
| 256 |  |  | if ((fp = popen(buf, "r")) == NULL) { | 
| 257 |  |  | error(COMMAND, "cannot run rad"); | 
| 258 |  |  | return; | 
| 259 |  |  | } | 
| 260 |  |  | buf[0] = '\0'; | 
| 261 |  |  | fgets(buf, sizeof(buf), fp); | 
| 262 |  |  | pclose(fp); | 
| 263 | schorsch | 2.42 | nv = stdview; | 
| 264 | greg | 2.18 | if (!sscanview(&nv, buf)) { | 
| 265 |  |  | error(COMMAND, "rad error -- no such view?"); | 
| 266 |  |  | return; | 
| 267 |  |  | } | 
| 268 |  |  | newview(&nv); | 
| 269 | greg | 1.1 | } | 
| 270 |  |  |  | 
| 271 |  |  |  | 
| 272 | greg | 2.55 | void | 
| 273 | schorsch | 2.43 | getaim(                         /* aim camera */ | 
| 274 |  |  | char  *s | 
| 275 |  |  | ) | 
| 276 | greg | 1.1 | { | 
| 277 | greg | 2.44 | VIEW  nv = ourview; | 
| 278 | greg | 1.1 | double  zfact; | 
| 279 |  |  |  | 
| 280 |  |  | if (getinterest(s, 1, nv.vdir, &zfact) < 0) | 
| 281 |  |  | return; | 
| 282 | greg | 1.21 | zoomview(&nv, zfact); | 
| 283 | greg | 1.1 | newview(&nv); | 
| 284 |  |  | } | 
| 285 |  |  |  | 
| 286 |  |  |  | 
| 287 | greg | 2.55 | void | 
| 288 | greg | 2.45 | getfocus(                               /* set focus distance */ | 
| 289 |  |  | char *s | 
| 290 |  |  | ) | 
| 291 |  |  | { | 
| 292 | greg | 2.53 | char  buf[64]; | 
| 293 | greg | 2.45 | double  dist; | 
| 294 |  |  |  | 
| 295 |  |  | if (sscanf(s, "%lf", &dist) < 1) { | 
| 296 |  |  | int     x, y; | 
| 297 |  |  | RAY     thisray; | 
| 298 |  |  | if (dev->getcur == NULL) | 
| 299 |  |  | return; | 
| 300 |  |  | (*dev->comout)("Pick focus point\n"); | 
| 301 |  |  | if ((*dev->getcur)(&x, &y) == ABORT) | 
| 302 |  |  | return; | 
| 303 |  |  | if ((thisray.rmax = viewray(thisray.rorg, thisray.rdir, | 
| 304 |  |  | &ourview, (x+.5)/hresolu, (y+.5)/vresolu)) < -FTINY) { | 
| 305 |  |  | error(COMMAND, "not on image"); | 
| 306 |  |  | return; | 
| 307 |  |  | } | 
| 308 | greg | 2.49 | rayorigin(&thisray, PRIMARY, NULL, NULL); | 
| 309 | greg | 2.45 | if (!localhit(&thisray, &thescene)) { | 
| 310 |  |  | error(COMMAND, "not a local object"); | 
| 311 |  |  | return; | 
| 312 |  |  | } | 
| 313 |  |  | dist = thisray.rot; | 
| 314 |  |  | } else if (dist <= .0) { | 
| 315 |  |  | error(COMMAND, "focus distance must be positive"); | 
| 316 |  |  | return; | 
| 317 |  |  | } | 
| 318 |  |  | ourview.vdist = dist; | 
| 319 | greg | 2.53 | sprintf(buf, "Focus distance set to %f\n", dist); | 
| 320 |  |  | (*dev->comout)(buf); | 
| 321 | greg | 2.45 | } | 
| 322 |  |  |  | 
| 323 |  |  |  | 
| 324 | greg | 2.55 | void | 
| 325 | schorsch | 2.43 | getmove(                                /* move camera */ | 
| 326 |  |  | char  *s | 
| 327 |  |  | ) | 
| 328 | greg | 1.1 | { | 
| 329 |  |  | FVECT  vc; | 
| 330 |  |  | double  mag; | 
| 331 |  |  |  | 
| 332 |  |  | if (getinterest(s, 0, vc, &mag) < 0) | 
| 333 |  |  | return; | 
| 334 | greg | 1.4 | moveview(0.0, 0.0, mag, vc); | 
| 335 | greg | 1.1 | } | 
| 336 |  |  |  | 
| 337 |  |  |  | 
| 338 | greg | 2.55 | void | 
| 339 | schorsch | 2.43 | getrotate(                              /* rotate camera */ | 
| 340 |  |  | char  *s | 
| 341 |  |  | ) | 
| 342 | greg | 1.1 | { | 
| 343 | greg | 2.44 | VIEW  nv = ourview; | 
| 344 | greg | 1.1 | FVECT  v1; | 
| 345 | greg | 1.4 | double  angle, elev, zfact; | 
| 346 | greg | 1.1 |  | 
| 347 | greg | 1.4 | elev = 0.0; zfact = 1.0; | 
| 348 |  |  | if (sscanf(s, "%lf %lf %lf", &angle, &elev, &zfact) < 1) { | 
| 349 | greg | 1.1 | error(COMMAND, "missing angle"); | 
| 350 |  |  | return; | 
| 351 |  |  | } | 
| 352 |  |  | spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.)); | 
| 353 |  |  | if (elev != 0.0) { | 
| 354 |  |  | fcross(v1, nv.vdir, ourview.vup); | 
| 355 |  |  | normalize(v1); | 
| 356 |  |  | spinvector(nv.vdir, nv.vdir, v1, elev*(PI/180.)); | 
| 357 |  |  | } | 
| 358 | greg | 1.21 | zoomview(&nv, zfact); | 
| 359 | greg | 1.1 | newview(&nv); | 
| 360 |  |  | } | 
| 361 |  |  |  | 
| 362 |  |  |  | 
| 363 | greg | 2.55 | void | 
| 364 | schorsch | 2.43 | getpivot(                               /* pivot viewpoint */ | 
| 365 | greg | 2.55 | char  *s | 
| 366 | schorsch | 2.43 | ) | 
| 367 | greg | 1.1 | { | 
| 368 |  |  | FVECT  vc; | 
| 369 | greg | 1.4 | double  angle, elev, mag; | 
| 370 | greg | 1.1 |  | 
| 371 | greg | 1.4 | elev = 0.0; | 
| 372 |  |  | if (sscanf(s, "%lf %lf", &angle, &elev) < 1) { | 
| 373 | greg | 1.1 | error(COMMAND, "missing angle"); | 
| 374 |  |  | return; | 
| 375 |  |  | } | 
| 376 | greg | 2.33 | if (getinterest(sskip2(s,2), 0, vc, &mag) < 0) | 
| 377 | greg | 1.1 | return; | 
| 378 | greg | 1.4 | moveview(angle, elev, mag, vc); | 
| 379 | greg | 1.1 | } | 
| 380 |  |  |  | 
| 381 |  |  |  | 
| 382 | greg | 2.55 | void | 
| 383 | schorsch | 2.43 | getexposure(                            /* get new exposure */ | 
| 384 |  |  | char  *s | 
| 385 |  |  | ) | 
| 386 | greg | 1.1 | { | 
| 387 |  |  | char  buf[128]; | 
| 388 | greg | 2.55 | char  *cp; | 
| 389 | greg | 1.1 | int  x, y; | 
| 390 | greg | 2.55 | PNODE  *p = &ptrunk; | 
| 391 | greg | 2.16 | int  adapt = 0; | 
| 392 |  |  | double  e = 1.0; | 
| 393 | greg | 2.25 |  | 
| 394 | greg | 1.1 | for (cp = s; isspace(*cp); cp++) | 
| 395 |  |  | ; | 
| 396 | greg | 2.16 | if (*cp == '@') { | 
| 397 |  |  | adapt++; | 
| 398 | greg | 2.25 | while (isspace(*++cp)) | 
| 399 |  |  | ; | 
| 400 | greg | 2.16 | } | 
| 401 | greg | 1.1 | if (*cp == '\0') {              /* normalize to point */ | 
| 402 |  |  | if (dev->getcur == NULL) | 
| 403 |  |  | return; | 
| 404 |  |  | (*dev->comout)("Pick point for exposure\n"); | 
| 405 |  |  | if ((*dev->getcur)(&x, &y) == ABORT) | 
| 406 |  |  | return; | 
| 407 | greg | 2.55 | p = findrect(x, y, &ptrunk, -1); | 
| 408 | greg | 1.1 | } else { | 
| 409 |  |  | if (*cp == '=') {       /* absolute setting */ | 
| 410 |  |  | p = NULL; | 
| 411 |  |  | e = 1.0/exposure; | 
| 412 |  |  | for (cp++; isspace(*cp); cp++) | 
| 413 |  |  | ; | 
| 414 |  |  | if (*cp == '\0') {      /* interactive */ | 
| 415 | greg | 2.15 | sprintf(buf, "exposure (%f): ", exposure); | 
| 416 | greg | 1.1 | (*dev->comout)(buf); | 
| 417 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 418 | greg | 1.1 | for (cp = buf; isspace(*cp); cp++) | 
| 419 |  |  | ; | 
| 420 |  |  | if (*cp == '\0') | 
| 421 |  |  | return; | 
| 422 |  |  | } | 
| 423 |  |  | } | 
| 424 |  |  | if (*cp == '+' || *cp == '-')   /* f-stops */ | 
| 425 |  |  | e *= pow(2.0, atof(cp)); | 
| 426 |  |  | else                            /* multiplier */ | 
| 427 |  |  | e *= atof(cp); | 
| 428 |  |  | } | 
| 429 |  |  | if (p != NULL) {                /* relative setting */ | 
| 430 | greg | 2.58 | compavg(p); | 
| 431 | greg | 2.14 | if (bright(p->v) < 1e-15) { | 
| 432 | greg | 1.1 | error(COMMAND, "cannot normalize to zero"); | 
| 433 |  |  | return; | 
| 434 |  |  | } | 
| 435 | greg | 2.16 | if (adapt) | 
| 436 |  |  | e *= 106./pow(1.219+pow(luminance(p->v)/exposure,.4),2.5)/exposure; | 
| 437 |  |  | else | 
| 438 |  |  | e *= 0.5 / bright(p->v); | 
| 439 | greg | 1.1 | } | 
| 440 |  |  | if (e <= FTINY || fabs(1.0 - e) <= FTINY) | 
| 441 |  |  | return; | 
| 442 |  |  | scalepict(&ptrunk, e); | 
| 443 |  |  | exposure *= e; | 
| 444 |  |  | redraw(); | 
| 445 |  |  | } | 
| 446 |  |  |  | 
| 447 | greg | 2.38 | typedef union {int i; double d; COLOR C;}       *MyUptr; | 
| 448 | greg | 2.57 | #define  FEQ(x,y)     (fabs((x)-(y)) <= FTINY) | 
| 449 | greg | 1.1 |  | 
| 450 | greg | 2.55 | int | 
| 451 | schorsch | 2.43 | getparam(               /* get variable from user */ | 
| 452 |  |  | char  *str, | 
| 453 |  |  | char  *dsc, | 
| 454 |  |  | int  typ, | 
| 455 |  |  | void  *p | 
| 456 |  |  | ) | 
| 457 | greg | 1.24 | { | 
| 458 | greg | 2.55 | MyUptr  ptr = (MyUptr)p; | 
| 459 | greg | 1.24 | int  i0; | 
| 460 |  |  | double  d0, d1, d2; | 
| 461 |  |  | char  buf[48]; | 
| 462 |  |  |  | 
| 463 |  |  | switch (typ) { | 
| 464 |  |  | case 'i':                       /* integer */ | 
| 465 |  |  | if (sscanf(str, "%d", &i0) != 1) { | 
| 466 |  |  | (*dev->comout)(dsc); | 
| 467 |  |  | sprintf(buf, " (%d): ", ptr->i); | 
| 468 |  |  | (*dev->comout)(buf); | 
| 469 |  |  | (*dev->comin)(buf, NULL); | 
| 470 |  |  | if (sscanf(buf, "%d", &i0) != 1) | 
| 471 | greg | 1.31 | return(0); | 
| 472 | greg | 1.24 | } | 
| 473 | greg | 2.57 | if (ptr->i == i0) | 
| 474 |  |  | return(0); | 
| 475 | greg | 1.24 | ptr->i = i0; | 
| 476 | greg | 2.55 | break; | 
| 477 | greg | 1.24 | case 'r':                       /* real */ | 
| 478 |  |  | if (sscanf(str, "%lf", &d0) != 1) { | 
| 479 |  |  | (*dev->comout)(dsc); | 
| 480 |  |  | sprintf(buf, " (%.6g): ", ptr->d); | 
| 481 |  |  | (*dev->comout)(buf); | 
| 482 |  |  | (*dev->comin)(buf, NULL); | 
| 483 |  |  | if (sscanf(buf, "%lf", &d0) != 1) | 
| 484 | greg | 1.31 | return(0); | 
| 485 | greg | 1.24 | } | 
| 486 | greg | 2.57 | if (FEQ(ptr->d, d0)) | 
| 487 |  |  | return(0); | 
| 488 | greg | 1.24 | ptr->d = d0; | 
| 489 | greg | 2.55 | break; | 
| 490 | greg | 1.24 | case 'b':                       /* boolean */ | 
| 491 |  |  | if (sscanf(str, "%1s", buf) != 1) { | 
| 492 |  |  | (*dev->comout)(dsc); | 
| 493 | greg | 2.2 | sprintf(buf, "? (%c): ", ptr->i ? 'y' : 'n'); | 
| 494 | greg | 1.24 | (*dev->comout)(buf); | 
| 495 |  |  | (*dev->comin)(buf, NULL); | 
| 496 | greg | 2.57 | if (buf[0] == '\0') | 
| 497 | greg | 1.31 | return(0); | 
| 498 | greg | 1.24 | } | 
| 499 | greg | 2.57 | if (strchr("yY+1tTnN-0fF", buf[0]) == NULL) | 
| 500 |  |  | return(0); | 
| 501 |  |  | i0 = strchr("yY+1tT", buf[0]) != NULL; | 
| 502 |  |  | if (ptr->i == i0) | 
| 503 |  |  | return(0); | 
| 504 |  |  | ptr->i = i0; | 
| 505 | greg | 2.55 | break; | 
| 506 | greg | 1.24 | case 'C':                       /* color */ | 
| 507 |  |  | if (sscanf(str, "%lf %lf %lf", &d0, &d1, &d2) != 3) { | 
| 508 |  |  | (*dev->comout)(dsc); | 
| 509 |  |  | sprintf(buf, " (%.6g %.6g %.6g): ", | 
| 510 |  |  | colval(ptr->C,RED), | 
| 511 |  |  | colval(ptr->C,GRN), | 
| 512 |  |  | colval(ptr->C,BLU)); | 
| 513 |  |  | (*dev->comout)(buf); | 
| 514 |  |  | (*dev->comin)(buf, NULL); | 
| 515 |  |  | if (sscanf(buf, "%lf %lf %lf", &d0, &d1, &d2) != 3) | 
| 516 | greg | 1.31 | return(0); | 
| 517 | greg | 1.24 | } | 
| 518 | greg | 2.57 | if (FEQ(colval(ptr->C,RED), d0) && | 
| 519 |  |  | FEQ(colval(ptr->C,GRN), d1) && | 
| 520 |  |  | FEQ(colval(ptr->C,BLU), d2)) | 
| 521 |  |  | return(0); | 
| 522 | greg | 1.24 | setcolor(ptr->C, d0, d1, d2); | 
| 523 | greg | 2.55 | break; | 
| 524 |  |  | default: | 
| 525 |  |  | return(0);              /* shouldn't happen */ | 
| 526 | greg | 1.24 | } | 
| 527 | greg | 2.55 | newparam++; | 
| 528 |  |  | return(1); | 
| 529 | greg | 1.24 | } | 
| 530 |  |  |  | 
| 531 |  |  |  | 
| 532 | greg | 2.55 | void | 
| 533 | schorsch | 2.43 | setparam(                               /* get/set program parameter */ | 
| 534 | greg | 2.55 | char  *s | 
| 535 | schorsch | 2.43 | ) | 
| 536 | greg | 1.1 | { | 
| 537 | greg | 1.24 | char  buf[128]; | 
| 538 | greg | 1.1 |  | 
| 539 |  |  | if (s[0] == '\0') { | 
| 540 | greg | 1.23 | (*dev->comout)( | 
| 541 | greg | 2.51 | "aa ab ad ar as av aw b bv dc dv dj ds dt i lr lw me ma mg ms ps pt sj st u: "); | 
| 542 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 543 | greg | 1.1 | s = buf; | 
| 544 |  |  | } | 
| 545 |  |  | switch (s[0]) { | 
| 546 | greg | 2.51 | case 'u':                       /* uncorrelated sampling */ | 
| 547 |  |  | getparam(s+1, "uncorrelated sampling", 'b', | 
| 548 | greg | 2.50 | (void *)&rand_samp); | 
| 549 |  |  | break; | 
| 550 | greg | 1.1 | case 'l':                       /* limit */ | 
| 551 |  |  | switch (s[1]) { | 
| 552 |  |  | case 'w':                       /* weight */ | 
| 553 | greg | 2.38 | getparam(s+2, "limit weight", 'r', | 
| 554 |  |  | (void *)&minweight); | 
| 555 | greg | 1.1 | break; | 
| 556 |  |  | case 'r':                       /* reflection */ | 
| 557 | greg | 2.38 | getparam(s+2, "limit reflection", 'i', | 
| 558 |  |  | (void *)&maxdepth); | 
| 559 | greg | 1.1 | break; | 
| 560 |  |  | default: | 
| 561 |  |  | goto badparam; | 
| 562 |  |  | } | 
| 563 |  |  | break; | 
| 564 | greg | 1.5 | case 'd':                       /* direct */ | 
| 565 |  |  | switch (s[1]) { | 
| 566 |  |  | case 'j':                       /* jitter */ | 
| 567 | greg | 2.38 | getparam(s+2, "direct jitter", 'r', | 
| 568 |  |  | (void *)&dstrsrc); | 
| 569 | greg | 1.9 | break; | 
| 570 |  |  | case 'c':                       /* certainty */ | 
| 571 | greg | 2.38 | getparam(s+2, "direct certainty", 'r', | 
| 572 |  |  | (void *)&shadcert); | 
| 573 | greg | 1.5 | break; | 
| 574 | greg | 1.8 | case 't':                       /* threshold */ | 
| 575 | greg | 2.38 | getparam(s+2, "direct threshold", 'r', | 
| 576 |  |  | (void *)&shadthresh); | 
| 577 | greg | 1.5 | break; | 
| 578 | greg | 2.13 | case 'v':                       /* visibility */ | 
| 579 | greg | 2.38 | getparam(s+2, "direct visibility", 'b', | 
| 580 |  |  | (void *)&directvis); | 
| 581 | greg | 1.30 | break; | 
| 582 |  |  | case 's':                       /* sampling */ | 
| 583 | greg | 2.38 | getparam(s+2, "direct sampling", 'r', | 
| 584 |  |  | (void *)&srcsizerat); | 
| 585 | greg | 1.23 | break; | 
| 586 | greg | 1.5 | default: | 
| 587 | greg | 1.1 | goto badparam; | 
| 588 |  |  | } | 
| 589 | greg | 1.23 | break; | 
| 590 | greg | 2.27 | case 'b':                       /* back faces or black and white */ | 
| 591 |  |  | switch (s[1]) { | 
| 592 |  |  | case 'v':                       /* back face visibility */ | 
| 593 | greg | 2.38 | getparam(s+2, "back face visibility", 'b', | 
| 594 |  |  | (void *)&backvis); | 
| 595 | greg | 2.27 | break; | 
| 596 |  |  | case '\0':                      /* black and white */ | 
| 597 | greg | 2.29 | case ' ': | 
| 598 | greg | 2.27 | case 'y': case 'Y': case 't': case 'T': case '1': case '+': | 
| 599 |  |  | case 'n': case 'N': case 'f': case 'F': case '0': case '-': | 
| 600 | greg | 2.38 | getparam(s+1, "black and white", 'b', | 
| 601 |  |  | (void *)&greyscale); | 
| 602 | greg | 2.27 | break; | 
| 603 |  |  | default: | 
| 604 |  |  | goto badparam; | 
| 605 |  |  | } | 
| 606 | greg | 1.23 | break; | 
| 607 |  |  | case 'i':                       /* irradiance */ | 
| 608 | greg | 2.38 | getparam(s+1, "irradiance", 'b', | 
| 609 |  |  | (void *)&do_irrad); | 
| 610 | greg | 1.1 | break; | 
| 611 |  |  | case 'a':                       /* ambient */ | 
| 612 |  |  | switch (s[1]) { | 
| 613 |  |  | case 'v':                       /* value */ | 
| 614 | greg | 2.38 | getparam(s+2, "ambient value", 'C', | 
| 615 |  |  | (void *)ambval); | 
| 616 | greg | 2.34 | break; | 
| 617 |  |  | case 'w':                       /* weight */ | 
| 618 | greg | 2.38 | getparam(s+2, "ambient value weight", 'i', | 
| 619 |  |  | (void *)&ambvwt); | 
| 620 | greg | 1.1 | break; | 
| 621 |  |  | case 'a':                       /* accuracy */ | 
| 622 | greg | 2.38 | if (getparam(s+2, "ambient accuracy", 'r', | 
| 623 |  |  | (void *)&ambacc)) | 
| 624 | greg | 2.17 | setambacc(ambacc); | 
| 625 | greg | 1.1 | break; | 
| 626 |  |  | case 'd':                       /* divisions */ | 
| 627 | greg | 2.38 | getparam(s+2, "ambient divisions", 'i', | 
| 628 |  |  | (void *)&ambdiv); | 
| 629 | greg | 1.1 | break; | 
| 630 |  |  | case 's':                       /* samples */ | 
| 631 | greg | 2.38 | getparam(s+2, "ambient super-samples", 'i', | 
| 632 |  |  | (void *)&ambssamp); | 
| 633 | greg | 1.1 | break; | 
| 634 |  |  | case 'b':                       /* bounces */ | 
| 635 | greg | 2.38 | getparam(s+2, "ambient bounces", 'i', | 
| 636 |  |  | (void *)&ambounce); | 
| 637 | greg | 1.1 | break; | 
| 638 |  |  | case 'r': | 
| 639 | greg | 2.38 | if (getparam(s+2, "ambient resolution", 'i', | 
| 640 |  |  | (void *)&ambres)) | 
| 641 | greg | 2.8 | setambres(ambres); | 
| 642 | greg | 2.32 | break; | 
| 643 |  |  | default: | 
| 644 |  |  | goto badparam; | 
| 645 |  |  | } | 
| 646 |  |  | break; | 
| 647 |  |  | case 'm':                       /* medium */ | 
| 648 |  |  | switch (s[1]) { | 
| 649 |  |  | case 'e':                       /* extinction coefficient */ | 
| 650 |  |  | getparam(s+2, "extinction coefficient", 'C', | 
| 651 | greg | 2.38 | (void *)cextinction); | 
| 652 | greg | 2.32 | break; | 
| 653 |  |  | case 'a':                       /* scattering albedo */ | 
| 654 | greg | 2.35 | getparam(s+2, "scattering albedo", 'C', | 
| 655 | greg | 2.38 | (void *)salbedo); | 
| 656 | greg | 2.32 | break; | 
| 657 |  |  | case 'g':                       /* scattering eccentricity */ | 
| 658 | greg | 2.38 | getparam(s+2, "scattering eccentricity", 'r', | 
| 659 |  |  | (void *)&seccg); | 
| 660 | greg | 2.32 | break; | 
| 661 |  |  | case 's':                       /* sampling distance */ | 
| 662 |  |  | getparam(s+2, "mist sampling distance", 'r', | 
| 663 | greg | 2.38 | (void *)&ssampdist); | 
| 664 | greg | 1.1 | break; | 
| 665 |  |  | default: | 
| 666 |  |  | goto badparam; | 
| 667 |  |  | } | 
| 668 |  |  | break; | 
| 669 | greg | 2.5 | case 'p':                       /* pixel */ | 
| 670 | greg | 1.1 | switch (s[1]) { | 
| 671 | greg | 2.5 | case 's':                       /* sample */ | 
| 672 | greg | 2.38 | if (getparam(s+2, "pixel sample", 'i', | 
| 673 |  |  | (void *)&psample)) | 
| 674 | greg | 1.31 | pdepth = 0; | 
| 675 | greg | 1.1 | break; | 
| 676 | greg | 1.8 | case 't':                       /* threshold */ | 
| 677 | greg | 2.38 | if (getparam(s+2, "pixel threshold", 'r', | 
| 678 |  |  | (void *)&maxdiff)) | 
| 679 | greg | 1.31 | pdepth = 0; | 
| 680 | greg | 2.5 | break; | 
| 681 |  |  | default: | 
| 682 |  |  | goto badparam; | 
| 683 |  |  | } | 
| 684 |  |  | break; | 
| 685 |  |  | case 's':                       /* specular */ | 
| 686 |  |  | switch (s[1]) { | 
| 687 |  |  | case 'j':                       /* jitter */ | 
| 688 | greg | 2.38 | getparam(s+2, "specular jitter", 'r', | 
| 689 |  |  | (void *)&specjitter); | 
| 690 | greg | 2.5 | break; | 
| 691 |  |  | case 't':                       /* threshold */ | 
| 692 | greg | 2.38 | getparam(s+2, "specular threshold", 'r', | 
| 693 |  |  | (void *)&specthresh); | 
| 694 | greg | 1.1 | break; | 
| 695 |  |  | default: | 
| 696 |  |  | goto badparam; | 
| 697 |  |  | } | 
| 698 |  |  | break; | 
| 699 |  |  | case '\0':                      /* nothing */ | 
| 700 |  |  | break; | 
| 701 |  |  | default:; | 
| 702 |  |  | badparam: | 
| 703 | greg | 1.26 | *sskip(s) = '\0'; | 
| 704 | greg | 1.1 | sprintf(errmsg, "%s: unknown variable", s); | 
| 705 |  |  | error(COMMAND, errmsg); | 
| 706 |  |  | break; | 
| 707 |  |  | } | 
| 708 |  |  | } | 
| 709 |  |  |  | 
| 710 |  |  |  | 
| 711 | greg | 2.55 | void | 
| 712 |  |  | traceray(                               /* trace a single ray */ | 
| 713 |  |  | char  *s | 
| 714 |  |  | ) | 
| 715 | greg | 1.1 | { | 
| 716 |  |  | char  buf[128]; | 
| 717 |  |  | int  x, y; | 
| 718 | greg | 2.37 | OBJREC  *ino; | 
| 719 | greg | 1.1 | RAY  thisray; | 
| 720 |  |  |  | 
| 721 | greg | 2.26 | thisray.rmax = 0.0; | 
| 722 |  |  |  | 
| 723 | greg | 2.6 | if (!sscanvec(s, thisray.rorg) || | 
| 724 | greg | 2.33 | !sscanvec(sskip2(s,3), thisray.rdir)) { | 
| 725 | greg | 1.1 |  | 
| 726 |  |  | if (dev->getcur == NULL) | 
| 727 |  |  | return; | 
| 728 |  |  | (*dev->comout)("Pick ray\n"); | 
| 729 |  |  | if ((*dev->getcur)(&x, &y) == ABORT) | 
| 730 |  |  | return; | 
| 731 |  |  |  | 
| 732 | greg | 2.26 | if ((thisray.rmax = viewray(thisray.rorg, thisray.rdir, | 
| 733 |  |  | &ourview, (x+.5)/hresolu, (y+.5)/vresolu)) < -FTINY) { | 
| 734 | greg | 1.21 | error(COMMAND, "not on image"); | 
| 735 |  |  | return; | 
| 736 |  |  | } | 
| 737 |  |  |  | 
| 738 | greg | 1.1 | } else if (normalize(thisray.rdir) == 0.0) { | 
| 739 |  |  | error(COMMAND, "zero ray direction"); | 
| 740 |  |  | return; | 
| 741 |  |  | } | 
| 742 |  |  |  | 
| 743 | greg | 2.56 | ray_trace(&thisray); | 
| 744 | greg | 1.1 |  | 
| 745 |  |  | if (thisray.ro == NULL) | 
| 746 |  |  | (*dev->comout)("ray hit nothing"); | 
| 747 |  |  | else { | 
| 748 | greg | 2.47 | OBJREC  *mat = NULL; | 
| 749 | greg | 2.52 | OBJREC  *mod = NULL; | 
| 750 | greg | 2.47 | char    matspec[256]; | 
| 751 |  |  | matspec[0] = '\0'; | 
| 752 |  |  | if (thisray.ro->omod != OVOID) { | 
| 753 |  |  | mod = objptr(thisray.ro->omod); | 
| 754 |  |  | mat = findmaterial(mod); | 
| 755 |  |  | } | 
| 756 |  |  | if (thisray.rod < 0.0) | 
| 757 |  |  | strcpy(matspec, "back of "); | 
| 758 |  |  | if (mod != NULL) { | 
| 759 |  |  | strcat(matspec, mod->oname); | 
| 760 | greg | 2.48 | if (mat != mod && mat != NULL) | 
| 761 |  |  | sprintf(matspec+strlen(matspec), | 
| 762 |  |  | " (%s)", mat->oname); | 
| 763 | greg | 2.47 | } else | 
| 764 |  |  | strcat(matspec, VOIDID); | 
| 765 |  |  | sprintf(buf, "ray hit %s %s \"%s\"", matspec, | 
| 766 | greg | 1.1 | ofun[thisray.ro->otype].funame, | 
| 767 |  |  | thisray.ro->oname); | 
| 768 | greg | 2.37 | if ((ino = objptr(thisray.robj)) != thisray.ro) | 
| 769 |  |  | sprintf(buf+strlen(buf), " in %s \"%s\"", | 
| 770 |  |  | ofun[ino->otype].funame, ino->oname); | 
| 771 | greg | 1.1 | (*dev->comout)(buf); | 
| 772 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 773 | greg | 1.1 | if (thisray.rot >= FHUGE) | 
| 774 |  |  | (*dev->comout)("at infinity"); | 
| 775 |  |  | else { | 
| 776 | greg | 2.30 | sprintf(buf, "at (%.6g %.6g %.6g) (%.6g)", | 
| 777 |  |  | thisray.rop[0], thisray.rop[1], | 
| 778 |  |  | thisray.rop[2], thisray.rt); | 
| 779 | greg | 1.1 | (*dev->comout)(buf); | 
| 780 |  |  | } | 
| 781 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 782 | greg | 2.38 | sprintf(buf, "value (%.5g %.5g %.5g) (%.3gL)", | 
| 783 | greg | 1.1 | colval(thisray.rcol,RED), | 
| 784 |  |  | colval(thisray.rcol,GRN), | 
| 785 | greg | 1.32 | colval(thisray.rcol,BLU), | 
| 786 |  |  | luminance(thisray.rcol)); | 
| 787 | greg | 1.1 | (*dev->comout)(buf); | 
| 788 |  |  | } | 
| 789 | greg | 1.18 | (*dev->comin)(buf, NULL); | 
| 790 | greg | 1.1 | } | 
| 791 |  |  |  | 
| 792 |  |  |  | 
| 793 | greg | 2.55 | void | 
| 794 |  |  | writepict(                              /* write the picture to a file */ | 
| 795 |  |  | char  *s | 
| 796 |  |  | ) | 
| 797 | greg | 1.1 | { | 
| 798 |  |  | static char  buf[128]; | 
| 799 | greg | 1.10 | char  *fname; | 
| 800 | greg | 1.1 | FILE  *fp; | 
| 801 |  |  | COLR  *scanline; | 
| 802 |  |  | int  y; | 
| 803 |  |  |  | 
| 804 | greg | 2.20 | while (isspace(*s)) | 
| 805 |  |  | s++; | 
| 806 |  |  | if (*s) | 
| 807 |  |  | atos(buf, sizeof(buf), s); | 
| 808 |  |  | else if (buf[0] == '\0') { | 
| 809 | greg | 1.1 | error(COMMAND, "no file"); | 
| 810 |  |  | return; | 
| 811 |  |  | } | 
| 812 | greg | 1.10 | if ((fname = getpath(buf, NULL, 0)) == NULL || | 
| 813 | greg | 1.1 | (fp = fopen(fname, "w")) == NULL) { | 
| 814 |  |  | sprintf(errmsg, "cannot open \"%s\"", buf); | 
| 815 |  |  | error(COMMAND, errmsg); | 
| 816 |  |  | return; | 
| 817 |  |  | } | 
| 818 | schorsch | 2.40 | SET_FILE_BINARY(fp); | 
| 819 | greg | 1.1 | (*dev->comout)("writing \""); | 
| 820 |  |  | (*dev->comout)(fname); | 
| 821 |  |  | (*dev->comout)("\"...\n"); | 
| 822 |  |  | /* write header */ | 
| 823 | greg | 2.24 | newheader("RADIANCE", fp); | 
| 824 | greg | 1.1 | fputs(progname, fp); | 
| 825 |  |  | fprintview(&ourview, fp); | 
| 826 | greg | 1.29 | if (octname != NULL) | 
| 827 |  |  | fprintf(fp, " %s\n", octname); | 
| 828 |  |  | else | 
| 829 |  |  | putc('\n', fp); | 
| 830 | greg | 1.28 | fprintf(fp, "SOFTWARE= %s\n", VersionID); | 
| 831 | greg | 2.38 | fputnow(fp); | 
| 832 | greg | 1.1 | if (exposure != 1.0) | 
| 833 | greg | 1.13 | fputexpos(exposure, fp); | 
| 834 | greg | 1.15 | if (dev->pixaspect != 1.0) | 
| 835 |  |  | fputaspect(dev->pixaspect, fp); | 
| 836 | greg | 1.22 | fputformat(COLRFMT, fp); | 
| 837 | greg | 1.11 | putc('\n', fp); | 
| 838 | greg | 1.33 | fprtresolu(hresolu, vresolu, fp); | 
| 839 | greg | 1.1 |  | 
| 840 | greg | 1.14 | scanline = (COLR *)malloc(hresolu*sizeof(COLR)); | 
| 841 | greg | 1.27 | if (scanline == NULL) { | 
| 842 |  |  | error(COMMAND, "not enough memory!"); | 
| 843 |  |  | fclose(fp); | 
| 844 |  |  | unlink(fname); | 
| 845 |  |  | return; | 
| 846 |  |  | } | 
| 847 | greg | 1.14 | for (y = vresolu-1; y >= 0; y--) { | 
| 848 |  |  | getpictcolrs(y, scanline, &ptrunk, hresolu, vresolu); | 
| 849 |  |  | if (fwritecolrs(scanline, hresolu, fp) < 0) | 
| 850 | greg | 1.1 | break; | 
| 851 |  |  | } | 
| 852 | greg | 2.38 | free((void *)scanline); | 
| 853 | greg | 1.1 | if (fclose(fp) < 0) | 
| 854 |  |  | error(COMMAND, "write error"); | 
| 855 |  |  | } |