(JEC) 14/4/15 fix bug in GetPixelPositions; add library making in Makefile

Makefile

@@ -46,6 +46,11 @@ CXXOBJ = $(OBJ)laguerreBuilder.o \
 	$(OBJ)lagSphericTransform.o \
 	$(OBJ)walltimer.o
 
+CXXSHOBJ = laguerreBuilder.o \
+	laguerreTransform.o \
+	lagSphericTransform.o \
+	walltimer.o
+
 
 #C++ common Headers
 CXXHDR = lagsht_exceptions.h \
@@ -68,6 +73,18 @@ $(OBJ)%.o: %.cc $(CXXHDR)
 	echo "compile... $<",
 	$(CXXCOMPILE) $< -o $@
 
+######################
+.PHONY: sharelib
+sharelib : $(CXXOBJ)
+	echo $(CXXSHOBJ) ; \
+	cd $(OBJ); \
+	$(CMDSHLCXX) -o $liblagsht.$(SLEXT) $(CXXSHOBJ) $(LDFLAGS)
+.PHONY: lib
+lib : $(OBJ)liblagsht.a
+
+$(OBJ)liblagsht.a : $(CXXOBJ)
+	$(AR) $(ARFLAGS) $@ $?
+
 ######################
 .PHONY:  lagsht_testsuite
 lagsht_testsuite: $(EXE)lagsht_testsuite

lagSphericTransform.cc

@@ -82,7 +82,7 @@ void BaseLagSphTransform::SetThetaPhiMap(string choice, int nrings, int nphi){
 }//SetThetaPhiMap
 
 BaseLagSphTransform::geom_t BaseLagSphTransform::GetPixelPositions() {
-  geom_t pixels(Npix_);
+  geom_t pixels;
   sharp_geom_info *ginfo; 
   ginfo = sharpjob_.get_geom_info();
   for (int i=0; i<ginfo->npairs; ++i) {
@@ -92,7 +92,7 @@ BaseLagSphTransform::geom_t BaseLagSphTransform::GetPixelPositions() {
       double phi0 = ginfo->pair[i].r1.phi0;
       for(int j=0; j<nph1; j++) {
 	double phi1 = phi0 + j*(2.0*M_PI/nph1);
-	cout << "pixels ( "<< theta1 << ", " << phi1 << ")" << endl;
+	//	cout << "pixels ( "<< theta1 << ", " << phi1 << ")" << endl;
 	pixels.push_back(make_pair(theta1,phi1));
       }
       if (ginfo->pair[i].r2.nph>0) {// the second ring in this pair exists
@@ -101,7 +101,7 @@ BaseLagSphTransform::geom_t BaseLagSphTransform::GetPixelPositions() {
 	double phi0 = ginfo->pair[i].r2.phi0;
 	for(int j=0; j<nph2; j++) {
 	  double phi2 = phi0 + j*(2.0*M_PI/nph2);
-	  cout << "pixels ( "<< theta2 << ", " << phi2 << ")" << endl;
+	  //	  cout << "pixels ( "<< theta2 << ", " << phi2 << ")" << endl;
 	  pixels.push_back(make_pair(theta2,phi2));
 	}
       }