Changeset 873 in Sophya
- Timestamp:
- Apr 11, 2000, 10:01:21 AM (25 years ago)
- Location:
- trunk/SophyaExt/FitsIOServer
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaExt/FitsIOServer/fitsntuple.cc
r860 r873 56 56 int nbcols, nbentries; 57 57 nbcols = fn.NbColsFromFits(); 58 cout << " nbcols= " << nbcols << endl;59 58 nbentries = 0; 60 59 for (int k=0; k<nbcols; k++) nbentries=max( nbentries, fn.NentriesFromFits(k) ); … … 99 98 } 100 99 delete [] ColName; 101 cout << " preparation du tabeau column" << endl;102 100 if (column_ != NULL) delete [] column_; 103 101 column_ = new float[nbentries]; … … 105 103 // j'initialise le NTuple a zero, pour le dimensionner 106 104 // (SetVal suppose que le ntuple est deja dimensionne) 107 cout << " nbcol= " << nbcols << endl;108 105 r_4* ligne = new r_4[nbcols]; 109 106 for (int k=0; k<nbcols; k++) ligne[k]=0.; 110 107 for (int k=0; k<nbentries;k++) dobj_->Fill(ligne); 111 cout << " ntuple initialise a zero" << endl;112 108 delete [] ligne; 113 cout << "debut boucle nbcol= " << nbcols << endl;114 109 for (int k=0; k<nbcols;k++) 115 110 { 116 cout << " avant lecture col. no " << k << endl;117 111 fn.GetBinTabFCol(column_, nbentries, k); 118 cout << " lecture effectuee " << endl;119 112 for (int nent=0; nent<nbentries; nent++) dobj_->SetVal(nent,k, column_[nent]); 120 cout << " valeurs inserre dans ntuple " << endl;121 113 } 122 114 dobj_->Info()=fn.DVListFromFits(); … … 147 139 Noms[k]= new char[LONNOM1]; 148 140 strncpy(Noms[k],dobj_->NomIndex(k),LONNOM1); 149 cout << " nom entre dans ntuple: " << Noms[k] << endl;150 141 } 151 142 // la librairie fitsio ecrit colonne par colonne … … 163 154 for (int k=0; k<ncols;k++) putColToFits(k, nentries, getColFromObj(k)); 164 155 165 cout << " fin d'ecriture "<< endl;166 156 } 167 157 -
trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc
r854 r873 96 96 strncpy(extname,dvl.GetS("EXTNAME").c_str(),15); 97 97 98 cout << " nom entre dans sphereHealpix: " << Noms[0] << endl;99 98 char Type[2]; 100 99 if (strncmp( typeid(T).name(),"r_8",3)==0 || strncmp( typeid(T).name(),"double",6)==0 ) Type[0]='D'; … … 133 132 // const DVList* dvl = &fn.DVListFromFits(); 134 133 DVList dvl=fn.DVListFromFits(); 135 dvl.Print();134 // dvl.Print(); 136 135 nbentries = fn.NentriesFromFits(0); 137 136 int lastpix=dvl.GetI("LASTPIX"); … … 174 173 if (nPix != nbentries) 175 174 { 176 cout << " le nombre de pixels relu " << nbentries << " est incoherent avec la pixelisation Healpix qui donne nPix= " << nPix << endl;175 cout << "WARNING: le nombre de pixels relu " << nbentries << " est incoherent avec la pixelisation Healpix qui donne nPix= " << nPix << endl; 177 176 } 178 177 double Omega= 4.0*Pi/nPix; 179 178 dobj_->setParameters(nside,nPix,Omega); 180 cout << " sphere Healpix : las parametres sont introduits" << endl;181 179 // On lit les DataBlocks; 182 cout << "nside= " << nside << "nPix= " << nPix << endl;183 180 dobj_->pixels_.ReSize(nPix); 184 cout << " sphere Healpix : resize fait" << endl;185 181 fn.GetSingleColumn(dobj_->pixels_.Data(),nPix); 186 182 -
trunk/SophyaExt/FitsIOServer/fitstarray.cc
r863 r873 99 99 throw PException(" FITS_TARRAY: number of dimensions of the array = 0?"); 100 100 } 101 cout << "FITS_TARRAY: nombre de dimension du tableau a ecrire: " << nbdim << endl;101 cout << "FITS_TARRAY: nombre de dimension du tableau a ecrire: " << nbdim << endl; 102 102 int* naxisn = new int[nbdim]; 103 103 for (int k=0; k< nbdim; k++) 104 104 { 105 105 naxisn[k] = dobj_->Size(k); 106 cout << " taille dela dimension " << k << " : " << naxisn[k] << endl;106 cout << " nombre de donnees dans la la dimension " << k << " : " << naxisn[k] << endl; 107 107 } 108 108 char type; -
trunk/SophyaExt/FitsIOServer/fitsxntuple.cc
r860 r873 73 73 int nbcols, nbentries; 74 74 nbcols = fn.NbColsFromFits(); 75 cout << " nbcols= " << nbcols << endl;76 75 nbentries = 0; 77 76 for (int k=0; k<nbcols; k++) nbentries=max( nbentries, fn.NentriesFromFits(k) ); 78 cout << " nbentries= " << nbentries << endl;79 77 80 78 // … … 148 146 // j'initialise le NTuple a zero, pour le dimensionner 149 147 // (SetXVal suppose que le ntuple est deja dimensionne) 150 cout << " nbcol= " << nbcols << endl;151 148 r_8* dligne; 152 149 r_4* fligne; … … 155 152 if (DfitsCol.size()>0) 156 153 { 157 cout << " preparation du tabeau dcolumn" << endl;158 154 if (dcolumn_ != NULL) delete [] dcolumn_; 159 155 dcolumn_ = new double[nbentries]; … … 164 160 if (FfitsCol.size()>0) 165 161 { 166 cout << " preparation du tabeau fcolumn" << endl;167 162 if (fcolumn_ != NULL) delete [] fcolumn_; 168 163 fcolumn_ = new float[nbentries]; … … 173 168 if (IfitsCol.size()>0) 174 169 { 175 cout << " preparation du tabeau icolumn" << endl;176 170 if (icolumn_ != NULL) delete [] icolumn_; 177 171 icolumn_ = new int[nbentries]; … … 182 176 if (SfitsCol.size()>0) 183 177 { 184 cout << " preparation du tabeau ccolumn" << endl;185 178 if (ccolumn_ != NULL) 186 179 { … … 192 185 int taille_des_chaines=0; 193 186 for (int k=0; k< SfitsCol.size(); k++) taille_des_chaines = max( taille_des_chaines, fn.ColStringLengthFromFits(SfitsCol[k]) ); 194 cout << " xntu: je dimensionne a taille_des_chaines= " << taille_des_chaines << endl;195 187 for (int k=0; k<nbentries;k++) ccolumn_[k]= new char[taille_des_chaines+1]; 196 188 cligne = new char*[SfitsCol.size()]; … … 200 192 201 193 for (int k=0; k<nbentries;k++) dobj_->Fill(dligne, fligne, iligne, cligne); 202 cout << " ntuple initialise a zero" << endl;203 194 delete [] dligne; 204 195 delete [] fligne; … … 232 223 { 233 224 fn.GetBinTabFCol(ccolumn_,nbentries, SfitsCol[k]); 234 cout << " chaine recue par xntuple " << endl;235 225 for (int nent=0; nent<nbentries; nent++) dobj_->SetSVal(nent,compt, ccolumn_[nent]); 236 226 compt++; 237 227 } 238 cout << " valeurs inserees dans xntuple " << endl;239 228 dobj_->Info()=fn.DVListFromFits(); 240 229 … … 262 251 Noms[k]= new char[LONNOM+1]; 263 252 strncpy(Noms[k],dobj_->NomIndex(k).c_str(),LONNOM+1); 264 cout << " nom entre dans ntuple: " << Noms[k] << endl;265 253 } 266 254 char* types=new char[ncols+1]; … … 314 302 compt++; 315 303 } 316 cout << " fin d'ecriture "<< endl; 317 } 318 /* 319 void* FITS_XNTuple::getColFromObj(int colNr) 320 { 321 if (colNr < dobj_->NDVar()) return (void*)getColDFromObj(colNr); 322 else 323 if (colNr < dobj_->NDVar()+dobj_->NFVar()) return (void*)getColFFromObj(colNr); 324 else 325 if (colNr < dobj_->NDVar()+dobj_->NFVar()+dobj_->NIVar()) return (void*)getColIFromObj(colNr); 326 else return (void*)getColSFromObj(colNr); 327 328 } 329 */ 304 } 330 305 double* FITS_XNTuple::getColDFromObj(int colNr) 331 306 { … … 377 352 ccolumn_[j] = new char[dobj_->mStrSz+1]; 378 353 strcpy(ccolumn_[j],s.c_str()); 379 // ccolumn_[j][strlen(ccolumn_[j])-1]='\0';380 cout << " chaine envoyee xntuple: " << ccolumn_[j] << " longueur " << strlen( ccolumn_[j]) << endl;381 354 } 382 355 return ccolumn_;
Note:
See TracChangeset
for help on using the changeset viewer.