Changeset 923 in Sophya for trunk/SophyaExt/FitsIOServer/fitsioserver.cc
- Timestamp:
- Apr 13, 2000, 7:34:37 PM (25 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaExt/FitsIOServer/fitsioserver.cc
r854 r923 329 329 long repeat; 330 330 int ntest= 0; 331 for(int ii = 0; ii < ncols; ii++) 331 int ii; 332 for(ii = 0; ii < ncols; ii++) 332 333 { 333 334 int DTYPE; … … 354 355 char **clname; 355 356 clname= new char*[tfields.size()]; 356 for(i nt ii = 0; ii < tfields.size(); ii++)357 for(ii = 0; ii < tfields.size(); ii++) 357 358 { 358 359 // ostringstream ne passe pas sur tous les compilos .. Reza 23/12/99 … … 367 368 // create a NTuple 368 369 NTuple nt0(tfields.size(),clname); 369 for(i nt ii = 0; ii < tfields.size(); ii++)370 for(ii = 0; ii < tfields.size(); ii++) 370 371 delete [] clname[ii]; 371 372 delete [] clname; … … 615 616 } 616 617 } 617 618 for( intj = 0; j < sph1.NbPixels(); j++) sph1(j)= (double)r_8tab_[j];619 for( intj = 0; j < sph2.NbPixels(); j++) sph2(j)= (double)r_8tab_[j+npixels];620 for( intj = 0; j < sph3.NbPixels(); j++) sph3(j)= (double)r_8tab_[j+2*npixels];618 int j; 619 for(j = 0; j < sph1.NbPixels(); j++) sph1(j)= (double)r_8tab_[j]; 620 for(j = 0; j < sph2.NbPixels(); j++) sph2(j)= (double)r_8tab_[j+npixels]; 621 for(j = 0; j < sph3.NbPixels(); j++) sph3(j)= (double)r_8tab_[j+2*npixels]; 621 622 } 622 623 … … 656 657 } 657 658 } 658 659 for( intj = 0; j < sph1.NbPixels(); j++) sph1(j)= (float)r_4tab_[j];660 for( intj = 0; j < sph2.NbPixels(); j++) sph2(j)= (float)r_4tab_[j+npixels];661 for( intj = 0; j < sph3.NbPixels(); j++) sph3(j)= (float)r_4tab_[j+2*npixels];659 int j; 660 for(j = 0; j < sph1.NbPixels(); j++) sph1(j)= (float)r_4tab_[j]; 661 for(j = 0; j < sph2.NbPixels(); j++) sph2(j)= (float)r_4tab_[j+npixels]; 662 for(j = 0; j < sph3.NbPixels(); j++) sph3(j)= (float)r_4tab_[j+2*npixels]; 662 663 } 663 664 … … 950 951 printerror( status ); 951 952 952 for( 953 for(int ii = 0; ii < tfields; ii++) 953 954 { 954 955 delete [] ttype[ii]; … … 1119 1120 r_8tab_=new r_8[3*npixels]; 1120 1121 1121 1122 for ( intj = 0; j < npixels; j++) r_8tab_[j]= (r_8)sph1(j);1123 for ( intj = 0; j < npixels; j++) r_8tab_[j+npixels]= (r_8)sph2(j);1124 for ( intj = 0; j < npixels; j++) r_8tab_[j+2*npixels]= (r_8)sph3(j);1122 int j; 1123 for (j = 0; j < npixels; j++) r_8tab_[j]= (r_8)sph1(j); 1124 for (j = 0; j < npixels; j++) r_8tab_[j+npixels]= (r_8)sph2(j); 1125 for (j = 0; j < npixels; j++) r_8tab_[j+2*npixels]= (r_8)sph3(j); 1125 1126 DVList dvl; 1126 1127 dvl.SetS("PIXTYPE","HEALPIX "); … … 1150 1151 if (r_4tab_ != NULL ) { delete[] r_4tab_; r_4tab_ = NULL; } 1151 1152 r_4tab_=new r_4[3*npixels]; 1152 1153 for (int j = 0; j < npixels; j++) r_4tab_[j]= (r_4)sph1(j); 1154 for (int j = 0; j < npixels; j++) r_4tab_[j+npixels]= (r_4)sph2(j); 1155 for (int j = 0; j < npixels; j++) r_4tab_[j+2*npixels]= (r_4)sph3(j); 1153 1154 int j; 1155 for (j = 0; j < npixels; j++) r_4tab_[j]= (r_4)sph1(j); 1156 for (j = 0; j < npixels; j++) r_4tab_[j+npixels]= (r_4)sph2(j); 1157 for (j = 0; j < npixels; j++) r_4tab_[j+2*npixels]= (r_4)sph3(j); 1156 1158 1157 1159 DVList dvl; … … 1534 1536 delete [] tform; 1535 1537 long nelements= npixels; 1538 int k; 1536 1539 switch (FITS_tab_typ_) 1537 1540 { 1538 1541 case TDOUBLE : 1539 for ( intk=1; k<=tfields; k++)1542 for (k=1; k<=tfields; k++) 1540 1543 { 1541 1544 fits_write_col(fptr, TDOUBLE, k, 1, 1, nelements, r_8tab_+(k-1)*nelements, &status); … … 1547 1550 case TFLOAT : 1548 1551 //DBG!!! $CHECK$ Reza for (int kk=0; kk<10; kk++) cout << r_4tab_[kk] << endl; 1549 for ( intk=1; k<=tfields; k++)1552 for (k=1; k<=tfields; k++) 1550 1553 { 1551 1554 fits_write_col(fptr, TFLOAT, k, 1, 1, nelements, r_4tab_+(k-1)*nelements, &status); … … 1554 1557 break; 1555 1558 case TINT : 1556 for ( intk=1; k<=tfields; k++)1559 for (k=1; k<=tfields; k++) 1557 1560 { 1558 1561 fits_write_col(fptr, TINT, k, 1, 1, nelements, i_4tab_+(k-1)*nelements, &status); … … 1869 1872 // le type suppose par l'utilisateur de fitsioserver 1870 1873 // 1874 int k; 1871 1875 switch ( FITS_tab_typ_) 1872 1876 { … … 1879 1883 if (r_8tab_ != NULL) { delete [] r_8tab_; r_8tab_ = NULL; } 1880 1884 r_8tab_=new r_8[ tfields*npixels]; 1881 for ( intk=1; k<=tfields; k++)1885 for (k=1; k<=tfields; k++) 1882 1886 { 1883 1887 fits_read_col(fptr, TDOUBLE, k, 1, 1, npixels, &dnullval, … … 1896 1900 if (r_4tab_ != NULL) { delete [] r_4tab_; r_4tab_ = NULL; } 1897 1901 r_4tab_=new r_4[ tfields*npixels]; 1898 for ( intk=1; k<=tfields; k++)1902 for (k=1; k<=tfields; k++) 1899 1903 { 1900 1904 fits_read_col(fptr, TFLOAT, k, 1, 1, npixels, &fnullval, … … 1913 1917 if (i_4tab_ != NULL) { delete [] i_4tab_; i_4tab_ = NULL; } 1914 1918 i_4tab_=new int_4[ tfields*npixels]; 1915 for ( intk=1; k<=tfields; k++)1919 for (k=1; k<=tfields; k++) 1916 1920 { 1917 1921 fits_read_col(fptr, TLONG, k, 1, 1, npixels, &inullval,
Note:
See TracChangeset
for help on using the changeset viewer.