Index: /trunk/psLib/src/astro/psCoord.c
===================================================================
--- /trunk/psLib/src/astro/psCoord.c	(revision 2601)
+++ /trunk/psLib/src/astro/psCoord.c	(revision 2602)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.37 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-12-02 21:12:51 $
+*  @version $Revision: 1.38 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2004-12-03 19:43:43 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -208,7 +208,9 @@
     y = coord->d;
     dx = x - transform->phiP;
-    sinY = cos(y) * sin(dx) * transform->sinDeltaP + sin(y) * transform->cosDeltaP;
+    sinY = cos(y) * sin(dx) * transform->sinDeltaP +
+           sin(y) * transform->cosDeltaP;
     cosY = sqrt(1.0 - sinY * sinY);
-    sinX = (cos(y) * sin(dx) * transform->cosDeltaP - sin(y) * transform->sinDeltaP) / cos(y);
+    sinX = (cos(y) * sin(dx) * transform->cosDeltaP -
+            sin(y) * transform->sinDeltaP) / cos(y);
     cosX = cos(y) * cos(dx) / cos(y);
 
Index: /trunk/psLib/src/astronomy/Makefile
===================================================================
--- /trunk/psLib/src/astronomy/Makefile	(revision 2601)
+++ /trunk/psLib/src/astronomy/Makefile	(revision 2602)
@@ -3,8 +3,12 @@
 endif
 
+# XXX: Must put appropiate conditional here:
+LIBXMLDIR = /usr/include/libxml2
+LIBXMLDIR = /usr/local/include/libxml2
+
 TARGET_STATIC  = libpslib.a
 include ../Makefile.Globals
 CFLAGS := $(CFLAGS_RELOC) -I. -I../sysUtils -I../collections -I../dataManip \
-    -I.. -I../image -I../fileUtils -I/usr/local/include/libxml2 \
+    -I.. -I../image -I../fileUtils -I$(LIBXMLDIR) \
     -DTIME_CONFIG_FILE="\"$(TIME_CONFIG_FILE)\""
 
Index: /trunk/psLib/src/astronomy/psCoord.c
===================================================================
--- /trunk/psLib/src/astronomy/psCoord.c	(revision 2601)
+++ /trunk/psLib/src/astronomy/psCoord.c	(revision 2602)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.37 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-12-02 21:12:51 $
+*  @version $Revision: 1.38 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2004-12-03 19:43:43 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -208,7 +208,9 @@
     y = coord->d;
     dx = x - transform->phiP;
-    sinY = cos(y) * sin(dx) * transform->sinDeltaP + sin(y) * transform->cosDeltaP;
+    sinY = cos(y) * sin(dx) * transform->sinDeltaP +
+           sin(y) * transform->cosDeltaP;
     cosY = sqrt(1.0 - sinY * sinY);
-    sinX = (cos(y) * sin(dx) * transform->cosDeltaP - sin(y) * transform->sinDeltaP) / cos(y);
+    sinX = (cos(y) * sin(dx) * transform->cosDeltaP -
+            sin(y) * transform->sinDeltaP) / cos(y);
     cosX = cos(y) * cos(dx) / cos(y);
 
Index: /trunk/psLib/src/image/psImageStats.c
===================================================================
--- /trunk/psLib/src/image/psImageStats.c	(revision 2601)
+++ /trunk/psLib/src/image/psImageStats.c	(revision 2602)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.55 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-11-24 22:14:39 $
+ *  @version $Revision: 1.56 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-12-03 19:43:43 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -127,4 +127,6 @@
 0:512 to -1:1.  This routine takes as input an integer N and produces as
 output a vector of evenly spaced floating point values between -1.0:1.0.
+ 
+XXX: Use the p_psNormalizeVector here?
  *****************************************************************************/
 double* calcScaleFactors(psS32 n)
@@ -253,13 +255,6 @@
                         pixel = input->data.F64[x][y];
                     }
-                    sums[i][j] += pixel *
-                                  psPolynomial1DEval(
-                                      chebPolys[i],
-                                      rScalingFactors[x]
-                                  ) *
-                                  psPolynomial1DEval(
-                                      chebPolys[j],
-                                      cScalingFactors[y]
-                                  );
+                    sums[i][j] += pixel * psPolynomial1DEval(chebPolys[i],rScalingFactors[x]) *
+                                  psPolynomial1DEval(chebPolys[j], cScalingFactors[y]);
                 }
             }
@@ -375,5 +370,5 @@
     // Free the Chebyshev polynomials that were
     // created in this routine.
-    // XXX: Use dtatic data structures here.
+    // XXX: Use static data structures here.
     for (i = 0; i < maxChebyPoly; i++) {
         psFree(chebPolys[i]);
Index: /trunk/psLib/src/imageops/psImageStats.c
===================================================================
--- /trunk/psLib/src/imageops/psImageStats.c	(revision 2601)
+++ /trunk/psLib/src/imageops/psImageStats.c	(revision 2602)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.55 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-11-24 22:14:39 $
+ *  @version $Revision: 1.56 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-12-03 19:43:43 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -127,4 +127,6 @@
 0:512 to -1:1.  This routine takes as input an integer N and produces as
 output a vector of evenly spaced floating point values between -1.0:1.0.
+ 
+XXX: Use the p_psNormalizeVector here?
  *****************************************************************************/
 double* calcScaleFactors(psS32 n)
@@ -253,13 +255,6 @@
                         pixel = input->data.F64[x][y];
                     }
-                    sums[i][j] += pixel *
-                                  psPolynomial1DEval(
-                                      chebPolys[i],
-                                      rScalingFactors[x]
-                                  ) *
-                                  psPolynomial1DEval(
-                                      chebPolys[j],
-                                      cScalingFactors[y]
-                                  );
+                    sums[i][j] += pixel * psPolynomial1DEval(chebPolys[i],rScalingFactors[x]) *
+                                  psPolynomial1DEval(chebPolys[j], cScalingFactors[y]);
                 }
             }
@@ -375,5 +370,5 @@
     // Free the Chebyshev polynomials that were
     // created in this routine.
-    // XXX: Use dtatic data structures here.
+    // XXX: Use static data structures here.
     for (i = 0; i < maxChebyPoly; i++) {
         psFree(chebPolys[i]);