| 125 |
|
error(SYSTEM, errmsg); |
| 126 |
|
return NULL; |
| 127 |
|
} |
| 128 |
< |
return new RdataShareFile(name, RSDOflags[op], siz); |
| 128 |
> |
RdataShare * rds = new RdataShareFile(name, RSDOflags[op], |
| 129 |
> |
siz*(op != RCOforce)); |
| 130 |
> |
|
| 131 |
> |
if (!rds || (op == RCOforce && rds->Resize(siz) < siz)) { |
| 132 |
> |
delete rds; |
| 133 |
> |
sprintf(errmsg, "cannot create %lu byte output file '%s'", |
| 134 |
> |
(unsigned long)siz, name); |
| 135 |
> |
error(SYSTEM, errmsg); |
| 136 |
> |
return NULL; |
| 137 |
> |
} |
| 138 |
> |
return rds; |
| 139 |
|
} |
| 140 |
|
|
| 141 |
|
// Memory-mapped data share function |
| 147 |
|
error(SYSTEM, errmsg); |
| 148 |
|
return NULL; |
| 149 |
|
} |
| 150 |
< |
return new RdataShareMap(name, RSDOflags[op], siz); |
| 150 |
> |
RdataShare * rds = new RdataShareMap(name, RSDOflags[op], |
| 151 |
> |
siz*(op != RCOforce)); |
| 152 |
> |
|
| 153 |
> |
if (!rds || (op == RCOforce && rds->Resize(siz) < siz)) { |
| 154 |
> |
delete rds; |
| 155 |
> |
sprintf(errmsg, "cannot create %lu byte output map '%s'", |
| 156 |
> |
(unsigned long)siz, name); |
| 157 |
> |
error(SYSTEM, errmsg); |
| 158 |
> |
return NULL; |
| 159 |
> |
} |
| 160 |
> |
return rds; |
| 161 |
|
} |
| 162 |
|
|
| 163 |
|
// Set output format ('f', 'd', or 'c') |
| 203 |
|
if (!mp) |
| 204 |
|
return 0; // not in our modifier list |
| 205 |
|
|
| 206 |
< |
int bi = 0; // get bin index |
| 206 |
> |
int i; // pre-emptive check for zero |
| 207 |
> |
if (rcp->HasFlag(RCcontrib)) { |
| 208 |
> |
for (i = NCSAMP; i--; ) |
| 209 |
> |
if (r->rcoef[i]*r->rcol[i] > FTINY) |
| 210 |
> |
break; |
| 211 |
> |
if (i < 0) |
| 212 |
> |
return 0; // zero contribution |
| 213 |
> |
} else if (sintens(r->rcoef) <= FTINY) |
| 214 |
> |
return 0; // zero coefficient |
| 215 |
> |
|
| 216 |
> |
int bi = 0; // get bin index |
| 217 |
|
if (mp->binv) { |
| 218 |
|
worldfunc(RCCONTEXT, r); |
| 219 |
|
set_eparams(mp->params); |
| 233 |
|
if (rcp->HasFlag(RCcontrib)) |
| 234 |
|
smultscolor(contr, r->rcol); // -> value contribution |
| 235 |
|
|
| 236 |
< |
for (int i = 0; i < NCSAMP; i++) |
| 236 |
> |
for (i = 0; i < NCSAMP; i++) |
| 237 |
|
*dvp++ += contr[i]; // accumulate color/spectrum |
| 238 |
|
return 1; |
| 239 |
|
} |
| 479 |
|
sprintf(hdr+begData, "%s%d\n", NCOMPSTR, NCSAMP); |
| 480 |
|
begData += strlen(hdr+begData); |
| 481 |
|
if (NCSAMP > 3) { |
| 482 |
< |
sprintf(hdr+begData, "%s %f %f %f %f\n", WLSPLTSTR, |
| 482 |
> |
sprintf(hdr+begData, "%s %g %g %g %g\n", WLSPLTSTR, |
| 483 |
|
WLPART[0], WLPART[1], WLPART[2], WLPART[3]); |
| 484 |
|
begData += strlen(hdr+begData); |
| 485 |
|
} |