Changeset 3572 in Sophya for trunk/SophyaExt
- Timestamp:
- Feb 7, 2009, 10:50:34 PM (17 years ago)
- Location:
- trunk/SophyaExt
- Files:
-
- 4 deleted
- 28 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaExt/FitsIOServer/Makefile
r3407 r3572 4 4 5 5 clean: 6 rm -f $(SOPHYAOBJP)datacirclefits.o $(SOPHYAOBJP)fabtcolread.o $(SOPHYAOBJP)fabtwriter.o $(SOPHYAOBJP)fboloread.o $(SOPHYAOBJP)fbtntintf.o $(SOPHYAOBJP)fiosinit.o $(SOPHYAOBJP)fits autoreader.o $(SOPHYAOBJP)fitsbntbllineRW.o $(SOPHYAOBJP)fitsfile.o $(SOPHYAOBJP)fitsgenedata.o $(SOPHYAOBJP)fitshdtable.o $(SOPHYAOBJP)fitshisterr.o $(SOPHYAOBJP)fitshistos.o $(SOPHYAOBJP)fitsinoutfile.o $(SOPHYAOBJP)fitslocalmap.o $(SOPHYAOBJP)fitsmanager.o $(SOPHYAOBJP)fitsntuple.o $(SOPHYAOBJP)fitsspherehealpix.o $(SOPHYAOBJP)fitsspherethetaphi.o $(SOPHYAOBJP)fitstarray.o $(SOPHYAOBJP)fitsxntuple.o $(SOPHYAOBJP)swfitsdtable.o6 rm -f $(SOPHYAOBJP)datacirclefits.o $(SOPHYAOBJP)fabtcolread.o $(SOPHYAOBJP)fabtwriter.o $(SOPHYAOBJP)fboloread.o $(SOPHYAOBJP)fbtntintf.o $(SOPHYAOBJP)fiosinit.o $(SOPHYAOBJP)fitsbntbllineRW.o $(SOPHYAOBJP)fitsfile.o $(SOPHYAOBJP)fitsgenedata.o $(SOPHYAOBJP)fitshdtable.o $(SOPHYAOBJP)fitshisterr.o $(SOPHYAOBJP)fitshistos.o $(SOPHYAOBJP)fitsinoutfile.o $(SOPHYAOBJP)fitslocalmap.o $(SOPHYAOBJP)fitsmanager.o $(SOPHYAOBJP)fitsntuple.o $(SOPHYAOBJP)fitsspherehealpix.o $(SOPHYAOBJP)fitsspherethetaphi.o $(SOPHYAOBJP)fitstarray.o $(SOPHYAOBJP)swfitsdtable.o 7 7 rm -f $(SOPHYALIBP)libFitsIOServer.a 8 8 9 $(SOPHYALIBP)libFitsIOServer.a : $(SOPHYAOBJP)datacirclefits.o $(SOPHYAOBJP)fabtcolread.o $(SOPHYAOBJP)fabtwriter.o $(SOPHYAOBJP)fboloread.o $(SOPHYAOBJP)fbtntintf.o $(SOPHYAOBJP)fiosinit.o $(SOPHYAOBJP)fits autoreader.o $(SOPHYAOBJP)fitsbntbllineRW.o $(SOPHYAOBJP)fitsfile.o $(SOPHYAOBJP)fitsgenedata.o $(SOPHYAOBJP)fitshdtable.o $(SOPHYAOBJP)fitshisterr.o $(SOPHYAOBJP)fitshistos.o $(SOPHYAOBJP)fitsinoutfile.o $(SOPHYAOBJP)fitslocalmap.o $(SOPHYAOBJP)fitsmanager.o $(SOPHYAOBJP)fitsntuple.o $(SOPHYAOBJP)fitsspherehealpix.o $(SOPHYAOBJP)fitsspherethetaphi.o $(SOPHYAOBJP)fitstarray.o $(SOPHYAOBJP)fitsxntuple.o $(SOPHYAOBJP)swfitsdtable.o9 $(SOPHYALIBP)libFitsIOServer.a : $(SOPHYAOBJP)datacirclefits.o $(SOPHYAOBJP)fabtcolread.o $(SOPHYAOBJP)fabtwriter.o $(SOPHYAOBJP)fboloread.o $(SOPHYAOBJP)fbtntintf.o $(SOPHYAOBJP)fiosinit.o $(SOPHYAOBJP)fitsbntbllineRW.o $(SOPHYAOBJP)fitsfile.o $(SOPHYAOBJP)fitsgenedata.o $(SOPHYAOBJP)fitshdtable.o $(SOPHYAOBJP)fitshisterr.o $(SOPHYAOBJP)fitshistos.o $(SOPHYAOBJP)fitsinoutfile.o $(SOPHYAOBJP)fitslocalmap.o $(SOPHYAOBJP)fitsmanager.o $(SOPHYAOBJP)fitsntuple.o $(SOPHYAOBJP)fitsspherehealpix.o $(SOPHYAOBJP)fitsspherethetaphi.o $(SOPHYAOBJP)fitstarray.o $(SOPHYAOBJP)swfitsdtable.o 10 10 $(AR) $(ARFLAGS) $@ $? 11 11 touch $(SOPHYAINCP)/SophyaConfInfo/libextsophya.objlist … … 184 184 $(CXXCOMPILE) $(CXXTEMPFLG) -o $@ fiosinit.cc 185 185 186 $(SOPHYAOBJP)fitsautoreader.o: fitsautoreader.cc \187 $(SOPHYAINCP)sopnamsp.h \188 $(SOPHYAINCP)pexceptions.h \189 $(SOPHYAINCP)machdefs.h fitsautoreader.h \190 $(SOPHYAINCP)machdefs.h \191 $(SOPHYAINCP)anydataobj.h fitsfile.h \192 $(SOPHYAINCP)ndatablock.h \193 $(SOPHYAINCP)anydataobj.h \194 $(SOPHYAINCP)dvlist.h \195 $(SOPHYAINCP)objfio.h \196 $(SOPHYAINCP)ppersist.h \197 $(SOPHYAINCP)pexceptions.h \198 $(SOPHYAINCP)gnumd5.h \199 $(SOPHYAINCP)ppfbinstream.h \200 $(SOPHYAINCP)rawstream.h \201 $(SOPHYAINCP)ppfnametag.h \202 $(SOPHYAINCP)mutyv.h \203 $(SOPHYAINCP)timestamp.h fitsinoutfile.h \204 $(SOPHYAINCP)FitsIO/fitsio.h \205 $(SOPHYAINCP)FitsIO/longnam.h fitshandler.h \206 fitstarray.h $(SOPHYAINCP)tarray.h \207 $(SOPHYAINCP)basarr.h \208 $(SOPHYAINCP)dvlist.h \209 $(SOPHYAINCP)ndatablock.h \210 $(SOPHYAINCP)utilarr.h fitsspherehealpix.h \211 $(SOPHYAINCP)spherehealpix.h \212 $(SOPHYAINCP)sphericalmap.h \213 $(SOPHYAINCP)smathconst.h \214 $(SOPHYAINCP)pixelmap.h \215 $(SOPHYAINCP)spherepos.h \216 $(SOPHYAINCP)unitvector.h \217 $(SOPHYAINCP)vector3d.h \218 $(SOPHYAINCP)longlat.h \219 $(SOPHYAINCP)utilgeom.h \220 $(SOPHYAINCP)datatype.h \221 $(SOPHYAINCP)tvector.h \222 $(SOPHYAINCP)tmatrix.h \223 $(SOPHYAINCP)tarray.h \224 $(SOPHYAINCP)tmatrix_tsnl.h \225 $(SOPHYAINCP)tvector_tsnl.h \226 $(SOPHYAINCP)HEALPixUtils.h fitsntuple.h \227 $(SOPHYAINCP)ntuple.h \228 $(SOPHYAINCP)ntupintf.h \229 $(SOPHYAINCP)ppersist.h fitsxntuple.h \230 $(SOPHYAINCP)xntuple.h fitslocalmap.h \231 $(SOPHYAINCP)localmap.h fbtntintf.h \232 $(SOPHYAINCP)ntupintf.h233 $(CXXCOMPILE) $(CXXTEMPFLG) -o $@ fitsautoreader.cc234 235 186 $(SOPHYAOBJP)fitsbntbllineRW.o: fitsbntbllineRW.cc \ 236 187 $(SOPHYAINCP)sopnamsp.h \ … … 633 584 $(CXXCOMPILE) $(CXXTEMPFLG) -o $@ fitstarray.cc 634 585 635 $(SOPHYAOBJP)fitsxntuple.o: fitsxntuple.cc \636 $(SOPHYAINCP)sopnamsp.h \637 $(SOPHYAINCP)pexceptions.h \638 $(SOPHYAINCP)machdefs.h fitsxntuple.h \639 $(SOPHYAINCP)machdefs.h \640 $(SOPHYAINCP)xntuple.h \641 $(SOPHYAINCP)ntupintf.h \642 $(SOPHYAINCP)ppersist.h \643 $(SOPHYAINCP)pexceptions.h \644 $(SOPHYAINCP)gnumd5.h \645 $(SOPHYAINCP)ppfbinstream.h \646 $(SOPHYAINCP)rawstream.h \647 $(SOPHYAINCP)ppfnametag.h \648 $(SOPHYAINCP)dvlist.h \649 $(SOPHYAINCP)objfio.h \650 $(SOPHYAINCP)anydataobj.h \651 $(SOPHYAINCP)mutyv.h \652 $(SOPHYAINCP)timestamp.h \653 $(SOPHYAINCP)anydataobj.h \654 $(SOPHYAINCP)ppersist.h fitsfile.h \655 $(SOPHYAINCP)ndatablock.h \656 $(SOPHYAINCP)dvlist.h fitsinoutfile.h \657 $(SOPHYAINCP)FitsIO/fitsio.h \658 $(SOPHYAINCP)FitsIO/longnam.h fitshandler.h659 $(CXXCOMPILE) $(CXXTEMPFLG) -o $@ fitsxntuple.cc660 661 586 $(SOPHYAOBJP)swfitsdtable.o: swfitsdtable.cc swfitsdtable.h \ 662 587 $(SOPHYAINCP)basedtable.h \ -
trunk/SophyaExt/FitsIOServer/fabtcolread.cc
r3188 r3572 189 189 \return value into double 190 190 */ 191 double FitsOpenFile::ReadKey(fitsfile *fitsptr,c har *keyname)191 double FitsOpenFile::ReadKey(fitsfile *fitsptr,const char *keyname) 192 192 { 193 193 if(keyname==NULL || fitsptr==NULL) return 0.; … … 204 204 \return value into long 205 205 */ 206 long FitsOpenFile::ReadKeyL(fitsfile *fitsptr,c har *keyname)206 long FitsOpenFile::ReadKeyL(fitsfile *fitsptr,const char *keyname) 207 207 { 208 208 if(keyname==NULL || fitsptr==NULL) return 0; … … 219 219 \return value into long long 220 220 */ 221 LONGLONG FitsOpenFile::ReadKeyLL(fitsfile *fitsptr,c har *keyname)221 LONGLONG FitsOpenFile::ReadKeyLL(fitsfile *fitsptr,const char *keyname) 222 222 { 223 223 if(keyname==NULL || fitsptr==NULL) return 0; … … 234 234 \return value into string 235 235 */ 236 string FitsOpenFile::ReadKeyS(fitsfile *fitsptr,c har *keyname)236 string FitsOpenFile::ReadKeyS(fitsfile *fitsptr,const char *keyname) 237 237 { 238 238 if(keyname==NULL || fitsptr==NULL) return (string)""; … … 529 529 \return value into double 530 530 */ 531 double FitsABTColRd::ReadKey(c har *keyname)531 double FitsABTColRd::ReadKey(const char *keyname) 532 532 { 533 533 return FitsOpenFile::ReadKey(GetFitsPtr(),keyname); … … 539 539 \return value into long 540 540 */ 541 long FitsABTColRd::ReadKeyL(c har *keyname)541 long FitsABTColRd::ReadKeyL(const char *keyname) 542 542 { 543 543 return FitsOpenFile::ReadKeyL(GetFitsPtr(),keyname); … … 549 549 \return value into long long 550 550 */ 551 LONGLONG FitsABTColRd::ReadKeyLL(c har *keyname)551 LONGLONG FitsABTColRd::ReadKeyLL(const char *keyname) 552 552 { 553 553 return FitsOpenFile::ReadKeyLL(GetFitsPtr(),keyname); … … 559 559 \return value into string 560 560 */ 561 string FitsABTColRd::ReadKeyS(c har *keyname)561 string FitsABTColRd::ReadKeyS(const char *keyname) 562 562 { 563 563 return FitsOpenFile::ReadKeyS(GetFitsPtr(),keyname); … … 1080 1080 \return value into double 1081 1081 */ 1082 double FitsImg2DRd::ReadKey(c har *keyname)1082 double FitsImg2DRd::ReadKey(const char *keyname) 1083 1083 { 1084 1084 return FitsOpenFile::ReadKey(GetFitsPtr(),keyname); … … 1090 1090 \return value into long 1091 1091 */ 1092 long FitsImg2DRd::ReadKeyL(c har *keyname)1092 long FitsImg2DRd::ReadKeyL(const char *keyname) 1093 1093 { 1094 1094 return FitsOpenFile::ReadKeyL(GetFitsPtr(),keyname); … … 1100 1100 \return value into long long 1101 1101 */ 1102 LONGLONG FitsImg2DRd::ReadKeyLL(c har *keyname)1102 LONGLONG FitsImg2DRd::ReadKeyLL(const char *keyname) 1103 1103 { 1104 1104 return FitsOpenFile::ReadKeyLL(GetFitsPtr(),keyname); … … 1110 1110 \return value into string 1111 1111 */ 1112 string FitsImg2DRd::ReadKeyS(c har *keyname)1112 string FitsImg2DRd::ReadKeyS(const char *keyname) 1113 1113 { 1114 1114 return FitsOpenFile::ReadKeyS(GetFitsPtr(),keyname); … … 1428 1428 \return value into double 1429 1429 */ 1430 double FitsImg3DRd::ReadKey(c har *keyname)1430 double FitsImg3DRd::ReadKey(const char *keyname) 1431 1431 { 1432 1432 return FitsOpenFile::ReadKey(GetFitsPtr(),keyname); … … 1438 1438 \return value into long 1439 1439 */ 1440 long FitsImg3DRd::ReadKeyL(c har *keyname)1440 long FitsImg3DRd::ReadKeyL(const char *keyname) 1441 1441 { 1442 1442 return FitsOpenFile::ReadKeyL(GetFitsPtr(),keyname); … … 1448 1448 \return value into long long 1449 1449 */ 1450 LONGLONG FitsImg3DRd::ReadKeyLL(c har *keyname)1450 LONGLONG FitsImg3DRd::ReadKeyLL(const char *keyname) 1451 1451 { 1452 1452 return FitsOpenFile::ReadKeyLL(GetFitsPtr(),keyname); … … 1458 1458 \return value into string 1459 1459 */ 1460 string FitsImg3DRd::ReadKeyS(c har *keyname)1460 string FitsImg3DRd::ReadKeyS(const char *keyname) 1461 1461 { 1462 1462 return FitsOpenFile::ReadKeyS(GetFitsPtr(),keyname); -
trunk/SophyaExt/FitsIOServer/fabtcolread.h
r3493 r3572 48 48 void Print(void); 49 49 50 static double ReadKey(fitsfile *fitsptr,c har *keyname);51 static long ReadKeyL(fitsfile *fitsptr,c har *keyname);52 static LONGLONG ReadKeyLL(fitsfile *fitsptr,c har *keyname);53 static string ReadKeyS(fitsfile *fitsptr,c har *keyname);50 static double ReadKey(fitsfile *fitsptr,const char *keyname); 51 static long ReadKeyL(fitsfile *fitsptr,const char *keyname); 52 static LONGLONG ReadKeyLL(fitsfile *fitsptr,const char *keyname); 53 static string ReadKeyS(fitsfile *fitsptr,const char *keyname); 54 54 static void printerror(int sta); 55 55 … … 78 78 void ChangeBuffer(long blen=100,long bsens=1); 79 79 80 double ReadKey(c har *keyname);81 long ReadKeyL(c har *keyname);82 LONGLONG ReadKeyLL(c har *keyname);83 string ReadKeyS(c har *keyname);80 double ReadKey(const char *keyname); 81 long ReadKeyL(const char *keyname); 82 LONGLONG ReadKeyLL(const char *keyname); 83 string ReadKeyS(const char *keyname); 84 84 85 85 double Read(LONGLONG n,bool usebuffer=true); … … 195 195 virtual ~FitsImg2DRd(); 196 196 197 double ReadKey(c har *keyname);198 long ReadKeyL(c har *keyname);199 LONGLONG ReadKeyLL(c har *keyname);200 string ReadKeyS(c har *keyname);197 double ReadKey(const char *keyname); 198 long ReadKeyL(const char *keyname); 199 LONGLONG ReadKeyLL(const char *keyname); 200 string ReadKeyS(const char *keyname); 201 201 202 202 LONGLONG Read(TMatrix<uint_2>& data); … … 260 260 virtual ~FitsImg3DRd(); 261 261 262 double ReadKey(c har *keyname);263 long ReadKeyL(c har *keyname);264 LONGLONG ReadKeyLL(c har *keyname);265 string ReadKeyS(c har *keyname);262 double ReadKey(const char *keyname); 263 long ReadKeyL(const char *keyname); 264 LONGLONG ReadKeyLL(const char *keyname); 265 string ReadKeyS(const char *keyname); 266 266 267 267 LONGLONG Read(TArray<uint_2>& data); -
trunk/SophyaExt/FitsIOServer/fabtwriter.cc
r3128 r3572 105 105 106 106 /*! Write a double value into Fits Header */ 107 void FitsWriter::WriteKey(const char *keyname,double val,c har* comment)107 void FitsWriter::WriteKey(const char *keyname,double val,const char* comment) 108 108 { 109 109 if(keyname==NULL || strlen(keyname)<=0) return; … … 116 116 117 117 /*! Write a long value into Fits Header */ 118 void FitsWriter::WriteKey(const char *keyname,long val,c har* comment)118 void FitsWriter::WriteKey(const char *keyname,long val,const char* comment) 119 119 { 120 120 if(keyname==NULL || strlen(keyname)<=0) return; … … 127 127 128 128 /*! Write a long long value into Fits Header */ 129 void FitsWriter::WriteKey(const char *keyname,LONGLONG val,c har* comment)129 void FitsWriter::WriteKey(const char *keyname,LONGLONG val,const char* comment) 130 130 { 131 131 if(keyname==NULL || strlen(keyname)<=0) return; … … 138 138 139 139 /*! Write a string value into Fits Header */ 140 void FitsWriter::WriteKey(const char *keyname,string val,c har* comment)140 void FitsWriter::WriteKey(const char *keyname,string val,const char* comment) 141 141 { 142 142 if(keyname==NULL || strlen(keyname)<=0) return; -
trunk/SophyaExt/FitsIOServer/fabtwriter.h
r3493 r3572 36 36 37 37 //! Write a double value in Fits header. 38 void WriteKey(const char *keyname,double val,c har* comment=NULL);38 void WriteKey(const char *keyname,double val,const char* comment=NULL); 39 39 //! Write a long value in Fits header. 40 void WriteKey(const char *keyname,long val,c har* comment=NULL);40 void WriteKey(const char *keyname,long val,const char* comment=NULL); 41 41 //! Write a long long value in Fits header. 42 void WriteKey(const char *keyname,LONGLONG val,c har* comment=NULL);42 void WriteKey(const char *keyname,LONGLONG val,const char* comment=NULL); 43 43 //! Write a string value in Fits header. 44 void WriteKey(const char *keyname,string val,c har* comment=NULL);44 void WriteKey(const char *keyname,string val,const char* comment=NULL); 45 45 //! Write a character string value in Fits header. 46 inline void WriteKey(const char *keyname,c har* val,char* comment=NULL)46 inline void WriteKey(const char *keyname,const char* val,const char* comment=NULL) 47 47 {string dum=val; WriteKey(keyname,dum,comment);} 48 48 //! Set debug level … … 93 93 inline void SetExtName(string extname=string("")) {ExtName = extname;} 94 94 //! Set the FITS table extension name 95 inline void SetExtName(c har* extname="") {ExtName = extname;}95 inline void SetExtName(const char* extname="") {ExtName = extname;} 96 96 97 97 //! Add a new column to the FITS table and return its number (see addcol). -
trunk/SophyaExt/FitsIOServer/fitsblkrw.h
r3493 r3572 159 159 string msg = "FitsBlockRW<std::string>::WriteColumnData() Error: " ; 160 160 msg += buff; 161 sprintf(buff," kk=%ld", kk); msg += buff;161 sprintf(buff," kk=%ld",(long)kk); msg += buff; 162 162 throw FitsIOException(msg); 163 163 } … … 187 187 string msg = "FitsBlockRW<std::string>::ReadColumnData() Error: " ; 188 188 msg += buff; 189 sprintf(buff," kk=%ld", kk); msg += buff;189 sprintf(buff," kk=%ld",(long)kk); msg += buff; 190 190 throw FitsIOException(msg); 191 191 } -
trunk/SophyaExt/FitsIOServer/fitsbntbllineRW.cc
r2615 r3572 64 64 else if (ss == FitsFile::FitsDataType_char) SfitsCol.push_back(k); 65 65 else { 66 cout << " FITS_XNTuple: colonne fits " << k << " type= " << (int) ss << endl;66 cout << "FITS_BntblLineReader: colonne fits " << k << " type= " << (int) ss << endl; 67 67 throw IOExc("type de champ inconnu"); 68 68 } … … 113 113 114 114 115 FITS_BntblLineWriter::FITS_BntblLineWriter(c har inputfile[],int dc, int fc,int ic, int lc, int bc, int cc, vector<string> names,DVList* ptr_dvl, FitsFile::WriteMode wrm)115 FITS_BntblLineWriter::FITS_BntblLineWriter(const char inputfile[],int dc, int fc,int ic, int lc, int bc, int cc, vector<string> names,DVList* ptr_dvl, FitsFile::WriteMode wrm) 116 116 { 117 117 int k; 118 118 int nbcols = dc+fc+ic+cc+lc+bc; 119 if (nbcols != names.size())119 if (nbcols != (int)names.size()) 120 120 { 121 121 cout << " WARNING: BnTblLineWriter:: length of vector of column names not equal to total number of columns" << endl; -
trunk/SophyaExt/FitsIOServer/fitsbntbllineRW.h
r2201 r3572 6 6 #include "dvlist.h" 7 7 #include "fitsfile.h" 8 //#include "xntuple.h"9 8 10 9 namespace SOPHYA { … … 66 65 67 66 public: 68 FITS_BntblLineWriter(c har inputfile[],int dc, int fc, int ic, int lc, int bc,int cc, vector<string> names, DVList* dvl=NULL, FitsFile::WriteMode wrm = FitsFile::clear);67 FITS_BntblLineWriter(const char inputfile[],int dc, int fc, int ic, int lc, int bc,int cc, vector<string> names, DVList* dvl=NULL, FitsFile::WriteMode wrm = FitsFile::clear); 69 68 ~FITS_BntblLineWriter(); 70 69 void WriteNextLine( BnTblLine& WorkLine); -
trunk/SophyaExt/FitsIOServer/fitsfile.cc
r2907 r3572 228 228 */ 229 229 230 /*! \fn void SOPHYA::FitsIOHandler::Read(c har flnm[],int hdunum)230 /*! \fn void SOPHYA::FitsIOHandler::Read(const char flnm[],int hdunum) 231 231 232 232 this method is called from inherited objects : … … 238 238 calls the method 'ReadFromFits' from the inherited object 239 239 */ 240 void FitsIOHandler::Read(c har flnm[],int hdunum)240 void FitsIOHandler::Read(const char flnm[],int hdunum) 241 241 { 242 242 FitsInFile ifts(flnm); … … 274 274 275 275 */ 276 void FitsIOHandler::Write(c har flnm[])276 void FitsIOHandler::Write(const char flnm[]) 277 277 278 278 { … … 327 327 } 328 328 329 void FitsFile::printerror(int& status, c har* texte)329 void FitsFile::printerror(int& status, const char* texte) 330 330 //*****************************************************/ 331 331 //* Print out cfitsio error messages and exit program */ … … 338 338 throw IOExc("FitsFile:: error FITSIO status"); 339 339 } 340 void FitsFile::printerrorAndContinue(int& status, c har* texte)340 void FitsFile::printerrorAndContinue(int& status, const char* texte) 341 341 //*****************************************************/ 342 342 //* Print out cfitsio error messages and exit program */ … … 1534 1534 char dtype; 1535 1535 char card[FLEN_CARD]; 1536 c har *comkey = "COMMENT";1536 const char *comkey = "COMMENT"; 1537 1537 char comment[FLEN_COMMENT]; 1538 1538 … … 2198 2198 char comment[FLEN_COMMENT]; 2199 2199 char strval[FLEN_VALUE]= ""; 2200 c har *comkey = "COMMENT";2200 const char *comkey = "COMMENT"; 2201 2201 // fits_read_keyword(fptr_, keyname, strval, NULL, &status); 2202 2202 // if (status != 0 || strncmp(keyname,comkey,LEN_KEYWORD-1) == 0 ) … … 2527 2527 default : 2528 2528 { 2529 c har *comkey = "COMMENT";2529 const char *comkey = "COMMENT"; 2530 2530 if(strncmp(keyname,comkey,LEN_KEYWORD-1) == 0) 2531 2531 { -
trunk/SophyaExt/FitsIOServer/fitsfile.h
r3493 r3572 64 64 65 65 virtual ~FitsIOHandler() {} 66 void Read(c har flnm[],int hdunum= 0);67 void Write(c har flnm[]) ;66 void Read(const char flnm[],int hdunum= 0); 67 void Write(const char flnm[]) ; 68 68 69 69 //Implementation par defaut de l'interface FitsHandlerInterface … … 181 181 void ResetStatus(int& status) ; 182 182 static void printerror(int&) ; 183 static void printerror(int&,c har* texte) ;184 static void printerrorAndContinue(int& status, char* texte);183 static void printerror(int&,const char* texte) ; 184 static void printerrorAndContinue(int& status,const char* texte); 185 185 inline void InitNull() 186 186 { -
trunk/SophyaExt/FitsIOServer/fitsinoutfile.cc
r3452 r3572 157 157 } 158 158 /* -- Fonction utilitaire pour verifier le code d'erreur fitsio -- */ 159 static inline bool FitsCheckStatus(int st, c har * emsg = NULL) {159 static inline bool FitsCheckStatus(int st, const char * emsg = NULL) { 160 160 if (st) { 161 161 fits_report_error(stderr, st); … … 272 272 int FitsInOutFile::CurrentHDU() const 273 273 { 274 int status = 0;274 //unused: int status = 0; 275 275 int curhdu = 0; 276 276 fits_get_hdu_num(FitsPtr() , &curhdu); … … 496 496 int colnum, typecode; 497 497 LONGLONG repeat, colw; // $CHECK$ LONGLONG ??? 498 int ncols = 0;498 //unused: int ncols = 0; 499 499 char colname[128]; // longueur max d'un nom de colonne 500 500 501 501 while (status != COL_NOT_FOUND) { 502 fits_get_colname(FitsPtr(), CASEINSEN, "*", colname, &colnum, &status); 502 char dum[2] = {'*','\0'}; 503 fits_get_colname(FitsPtr(), CASEINSEN, dum, colname, &colnum, &status); 503 504 if (status == COL_NOT_FOUND) break; 504 505 if ( (status != COL_NOT_UNIQUE) && (status != 0) ) { … … 665 666 666 667 CheckFitsPtr(FitsPtr()); 667 int status = 0;668 //unused: int status = 0; 668 669 DVList::ValList::const_iterator it; 669 670 for(it = dvl.Begin(); it != dvl.End(); it++) -
trunk/SophyaExt/FitsIOServer/fitsioserver.h
r3493 r3572 19 19 #include "fitsspherehealpix.h" 20 20 #include "fitsntuple.h" 21 // #include "fitsxntuple.h"22 21 #include "fitslocalmap.h" 23 22 -
trunk/SophyaExt/FitsIOServer/fitsmanager.cc
r3167 r3572 80 80 int hfg = 0; 81 81 int bhfg = 0; 82 int clev = 0;82 //unused: int clev = 0; 83 83 int blev = 0; 84 84 for(it = hlistp->begin(); it != hlistp->end(); it++) { … … 196 196 LONGLONG naxes[5] = {0,0,0,0,0}; 197 197 int naxis=5; 198 int imgtyp = is.GetImageHDUInfo(naxis, naxes);198 //unused: int imgtyp = is.GetImageHDUInfo(naxis, naxes); 199 199 os << ">> IMAGE_HDU: naxis= " << naxis << " : "; 200 200 for(int i=0; i<naxis; i++) { … … 212 212 else os << ">> ASCII_TBL : NRows= " << is.GetNbRows(); 213 213 os << " x NCols= " << ncols << endl; 214 for( int kk=0; kk<colnames.size(); kk++) {214 for(size_t kk=0; kk<colnames.size(); kk++) { 215 215 os << "Col[" << kk+1 << "] Name= " << colnames[kk] 216 216 << " Type= " << FitsTypes::DataTypeToTypeString(coltypes[kk]) -
trunk/SophyaExt/FitsIOServer/fitsntuple.h
r3493 r3572 53 53 54 54 /*! 55 fill the XNTuple only with lines from firstLine-th to (firstLine+numberOfLines-1)-th of the FITS-file inputfile.55 fill the NTuple only with lines from firstLine-th to (firstLine+numberOfLines-1)-th of the FITS-file inputfile. 56 56 \param <firstLine> first line to be read (the first line of the file is numbered 0) 57 57 \param <numberOfLines> number of lines to be read -
trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc
r3167 r3572 27 27 28 28 template <class T> 29 FITS_SphereHEALPix<T>::FITS_SphereHEALPix(c har inputfile[],int hdunum)29 FITS_SphereHEALPix<T>::FITS_SphereHEALPix(const char inputfile[],int hdunum) 30 30 { 31 31 dobj_= new SphereHEALPix<T>; … … 266 266 267 267 template <class T> 268 void FITS_SphereHEALPix<T>::Mollweide_picture_projection(c har filename[])268 void FITS_SphereHEALPix<T>::Mollweide_picture_projection(const char filename[]) 269 269 { 270 270 int ni = 300; … … 297 297 } 298 298 template <class T> 299 void FITS_SphereHEALPix<T>::sinus_picture_projection(c har filename[])299 void FITS_SphereHEALPix<T>::sinus_picture_projection(const char filename[]) 300 300 { 301 301 int ni = 300; -
trunk/SophyaExt/FitsIOServer/fitsspherehealpix.h
r3493 r3572 24 24 25 25 FITS_SphereHEALPix(); 26 FITS_SphereHEALPix(c har inputfile[],int hdunum=0);26 FITS_SphereHEALPix(const char inputfile[],int hdunum=0); 27 27 FITS_SphereHEALPix(const SphereHEALPix<T>& obj); 28 28 FITS_SphereHEALPix(SphereHEALPix<T>* obj); … … 35 35 virtual FitsHandlerInterface* Clone(); 36 36 37 void Mollweide_picture_projection(c har flnm[]);38 void sinus_picture_projection(c har flnm[]);37 void Mollweide_picture_projection(const char flnm[]); 38 void sinus_picture_projection(const char flnm[]); 39 39 //void ReWriteFromFitsToFits(FitsInFile& headerin, FitsOutFile& os); 40 40 -
trunk/SophyaExt/FitsIOServer/fitsspherethetaphi.cc
r3167 r3572 26 26 27 27 template <class T> 28 FITS_SphereThetaPhi<T>::FITS_SphereThetaPhi(c har inputfile[],int hdunum)28 FITS_SphereThetaPhi<T>::FITS_SphereThetaPhi(const char inputfile[],int hdunum) 29 29 { 30 30 dobj_= new SphereThetaPhi<T>; … … 255 255 256 256 template <class T> 257 void FITS_SphereThetaPhi<T>::Mollweide_picture_projection(c har filename[])257 void FITS_SphereThetaPhi<T>::Mollweide_picture_projection(const char filename[]) 258 258 { 259 259 int ni = 300; … … 286 286 } 287 287 template <class T> 288 void FITS_SphereThetaPhi<T>::sinus_picture_projection(c har filename[])288 void FITS_SphereThetaPhi<T>::sinus_picture_projection(const char filename[]) 289 289 { 290 290 int ni = 300; -
trunk/SophyaExt/FitsIOServer/fitsspherethetaphi.h
r3493 r3572 22 22 23 23 FITS_SphereThetaPhi(); 24 FITS_SphereThetaPhi(c har inputfile[],int hdunum=0);24 FITS_SphereThetaPhi(const char inputfile[],int hdunum=0); 25 25 FITS_SphereThetaPhi(const SphereThetaPhi<T>& obj); 26 26 FITS_SphereThetaPhi(SphereThetaPhi<T>* obj); … … 33 33 virtual FitsHandlerInterface* Clone(); 34 34 35 void Mollweide_picture_projection(c har flnm[]);36 void sinus_picture_projection(c har flnm[]);35 void Mollweide_picture_projection(const char flnm[]); 36 void sinus_picture_projection(const char flnm[]); 37 37 38 38 -
trunk/SophyaExt/FitsIOServer/fitstarray.cc
r2874 r3572 24 24 25 25 template <class T> 26 FITS_TArray<T>::FITS_TArray(c har inputfile[],int hdunum)26 FITS_TArray<T>::FITS_TArray(const char inputfile[],int hdunum) 27 27 { 28 28 dobj_=NULL; -
trunk/SophyaExt/FitsIOServer/fitstarray.h
r1231 r3572 25 25 26 26 FITS_TArray(); 27 FITS_TArray(c har inputfile[],int hdunum=0);27 FITS_TArray(const char inputfile[],int hdunum=0); 28 28 FITS_TArray(const TArray<T> & obj); 29 29 FITS_TArray(TArray<T> *obj); -
trunk/SophyaExt/FitsIOServer/objlist.list
r3123 r3572 5 5 fbtntintf.o 6 6 fiosinit.o 7 fitsautoreader.o8 7 fitsbntbllineRW.o 9 8 fitsfile.o … … 19 18 fitsspherethetaphi.o 20 19 fitstarray.o 21 fitsxntuple.o22 20 swfitsdtable.o -
trunk/SophyaExt/FitsIOServer/smakefile
r3407 r3572 4 4 5 5 clean: 6 rm -f $(SOPHYAOBJP)datacirclefits.o $(SOPHYAOBJP)fabtcolread.o $(SOPHYAOBJP)fabtwriter.o $(SOPHYAOBJP)fboloread.o $(SOPHYAOBJP)fbtntintf.o $(SOPHYAOBJP)fiosinit.o $(SOPHYAOBJP)fits autoreader.o $(SOPHYAOBJP)fitsbntbllineRW.o $(SOPHYAOBJP)fitsfile.o $(SOPHYAOBJP)fitsgenedata.o $(SOPHYAOBJP)fitshdtable.o $(SOPHYAOBJP)fitshisterr.o $(SOPHYAOBJP)fitshistos.o $(SOPHYAOBJP)fitsinoutfile.o $(SOPHYAOBJP)fitslocalmap.o $(SOPHYAOBJP)fitsmanager.o $(SOPHYAOBJP)fitsntuple.o $(SOPHYAOBJP)fitsspherehealpix.o $(SOPHYAOBJP)fitsspherethetaphi.o $(SOPHYAOBJP)fitstarray.o $(SOPHYAOBJP)fitsxntuple.o $(SOPHYAOBJP)swfitsdtable.o6 rm -f $(SOPHYAOBJP)datacirclefits.o $(SOPHYAOBJP)fabtcolread.o $(SOPHYAOBJP)fabtwriter.o $(SOPHYAOBJP)fboloread.o $(SOPHYAOBJP)fbtntintf.o $(SOPHYAOBJP)fiosinit.o $(SOPHYAOBJP)fitsbntbllineRW.o $(SOPHYAOBJP)fitsfile.o $(SOPHYAOBJP)fitsgenedata.o $(SOPHYAOBJP)fitshdtable.o $(SOPHYAOBJP)fitshisterr.o $(SOPHYAOBJP)fitshistos.o $(SOPHYAOBJP)fitsinoutfile.o $(SOPHYAOBJP)fitslocalmap.o $(SOPHYAOBJP)fitsmanager.o $(SOPHYAOBJP)fitsntuple.o $(SOPHYAOBJP)fitsspherehealpix.o $(SOPHYAOBJP)fitsspherethetaphi.o $(SOPHYAOBJP)fitstarray.o $(SOPHYAOBJP)swfitsdtable.o 7 7 rm -f $(SOPHYALIBP)libFitsIOServer.a 8 8 9 $(SOPHYALIBP)libFitsIOServer.a : $(SOPHYAOBJP)datacirclefits.o $(SOPHYAOBJP)fabtcolread.o $(SOPHYAOBJP)fabtwriter.o $(SOPHYAOBJP)fboloread.o $(SOPHYAOBJP)fbtntintf.o $(SOPHYAOBJP)fiosinit.o $(SOPHYAOBJP)fits autoreader.o $(SOPHYAOBJP)fitsbntbllineRW.o $(SOPHYAOBJP)fitsfile.o $(SOPHYAOBJP)fitsgenedata.o $(SOPHYAOBJP)fitshdtable.o $(SOPHYAOBJP)fitshisterr.o $(SOPHYAOBJP)fitshistos.o $(SOPHYAOBJP)fitsinoutfile.o $(SOPHYAOBJP)fitslocalmap.o $(SOPHYAOBJP)fitsmanager.o $(SOPHYAOBJP)fitsntuple.o $(SOPHYAOBJP)fitsspherehealpix.o $(SOPHYAOBJP)fitsspherethetaphi.o $(SOPHYAOBJP)fitstarray.o $(SOPHYAOBJP)fitsxntuple.o $(SOPHYAOBJP)swfitsdtable.o9 $(SOPHYALIBP)libFitsIOServer.a : $(SOPHYAOBJP)datacirclefits.o $(SOPHYAOBJP)fabtcolread.o $(SOPHYAOBJP)fabtwriter.o $(SOPHYAOBJP)fboloread.o $(SOPHYAOBJP)fbtntintf.o $(SOPHYAOBJP)fiosinit.o $(SOPHYAOBJP)fitsbntbllineRW.o $(SOPHYAOBJP)fitsfile.o $(SOPHYAOBJP)fitsgenedata.o $(SOPHYAOBJP)fitshdtable.o $(SOPHYAOBJP)fitshisterr.o $(SOPHYAOBJP)fitshistos.o $(SOPHYAOBJP)fitsinoutfile.o $(SOPHYAOBJP)fitslocalmap.o $(SOPHYAOBJP)fitsmanager.o $(SOPHYAOBJP)fitsntuple.o $(SOPHYAOBJP)fitsspherehealpix.o $(SOPHYAOBJP)fitsspherethetaphi.o $(SOPHYAOBJP)fitstarray.o $(SOPHYAOBJP)swfitsdtable.o 10 10 $(AR) $(ARFLAGS) $@ $? 11 11 touch $(SOPHYAINCP)/SophyaConfInfo/libextsophya.objlist … … 184 184 $(CXXCOMPILE) $(CXXTEMPFLG) -o $@ fiosinit.cc 185 185 186 $(SOPHYAOBJP)fitsautoreader.o: fitsautoreader.cc \187 $(SOPHYAINCP)sopnamsp.h \188 $(SOPHYAINCP)pexceptions.h \189 $(SOPHYAINCP)machdefs.h fitsautoreader.h \190 $(SOPHYAINCP)machdefs.h \191 $(SOPHYAINCP)anydataobj.h fitsfile.h \192 $(SOPHYAINCP)ndatablock.h \193 $(SOPHYAINCP)anydataobj.h \194 $(SOPHYAINCP)dvlist.h \195 $(SOPHYAINCP)objfio.h \196 $(SOPHYAINCP)ppersist.h \197 $(SOPHYAINCP)pexceptions.h \198 $(SOPHYAINCP)gnumd5.h \199 $(SOPHYAINCP)ppfbinstream.h \200 $(SOPHYAINCP)rawstream.h \201 $(SOPHYAINCP)ppfnametag.h \202 $(SOPHYAINCP)mutyv.h \203 $(SOPHYAINCP)timestamp.h fitsinoutfile.h \204 $(SOPHYAINCP)FitsIO/fitsio.h \205 $(SOPHYAINCP)FitsIO/longnam.h fitshandler.h \206 fitstarray.h $(SOPHYAINCP)tarray.h \207 $(SOPHYAINCP)basarr.h \208 $(SOPHYAINCP)dvlist.h \209 $(SOPHYAINCP)ndatablock.h \210 $(SOPHYAINCP)utilarr.h fitsspherehealpix.h \211 $(SOPHYAINCP)spherehealpix.h \212 $(SOPHYAINCP)sphericalmap.h \213 $(SOPHYAINCP)smathconst.h \214 $(SOPHYAINCP)pixelmap.h \215 $(SOPHYAINCP)spherepos.h \216 $(SOPHYAINCP)unitvector.h \217 $(SOPHYAINCP)vector3d.h \218 $(SOPHYAINCP)longlat.h \219 $(SOPHYAINCP)utilgeom.h \220 $(SOPHYAINCP)datatype.h \221 $(SOPHYAINCP)tvector.h \222 $(SOPHYAINCP)tmatrix.h \223 $(SOPHYAINCP)tarray.h \224 $(SOPHYAINCP)tmatrix_tsnl.h \225 $(SOPHYAINCP)tvector_tsnl.h \226 $(SOPHYAINCP)HEALPixUtils.h fitsntuple.h \227 $(SOPHYAINCP)ntuple.h \228 $(SOPHYAINCP)ntupintf.h \229 $(SOPHYAINCP)ppersist.h fitsxntuple.h \230 $(SOPHYAINCP)xntuple.h fitslocalmap.h \231 $(SOPHYAINCP)localmap.h fbtntintf.h \232 $(SOPHYAINCP)ntupintf.h233 $(CXXCOMPILE) $(CXXTEMPFLG) -o $@ fitsautoreader.cc234 235 186 $(SOPHYAOBJP)fitsbntbllineRW.o: fitsbntbllineRW.cc \ 236 187 $(SOPHYAINCP)sopnamsp.h \ … … 633 584 $(CXXCOMPILE) $(CXXTEMPFLG) -o $@ fitstarray.cc 634 585 635 $(SOPHYAOBJP)fitsxntuple.o: fitsxntuple.cc \636 $(SOPHYAINCP)sopnamsp.h \637 $(SOPHYAINCP)pexceptions.h \638 $(SOPHYAINCP)machdefs.h fitsxntuple.h \639 $(SOPHYAINCP)machdefs.h \640 $(SOPHYAINCP)xntuple.h \641 $(SOPHYAINCP)ntupintf.h \642 $(SOPHYAINCP)ppersist.h \643 $(SOPHYAINCP)pexceptions.h \644 $(SOPHYAINCP)gnumd5.h \645 $(SOPHYAINCP)ppfbinstream.h \646 $(SOPHYAINCP)rawstream.h \647 $(SOPHYAINCP)ppfnametag.h \648 $(SOPHYAINCP)dvlist.h \649 $(SOPHYAINCP)objfio.h \650 $(SOPHYAINCP)anydataobj.h \651 $(SOPHYAINCP)mutyv.h \652 $(SOPHYAINCP)timestamp.h \653 $(SOPHYAINCP)anydataobj.h \654 $(SOPHYAINCP)ppersist.h fitsfile.h \655 $(SOPHYAINCP)ndatablock.h \656 $(SOPHYAINCP)dvlist.h fitsinoutfile.h \657 $(SOPHYAINCP)FitsIO/fitsio.h \658 $(SOPHYAINCP)FitsIO/longnam.h fitshandler.h659 $(CXXCOMPILE) $(CXXTEMPFLG) -o $@ fitsxntuple.cc660 661 586 $(SOPHYAOBJP)swfitsdtable.o: swfitsdtable.cc swfitsdtable.h \ 662 587 $(SOPHYAINCP)basedtable.h \ -
trunk/SophyaExt/IFFTW/fftw3server.cc
r3412 r3572 209 209 sa_size_t n = in.NElts(); 210 210 sa_size_t ncs = n/2+1; 211 sa_size_t nc = (n%2 != 0) ? n/2+1 : n/2;211 //unused: sa_size_t nc = (n%2 != 0) ? n/2+1 : n/2; 212 212 if (out.NElts() != ncs) 213 213 throw SzMismatchError("FFTWServer::ReShapetoCompl() - Wrong output array size !"); … … 398 398 sa_size_t n = in.NElts(); 399 399 sa_size_t ncs = n/2+1; 400 sa_size_t nc = (n%2 != 0) ? n/2+1 : n/2;400 //unused: sa_size_t nc = (n%2 != 0) ? n/2+1 : n/2; 401 401 if (out.NElts() != ncs) 402 402 throw SzMismatchError("FFTWServer::ReShapetoCompl() - Wrong output array size !"); -
trunk/SophyaExt/LinAlg/intflapack.cc
r3534 r3572 178 178 extern "C" { 179 179 // Le calculateur de workingspace 180 int_4 ilaenv(int_4 *ispec,c har *name,char *opts,int_4 *n1,int_4 *n2,int_4 *n3,int_4 *n4,180 int_4 ilaenv(int_4 *ispec,const char *name,const char *opts,int_4 *n1,int_4 *n2,int_4 *n3,int_4 *n4, 181 181 int_4 nc1,int_4 nc2); 182 182 … … 300 300 //////////////////////////////////////////////////////////////////////////////////// 301 301 template <class T> 302 int_4 LapackServer<T>::ilaenv_en_C(int_4 ispec,c har *name,char *opts,int_4 n1,int_4 n2,int_4 n3,int_4 n4)302 int_4 LapackServer<T>::ilaenv_en_C(int_4 ispec,const char *name,const char *opts,int_4 n1,int_4 n2,int_4 n3,int_4 n4) 303 303 { 304 304 int_4 nc1 = strlen(name), nc2 = strlen(opts), rc=0; … … 519 519 int_4 info; 520 520 521 int_4 minmn = (m < n) ? m : n;521 //unused: int_4 minmn = (m < n) ? m : n; 522 522 int_4 maxmn = (m > n) ? m : n; 523 523 int_4 maxmnrhs = (nrhs > maxmn) ? nrhs : maxmn; … … 788 788 int_4 m = a.Size(rowa); 789 789 int_4 n = a.Size(cola); 790 int_4 maxmn = (m > n) ? m : n;790 //unused: int_4 maxmn = (m > n) ? m : n; 791 791 int_4 minmn = (m < n) ? m : n; 792 792 -
trunk/SophyaExt/LinAlg/intflapack.h
r2906 r3572 40 40 int SVDDriver(TArray<T>& a, TArray<T> & s, 41 41 TArray<T>* up=NULL, TArray<T> * vtp=NULL); 42 int_4 ilaenv_en_C(int_4 ispec,c har *name,char *opts,int_4 n1,int_4 n2,int_4 n3,int_4 n4);42 int_4 ilaenv_en_C(int_4 ispec,const char *name,const char *opts,int_4 n1,int_4 n2,int_4 n3,int_4 n4); 43 43 int_4 type2i4(void *val,int nbytes); 44 44 -
trunk/SophyaExt/XAstroPack/mollweide.cc
r2730 r3572 48 48 \endverbatim 49 49 */ 50 double MollWeide_XY(double longitude,double latitude,double * xmoll,double*ymoll)50 double MollWeide_XY(double longitude,double latitude,double& xmoll,double& ymoll) 51 51 { 52 52 int i,niter=2; … … 90 90 91 91 t0 *= asgn; 92 *xmoll = 2.*M_SQRT2*(longitude-M_PI)/M_PI*cos(t0);93 *ymoll = M_SQRT2*sin(t0);92 xmoll = 2.*M_SQRT2*(longitude-M_PI)/M_PI*cos(t0); 93 ymoll = M_SQRT2*sin(t0); 94 94 95 95 return t0; … … 118 118 \endverbatim 119 119 */ 120 double MollWeide_LL(double xmoll,double ymoll,double * longitude,double*latitude)120 double MollWeide_LL(double xmoll,double ymoll,double& longitude,double& latitude) 121 121 { 122 122 double t,a,eps=1.e-20; 123 123 124 *longitude = *latitude = -999.;124 longitude = latitude = -999.; 125 125 if(xmoll<-2*M_SQRT2 || xmoll>2*M_SQRT2) return(-901.); 126 126 if(ymoll<-M_SQRT2 || ymoll>M_SQRT2) return(-902.); … … 130 130 131 131 a = (2.*t+sin(2.*t))/M_PI; if(a<-1.) a=-1.; else if(a>1.) a=1.; 132 *latitude = asin(a);132 latitude = asin(a); 133 133 134 134 a= cos(t); 135 135 if(fabs(a)<eps) return(-903.); 136 *longitude = M_PI*(xmoll/(2.*M_SQRT2*a)+1.);136 longitude = M_PI*(xmoll/(2.*M_SQRT2*a)+1.); 137 137 138 138 return(t); -
trunk/SophyaExt/XAstroPack/mollweide.h
r1811 r3572 1 1 #ifndef MOLLWEIDE_H_SEEN 2 2 #define MOLLWEIDE_H_SEEN 3 double MollWeide_XY(double longitude,double latitude,double * xmoll,double*ymoll);4 double MollWeide_LL(double xmoll,double ymoll,double * longitude,double*latitude);3 double MollWeide_XY(double longitude,double latitude,double& xmoll,double& ymoll); 4 double MollWeide_LL(double xmoll,double ymoll,double& longitude,double& latitude); 5 5 #endif -
trunk/SophyaExt/XAstroPack/tstmollweide.cc
r2743 r3572 20 20 a = 2.*tvrai+sin(2.*tvrai); 21 21 lat = asin(a/M_PI); 22 tcalc = MollWeide_XY(0.,lat, &x,&y);22 tcalc = MollWeide_XY(0.,lat,x,y); 23 23 if(fabs(tcalc-tvrai)>tll) {tll = fabs(tcalc-tvrai); txy=tvrai;} 24 24 if(file) fprintf(file,"%e %e %e %e\n",tvrai,tcalc,a,tcalc-tvrai); … … 28 28 29 29 for(lon=0.;lon<360.;lon+=2.5) for(lat=-90.;lat<=90.1;lat+=2.5) { 30 txy = MollWeide_XY(lon*M_PI/180.,lat*M_PI/180., &x,&y);31 tll = MollWeide_LL(x,y, &longitude,&latitude);30 txy = MollWeide_XY(lon*M_PI/180.,lat*M_PI/180.,x,y); 31 tll = MollWeide_LL(x,y,longitude,latitude); 32 32 longitude *= 180./M_PI; 33 33 latitude *= 180./M_PI;
Note:
See TracChangeset
for help on using the changeset viewer.