Index: /trunk/psLib/configure.ac
===================================================================
--- /trunk/psLib/configure.ac	(revision 4143)
+++ /trunk/psLib/configure.ac	(revision 4144)
@@ -38,9 +38,4 @@
 dnl enable POSIX/XSI and C99 compliance
 CFLAGS="${CFLAGS=} -D_XOPEN_SOURCE=600 -D_POSIX_C_SOURCE=200112L -std=c99"
-
-AC_LANG_PUSH(Fortran 77)
-AC_PROG_F77([gfortran g77])
-AC_F77_WRAPPERS
-AC_LANG_POP(Fortran 77)
 
 AC_PROG_INSTALL
Index: /trunk/psLib/src/dataManip/psMatrix.c
===================================================================
--- /trunk/psLib/src/dataManip/psMatrix.c	(revision 4143)
+++ /trunk/psLib/src/dataManip/psMatrix.c	(revision 4144)
@@ -21,6 +21,6 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.29 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-05-25 20:26:55 $
+ *  @version $Revision: 1.30 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 02:38:02 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -195,5 +195,5 @@
     outImage = psImageRecycle(outImage, inImage->numCols, inImage->numRows, inImage->type.type);
 
-    PS_CHECK_SQUARE(inImage, psMatrixLUD_EXIT);
+    PS_CHECK_SQUARE(inImage, psMatrixLUD_EXIT); // gsl_linalg_LU_decomp would fail on non-square input.
     PS_CHECK_SQUARE(outImage, psMatrixLUD_EXIT);
 
@@ -223,5 +223,5 @@
 
     // Calculate LU decomposition
-    gsl_linalg_LU_decomp(lu, &perm, &signum);
+    gsl_linalg_LU_decomp(lu, &perm, &signum); // N.B., uses Gaussian Elimination with partial pivoting.
 
     // Copy GSL matrix data to psImage data
@@ -412,4 +412,6 @@
     // Perform multiplication
     gsl_linalg_matmult(m1, m2, m3);
+    // N.B., gsl_blas_sgemm/dgemm could have been used instead, e.g.,
+    // gsl_blas_dgemm(CblasNoTrans, CblasNoTrans, 1.0, m1, m2, m3)
 
     // Copy GSL matrix data to psImage data
Index: /trunk/psLib/src/math/psMatrix.c
===================================================================
--- /trunk/psLib/src/math/psMatrix.c	(revision 4143)
+++ /trunk/psLib/src/math/psMatrix.c	(revision 4144)
@@ -21,6 +21,6 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.29 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-05-25 20:26:55 $
+ *  @version $Revision: 1.30 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 02:38:02 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -195,5 +195,5 @@
     outImage = psImageRecycle(outImage, inImage->numCols, inImage->numRows, inImage->type.type);
 
-    PS_CHECK_SQUARE(inImage, psMatrixLUD_EXIT);
+    PS_CHECK_SQUARE(inImage, psMatrixLUD_EXIT); // gsl_linalg_LU_decomp would fail on non-square input.
     PS_CHECK_SQUARE(outImage, psMatrixLUD_EXIT);
 
@@ -223,5 +223,5 @@
 
     // Calculate LU decomposition
-    gsl_linalg_LU_decomp(lu, &perm, &signum);
+    gsl_linalg_LU_decomp(lu, &perm, &signum); // N.B., uses Gaussian Elimination with partial pivoting.
 
     // Copy GSL matrix data to psImage data
@@ -412,4 +412,6 @@
     // Perform multiplication
     gsl_linalg_matmult(m1, m2, m3);
+    // N.B., gsl_blas_sgemm/dgemm could have been used instead, e.g.,
+    // gsl_blas_dgemm(CblasNoTrans, CblasNoTrans, 1.0, m1, m2, m3)
 
     // Copy GSL matrix data to psImage data
Index: /trunk/psLib/test/image/Makefile.am
===================================================================
--- /trunk/psLib/test/image/Makefile.am	(revision 4143)
+++ /trunk/psLib/test/image/Makefile.am	(revision 4144)
@@ -48,5 +48,4 @@
 tst_psImageStructManip_SOURCES =  tst_psImageStructManip.c
 tst_psImageConvolve_SOURCES =  tst_psImageConvolve.c
-tst_psImageIO_SOURCES =  tst_psImageIO.c
 tst_psImageInterpolate_SOURCES =  tst_psImageInterpolate.c