Changeset 934 for trunk/src/STMath.cpp
- Timestamp:
- 03/27/06 19:32:14 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/STMath.cpp
r926 r934 70 70 WeightType wtype = stringToWeight(weight); 71 71 72 std::vector<CountedPtr<Scantable> >* in = const_cast<std::vector<CountedPtr<Scantable> >*>(&scantbls); 72 std::vector<CountedPtr<Scantable> >* in = 73 const_cast<std::vector<CountedPtr<Scantable> >*>(&scantbls); 73 74 std::vector<CountedPtr<Scantable> > inaligned; 74 75 if ( alignfreq ) { … … 940 941 ArrayColumn<Float> speccol(t, "SPECTRA"); 941 942 Matrix<Float> pols = speccol.getColumn(); 942 Matrix<Float> out;943 943 try { 944 944 stpol->setSpectra(pols); 945 945 (stpol->*fptr)(phase); 946 946 speccol.putColumn(stpol->getSpectra()); 947 delete stpol;stpol=0;948 947 } catch (AipsError& e) { 949 948 delete stpol;stpol=0; … … 952 951 ++iter; 953 952 } 953 delete stpol;stpol=0; 954 954 return out; 955 955 } … … 1016 1016 } 1017 1017 MFrequency::Types system = in->frequencies().getFrame(); 1018 out->frequencies().setFrame(system );1018 out->frequencies().setFrame(system, true); 1019 1019 // set up the iterator 1020 1020 Block<String> cols(4); … … 1051 1051 SpectralCoordinate sc2 = fa.alignedSpectralCoordinate(linear); 1052 1052 sc2.setWorldAxisUnits(units); 1053 out->frequencies().addEntry(sc2.referencePixel()[0], 1054 sc2.referenceValue()[0], 1055 sc2.increment()[0]); 1053 uInt id = out->frequencies().addEntry(sc2.referencePixel()[0], 1054 sc2.referenceValue()[0], 1055 sc2.increment()[0]); 1056 TableVector<uInt> tvec(ftab, "FREQ_ID"); 1057 tvec = id; 1056 1058 // create the "global" abcissa for alignment with same FREQ_ID 1057 1059 Vector<Double> abc(nchan);
Note: See TracChangeset
for help on using the changeset viewer.