Changeset 604 in Sophya for trunk/SophyaLib/Samba/spheregorski.cc
- Timestamp:
- Nov 20, 1999, 9:39:49 PM (26 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/Samba/spheregorski.cc
r574 r604 1 #include "machdefs.h" 2 #include <math.h> 3 #include <complex> 4 5 #include "pexceptions.h" 1 6 #include "spheregorski.h" 2 7 #include "strutil.h" 3 #include <math.h>4 #include <complex>5 #include "piocmplx.h"6 8 7 9 extern "C" … … 179 181 { 180 182 InitNul(); 183 pixels_.Reset(); 181 184 } 182 185 … … 196 199 { 197 200 cout << "SphereGorski : m hors bornes [0,8192], m= " << m << endl; 198 exit(1);201 throw RangeCheckError("SphereGorski<T>::SphereGorski() - Out of bound nside (< 8192)!"); 199 202 } 200 203 // verifier que m est une puissance de deux … … 204 207 { 205 208 cout<<"SphereGorski: m doit etre une puissance de deux, m= "<<m<<endl; 206 exit(1);209 throw ParmError("SphereGorski<T>::SphereGorski() - nside != 2^n !"); 207 210 } 208 211 InitNul(); … … 233 236 //-- 234 237 { 235 InitNul(); 238 // InitNul(); - Non , le destructeur ne doit pas mettre les pixels a zero ! 239 // Reza 20/11/99 236 240 } 237 241 … … 303 307 nPix_ = 0; 304 308 omeg_ = 0.; 305 pixels_.Reset(); 306 309 // pixels_.Reset(); - Il ne faut pas mettre les pixels a zero si share ! 307 310 } 308 311 … … 1114 1117 { 1115 1118 dobj= new SphereGorski<T>; 1119 dobj->DataBlock().SetTemp(true); 1116 1120 ownobj= true; 1117 1121 Read(filename); … … 1122 1126 { 1123 1127 dobj= new SphereGorski<T>(obj, true); 1128 dobj->DataBlock().SetTemp(true); 1124 1129 ownobj= true; 1125 1130 } … … 1151 1156 { 1152 1157 dobj= new SphereGorski<T>; 1158 dobj->DataBlock().SetTemp(true); 1159 ownobj= true; 1153 1160 } 1154 1161 … … 1179 1186 dobj->setPixSolAngle(Omega); 1180 1187 1181 T* pixels= new T[nPix]; 1182 PIOSReadArray(is, pixels, nPix); 1183 dobj->setDataBlock(pixels, nPix); 1184 delete [] pixels; 1185 1188 // On lit le DataBlock; 1189 is >> dobj->DataBlock(); 1186 1190 } 1187 1191 … … 1218 1222 os.PutI4(nPix); 1219 1223 os.PutR8(dobj->PixSolAngle(0)); 1220 1221 PIOSWriteArray(os,(dobj->getDataBlock())->Data(), nPix); 1222 1223 1224 // On ecrit le dataBlock 1225 os << dobj->DataBlock(); 1224 1226 } 1225 1227
Note:
See TracChangeset
for help on using the changeset viewer.