- Timestamp:
- 02/11/05 17:50:42 (20 years ago)
- Location:
- trunk/src
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/Makefile
r356 r414 8 8 9 9 CXXFLAGS := -O -fPIC 10 #CXXFLAGS += -fpermissive # This should be removed. 10 11 11 #CXXFLAGS += -ansi 12 12 CXXFLAGS += -pedantic -Wno-long-long -
trunk/src/MathUtils.h
r384 r414 39 39 // Add new entry to list 40 40 template <class T> 41 casa::uInt addEntry 41 casa::uInt addEntry(casa::Vector<T>& list, T val); 42 42 43 43 // Find the scan boundaries from a list of ScanIDs … … 53 53 54 54 // Generate specified statistic 55 float statistics (const casa::String& which, const casa::MaskedArray<casa::Float>& data); 55 float statistics(const casa::String& which, 56 const casa::MaskedArray<casa::Float>& data); 56 57 57 58 // Replace masked value by zero 58 void replaceMaskByZero (casa::Vector<casa::Float>& data, const casa::Vector<casa::Bool>& mask); 59 void replaceMaskByZero(casa::Vector<casa::Float>& data, 60 const casa::Vector<casa::Bool>& mask); 59 61 }; 60 62 -
trunk/src/SDAsciiWriter.cc
r215 r414 175 175 // 176 176 of.close(); 177 c err<< "Wrote " << nRows << " rows into file " << fileName << endl;177 cout << "Wrote " << nRows << " rows into file " << fileName << endl; 178 178 // 179 179 return True; -
trunk/src/SDContainer.cc
r396 r414 520 520 cerr << "Source ID" << endl; 521 521 for (uInt i=0; i<n_; i++) { 522 c err<< setw(11) << source_(i) << ID_(i) << endl;522 cout << setw(11) << source_(i) << ID_(i) << endl; 523 523 } 524 524 } -
trunk/src/SDFITSImageWriter.cc
r389 r414 117 117 Directory dir(dirName2); 118 118 dir.create(True); 119 c err<< "Created directory '" << dirName2 << "' for output files" << endl;119 cout << "Created directory '" << dirName2 << "' for output files" << endl; 120 120 121 121 // Temps … … 229 229 String tS(oss); 230 230 String fileName = dirName2 + String("/") + tS; 231 if (verbose) c err<< "Writing row " << iRow231 if (verbose) cout << "Writing row " << iRow 232 232 << " into file '" << tS << "'" << endl; 233 233 // … … 249 249 // 250 250 if (!verbose) { 251 c err<< "Wrote " << nRows << " into individual FITS files" << endl;251 cout << "Wrote " << nRows << " into individual FITS files" << endl; 252 252 } 253 253 // -
trunk/src/SDMath.cc
r410 r414 1089 1089 Double cFac = 1.0; 1090 1090 if (fluxUnit==JY) { 1091 c err<< "Converting to K" << endl;1091 cout << "Converting to K" << endl; 1092 1092 // 1093 1093 Quantum<Double> t(1.0,fluxUnit); … … 1098 1098 sh.fluxunit = "K"; 1099 1099 } else if (fluxUnit==K) { 1100 c err<< "Converting to Jy" << endl;1100 cout << "Converting to Jy" << endl; 1101 1101 // 1102 1102 Quantum<Double> t(1.0,fluxUnit); … … 1121 1121 if (toKelvin) factor = 1.0 / JyPerK; 1122 1122 // 1123 c err<< "Applying supplied conversion factor = " << factor << endl;1123 cout << "Applying supplied conversion factor = " << factor << endl; 1124 1124 Vector<Float> factors(in.nRow(), factor); 1125 1125 correctFromVector (pTabOut, in, doAll, factors); … … 1130 1130 } 1131 1131 // 1132 c err<< "Applying supplied conversion factor = " << factor << endl;1132 cout << "Applying supplied conversion factor = " << factor << endl; 1133 1133 Vector<Float> factors(in.nRow(), factor); 1134 1134 correctFromVector (pTabOut, in, doAll, factors); … … 1140 1140 // change per integration. 1141 1141 1142 c err<< "Looking up conversion factors" << endl;1142 cout << "Looking up conversion factors" << endl; 1143 1143 convertBrightnessUnits (pTabOut, in, toKelvin, cFac, doAll); 1144 1144 } … … 1204 1204 } 1205 1205 // 1206 c err<< "Making polynomial correction with " << msg << " coefficients" << endl;1206 cout << "Making polynomial correction with " << msg << " coefficients" << endl; 1207 1207 const uInt nRow = in.nRow(); 1208 1208 Vector<Float> factor(nRow); … … 1222 1222 // Read and correct 1223 1223 1224 c err<< "Making correction from ascii Table" << endl;1224 cout << "Making correction from ascii Table" << endl; 1225 1225 correctFromAsciiTable (pTabOut, in, fileName, col0, col1, 1226 1226 methodStr, doAll, x); … … 1301 1301 // 1302 1302 Vector<Float> JyPerK = sdAtt.JyPerK(inst, dateObs, freqs); 1303 c err<< "Jy/K = " << JyPerK << endl;1303 cout << "Jy/K = " << JyPerK << endl; 1304 1304 Vector<Float> factors = cFac * JyPerK; 1305 1305 if (toKelvin) factors = Float(1.0) / factors; … … 1405 1405 refEpoch = in.getEpoch(0); 1406 1406 } 1407 c err << "Aligning at reference Epoch " << formatEpoch(refEpoch) <<1408 1409 1407 cout << "Aligning at reference Epoch " << formatEpoch(refEpoch) 1408 << " in frame " << MFrequency::showType(freqSystem) << endl; 1409 1410 1410 // Get Reference Position 1411 1411 … … 1452 1452 // New output Table 1453 1453 1454 c err<< "Create output table" << endl;1454 cout << "Create output table" << endl; 1455 1455 SDMemTable* pTabOut = new SDMemTable(in,True); 1456 1456 pTabOut->putSDFreqTable(freqTabOut); … … 1471 1471 for (uInt iRow=0; iRow<nRows; ++iRow) { 1472 1472 if (iRow%10==0) { 1473 c err<< "Processing row " << iRow << endl;1473 cout << "Processing row " << iRow << endl; 1474 1474 } 1475 1475 -
trunk/src/SDMemTable.cc
r401 r414 1381 1381 throw (AipsError("Number of rest frequencies must be equal to the number of IFs")); 1382 1382 } 1383 c err<< "Replacing rest frequencies with given list, one per IF" << endl;1383 cout << "Replacing rest frequencies with given list, one per IF" << endl; 1384 1384 // 1385 1385 sdft.deleteRestFrequencies(); … … 1394 1394 Quantum<Double> rf(restFreqs[0], unit); 1395 1395 idx = sdft.addRestFrequency(rf.getValue("Hz")); 1396 c err<< "Selecting given rest frequency" << endl;1396 cout << "Selecting given rest frequency" << endl; 1397 1397 } 1398 1398 … … 1438 1438 } 1439 1439 1440 void SDMemTable::spectralLines 1440 void SDMemTable::spectralLines() const 1441 1441 { 1442 1442 Vector<String> lines = MeasTable::Lines(); … … 1444 1444 Double freq; 1445 1445 // 1446 c err.flags(std::ios_base::left);1447 c err<< "Line Frequency (Hz)" << endl;1448 c err<< "-----------------------" << endl;1446 cout.flags(std::ios_base::left); 1447 cout << "Line Frequency (Hz)" << endl; 1448 cout << "-----------------------" << endl; 1449 1449 for (uInt i=0; i<lines.nelements(); i++) { 1450 1450 MeasTable::Line(lineFreq, lines[i]); 1451 1451 freq = lineFreq.getValue().getValue(); // Hz 1452 1452 // 1453 c err<< setw(11) << lines[i] << setprecision(10) << freq << endl;1453 cout << setw(11) << lines[i] << setprecision(10) << freq << endl; 1454 1454 } 1455 1455 } -
trunk/src/SDReader.cc
r405 r414 219 219 Vector<Complex> xPol; 220 220 uInt n = seq.size(); 221 221 222 222 223 uInt stepsize = header_->nif*header_->nbeam;
Note:
See TracChangeset
for help on using the changeset viewer.