Changeset 3572 in Sophya for trunk/SophyaExt/XAstroPack/mollweide.cc

Timestamp:

Feb 7, 2009, 10:50:34 PM (17 years ago)

Author:

cmv

Message:

char* -> const char* pour regler les problemes de deprecated string const... + comparaison unsigned signed + suppression EVOL_PLANCK rz+cmv 07/02/2009

File:

• trunk/SophyaExt/XAstroPack/mollweide.cc (modified) (4 diffs)

trunk/SophyaExt/XAstroPack/mollweide.cc

@@ -48 +48 @@
  \endverbatim
 */
-double MollWeide_XY(double longitude,double latitude,double* xmoll,double*ymoll)
+double MollWeide_XY(double longitude,double latitude,double& xmoll,double& ymoll)
 {
  int i,niter=2;
@@ -90 +90 @@
 
  t0 *= asgn;
- *xmoll = 2.*M_SQRT2*(longitude-M_PI)/M_PI*cos(t0);
- *ymoll = M_SQRT2*sin(t0);
+ xmoll = 2.*M_SQRT2*(longitude-M_PI)/M_PI*cos(t0);
+ ymoll = M_SQRT2*sin(t0);
 
  return t0;
@@ -118 +118 @@
  \endverbatim
 */
-double MollWeide_LL(double xmoll,double ymoll,double* longitude,double* latitude)
+double MollWeide_LL(double xmoll,double ymoll,double& longitude,double& latitude)
 {
  double t,a,eps=1.e-20;
 
- *longitude = *latitude = -999.;
+ longitude = latitude = -999.;
  if(xmoll<-2*M_SQRT2 || xmoll>2*M_SQRT2) return(-901.);
  if(ymoll<-M_SQRT2 || ymoll>M_SQRT2) return(-902.);
@@ -130 +130 @@
 
  a = (2.*t+sin(2.*t))/M_PI; if(a<-1.) a=-1.; else if(a>1.) a=1.;
- *latitude = asin(a);
+ latitude = asin(a);
 
  a= cos(t);
  if(fabs(a)<eps) return(-903.);
- *longitude = M_PI*(xmoll/(2.*M_SQRT2*a)+1.);
+ longitude = M_PI*(xmoll/(2.*M_SQRT2*a)+1.);
 
  return(t);