Changeset 923 in Sophya for trunk/SophyaExt
- Timestamp:
- Apr 13, 2000, 7:34:37 PM (25 years ago)
- Location:
- trunk/SophyaExt/FitsIOServer
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaExt/FitsIOServer/fitsfile.cc
r903 r923 568 568 taille_des_chaines.clear(); 569 569 char **ttype = new char*[nbcols]; 570 for (int ii=0; ii < nbcols; ii++) ttype[ii]=new char[FLEN_VALUE]; 570 int ii; 571 for (ii=0; ii < nbcols; ii++) ttype[ii]=new char[FLEN_VALUE]; 571 572 int nfound; 572 573 fits_read_keys_str(fileptr, "TTYPE",1,nbcols,ttype,&nfound, &status); 573 574 if( status ) printerror( status,"erreur lecture des noms de colonne"); 574 575 int rept=0; 575 for(i nt ii = 0; ii < nbcols; ii++)576 for(ii = 0; ii < nbcols; ii++) 576 577 { 577 578 int DTYPE; -
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, -
trunk/SophyaExt/FitsIOServer/fitsntuple.cc
r873 r923 57 57 nbcols = fn.NbColsFromFits(); 58 58 nbentries = 0; 59 for (int k=0; k<nbcols; k++) nbentries=max( nbentries, fn.NentriesFromFits(k) ); 59 int k; 60 for (k=0; k<nbcols; k++) nbentries=max( nbentries, fn.NentriesFromFits(k) ); 60 61 61 62 char ** ColName = new char*[nbcols]; 62 63 63 for ( intk=0; k<nbcols;k++)64 for (k=0; k<nbcols;k++) 64 65 { 65 66 ColName[k] = new char[LONNOM1]; … … 67 68 ColName[k][LONNOM] = '\0'; 68 69 } 69 for ( intk=0; k<nbcols;k++)70 for (k=0; k<nbcols;k++) 70 71 { 71 72 char ss= fn.ColTypeFromFits(k); … … 93 94 (*dobj_) = NTuple(nbcols,ColName); 94 95 } 95 for ( intk=0; k<nbcols;k++)96 for (k=0; k<nbcols;k++) 96 97 { 97 98 delete [] ColName[k]; … … 104 105 // (SetVal suppose que le ntuple est deja dimensionne) 105 106 r_4* ligne = new r_4[nbcols]; 106 for ( intk=0; k<nbcols; k++) ligne[k]=0.;107 for ( intk=0; k<nbentries;k++) dobj_->Fill(ligne);107 for (k=0; k<nbcols; k++) ligne[k]=0.; 108 for (k=0; k<nbentries;k++) dobj_->Fill(ligne); 108 109 delete [] ligne; 109 for ( intk=0; k<nbcols;k++)110 for (k=0; k<nbcols;k++) 110 111 { 111 112 fn.GetBinTabFCol(column_, nbentries, k); … … 135 136 136 137 char** Noms = new char*[ncols]; 137 for (int k=0; k< ncols; k++) 138 int k; 139 for (k=0; k< ncols; k++) 138 140 { 139 141 Noms[k]= new char[LONNOM1]; … … 142 144 // la librairie fitsio ecrit colonne par colonne 143 145 char* type= new char[ncols+1]; 144 for ( intk=0;k<ncols+1;k++) type[k]='E';146 for (k=0;k<ncols+1;k++) type[k]='E'; 145 147 type[ncols]='\0'; 146 148 vector<int> dummy; 147 149 fn.makeHeaderBntblOnFits(type,Noms, nentries, ncols, dvl, extname, dummy); 148 for ( intk=0; k< ncols; k++)150 for (k=0; k< ncols; k++) 149 151 { 150 152 delete [] Noms[k]; … … 152 154 delete [] Noms; 153 155 delete [] type; 154 for ( intk=0; k<ncols;k++) putColToFits(k, nentries, getColFromObj(k));156 for (k=0; k<ncols;k++) putColToFits(k, nentries, getColFromObj(k)); 155 157 156 158 } -
trunk/SophyaExt/FitsIOServer/fitstarray.cc
r874 r923 101 101 cout << "FITS_TARRAY: nombre de dimension du tableau a ecrire: " << nbdim << endl; 102 102 int* naxisn = new int[nbdim]; 103 for (int k=0; k< nbdim; k++) 103 int k; 104 for (k=0; k< nbdim; k++) 104 105 { 105 106 naxisn[k] = dobj_->Size(k); … … 119 120 120 121 int nbels = 1; 121 for ( intk=0; k< nbdim; k++)122 for (k=0; k< nbdim; k++) 122 123 { 123 124 if (naxisn[k] > 0) nbels *= naxisn[k]; -
trunk/SophyaExt/FitsIOServer/fitsxntuple.cc
r873 r923 74 74 nbcols = fn.NbColsFromFits(); 75 75 nbentries = 0; 76 for (int k=0; k<nbcols; k++) nbentries=max( nbentries, fn.NentriesFromFits(k) ); 76 int k; 77 for (k=0; k<nbcols; k++) nbentries=max( nbentries, fn.NentriesFromFits(k) ); 77 78 78 79 // … … 88 89 vector<int> IfitsCol; 89 90 vector<int> SfitsCol; 90 for ( intk=0; k<nbcols;k++)91 for (k=0; k<nbcols;k++) 91 92 { 92 93 char ss= fn.ColTypeFromFits(k); … … 102 103 char ** ColName = new char*[nbcols]; 103 104 int compt=0; 104 for ( intk=0; k<DfitsCol.size(); k++)105 for (k=0; k<DfitsCol.size(); k++) 105 106 { 106 107 ColName[compt] = new char[LONNOM+1]; … … 108 109 ColName[compt++][ LONNOM] = '\0'; 109 110 } 110 for ( intk=0; k<FfitsCol.size(); k++)111 for (k=0; k<FfitsCol.size(); k++) 111 112 { 112 113 ColName[compt] = new char[LONNOM+1]; … … 114 115 ColName[compt++][ LONNOM] = '\0'; 115 116 } 116 for ( intk=0; k<IfitsCol.size(); k++)117 for (k=0; k<IfitsCol.size(); k++) 117 118 { 118 119 ColName[compt] = new char[LONNOM+1]; … … 120 121 ColName[compt++][ LONNOM] = '\0'; 121 122 } 122 for ( intk=0; k<SfitsCol.size(); k++)123 for (k=0; k<SfitsCol.size(); k++) 123 124 { 124 125 ColName[compt] = new char[LONNOM+1]; … … 138 139 (*dobj_)= XNTuple(DfitsCol.size(), FfitsCol.size(), IfitsCol.size(), SfitsCol.size(),ColName); 139 140 } 140 for ( intk=0; k<nbcols;k++)141 for (k=0; k<nbcols;k++) 141 142 { 142 143 delete [] ColName[k]; … … 155 156 dcolumn_ = new double[nbentries]; 156 157 dligne = new r_8[DfitsCol.size()]; 157 for ( intk=0; k<DfitsCol.size(); k++) dligne[k]=0.;158 for (k=0; k<DfitsCol.size(); k++) dligne[k]=0.; 158 159 } 159 160 else dligne=NULL; … … 163 164 fcolumn_ = new float[nbentries]; 164 165 fligne = new r_4[FfitsCol.size()]; 165 for ( intk=0; k<FfitsCol.size(); k++) fligne[k]=0.;166 for (k=0; k<FfitsCol.size(); k++) fligne[k]=0.; 166 167 } 167 168 else fligne=NULL; … … 171 172 icolumn_ = new int[nbentries]; 172 173 iligne = new int_4[IfitsCol.size()]; 173 for ( intk=0; k<IfitsCol.size(); k++) iligne[k]=0;174 for (k=0; k<IfitsCol.size(); k++) iligne[k]=0; 174 175 } 175 176 else iligne=NULL; … … 178 179 if (ccolumn_ != NULL) 179 180 { 180 for ( intk=0; k<dobj_->NEntry(); k++) delete [] ccolumn_[k];181 for (k=0; k<dobj_->NEntry(); k++) delete [] ccolumn_[k]; 181 182 delete [] ccolumn_; 182 183 ccolumn_ = NULL; … … 184 185 ccolumn_ = new char*[nbentries]; 185 186 int taille_des_chaines=0; 186 for ( intk=0; k< SfitsCol.size(); k++) taille_des_chaines = max( taille_des_chaines, fn.ColStringLengthFromFits(SfitsCol[k]) );187 for ( intk=0; k<nbentries;k++) ccolumn_[k]= new char[taille_des_chaines+1];187 for (k=0; k< SfitsCol.size(); k++) taille_des_chaines = max( taille_des_chaines, fn.ColStringLengthFromFits(SfitsCol[k]) ); 188 for (k=0; k<nbentries;k++) ccolumn_[k]= new char[taille_des_chaines+1]; 188 189 cligne = new char*[SfitsCol.size()]; 189 for ( intk=0; k<SfitsCol.size(); k++) cligne[k]=" ";190 for (k=0; k<SfitsCol.size(); k++) cligne[k]=" "; 190 191 } 191 192 else cligne=NULL; 192 193 193 for ( intk=0; k<nbentries;k++) dobj_->Fill(dligne, fligne, iligne, cligne);194 for (k=0; k<nbentries;k++) dobj_->Fill(dligne, fligne, iligne, cligne); 194 195 delete [] dligne; 195 196 delete [] fligne; 196 197 delete [] iligne; 197 for ( intk=0; k< SfitsCol.size(); k++) delete [] cligne[k];198 for (k=0; k< SfitsCol.size(); k++) delete [] cligne[k]; 198 199 delete [] cligne; 199 200 200 201 compt=0; 201 for ( intk=0; k<DfitsCol.size(); k++)202 for (k=0; k<DfitsCol.size(); k++) 202 203 { 203 204 fn.GetBinTabFCol(dcolumn_, nbentries, DfitsCol[k]); … … 206 207 } 207 208 208 for ( intk=0; k<FfitsCol.size(); k++)209 for (k=0; k<FfitsCol.size(); k++) 209 210 { 210 211 fn.GetBinTabFCol(fcolumn_,nbentries, FfitsCol[k]); … … 213 214 } 214 215 215 for ( intk=0; k<IfitsCol.size(); k++)216 for (k=0; k<IfitsCol.size(); k++) 216 217 { 217 218 fn.GetBinTabFCol(icolumn_,nbentries, IfitsCol[k]); … … 220 221 } 221 222 222 for ( intk=0; k<SfitsCol.size(); k++)223 for (k=0; k<SfitsCol.size(); k++) 223 224 { 224 225 fn.GetBinTabFCol(ccolumn_,nbentries, SfitsCol[k]); … … 247 248 dvl.Print(); 248 249 char** Noms = new char*[ncols]; 249 for (int k=0; k< ncols; k++) 250 int k; 251 for (k=0; k< ncols; k++) 250 252 { 251 253 Noms[k]= new char[LONNOM+1]; … … 254 256 char* types=new char[ncols+1]; 255 257 int compt=0; 256 for ( intk=0; k<dobj_->NDVar();k++)258 for (k=0; k<dobj_->NDVar();k++) 257 259 { 258 260 types[compt++]='D'; 259 261 } 260 for ( intk=0; k<dobj_->NFVar();k++)262 for (k=0; k<dobj_->NFVar();k++) 261 263 { 262 264 types[compt++]='E'; 263 265 } 264 for ( intk=0; k<dobj_->NIVar();k++)266 for (k=0; k<dobj_->NIVar();k++) 265 267 { 266 268 types[compt++]='I'; 267 269 } 268 for ( intk=0; k<dobj_->NSVar();k++)270 for (k=0; k<dobj_->NSVar();k++) 269 271 { 270 272 types[compt++]='A'; … … 272 274 types[ncols]='\0'; 273 275 vector<int> StringSizes(dobj_->NSVar()); 274 for ( intk=0; k< StringSizes.size(); k++) StringSizes[k]=dobj_->mStrSz;276 for (k=0; k< StringSizes.size(); k++) StringSizes[k]=dobj_->mStrSz; 275 277 // la librairie fitsio ecrit colonne par colonne 276 278 fn.makeHeaderBntblOnFits(types, Noms, nrows, ncols, dvl, extname,StringSizes); 277 for ( intk=0; k< ncols; k++)279 for (k=0; k< ncols; k++) 278 280 { 279 281 delete [] Noms[k]; … … 282 284 delete [] types; 283 285 compt=0; 284 for ( intk=0; k<dobj_->NDVar();k++)286 for (k=0; k<dobj_->NDVar();k++) 285 287 { 286 288 putColToFits(compt, dobj_->NEntry(), getColDFromObj(compt)); 287 289 compt++; 288 290 } 289 for ( intk=0; k<dobj_->NFVar();k++)291 for (k=0; k<dobj_->NFVar();k++) 290 292 { 291 293 putColToFits(compt, dobj_->NEntry(), getColFFromObj(compt)); 292 294 compt++; 293 295 } 294 for ( intk=0; k<dobj_->NIVar();k++)296 for (k=0; k<dobj_->NIVar();k++) 295 297 { 296 298 putColToFits(compt, dobj_->NEntry(), getColIFromObj(compt)); 297 299 compt++; 298 300 } 299 for ( intk=0; k<dobj_->NSVar();k++)301 for (k=0; k<dobj_->NSVar();k++) 300 302 { 301 303 putColToFits(compt, dobj_->NEntry(), getColSFromObj(compt));
Note:
See TracChangeset
for help on using the changeset viewer.