122 |
|
il->matname, DSTSUF); |
123 |
|
printf("\n5 noneg"); |
124 |
|
dfp = dfopen(il, 0); |
125 |
< |
fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", 1.+.5/n, .5/n, n+1, |
125 |
> |
fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", |
126 |
> |
1.+.5/n, .5/n, n+1, |
127 |
|
0., 2.*PI, m+1); |
128 |
|
brightout(Ninv, 1, m, 1./il->nsamps/brt(il->col), dfp); |
129 |
|
brightout(da, n, m, 1./il->nsamps/brt(il->col), dfp); |
152 |
|
|
153 |
|
if ((Ninv = (float *)malloc(3*m*sizeof(float))) == NULL || |
154 |
|
(Sinv = (float *)malloc(3*m*sizeof(float))) == NULL) |
155 |
< |
error(SYSTEM, "out of memory in flatout"); |
155 |
> |
error(SYSTEM, "out of memory in roundout"); |
156 |
|
compinv(Ninv, da, m); |
157 |
|
compinv(Sinv, da+3*m*(n-1), m); |
158 |
|
if (il->flags & IL_COLDST) { |
176 |
|
il->matname, DSTSUF); |
177 |
|
printf("\n5 noneg"); |
178 |
|
dfp = dfopen(il, 0); |
179 |
< |
fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", 1.-1./n, -1.+1./n, n, |
179 |
> |
fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", |
180 |
> |
1.+1./n, -1.-1./n, n+2, |
181 |
|
0., 2.*PI, m+1); |
182 |
|
brightout(Ninv, 1, m, 1./il->nsamps/brt(il->col), dfp); |
183 |
|
brightout(da, n, m, 1./il->nsamps/brt(il->col), dfp); |