Commit 4ee8a49c authored by Johan Richard's avatar Johan Richard
Browse files

Corrected identations and cleanset 4 output

parent 86c07ad3
......@@ -46,8 +46,7 @@ void do_itos(double **im, struct pixlist *pl, int npl, double dlsds, double z
/*According to the subsampling in the image plane
*for each subpixel in the image plane*/
for (k = 0; k < npl; k++)
{
{
for (i = 0; i < imFrame.ech; i++)
for (j = 0; j < imFrame.ech; j++)
{
......
......@@ -161,7 +161,6 @@ void cubetosou(double zimage, char *sname )
*/
dlsds = dratio(lens[0].z, zimage);
// NPRINTF("%d%s", &cubeFrame.format, cubeFrame.pixfile,);
NPRINTF(stderr, "READ: cubeframe\n");
/* Allocate the maximum number of points to send to source plane */
pl = (struct pixlist*)malloc((unsigned int)
......@@ -215,21 +214,22 @@ void cubetosou(double zimage, char *sname )
cubeFrame.nx, cubeFrame.ny, ps.nx, ps.ny);
for (kk=0 ; kk < cubeFrame.nz; kk ++)
{ fprintf(stderr, "cubeFrame.nz= %d\n", kk); /* printing the lambda-frame for which source is being */
{
NPRINTF(stderr, "cubeFrame.nz= %d\n", kk); /* printing the lambda-frame for which source is being */
for (j = 0; j < cubeFrame.ny; j++) /* reconstructed */
for (j = 0; j < cubeFrame.ny; j++) /* reconstructed */
{
for (k = 0; k < cubeFrame.nx; k++)
{
ima_temp[j][k]=cubeo[j][k][kk]; /* writing 3d array into 2D array per wavelength frame */
}
for (k = 0; k < cubeFrame.nx; k++)
ima_temp[j][k]=cubeo[j][k][kk]; /* writing 3d array into 2D array per wavelength frame */
}
/* keep only the relevant points */
if (imFrame.ncont == 0)
s_pixlist(ima_temp, pl, &npl);
s_pixlist(ima_temp, pl, &npl);
else
keep_cl(ima_temp, pl, &npl);
keep_cl(ima_temp, pl, &npl);
do_itos(ima_temp, pl, npl, dlsds, zimage, source, erreur, imult);
imFrame.wcsinfo = cubeFrame.wcsinfo;
......
......@@ -301,12 +301,14 @@ void setBayesModel( long int iVal, long int nVal, double **array)
// cleanset = 4 case
if (M.iclean == 4)
{
{
extern struct g_pixel imFrame;
int ns = ps.nx*ps.ny;
imethod = -2;
imethod = -1;
fprintf(stderr,"imethod=%d\n",imethod);
for ( ilens = 0; ilens < ns; ilens++ )
sflux[ilens] = method(nVal, iParam++, array);
testpixel("testpixel_method-2.fits","sflux_method-2.fits",sflux,ero);
testpixel("image_avg.fits","source_avg.fits",sflux,ero);
}
// Set the cosmological parameters
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment