Index: /branches/rel5/psModules/.cvsignore
===================================================================
--- /branches/rel5/psModules/.cvsignore (revision 3634)
+++ /branches/rel5/psModules/.cvsignore (revision 3634)
@@ -0,0 +1,9 @@
+autom4te.cache
+Makefile
+config.log
+config.status
+libtool
+Makefile.in
+aclocal.m4
+configure
+
Index: /branches/rel5/psModules/psmodule.kdevses
===================================================================
--- /branches/rel5/psModules/psmodule.kdevses (revision 3634)
+++ /branches/rel5/psModules/psmodule.kdevses (revision 3634)
@@ -0,0 +1,35 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Index: /branches/rel5/psModules/src/.cvsignore
===================================================================
--- /branches/rel5/psModules/src/.cvsignore (revision 3634)
+++ /branches/rel5/psModules/src/.cvsignore (revision 3634)
@@ -0,0 +1,7 @@
+.deps
+Makefile
+config.h
+stamp-h1
+Makefile.in
+config.h.in
+
Index: /branches/rel5/psModules/src/Makefile.am
===================================================================
--- /branches/rel5/psModules/src/Makefile.am (revision 3634)
+++ /branches/rel5/psModules/src/Makefile.am (revision 3634)
@@ -0,0 +1,12 @@
+lib_LIBRARIES = libpsmodule.a
+
+libpsmodule_a_SOURCES = pmFlatField.c pmMaskBadPixels.c pmSubtractBias.c \
+ pmSubtractBias.c pmNonLinear.c pmReadoutCombine.c pmSubtractSky.c pmObjects.c
+
+libpsmodule_a_HEADERS = pmFlatFieldErrors.h pmMaskBadPixelsErrors.h \
+ pmFlatField.h pmMaskBadPixels.h pmSubtractBias.h pmSubtractBias.h \
+ pmNonLinear.h pmReadoutCombine.h pmSubtractSky.h pmObjects.h
+
+EXTRA_DIST = psErrorCodes.dat
+
+libpsmodule_adir = $(includedir)
Index: /branches/rel5/psModules/src/pmObjects.h
===================================================================
--- /branches/rel5/psModules/src/pmObjects.h (revision 3634)
+++ /branches/rel5/psModules/src/pmObjects.h (revision 3634)
@@ -0,0 +1,263 @@
+/** @file pmObjects.h
+ *
+ * This file will ...
+ *
+ * @author GLG, MHPCC
+ *
+ * @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
+ * @date $Date: 2005-04-01 21:38:27 $
+ *
+ * Copyright 2004 Maui High Performance Computing Center, University of Hawaii
+ *
+ */
+
+#if !defined(PM_OBJECTS_H)
+#define PM_OBJECTS_H
+
+#if HAVE_CONFIG_H
+#include
+#endif
+
+#include
+#include
+#include "pslib.h"
+// XXX: have Desonia do this correctly.
+#include "../../psLib/src/dataManip/psConstants.h"
+
+typedef enum {
+ PM_PEAK_LONE, ///< Isolated peak.
+ PM_PEAK_EDGE, ///< Peak on edge.
+ PM_PEAK_FLAT, ///< Peak has equal-value neighbors.
+ PM_PEAK_UNDEF ///< Undefined.
+} psPeakType;
+
+typedef struct
+{
+ psS32 x; ///< X-coordinate of peak pixel.
+ psS32 y; ///< Y-coordinate of peak pixel.
+ psF32 counts; ///< Value of peak pixel (above sky?).
+ psPeakType class; ///< Description of peak.
+}
+psPeak;
+
+typedef struct
+{
+ psF32 x; ///< X-coord of centroid.
+ psF32 y; ///< Y-coord of centroid.
+ psF32 Sx; ///< x-second moment.
+ psF32 Sy; ///< y-second moment.
+ psF32 Sxy; ///< xy cross moment.
+ psF32 Sum; ///< Pixel sum above sky (background).
+ psF32 Peak; ///< Peak counts above sky.
+ psF32 Sky; ///< Sky level (background).
+ psS32 nPixels; ///< Number of pixels used.
+}
+psMoments;
+
+typedef enum {
+ PS_MODEL_GAUSS, ///< Regular 2-D Gaussian
+ PS_MODEL_PGAUSS, ///< Psuedo 2-D Gaussian
+ PS_MODEL_TWIST_GAUSS, ///< 2-D Twisted Gaussian
+ PS_MODEL_WAUSS, ///< 2-D Waussian
+ PS_MODEL_SERSIC, ///< Sersic
+ PS_MODEL_SERSIC_CORE, ///< Sersic Core
+ PS_MODEL_UNDEFINED ///< Undefined
+} psModelType;
+
+// XXX: It will be better if params, and dparams, were psVectors.
+typedef struct
+{
+ psModelType type; ///< Model to be used.
+ psS32 Nparams; ///< Number of parameters.
+ psF32 *params; ///< Paramater values.
+ psF32 *dparams; ///< Parameter errors.
+ psF32 chisq; ///< Fit chi-squared.
+}
+psModel;
+
+// XXX: What is this enum?
+typedef enum {
+ PS_SOURCE_PSFSTAR,
+ PS_SOURCE_GALAXY,
+ PS_SOURCE_DEFECT,
+ PS_SOURCE_SATURATED,
+ PS_SOURCE_SATSTAR,
+ PS_SOURCE_FAINTSTAR,
+ PS_SOURCE_BRIGHTSTAR,
+ PS_SOURCE_OTHER
+} psSourceType;
+
+typedef struct
+{
+ psPeak *peak; ///< Description of peak pixel.
+ psImage *pixels; ///< Rectangular region including object pixels.
+ psImage *mask; ///< Mask which marks pixels associated with objects.
+ psMoments *moments; ///< Basic moments measure for the object.
+ psModel *models; ///< Model parameters and type.
+ psSourceType type; ///< Best identification of object.
+}
+psSource;
+
+psPeak *pmPeakAlloc(psS32 x, ///< Row-coordinate in image space
+ psS32 y, ///< Col-coordinate in image space
+ psF32 counts, ///< The value of the peak pixel
+ psPeakType class ///< The type of peak pixel
+ );
+
+psMoments *pmMomentsAlloc();
+psModel *pmModelAlloc(psModelType type);
+psSource *pmSourceAlloc();
+
+/******************************************************************************
+pmFindVectorPeaks(vector, threshold): Find all local peaks in the given vector
+above the given threshold. Returns a vector of type PS_TYPE_U32 containing
+the location (x value) of all peaks.
+ *****************************************************************************/
+psVector *pmFindVectorPeaks(const psVector *vector, ///< The input vector (psF32)
+ psF32 threshold ///< Threshold above which to find a peak
+ );
+
+/******************************************************************************
+pmFindImagePeaks(image, threshold): Find all local peaks in the given psImage
+above the given threshold. Returns a psList containing the location (x/y
+value) of all peaks.
+ *****************************************************************************/
+psList *pmFindImagePeaks(const psImage *image, ///< The input image where peaks will be found (psF32)
+ psF32 threshold ///< Threshold above which to find a peak
+ );
+
+/******************************************************************************
+psCullPeaks(peaks, maxValue, valid): eliminate peaks from the psList that have
+a peak value above the given maximum, or fall outside the valid region.
+
+XXX: Do we free the psList elements of those culled peaks?
+ *****************************************************************************/
+psList *pmCullPeaks(psList *peaks, ///< The psList of peaks to be culled
+ psF32 maxValue, ///< Cull peaks above this value
+ const psRegion *valid ///< Cull peaks otside this psRegion
+ );
+
+/******************************************************************************
+psSource *pmSourceLocalSky(image, peak, innerRadius, outerRadius):
+
+ *****************************************************************************/
+psSource *pmSourceLocalSky(const psImage *image, ///< The input image (psF32)
+ const psPeak *peak, ///< The peak for which the psSource struct is created.
+ psStatsOptions statsOptions, ///< The statistic used in calculating the background sky
+ psF32 innerRadius, ///< The inner radius of the suqare annulus for calculating sky
+ psF32 outerRadius ///< The outer radius of the suqare annulus for calculating sky
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+psSource *pmSourceMoments(psSource *source, ///< The input psSource for which moments will be computed
+ psF32 radius ///< Use a circle of pixels around the peak
+ );
+
+/******************************************************************************
+pmSourceRoughClass(pmArray *source, psMetaDeta *metadata): make a guess at the
+source classification.
+ *****************************************************************************/
+bool pmSourceRoughClass(psArray *source, ///< The input psSource
+ psMetadata *metadata ///< Contains classification parameters
+ );
+/******************************************************************************
+pmSourceSetPixelCircle(source, image, radius)
+ *****************************************************************************/
+bool pmSourceSetPixelCircle(psSource *source, ///< The input psSource
+ const psImage *image, ///< The input image (psF32)
+ psF32 radius ///< The radius of the circle
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+bool pmSourceModelGuess(psSource *source, ///< The input psSource
+ const psImage *image, ///< The input image (psF32)
+ psModelType model ///< The type of model to be created.
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+typedef enum {
+ PS_CONTOUR_CRUDE,
+} pmContourType;
+
+psArray *pmSourceContour(psSource *source, ///< The input psSource
+ const psImage *image, ///< The input image (psF32) (this arg should be removed)
+ psF32 level, ///< The level of the contour
+ pmContourType mode ///< Currently this must be PS_CONTOUR_CRUDE
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+bool pmSourceFitModel(psSource *source, ///< The input psSource
+ const psImage *image ///< The input image (psF32)
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+bool pmSourceAddModel(psImage *image, ///< The opuut image (psF32)
+ psSource *source, ///< The input psSource
+ bool center ///< A boolean flag that determines whether pixels are centered
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+bool pmSourceSubModel(psImage *image, ///< The output image (psF32)
+ psSource *source, ///< The input psSource
+ bool center ///< A boolean flag that determines whether pixels are centered
+ );
+
+/******************************************************************************
+XXX: Why only *x argument?
+ *****************************************************************************/
+psF32 pmMinLM_Gauss2D(psVector *deriv, ///< A possibly-NULL structure for the output derivatives
+ psVector *params, ///< A psVector which holds the parameters of this function
+ psVector *x ///< A psVector which holds the row/col coordinate
+ );
+
+/******************************************************************************
+XXX: Why only *x argument?
+ *****************************************************************************/
+psF32 pmMinLM_PsuedoGauss2D(psVector *deriv, ///< A possibly-NULL structure for the output derivatives
+ psVector *params, ///< A psVector which holds the parameters of this function
+ psVector *x ///< A psVector which holds the row/col coordinate
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+psF32 pmMinLM_Wauss2D(psVector *deriv, ///< A possibly-NULL structure for the output derivatives
+ psVector *params, ///< A psVector which holds the parameters of this function
+ psVector *x ///< A psVector which holds the row/col coordinate
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+psF32 pmMinLM_TwistGauss2D(psVector *deriv, ///< A possibly-NULL structure for the output derivatives
+ psVector *params, ///< A psVector which holds the parameters of this function
+ psVector *x ///< A psVector which holds the row/col coordinate
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+psF32 pmMinLM_Sersic(psVector *deriv, ///< A possibly-NULL structure for the output derivatives
+ psVector *params, ///< A psVector which holds the parameters of this function
+ psVector *x ///< A psVector which holds the row/col coordinate
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+psF32 pmMinLM_SersicCore(psVector *deriv, ///< A possibly-NULL structure for the output derivatives
+ psVector *params, ///< A psVector which holds the parameters of this function
+ psVector *x ///< A psVector which holds the row/col coordinate
+ );
+
+/******************************************************************************
+ *****************************************************************************/
+psF32 pmMinLM_PsuedoSersic(psVector *deriv, ///< A possibly-NULL structure for the output derivatives
+ psVector *params, ///< A psVector which holds the parameters of this function
+ psVector *x ///< A psVector which holds the row/col coordinate
+ );
+
+
+#endif
Index: /branches/rel5/psModules/test/.cvsignore
===================================================================
--- /branches/rel5/psModules/test/.cvsignore (revision 3634)
+++ /branches/rel5/psModules/test/.cvsignore (revision 3634)
@@ -0,0 +1,13 @@
+.deps
+.libs
+Makefile
+tst_pmFlatField
+tst_pmMaskBadPixels
+tst_pmNonLinear
+tst_pmSubtractBias
+temp
+Makefile.in
+tst_pmObjects01
+tst_pmReadoutCombine
+tst_pmSubtractSky
+
Index: /branches/rel5/psModules/test/Makefile.am
===================================================================
--- /branches/rel5/psModules/test/Makefile.am (revision 3634)
+++ /branches/rel5/psModules/test/Makefile.am (revision 3634)
@@ -0,0 +1,32 @@
+# Makefile for psModule tests
+
+AM_LDFLAGS = -L$(top_builddir)/src -lpsmodule $(LDFLAGS)
+AM_CFLAGS = -I$(top_srcdir)/src
+
+TESTS = tst_pmFlatField \
+ tst_pmMaskBadPixels \
+ tst_pmNonLinear \
+ tst_pmSubtractBias \
+ tst_pmReadoutCombine \
+ tst_pmSubtractSky \
+ tst_pmObjects01
+
+check_PROGRAMS = $(TESTS)
+
+TESTS_ENVIRONMENT = perl $(top_srcdir)/test/runTest --verified=$(srcdir)/verified
+
+tests: $(TESTS)
+
+EXTRA_DIST = runTest verified
+
+CLEANFILES = $(TESTS) temp/*
+
+tst_pmFlatField_SOURCES = tst_pmFlatField.c
+tst_pmMaskBadPixels_SOURCES = tst_pmMaskBadPixels.c
+tst_pmNonLinear_SOURCES = tst_pmNonLinear.c
+tst_pmSubtractBias_SOURCES = tst_pmSubtractBias.c
+tst_pmReadoutCombine_SOURCES = tst_pmReadoutCombine.c
+tst_pmSubtractSky_SOURCES = tst_pmSubtractSky.c
+tst_pmObjects01_SOURCES = tst_pmObjects01.c
+
+test: check
Index: /branches/rel5/psModules/test/runTest
===================================================================
--- /branches/rel5/psModules/test/runTest (revision 3634)
+++ /branches/rel5/psModules/test/runTest (revision 3634)
@@ -0,0 +1,358 @@
+#!/usr/bin/perl
+#
+# This is a perl script to execute a single test driver program. The script
+# will examine the return value of the test driver and capture STDOUT, STDERR
+# to files. If the return value of the test driver is not zero, the test is
+# deemed a failure. STDOUT and STDERR files will be place in a subdirectory
+# named temp and the files will be compared with verified STDOUT, STDERR files
+# in a subdirectory named verified. If there are no STDOUT, STDERR files in
+# the verified directory, the captured STDOUT, STDERR files will be scanned
+# for words which indicate a failure.
+#
+# Assumptions: A verified subdirectory with STDOUT, STDERR files exists for
+# the test driver specified in the arguements if an exact
+# match of the streams STDOUT, STDERR is necessary.
+#
+# Return values:
+# Bit mapped values
+# 0 Test run successfull and all tests passed
+# 1 Verified directory did not exist
+# 8 STDOUT files did not compare
+# 16 STDERR files did not compare
+# 32 Test driver doesn't exist or is not executable
+# 64 Test driver did not return zero
+#
+# SYNOSIS: runTest testDriverName
+#
+# $Revison: $ $Name: not supported by cvs2svn $
+# $Date: 2005-04-01 21:52:22 $
+#
+# Copyright 2004 Maui High Performance Computering Center, University of Hawaii
+#
+################################################################################
+
+# Provides functions for handling psS64 command line options
+use Getopt::Long;
+
+$verifiedDir = "verified";
+
+# Assign variables based on the presence of command line options to the script
+GetOptions(
+ "reset!" => \$reset,
+ "resetStderr!" => \$resetStderr,
+ "resetStdout!" => \$resetStdout,
+ "verified=s" => \$verifiedDir,
+ "help!" => \$help
+);
+
+if ($help || $#ARGV < 0) {
+ print
+"Usage: runTest [--reset] [--resetStderr] [--resetStdout] [--verified=DIR] testfile(s)\n\n";
+ exit(0);
+}
+
+if ($reset) {
+ $resetStderr = 1;
+ $resetStdout = 1;
+}
+
+# Initialize exit value
+$exitValue = 0;
+
+# Set up the PSLIB_ROOT environment variable if the user doesn't have it set
+if ( $ENV{'PSLIB_ROOT'} ) {
+
+ # Add PSLIB_ROOT/lib to LD_LIBRARY_PATH and DYLD_LIBRARY_PATH environment vars
+ $ENV{'LD_LIBRARY_PATH'} = "$ENV{'PSLIB_ROOT'}/src/.libs:$ENV{'LD_LIBRARY_PATH'}";
+ $ENV{'DYLD_LIBRARY_PATH'} = "$ENV{'PSLIB_ROOT'}/src/.libs:$ENV{'DYLD_LIBRARY_PATH'}";
+}
+
+# Loop through the arguements passed to the script
+foreach (@ARGV) {
+ chomp;
+ $testFile = $_;
+
+ # Check if a temp directory already exists in the current work directory
+ if ( !( -e "temp" ) ) {
+
+ # Create temp directory to store STDOUT, STDERR files during test run
+ `mkdir temp`;
+
+ # Display message of new directory created
+ print("Creating temp directory\n");
+ }
+
+ # Check if a verified directory exists in the current work directory
+ if ( !( -e $verifiedDir ) ) {
+
+ # Display message that verified subdirectory doesn't exist
+ print(" Verified directory doesn't exist.\n");
+
+ # Exit script since the test cannot be run with verified files
+ $exitValue = 1;
+ }
+
+ # Check if the test driver file exists and is executable
+ if ( ( -e $testFile ) && ( -x $testFile ) ) {
+
+ # Display message the test driver is being executed
+ print("--- Executing test driver: $testFile\n");
+
+ # Invoke the test driver
+ system("./$testFile 1> temp/$testFile.stdout 2> temp/$testFile.stderr");
+
+ # Check the return value of the test driver. A value other than 0
+ # indicates failure of the test driver unless the test driver is name
+ # with a leading 'a' which indicates the test driver is expecting an
+ # abort condition to terminate the driver.
+ if ( ( $? != 0 && ( $testFile !~ /^A/i ) )
+ || ( $? == 0 && ( $testFile =~ /^A/i ) ) )
+ {
+
+ # Display failure message with return value to user
+ if ( $? != 0 && ( $testFile !~ /^A/i ) ) {
+ print("Failed - Test Driver returned $?, expected 0.\n");
+ }
+ elsif ( $? == 0 && ( $testFile =~ /^A/i ) ) {
+ print("Failed - Test Driver returned $?, expected abort.\n");
+ }
+ $exitValue |= 64;
+ }
+
+ # Test driver succeeded.
+
+ # Create filtered version of stdout and stderr
+
+ # Open the STDOUT file for reading
+ open( OUTFILE, "< temp/$testFile.stdout" );
+
+ # Open mod file to place filtered STDOUT
+ open( MODFILE, "> temp/$testFile.stdout.mod" );
+ open( MODFILE2, "> $verifiedDir/$testFile.stdout" ) if $resetStdout;
+
+ # Replace the variable date, time and host information with constants
+ $hostname = `hostname`;
+ chop $hostname;
+ while () {
+ s/\s+\d+:\d+:\d+\w//g;
+ s/\d+:\d+:\d+//g;
+ s/$hostname\s*//g;
+ s/: Line \d+/: Line /g;
+ s/\(.*\:\d+\)/\(FILE\:LINENO\)/g;
+ s/\s+[\_\-\/\.\w]+\:\d+/ FILE\:LINENO/g;
+ s/allocate \d+ bytes at/allocate bytes at/g;
+ s/v\d+.\d+.\d+/vX.X.X/g;
+
+ # Filter lines with malloc. This is an artifact of memory testing
+ # with the Mac testbed
+ if ( !m/\*\*\*\smalloc/ ) {
+ print MODFILE ($_);
+ print MODFILE2 ($_) if $resetStdout;
+ }
+ }
+
+ # Close mod file
+ close(MODFILE);
+ close(MODFILE2) if $resetStdout;
+
+ # Close STDERR file
+ close(OUTFILE);
+
+ # Open the STDERR file for reading
+ open( OUTFILE, "< temp/$testFile.stderr" );
+
+ # Open mod file to place filtered STDERR
+ open( MODFILE, "> temp/$testFile.stderr.mod" );
+ open( MODFILE2, "> $verifiedDir/$testFile.stderr" ) if $resetStderr;
+
+ # Replace the variable date, time and host information with constants
+ while () {
+ s/\s+\d+:\d+:\d+\w//g;
+ s/\d+:\d+:\d+//g;
+ s/$hostname\s*//g;
+ s/: Line \d+/: Line /g;
+ s/\(.*\:\d+\)/\(FILE\:LINENO\)/g;
+ s/\s+[\_\-\/\.\w]+\:\d+/ FILE\:LINENO/g;
+ s/allocate \d+ bytes at/allocate bytes at/g;
+ s/v\d+.\d+.\d+/vX.X.X/g;
+
+ # Filter lines with malloc. This is an artifact of memory testing
+ # with the Mac testbed
+ if ( !m/\*\*\*\smalloc/ ) {
+ print MODFILE ($_);
+ print MODFILE2 ($_) if $resetStderr;
+ }
+ }
+
+ # Close mod file
+ close(MODFILE);
+ close(MODFILE2) if $resetStderr;
+
+ # Close STDERR file
+ close(OUTFILE);
+
+ # Compare STDOUT capture with verified file
+ $exitValue |= &compareStream("$verifiedDir/$testFile.stdout");
+
+ # Check exit value to determine if verified file doesn't exist
+ if ( $exitValue & 2 ) {
+
+ # STDOUT verified doesn't exist. Search STDOUT capture
+ # for strings indicating error or failure
+ $exitValue |= &errorStrSearch("$testFile.stdout");
+ }
+
+ # Compare STDERR capture with verified file
+ $exitValue |= &compareStream("$verifiedDir/$testFile.stderr");
+
+ # Check exit value to determine if verified file doesn't exist
+ if ( $exitValue & 4 ) {
+
+ # STDERR verified doesn't exist. Search STDERR capture
+ # for strings indicating error or failure
+ $exitValue |= &errorStrSearch("$testFile.stderr");
+ }
+ }
+ else {
+
+ # Since test driver doesn't exist or is not executable then display
+ # message to user.
+ print(
+" Need to specify an executable test file within directory.\n"
+ );
+
+ # Exit value set to indicate test driver doesn't exist or not executable
+ $exitValue |= 32;
+ }
+}
+
+# Exit for the script with exit value
+# Ignore the first three bit flags in exitValue since there are not indicators
+# of a failed test
+if ( ( $exitValue >> 3 ) != 0 ) {
+ print("Test failed - return status = $exitValue\n\n");
+}
+else {
+ print("Test successful\n\n");
+}
+exit($exitValue);
+
+################################################################################
+#
+# SUBROUTINE: errorStrSearch
+#
+# Description: This subroutine will search the file specified in its
+# parameter for error strings and if they do exists
+# the appropriate value will be returned.
+#
+# Parameter(s): fileName - input file to search for strings
+#
+# Return: 0 - No error strings found
+# 64 - Error strings found in STDOUT
+# 128 - Error strings found in STDERR
+#
+###############################################################################
+
+sub errorStrSearch {
+ local ($fileName) = @_;
+ local ($returnVal) = 0;
+
+ # Open the captured file
+ open( CAPT_FILE, ") {
+ if ( m/FAIL/i
+ || m/FAULT/i
+ || m/ERROR/i
+ || m/Not Found/i
+ || m/SIGNAL/i
+ || m/NO SUCH FILE/i
+ || m/\|E\|/i
+ || m/\|A\|/i )
+ {
+ if ( $fileName =~ m/out/ ) {
+ print(
+ " Failed - File $fileName contains error strings.\n"
+ );
+ $returnVal = 64;
+ }
+ elsif ( $fileName =~ m/err/ ) {
+ print(
+ " Failed - File $fileName contains error strings.\n"
+ );
+ $returnVal = 128;
+ }
+ last;
+ }
+ }
+
+ # Close the capture file
+ close(CAPT_FILE);
+
+ # Return the return value
+ return ($returnVal);
+}
+
+################################################################################
+#
+# SUBROUTINE: compareStream
+#
+# Description: This subroutine will compare the captured stream file with
+# a file in the verified directory if one exists.
+#
+# Parameter(s): streamFile - input file to compare
+#
+# Return: 0 - Compare successful
+# 8 - STDOUT verified file doesn't compare
+# 16 - STDERR verified file doesn't compare
+#
+###############################################################################
+
+sub compareStream {
+ local ($streamFile) = @_;
+ local ($returnVal) = 0;
+ local ($tempFile) = "";
+
+ # Check for existence of verified STD stream files
+ if ( -e $streamFile ) {
+
+ # Verified STD stream file exists
+
+ # Create name of the temp file to compare
+ $tempFile = $streamFile;
+ $tempFile =~ s/$verifiedDir/temp/;
+ $tempFile = $tempFile . ".mod";
+
+ # Perform difference on the STD stream files
+ $diffstdout = `diff $streamFile $tempFile`;
+
+ # Check the return value of the difference
+ if ( $? != 0 ) {
+
+ # Difference of STD stream files failed
+
+ # Check for STDOUT in file name to convey appropirate return value
+ if ( $streamFile =~ /out$/ ) {
+
+ # Display message of the failure of difference to user
+ print(" Failed - STDOUT difference\n$diffstdout\n");
+
+ # Exit value to indicate STDOUT did not compare
+ $returnVal |= 8;
+ }
+ elsif ( $streamFile =~ /err$/ ) {
+
+ # Display message of the failure of difference to user
+ print(" Failed - STDERR difference\n$diffstdout\n");
+
+ # Exit value to indicate STDERR did not compare
+ $returnVal |= 16;
+ }
+ }
+ }
+
+ # Return the result of the compare
+ return ($returnVal);
+}
+
Index: /branches/rel5/psModules/test/tst_pmObjects01.c
===================================================================
--- /branches/rel5/psModules/test/tst_pmObjects01.c (revision 3634)
+++ /branches/rel5/psModules/test/tst_pmObjects01.c (revision 3634)
@@ -0,0 +1,1757 @@
+/** @file tst_pmFindObjects.c
+ *
+ * @brief Contains the tests for pmSubtractSky.c:
+ *
+ * test00: This code will ...
+ *
+ * @author GLG, MHPCC
+ *
+ * XXX: Must test
+ * pmSourceRoughClass
+ *
+ * XXX: Must test output results for many other functions.
+ *
+ * XXX: There are many cases where row/col can be switched in the code.
+ * We must test that here by using non-square images. All tests
+ * in this file should be run with non-square images.
+ *
+ * XXX: Memory leaks are not being caught. If I allocated a psVector in these functions
+ * abd never deallocate, no error is generated.
+ *
+ * @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
+ * @date $Date: 2005-04-01 21:39:06 $
+ *
+ * Copyright 2004 Maui High Performance Computing Center, University of Hawaii
+ */
+#include "psTest.h"
+#include "pslib.h"
+#include "pmObjects.h"
+#define NUM_ROWS 10
+#define NUM_COLS 10
+#define ERROR_TOLERANCE 1.0
+static int test00(void);
+static int test01(void);
+static int test02(void);
+static int test03(void);
+static int test04(void);
+static int test05(void);
+static int test06(void);
+static int test07(void);
+static int test08(void);
+static int test09(void);
+static int test15(void);
+static int test16(void);
+static int test20(void);
+testDescription tests[] = {
+ {test00, 000, "pmObjects: structure allocators and deallocators", true, false},
+ {test01, 001, "pmObjects: psFindVectorPeaks()", true, false},
+ {test02, 001, "pmObjects: psFindImagePeaks()", true, false},
+ {test03, 001, "pmObjects: pmCullPeaks()", true, false},
+ {test04, 001, "pmObjects: pmSourceLocalSky()", true, false},
+ {test06, 001, "pmObjects: pmSourceSetPixelCircle()", true, false},
+ {test05, 001, "pmObjects: pmSourceMoments()", true, false},
+ {test07, 001, "pmObjects: pmMin()", true, false},
+ {test08, 001, "pmObjects: pmSourceModelGuess()", true, false},
+ {test09, 001, "pmObjects: pmSourceContour()", true, false},
+ {test15, 001, "pmObjects: pmSourceAddModel()", true, false},
+ {test16, 001, "pmObjects: pmSourceSubModel()", true, false},
+ {test20, 001, "pmObjects: pmSourceSubModel()", true, false},
+ {NULL}
+ };
+
+int main(int argc, char* argv[])
+{
+ return !runTestSuite(stderr, "Test Point Driver", tests, argc, argv);
+ test02();
+}
+
+/******************************************************************************
+test00(): Test the various allocators and deallocators.
+ *****************************************************************************/
+int test00( void )
+{
+ bool testStatus = true;
+ psTraceSetLevel(".", 0);
+
+ psPeak *tmpPeak = pmPeakAlloc(0, 0, 0.0, PM_PEAK_UNDEF);
+ if (tmpPeak == NULL) {
+ printf("TEST ERROR: pmPeakAlloc() returned a NULL psPeak\n");
+ testStatus = false;
+ }
+ // XXX: Test that values were set properly.
+ psFree(tmpPeak);
+
+ psMoments *tmpMoments = pmMomentsAlloc();
+ if (tmpMoments == NULL) {
+ printf("TEST ERROR: pmMomentsAlloc() returned a NULL psMoments\n");
+ testStatus = false;
+ }
+ psFree(tmpMoments);
+
+ // XXX: Must test this with variety of models, and verify param length.
+ psModel *tmpModel = pmModelAlloc(PS_MODEL_GAUSS);
+ if (tmpModel == NULL) {
+ printf("TEST ERROR: pmModelAlloc() returned a NULL psModel\n");
+ testStatus = false;
+ }
+ // XXX: Test that values were set properly.
+ psFree(tmpModel);
+
+ // XXX: Must test this with variety of models, and verify param length.
+ tmpModel = pmModelAlloc(PS_MODEL_PGAUSS);
+ if (tmpModel == NULL) {
+ printf("TEST ERROR: pmModelAlloc() returned a NULL psModel\n");
+ testStatus = false;
+ }
+ // XXX: Test that values were set properly.
+ psFree(tmpModel);
+
+ // XXX: Must test this with variety of models, and verify param length.
+ tmpModel = pmModelAlloc(PS_MODEL_TWIST_GAUSS);
+ if (tmpModel == NULL) {
+ printf("TEST ERROR: pmModelAlloc() returned a NULL psModel\n");
+ testStatus = false;
+ }
+ // XXX: Test that values were set properly.
+ psFree(tmpModel);
+
+ // XXX: Must test this with variety of models, and verify param length.
+ tmpModel = pmModelAlloc(PS_MODEL_WAUSS);
+ if (tmpModel == NULL) {
+ printf("TEST ERROR: pmModelAlloc() returned a NULL psModel\n");
+ testStatus = false;
+ }
+ // XXX: Test that values were set properly.
+ psFree(tmpModel);
+
+ // XXX: Must test this with variety of models, and verify param length.
+ tmpModel = pmModelAlloc(PS_MODEL_SERSIC);
+ if (tmpModel == NULL) {
+ printf("TEST ERROR: pmModelAlloc() returned a NULL psModel\n");
+ testStatus = false;
+ }
+ // XXX: Test that values were set properly.
+ psFree(tmpModel);
+
+ // XXX: Must test this with variety of models, and verify param length.
+ tmpModel = pmModelAlloc(PS_MODEL_SERSIC_CORE);
+ if (tmpModel == NULL) {
+ printf("TEST ERROR: pmModelAlloc() returned a NULL psModel\n");
+ testStatus = false;
+ }
+ // XXX: Test that values were set properly.
+ psFree(tmpModel);
+
+ psSource *tmpSource = pmSourceAlloc();
+ if (tmpSource == NULL) {
+ printf("TEST ERROR: pmSourceAlloc() returned a NULL psSource\n");
+ testStatus = false;
+ }
+ psFree(tmpSource);
+
+ return(testStatus);
+}
+
+#define VECTOR_SIZE 10
+bool test_pmFindVectorPeaks(int n)
+{
+ bool testStatus = true;
+ psVector *inData = psVectorAlloc(n, PS_TYPE_F32);
+ psVector *outData = NULL;
+
+ printf("-------------- Calling test_pmFindVectorPeaks on an %d size vector. --------------\n", n);
+ //
+ // Test first pixel peak.
+ //
+ printf("Test pmFindVectorPeaks() with a first-element peak.\n");
+ for (psS32 i = 0 ; i < n ; i++) {
+ inData->data.F32[i] = (float) (n-i);
+ }
+ inData->data.F32[0] = (float) n;
+
+ outData= pmFindVectorPeaks(inData, 0.0);
+ if (outData == NULL) {
+ printf("TEST ERROR: pmFindVectorPeaks returned a NULL psVector.\n");
+ testStatus = false;
+ } else {
+
+ if (outData->n != 1) {
+ printf("TEST ERROR: outData->n is %d\n", outData->n);
+ testStatus = false;
+ }
+ if (outData->data.U32[0] != 0) {
+ printf("TEST ERROR: Did not find peak at element 0.\n");
+ testStatus = false;
+ }
+ psFree(outData);
+ }
+
+ //
+ // Test first pixel peak, large threshold
+ //
+ printf("Test pmFindVectorPeaks() with a first-element peak, large threshold.\n");
+ for (psS32 i = 0 ; i < n ; i++) {
+ inData->data.F32[i] = (float) (n-i);
+ }
+ inData->data.F32[0] = (float) n;
+
+ outData= pmFindVectorPeaks(inData, (float) (n*n));
+ if (outData == NULL) {
+ printf("TEST ERROR: pmFindVectorPeaks returned a NULL psVector.\n");
+ testStatus = false;
+ } else {
+
+ if (outData->n != 0) {
+ printf("TEST ERROR: outData->n is %d\n", outData->n);
+ testStatus = false;
+ }
+ psFree(outData);
+
+ // Skip remaining tests if the input vector has length 1.
+ if (n == 1) {
+ psFree(inData);
+ return(testStatus);
+ }
+ }
+
+ //
+ // Test last pixel peak.
+ //
+ printf("Test pmFindVectorPeaks() with a last-element peak.\n");
+ for (psS32 i = 0 ; i < n ; i++) {
+ inData->data.F32[i] = (float) (i);
+ }
+ inData->data.F32[n-1] = (float) n;
+
+ outData= pmFindVectorPeaks(inData, 0.0);
+ if (outData == NULL) {
+ printf("TEST ERROR: pmFindVectorPeaks returned a NULL psVector.\n");
+ testStatus = false;
+ } else {
+
+ if (outData->n != 1) {
+ printf("TEST ERROR: outData->n is %d\n", outData->n);
+ testStatus = false;
+ }
+ if (outData->data.U32[0] != n-1) {
+ printf("TEST ERROR: Did not find peak at element %d.\n", n-1);
+ testStatus = false;
+ }
+ psFree(outData);
+ }
+
+ //
+ // Test last pixel peak, large threshold.
+ //
+ printf("Test pmFindVectorPeaks() with a last-element peak, large threshold.\n");
+ for (psS32 i = 0 ; i < n ; i++) {
+ inData->data.F32[i] = (float) (i);
+ }
+ inData->data.F32[n-1] = (float) n;
+
+ outData= pmFindVectorPeaks(inData, (float) (n*n));
+ if (outData == NULL) {
+ printf("TEST ERROR: pmFindVectorPeaks returned a NULL psVector.\n");
+ testStatus = false;
+ } else {
+
+ if (outData->n != 0) {
+ printf("TEST ERROR: outData->n is %d\n", outData->n);
+ testStatus = false;
+ }
+ psFree(outData);
+ }
+
+ //
+ // Test interior peaks.
+ // Set all even number elements to be peaks.
+ //
+ printf("Test pmFindVectorPeaks() with all even-numbered elements peak.\n");
+ for (psS32 i = 0 ; i < n ; i++) {
+ if (0 == i%2) {
+ inData->data.F32[i] = (float) (2 * i);
+ } else {
+ inData->data.F32[i] = (float) (i);
+ }
+ }
+ inData->data.F32[0] = (float) n;
+
+
+ outData= pmFindVectorPeaks(inData, 0.0);
+ if (outData == NULL) {
+ printf("TEST ERROR: pmFindVectorPeaks returned a NULL psVector.\n");
+ testStatus = false;
+ } else {
+
+ if (outData->n != n/2) {
+ printf("TEST ERROR: outData->n is %d\n", outData->n);
+ testStatus = false;
+ }
+
+ for (psS32 i = 0 ; i < outData->n ; i++) {
+ if (outData->data.U32[i] != (2 * i)) {
+ printf("TEST ERROR: the %d-th peak is element number %d\n", i, outData->data.U32[i]);
+ testStatus = false;
+ }
+ }
+ psFree(outData);
+ }
+
+ //
+ // Test interior peaks, with threshold = n*n.
+ // Should generate an empty output psVector.
+ //
+ printf("Test pmFindVectorPeaks() with all even-numbered elements peak, large threshold.\n");
+ outData= pmFindVectorPeaks(inData, (float) (n*n));
+ if (outData == NULL) {
+ printf("TEST ERROR: pmFindVectorPeaks returned a NULL psVector.\n");
+ testStatus = false;
+ } else {
+
+ if (outData->n != 0) {
+ printf("TEST ERROR: outData->n is %d\n", outData->n);
+ testStatus = false;
+ }
+ psFree(outData);
+ }
+
+
+ psFree(inData);
+ return(testStatus);
+}
+
+
+int test01( void )
+{
+ bool testStatus = true;
+ psVector *tmpVec = NULL;
+ psVector *tmpVecF64 = psVectorAlloc(10, PS_TYPE_F64);
+ psVector *tmpVecEmpty = psVectorAlloc(0, PS_TYPE_F32);
+
+ psTraceSetLevel(".", 0);
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmFindVectorPeaks with NULL psVector. Should generate error and return NULL.\n");
+ tmpVec = pmFindVectorPeaks(NULL, 0.0);
+ if (tmpVec != NULL) {
+ printf("TEST ERROR: pmFindVectorPeaks() returned a non-NULL psVector.\n");
+ testStatus = false;
+ psFree(tmpVec);
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmFindVectorPeaks with empty psVector. Should generate error and return NULL.\n");
+ tmpVec = pmFindVectorPeaks(tmpVecEmpty, 0.0);
+ if (tmpVec != NULL) {
+ printf("TEST ERROR: pmFindVectorPeaks() returned a non-NULL psVector.\n");
+ testStatus = false;
+ psFree(tmpVec);
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmFindVectorPeaks with PS_TYPE_F64 psVector. Should generate error and return NULL.\n");
+ tmpVec = pmFindVectorPeaks(tmpVecF64, 0.0);
+ if (tmpVec != NULL) {
+ printf("TEST ERROR: pmFindVectorPeaks() returned a non-NULL psVector.\n");
+ testStatus = false;
+ psFree(tmpVec);
+ }
+ testStatus&= test_pmFindVectorPeaks(1);
+ testStatus&= test_pmFindVectorPeaks(VECTOR_SIZE);
+
+ psFree(tmpVecF64);
+ psFree(tmpVecEmpty);
+ return(testStatus);
+}
+
+// XXX: Add tests for the PEAK_TYPE.
+// XXX: Add interior peaks, edge peaks, and flat peaks.
+bool test_pmFindImagePeaks(int numRows, int numCols)
+{
+ printf("-------------- Calling test_pmFindVectorPeaks on an %d-by-%d image. --------------\n", numRows, numCols);
+ bool testStatus = true;
+ psImage *inData = psImageAlloc(numCols, numRows, PS_TYPE_F32);
+ psList *outData = NULL;
+
+ //
+ // Test first pixel peak.
+ //
+ for (psS32 i = 0 ; i < numRows ; i++) {
+ for (psS32 j = 0 ; j < numCols ; j++) {
+ inData->data.F32[i][j] = PS_SQR(i - numRows/2) + PS_SQR(j-numCols/2);
+ }
+ }
+ inData->data.F32[0][0] = PS_SQR(numRows) + PS_SQR(numCols);
+ inData->data.F32[0][numCols-1] = PS_SQR(numRows) + PS_SQR(numCols);
+ inData->data.F32[numRows-1][0] = PS_SQR(numRows) + PS_SQR(numCols);
+ inData->data.F32[numRows-1][numCols-1] = PS_SQR(numRows) + PS_SQR(numCols);
+ for (psS32 i = 0 ; i < numRows ; i++) {
+ for (psS32 j = 0 ; j < numCols ; j++) {
+ printf("(%.1f) ", inData->data.F32[i][j]);
+ }
+ printf("\n");
+ }
+
+ outData = pmFindImagePeaks(inData, 0.0);
+
+ if (outData == NULL) {
+ printf("TEST ERROR: pmFindImagePeaks returned a NULL psList.\n");
+ testStatus = false;
+ } else {
+ if (outData->size != 4) {
+ printf("TEST ERROR: pmFindImagePeaks found only %d peaks (should be 4)\n", outData->size);
+ testStatus = false;
+ }
+
+ psListElem *tmpPeakLE = (psListElem *) outData->head;
+ while (tmpPeakLE != NULL) {
+ psPeak *tmpPeak = (psPeak *) tmpPeakLE->data;
+ if (((tmpPeak->x == 0) || (tmpPeak->x == numRows-1)) &&
+ ((tmpPeak->y == 0) || (tmpPeak->y == numCols-1))) {}
+ else {
+ printf("TEST ERROR: Peak at (%d, %d) (%f)\n", tmpPeak->x, tmpPeak->y, tmpPeak->counts);
+ testStatus = false;
+ }
+ tmpPeakLE = tmpPeakLE->next;
+ }
+ }
+
+ psFree(inData);
+ psFree(outData);
+ return(testStatus);
+}
+
+
+int test02( void )
+{
+ bool testStatus = true;
+ psList *tmpList = NULL;
+ psImage *tmpImageF64 = psImageAlloc(10, 10, PS_TYPE_F64);
+ psImage *tmpImageEmpty = psImageAlloc(0, 0, PS_TYPE_F32);
+
+ psTraceSetLevel(".", 0);
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmFindImagePeaks with NULL psImage. Should generate error and return NULL.\n");
+ tmpList = pmFindImagePeaks(NULL, 0.0);
+ if (tmpList != NULL) {
+ printf("TEST ERROR: pmFindImagePeaks() returned a non-NULL psImage.\n");
+ testStatus = false;
+ psFree(tmpList);
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmFindImagePeaks with empty psImage. Should generate error and return NULL.\n");
+ tmpList = pmFindImagePeaks(tmpImageEmpty, 0.0);
+ if (tmpList != NULL) {
+ printf("TEST ERROR: pmFindImagePeaks() returned a non-NULL psImage.\n");
+ testStatus = false;
+ psFree(tmpList);
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmFindImagePeaks with PS_TYPE_F64 psImage. Should generate error and return NULL.\n");
+ tmpList = pmFindImagePeaks(tmpImageF64, 0.0);
+ if (tmpList != NULL) {
+ printf("TEST ERROR: pmFindImagePeaks() returned a non-NULL psImage.\n");
+ testStatus = false;
+ psFree(tmpList);
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ // testStatus&= test_pmFindImagePeaks(1, 1);
+ // testStatus&= test_pmFindImagePeaks(2, 5);
+ // testStatus&= test_pmFindImagePeaks(5, 2);
+ testStatus&= test_pmFindImagePeaks(10, 10);
+
+
+ psFree(tmpImageF64);
+ psFree(tmpImageEmpty);
+ return(testStatus);
+}
+
+/******************************************************************************
+test03(): We first test pmCullPeaks() with various NULL and unallowable input
+parameters. Then we generate a list of peaks and test that pmCullPeaks()
+removes them correctly.
+ *****************************************************************************/
+int test03( void )
+{
+ bool testStatus = true;
+ psImage *imgData = psImageAlloc(NUM_COLS, NUM_ROWS, PS_TYPE_F32);
+ psList *outData = NULL;
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmCullPeaks with NULL psList. Should generate error and return NULL.\n");
+ outData = pmCullPeaks(NULL, 0.0, NULL);
+ if (outData != NULL) {
+ printf("TEST ERROR: pmCulPeaks() returned a non-NULL psList.\n");
+ testStatus = false;
+ }
+
+ //
+ // Set peaks in input image. All even-column and even-row pixels are
+ // set non-zero, all other pixels are set to zero.
+ //
+ psS32 numPeaksOrig = 0;
+ for (psS32 i = 0 ; i < NUM_ROWS ; i++) {
+ for (psS32 j = 0 ; j < NUM_COLS ; j++) {
+ if ((0 == i%2) && (0 == j%2)) {
+ imgData->data.F32[i][j] = (float) (i + 10);
+ numPeaksOrig++;
+ } else {
+ imgData->data.F32[i][j] = 0.0;
+ }
+ }
+ }
+ for (psS32 i = 0 ; i < NUM_ROWS ; i++) {
+ for (psS32 j = 0 ; j < NUM_COLS ; j++) {
+ printf("(%.1f) ", imgData->data.F32[i][j]);
+ }
+ printf("\n");
+ }
+ printf("Set %d peaks\n", numPeaksOrig);
+
+ //
+ // Call pmCullPeaks() with HUGE maxValue and NULL psRegion. Should not
+ // remove any peaks.
+ //
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmCullPeaks with large maxValue and NULL psRegion.\n");
+ outData = pmFindImagePeaks(imgData, 0.0);
+ outData = pmCullPeaks(outData, PS_MAX_F32, NULL);
+
+ if (outData == NULL) {
+ printf("TEST ERROR: pmCullPeaks() returned a non-NULL psList.\n");
+ testStatus = false;
+ return(testStatus);
+ }
+ if (outData->size != numPeaksOrig) {
+ printf("TEST ERROR (0): pmCullPeaks incorrectly removed peaks\n");
+ printf("The pmCullPeaks() output had %d peaks, should have had %d peaks.\n", outData->size, numPeaksOrig);
+ testStatus = false;
+ }
+ psFree(outData);
+
+ //
+ // Call pmCullPeaks() with TINY maxValue and NULL psRegion. Should
+ // remove all peaks.
+ //
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmCullPeaks with tiny maxValue and NULL psRegion.\n");
+ outData = pmFindImagePeaks(imgData, 0.0);
+ printf("pmFindImagePeaks found %d peaks\n", outData->size);
+ outData = pmCullPeaks(outData, 0.0, NULL);
+
+ if (outData == NULL) {
+ printf("TEST ERROR: pmCullPeaks() returned a non-NULL psList.\n");
+ testStatus = false;
+ return(testStatus);
+ }
+ if (outData->size != 0) {
+ printf("TEST ERROR (1): pmCullPeaks incorrectly removed peaks\n");
+ printf("The pmCullPeaks() output had %d peaks, should have had %d peaks.\n", outData->size, 0);
+ testStatus = false;
+ }
+ psFree(outData);
+
+ //
+ // Call pmCullPeaks() with HUGE maxValue and disjoint psRegion. Should
+ // not remove any peaks.
+ //
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmCullPeaks with large maxValue and disjoint psRegion.\n");
+ outData = pmFindImagePeaks(imgData, 0.0);
+ printf("pmFindImagePeaks found %d peaks\n", outData->size);
+ psRegion *tmpRegion = psRegionAlloc(10000.0, 20000.0, 10000.0, 20000.0);
+ outData = pmCullPeaks(outData, PS_MAX_F32, tmpRegion);
+
+ if (outData == NULL) {
+ printf("TEST ERROR: pmCullPeaks() returned a non-NULL psList.\n");
+ testStatus = false;
+ return(testStatus);
+ }
+ if (outData->size != numPeaksOrig) {
+ printf("TEST ERROR (2): pmCullPeaks incorrectly removed peaks\n");
+ printf("The pmCullPeaks() output had %d peaks, should have had %d peaks.\n", outData->size, numPeaksOrig);
+ testStatus = false;
+ }
+ psFree(outData);
+ psFree(tmpRegion);
+
+ //
+ // Call pmCullPeaks() with HUGE maxValue and non-disjoint psRegion. Should
+ // remove all peaks.
+ //
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmCullPeaks with large maxValue and non-disjoint psRegion.\n");
+ outData = pmFindImagePeaks(imgData, 0.0);
+ printf("pmFindImagePeaks found %d peaks\n", outData->size);
+ tmpRegion = psRegionAlloc(-PS_MAX_F32, PS_MAX_F32, -PS_MAX_F32, PS_MAX_F32);
+ outData = pmCullPeaks(outData, PS_MAX_F32, tmpRegion);
+
+ if (outData == NULL) {
+ printf("TEST ERROR: pmCullPeaks() returned a non-NULL psList.\n");
+ testStatus = false;
+ return(testStatus);
+ }
+ if (outData->size != 0) {
+ printf("TEST ERROR (3): pmCullPeaks incorrectly removed peaks\n");
+ printf("The pmCullPeaks() output had %d peaks, should have had %d peaks.\n", outData->size, 0);
+ testStatus = false;
+ }
+ psFree(outData);
+ psFree(tmpRegion);
+
+ printf("----------------------------------------------------------------------------------\n");
+ psFree(imgData);
+ return(testStatus);
+}
+
+#define TST04_NUM_ROWS 100
+#define TST04_NUM_COLS 100
+#define TST04_SKY 20.0
+#define TST04_INNER_RADIUS 3
+#define TST04_OUTER_RADIUS 5
+/******************************************************************************
+test04(): We first test pmSourceLocalSky() with various NULL and unallowable
+input parameters.
+
+XXX: Should we produce tests with boundary numbers for the inner/outer radius?
+
+XXX: Call this with varying sizes for the image.
+ *****************************************************************************/
+int test04( void )
+{
+ bool testStatus = true;
+ psImage *imgData = psImageAlloc(TST04_NUM_COLS, TST04_NUM_ROWS, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < imgData->numRows; i++) {
+ for (psS32 j = 0 ; j < imgData->numCols; j++) {
+ imgData->data.F32[i][j] = TST04_SKY;
+ }
+ }
+ psImage *imgDataF64 = psImageAlloc(TST04_NUM_COLS, TST04_NUM_ROWS, PS_TYPE_F64);
+ for (psS32 i = 0 ; i < imgDataF64->numRows; i++) {
+ for (psS32 j = 0 ; j < imgDataF64->numCols; j++) {
+ imgDataF64->data.F64[i][j] = 0.0;
+ }
+ }
+ psSource *rc = NULL;
+ psPeak *tmpPeak = pmPeakAlloc((psF32) (TST04_NUM_ROWS / 2),
+ (psF32) (TST04_NUM_COLS / 2),
+ 200.0,
+ PM_PEAK_LONE);
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with NULL psImage. Should generate error and return NULL.\n");
+ rc = pmSourceLocalSky(NULL, tmpPeak, PS_STAT_SAMPLE_MEAN, 10.0, 20.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a non-NULL psSource.\n");
+ psFree(rc);
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with wrong-type psImage. Should generate error and return NULL.\n");
+ rc = pmSourceLocalSky(imgDataF64, tmpPeak, PS_STAT_SAMPLE_MEAN, 10.0, 20.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a non-NULL psSource.\n");
+ psFree(rc);
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with NULL psPeak. Should generate error and return NULL.\n");
+ rc = pmSourceLocalSky(imgData, NULL, PS_STAT_SAMPLE_MEAN, 10.0, 20.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a non-NULL psSource.\n");
+ psFree(rc);
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with innerRadius<0.0. Should generate error and return NULL.\n");
+ rc = pmSourceLocalSky(imgData, tmpPeak, PS_STAT_SAMPLE_MEAN, -10.0, 20.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a non-NULL psSource.\n");
+ psFree(rc);
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with innerRadius>outerRadius. Should generate error and return NULL.\n");
+ rc = pmSourceLocalSky(imgData, tmpPeak, PS_STAT_SAMPLE_MEAN, 10.0, 5.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a non-NULL psSource.\n");
+ psFree(rc);
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with subImage startRow < 0. Should generate error and return NULL.\n");
+ tmpPeak->x = (psF32) (TST04_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST04_NUM_COLS / 2);
+ tmpPeak->y = 1;
+ rc = pmSourceLocalSky(imgData, tmpPeak, PS_STAT_SAMPLE_MEAN, 10.0, 20.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a non-NULL psSource.\n");
+ psFree(rc);
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with subImage endRow > numRows. Should generate error and return NULL.\n");
+ tmpPeak->x = (psF32) (TST04_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST04_NUM_COLS / 2);
+ tmpPeak->y = TST04_NUM_ROWS;
+ rc = pmSourceLocalSky(imgData, tmpPeak, PS_STAT_SAMPLE_MEAN, 10.0, 20.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a non-NULL psSource.\n");
+ psFree(rc);
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with subImage startCol < 0. Should generate error and return NULL.\n");
+ tmpPeak->x = (psF32) (TST04_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST04_NUM_COLS / 2);
+ tmpPeak->x = 1;
+ rc = pmSourceLocalSky(imgData, tmpPeak, PS_STAT_SAMPLE_MEAN, 10.0, 20.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a non-NULL psSource.\n");
+ psFree(rc);
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with subImage endCol > numCols. Should generate error and return NULL.\n");
+ tmpPeak->x = (psF32) (TST04_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST04_NUM_COLS / 2);
+ tmpPeak->x = TST04_NUM_COLS;
+ rc = pmSourceLocalSky(imgData, tmpPeak, PS_STAT_SAMPLE_MEAN, (psF32) TST04_INNER_RADIUS, (psF32) TST04_OUTER_RADIUS);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a non-NULL psSource.\n");
+ psFree(rc);
+ testStatus = false;
+ }
+
+
+ //
+ // XXX: The following code should be a separate function, and we should call it
+ // with a variety of image sizes, peaks.
+ //
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceLocalSky with valid data.\n");
+ tmpPeak->x = (psF32) (TST04_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST04_NUM_COLS / 2);
+ rc = pmSourceLocalSky(imgData,
+ tmpPeak,
+ PS_STAT_SAMPLE_MEAN,
+ (psF32) TST04_INNER_RADIUS,
+ (psF32) TST04_OUTER_RADIUS);
+
+ if (rc == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource.\n");
+ testStatus = false;
+ } else {
+ if (rc->peak == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource->peak.\n");
+ testStatus = false;
+ } else {
+ if (rc->peak->x != tmpPeak->x) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->peak->x was %d, should have been %d.\n", rc->peak->x, tmpPeak->x);
+ testStatus = false;
+ }
+
+ if (rc->peak->y != tmpPeak->y) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->peak->y was %d, should have been %d.\n", rc->peak->y, tmpPeak->y);
+ testStatus = false;
+ }
+
+ if (rc->peak->counts != tmpPeak->counts) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->peak->counts was %f, should have been %f.\n", rc->peak->counts, tmpPeak->counts);
+ testStatus = false;
+ }
+
+ if (rc->peak->class != tmpPeak->class) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->peak->class was %d, should have been %d.\n", rc->peak->class, tmpPeak->class);
+ testStatus = false;
+ }
+ }
+
+ if (rc->pixels == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->pixels was NULL.\n");
+ testStatus = false;
+ } else {
+ if (rc->pixels->numRows != (2 * TST04_OUTER_RADIUS)) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->pixels->numRows was %d, should have been %d.\n",
+ rc->pixels->numRows, (2 * TST04_OUTER_RADIUS));
+ testStatus = false;
+ }
+
+ if (rc->pixels->numCols != (2 * TST04_OUTER_RADIUS)) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->pixels->numCols was %d, should have been %d.\n",
+ rc->pixels->numCols, (2 * TST04_OUTER_RADIUS));
+ testStatus = false;
+ }
+
+ if (rc->pixels->col0 != (tmpPeak->x - TST04_OUTER_RADIUS)) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->pixels->col0 was %d, should have been %d.\n",
+ rc->pixels->col0, (tmpPeak->x - TST04_OUTER_RADIUS));
+ testStatus = false;
+ }
+
+ if (rc->pixels->row0 != (tmpPeak->y - TST04_OUTER_RADIUS)) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->pixels->row0 was %d, should have been %d.\n",
+ rc->pixels->row0, (tmpPeak->y - TST04_OUTER_RADIUS));
+ testStatus = false;
+ }
+
+ if (rc->pixels->type.type != PS_TYPE_F32) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->pixels->type was %d, should have been %d.\n",
+ rc->pixels->type.type, PS_TYPE_F32);
+ testStatus = false;
+ }
+
+ // XXX: Test the rc->pixels-> row/col offsets.
+ // XXX: Test that the pixels corresponds to the source image pixels.
+ }
+
+ if (rc->mask == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->mask was NULL.\n");
+ testStatus = false;
+ } else {
+ if (rc->mask->numRows != (2 * TST04_OUTER_RADIUS)) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->mask->numRows was %d, should have been %d.\n",
+ rc->mask->numRows, (2 * TST04_OUTER_RADIUS));
+ testStatus = false;
+ }
+
+ if (rc->mask->numCols != (2 * TST04_OUTER_RADIUS)) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->mask->numCols was %d, should have been %d.\n",
+ rc->mask->numCols, (2 * TST04_OUTER_RADIUS));
+ testStatus = false;
+ }
+
+ if (rc->mask->type.type != PS_TYPE_U8) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->mask->type was %d, should have been %d.\n",
+ rc->mask->type.type, PS_TYPE_U8);
+ testStatus = false;
+ }
+
+ // XXX: Test the rc->mask-> row/col offsets.
+ // XXX: Test that the correct pixels were masked, not merely the number of masked pixels.
+ psS32 unmaskedPixels = 0;
+ psS32 maskedPixels = 0;
+
+ for (psS32 row = 0 ; row < rc->mask->numRows; row++) {
+ for (psS32 col = 0 ; col < rc->mask->numCols; col++) {
+ if (rc->mask->data.U8[row][col] == 0) {
+ unmaskedPixels++;
+ } else {
+ maskedPixels++;
+ }
+ }
+ }
+ if (maskedPixels != PS_SQR(2*(TST04_INNER_RADIUS-1))) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->mask had %d masked pixels, should have been %d.\n",
+ maskedPixels, PS_SQR(2*(TST04_INNER_RADIUS-1)));
+ testStatus = false;
+ }
+ if (unmaskedPixels != (PS_SQR(2*TST04_OUTER_RADIUS) - PS_SQR(2*(TST04_INNER_RADIUS-1)))) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->mask had %d masked pixels, should have been %d.\n",
+ unmaskedPixels, (PS_SQR(2*TST04_OUTER_RADIUS) - PS_SQR(2*(TST04_INNER_RADIUS-1))));
+ testStatus = false;
+ }
+ }
+
+ if (rc->moments == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource->moments.\n");
+ testStatus = false;
+ } else {
+ if (rc->moments->Sky != TST04_SKY) {
+ printf("TEST ERROR: pmSourceLocalSky() psSource->moments->Sky was %f, should have been %f.\n", rc->moments->Sky, TST04_SKY);
+ testStatus = false;
+ }
+ }
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ psFree(rc);
+ psFree(imgData);
+ psFree(imgDataF64);
+ return(testStatus);
+}
+
+#define TST06_NUM_ROWS 100
+#define TST06_NUM_COLS 100
+#define TST06_SKY 20.0
+#define TST06_INNER_RADIUS 3
+#define TST06_OUTER_RADIUS 5
+/******************************************************************************
+test06(): We first test pmSourceLocalSky() with various NULL and unallowable
+input parameters.
+
+XXX: Should we produce tests with boundary numbers for the inner/outer radius?
+
+XXX: Call this with varying sizes for the image.
+ *****************************************************************************/
+int test06( void )
+{
+ bool testStatus = true;
+ psSource *tmpSource = NULL;
+ bool rc = false;
+ // Create the image used in this test.
+ psImage *imgData = psImageAlloc(TST06_NUM_COLS, TST06_NUM_ROWS, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < imgData->numRows; i++) {
+ for (psS32 j = 0 ; j < imgData->numCols; j++) {
+ imgData->data.F32[i][j] = TST06_SKY;
+ }
+ }
+ psImage *imgDataF64 = psImageAlloc(TST06_NUM_COLS, TST06_NUM_ROWS, PS_TYPE_F64);
+ for (psS32 i = 0 ; i < imgDataF64->numRows; i++) {
+ for (psS32 j = 0 ; j < imgDataF64->numCols; j++) {
+ imgDataF64->data.F64[i][j] = 0.0;
+ }
+ }
+
+ //
+ // Create a psPeak with the center pixel set to the peak.
+ //
+ psPeak *tmpPeak = pmPeakAlloc((psF32) (TST06_NUM_ROWS / 2),
+ (psF32) (TST06_NUM_COLS / 2),
+ 200.0,
+ PM_PEAK_LONE);
+
+ tmpSource = pmSourceLocalSky(imgData,
+ tmpPeak,
+ PS_STAT_SAMPLE_MEAN,
+ (psF32) TST06_INNER_RADIUS,
+ (psF32) TST06_OUTER_RADIUS);
+ if (tmpSource == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource.\n");
+ psFree(imgData);
+ psFree(imgDataF64);
+ psFree(tmpSource);
+ testStatus = false;
+ return(testStatus);
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSetPixelCircle with NULL psSource. Should generate error and return NULL.\n");
+ rc = pmSourceSetPixelCircle(NULL, imgData, 10.0);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSetPixelCircle() returned TRUE.\n");
+ testStatus = false;
+ }
+ // XXX: test with psSource->peaks NULL
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSetPixelCircle with NULL psImage. Should generate error and return NULL.\n");
+ rc = pmSourceSetPixelCircle(tmpSource, NULL, 10.0);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSetPixelCircle() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSetPixelCircle with wrong type psImage. Should generate error and return NULL.\n");
+ rc = pmSourceSetPixelCircle(tmpSource, imgDataF64, 10.0);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSetPixelCircle() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSetPixelCircle with radius < 0.0. Should generate error and return NULL.\n");
+ rc = pmSourceSetPixelCircle(tmpSource, imgData, -10.0);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSetPixelCircle() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSetPixelCircle with subImage startCol < 0. Should generate error and return NULL.\n");
+ tmpSource->peak->x = (psF32) (TST06_NUM_ROWS / 2);
+ tmpSource->peak->y = (psF32) (TST06_NUM_COLS / 2);
+ tmpSource->peak->x = 1;
+ rc = pmSourceSetPixelCircle(tmpSource, imgData, 10.0);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSetPixelCircle() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSetPixelCircle with subImage endCol > numCols. Should generate error and return NULL.\n");
+ tmpSource->peak->x = (psF32) (TST06_NUM_ROWS / 2);
+ tmpSource->peak->y = (psF32) (TST06_NUM_COLS / 2);
+ tmpSource->peak->x = TST06_NUM_COLS;
+ rc = pmSourceSetPixelCircle(tmpSource, imgData, 10.0);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSetPixelCircle() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSetPixelCircle with subImage startRow < 0. Should generate error and return NULL.\n");
+ tmpSource->peak->x = (psF32) (TST06_NUM_ROWS / 2);
+ tmpSource->peak->y = (psF32) (TST06_NUM_COLS / 2);
+ tmpSource->peak->y = 1;
+ rc = pmSourceSetPixelCircle(tmpSource, imgData, 10.0);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSetPixelCircle() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSetPixelCircle with subImage endRow > numRows. Should generate error and return NULL.\n");
+ tmpSource->peak->x = (psF32) (TST06_NUM_ROWS / 2);
+ tmpSource->peak->y = (psF32) (TST06_NUM_COLS / 2);
+ tmpSource->peak->y = TST06_NUM_ROWS;
+ rc = pmSourceSetPixelCircle(tmpSource, imgData, 10.0);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSetPixelCircle() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSetPixelCircle with valid data.\n");
+ tmpSource->peak->x = (psF32) (TST06_NUM_ROWS / 2);
+ tmpSource->peak->y = (psF32) (TST06_NUM_COLS / 2);
+ rc = pmSourceSetPixelCircle(tmpSource, imgData, 10.0);
+
+ if (rc == false) {
+ printf("TEST ERROR: pmSourceSetPixelCircle() returned FALSE.\n");
+ testStatus = false;
+ } else {
+ // XXX: Test correctness of the various psSetPixelCircle members.
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ psFree(imgData);
+ psFree(imgDataF64);
+ psFree(tmpSource);
+ return(testStatus);
+
+}
+
+#define TST05_NUM_ROWS 100
+#define TST05_NUM_COLS 100
+#define TST05_SKY 20.0
+#define TST05_INNER_RADIUS 3
+#define TST05_OUTER_RADIUS 5
+/******************************************************************************
+test05(): We first test pmSourceMoments() with various NULL and unallowable
+input parameters.
+
+XXX: Should we produce tests with boundary numbers for the inner/outer radius?
+
+XXX: Call this with varying sizes for the image.
+
+XXX: The actual values of the moments are not tested.
+ *****************************************************************************/
+int test05( void )
+{
+ bool testStatus = true;
+ psSource *tmpSource = NULL;
+ psSource *rc = NULL;
+ // Create the image used in this test.
+ psImage *imgData = psImageAlloc(TST05_NUM_COLS, TST05_NUM_ROWS, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < imgData->numRows; i++) {
+ for (psS32 j = 0 ; j < imgData->numCols; j++) {
+ imgData->data.F32[i][j] = TST05_SKY;
+ }
+ }
+
+ //
+ // Create a psPeak with the center pixel set to the peak.
+ //
+ psPeak *tmpPeak = pmPeakAlloc((psF32) (TST05_NUM_ROWS / 2),
+ (psF32) (TST05_NUM_COLS / 2),
+ 200.0,
+ PM_PEAK_LONE);
+
+ tmpSource = pmSourceLocalSky(imgData,
+ tmpPeak,
+ PS_STAT_SAMPLE_MEAN,
+ (psF32) TST05_INNER_RADIUS,
+ (psF32) TST05_OUTER_RADIUS);
+ if (tmpSource == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource.\n");
+ psFree(imgData);
+ psFree(tmpSource);
+ testStatus = false;
+ return(testStatus);
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceMoments with NULL psSource. Should generate error and return NULL.\n");
+ rc = pmSourceMoments(NULL, 10.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceMoments() returned a non-NULL psSource.\n");
+ testStatus = false;
+ }
+ // XXX: test with psSource->peaks NULL
+ // XXX: test with psSource->pixels NULL
+ // XXX: test with psSource->mask NULL
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceMoments with radius < 0.0. Should generate error and return NULL.\n");
+ rc = pmSourceMoments(tmpSource, -10.0);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceMoments() returned a non-NULL psSource.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ psFree(imgData);
+ psFree(tmpSource);
+ return(testStatus);
+
+}
+
+/******************************************************************************
+test07(): We first test the various psMinLM_... routines with various NULL and
+unallowable input parameters.
+
+XXX: We don't verify the numbers. Must do this.
+ *****************************************************************************/
+int test07( void )
+{
+ bool testStatus = true;
+ psF32 rc;
+ psVector *deriv = psVectorAlloc(7, PS_TYPE_F32);
+ psVector *params = psVectorAlloc(7, PS_TYPE_F32);
+ psVector *x = psVectorAlloc(2, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < 7 ; i++) {
+ deriv->data.F32[i] = 1.0;
+ params->data.F32[i] = 1.0;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_Gauss2D with NULL deriv vector. Should not generate error.\n");
+ rc = pmMinLM_Gauss2D(NULL, params, x);
+ if (isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_Gauss2D() returned a NAN psF32.\n");
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_Gauss2D with NULL params vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_Gauss2D(deriv, NULL, x);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_Gauss2D() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_Gauss2D with NULL x vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_Gauss2D(deriv, params, NULL);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_Gauss2D() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ psFree(deriv);
+ psFree(params);
+
+ deriv = psVectorAlloc(7, PS_TYPE_F32);
+ params = psVectorAlloc(7, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < 7 ; i++) {
+ deriv->data.F32[i] = 1.0;
+ params->data.F32[i] = 1.0;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_PsuedoGauss2D with NULL deriv vector. Should not generate error.\n");
+ rc = pmMinLM_PsuedoGauss2D(NULL, params, x);
+ if (isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_PsuedoGauss2D() returned a NAN psF32.\n");
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_PsuedoGauss2D with NULL params vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_PsuedoGauss2D(deriv, NULL, x);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_PsuedoGauss2D() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_PsuedoGauss2D with NULL x vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_PsuedoGauss2D(deriv, params, NULL);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_PsuedoGauss2D() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ psFree(deriv);
+ psFree(params);
+
+
+ deriv = psVectorAlloc(9, PS_TYPE_F32);
+ params = psVectorAlloc(9, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < 9 ; i++) {
+ deriv->data.F32[i] = 1.0;
+ params->data.F32[i] = 1.0;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_Wauss2D with NULL deriv vector. Should not generate error.\n");
+ rc = pmMinLM_Wauss2D(NULL, params, x);
+ if (isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_Wauss2D() returned a NAN psF32.\n");
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_Wauss2D with NULL params vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_Wauss2D(deriv, NULL, x);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_Wauss2D() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_Wauss2D with NULL x vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_Wauss2D(deriv, params, NULL);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_Wauss2D() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ psFree(deriv);
+ psFree(params);
+
+ deriv = psVectorAlloc(11, PS_TYPE_F32);
+ params = psVectorAlloc(11, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < 11 ; i++) {
+ deriv->data.F32[i] = 1.0;
+ params->data.F32[i] = 1.0;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_TwistGauss2D with NULL deriv vector. Should not generate error.\n");
+ rc = pmMinLM_TwistGauss2D(NULL, params, x);
+ if (isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_TwistGauss2D() returned a NAN psF32.\n");
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_TwistGauss2D with NULL params vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_TwistGauss2D(deriv, NULL, x);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_TwistGauss2D() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_TwistGauss2D with NULL x vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_TwistGauss2D(deriv, params, NULL);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_TwistGauss2D() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ psFree(deriv);
+ psFree(params);
+
+ deriv = psVectorAlloc(8, PS_TYPE_F32);
+ params = psVectorAlloc(8, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < 8 ; i++) {
+ deriv->data.F32[i] = 1.0;
+ params->data.F32[i] = 1.0;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_Sersic with NULL deriv vector. Should not generate error.\n");
+ rc = pmMinLM_Sersic(NULL, params, x);
+ if (isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_Sersic() returned a NAN psF32.\n");
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_Sersic with NULL params vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_Sersic(deriv, NULL, x);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_Sersic() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_Sersic with NULL x vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_Sersic(deriv, params, NULL);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_Sersic() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ psFree(deriv);
+ psFree(params);
+
+ deriv = psVectorAlloc(12, PS_TYPE_F32);
+ params = psVectorAlloc(12, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < 12 ; i++) {
+ deriv->data.F32[i] = 1.0;
+ params->data.F32[i] = 1.0;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_SersicCore with NULL deriv vector. Should not generate error.\n");
+ rc = pmMinLM_SersicCore(NULL, params, x);
+ if (isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_SersicCore() returned a NAN psF32.\n");
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_SersicCore with NULL params vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_SersicCore(deriv, NULL, x);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_SersicCore() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmMinLM_SersicCore with NULL x vector. Should generate error and return NAN.\n");
+ rc = pmMinLM_SersicCore(deriv, params, NULL);
+ if (!isnan(rc)) {
+ printf("TEST ERROR: pmMinLM_SersicCore() returned a non-NAN psF32 (%f).\n", rc);
+ testStatus = false;
+ }
+
+ psFree(deriv);
+ psFree(params);
+ psFree(x);
+ return(testStatus);
+}
+
+/******************************************************************************
+test08(): We first test pmSourceModelGuess() with various NULL and unallowable
+input parameters.
+
+XXX: We don't verify the numbers.
+ *****************************************************************************/
+int test08( void )
+{
+ bool testStatus = true;
+ psImage *imgData = psImageAlloc(TST04_NUM_COLS, TST04_NUM_ROWS, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < imgData->numRows; i++) {
+ for (psS32 j = 0 ; j < imgData->numCols; j++) {
+ imgData->data.F32[i][j] = TST04_SKY;
+ }
+ }
+ psSource *mySrc = NULL;
+ bool rc = false;
+ psPeak *tmpPeak = pmPeakAlloc((psF32) (TST04_NUM_ROWS / 2),
+ (psF32) (TST04_NUM_COLS / 2),
+ 200.0,
+ PM_PEAK_LONE);
+
+ printf("Calling pmSourceLocalSky with valid data.\n");
+ tmpPeak->x = (psF32) (TST04_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST04_NUM_COLS / 2);
+ mySrc = pmSourceLocalSky(imgData,
+ tmpPeak,
+ PS_STAT_SAMPLE_MEAN,
+ (psF32) TST04_INNER_RADIUS,
+ (psF32) TST04_OUTER_RADIUS);
+
+ if (mySrc == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceModelGuess with NULL psSource. Should generate error, return FALSE.\n");
+ rc = pmSourceModelGuess(NULL, imgData, PS_MODEL_GAUSS);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceModelGuess with NULL psImage. Should generate error, return FALSE.\n");
+ rc = pmSourceModelGuess(mySrc, NULL, PS_MODEL_GAUSS);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceModelGuess with bad model type. Should generate error, return FALSE.\n");
+ rc = pmSourceModelGuess(mySrc, imgData, PS_MODEL_UNDEFINED);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceModelGuess with PS_MODEL_GAUSS\n");
+ rc = pmSourceModelGuess(mySrc, imgData, PS_MODEL_GAUSS);
+ if (rc != true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned FALSE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceModelGuess with PS_MODEL_PGAUSS\n");
+ rc = pmSourceModelGuess(mySrc, imgData, PS_MODEL_PGAUSS);
+ if (rc != true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned FALSE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceModelGuess with PS_MODEL_TWIST_GAUSS\n");
+ rc = pmSourceModelGuess(mySrc, imgData, PS_MODEL_TWIST_GAUSS);
+ if (rc != true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned FALSE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceModelGuess with PS_MODEL_WAUSS\n");
+ rc = pmSourceModelGuess(mySrc, imgData, PS_MODEL_WAUSS);
+ if (rc != true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned FALSE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceModelGuess with PS_MODEL_SERSIC\n");
+ rc = pmSourceModelGuess(mySrc, imgData, PS_MODEL_SERSIC);
+ if (rc != true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned FALSE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceModelGuess with PS_MODEL_SERSIC_CORE\n");
+ rc = pmSourceModelGuess(mySrc, imgData, PS_MODEL_SERSIC_CORE);
+ if (rc != true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned FALSE.\n");
+ testStatus = false;
+ }
+
+ psFree(mySrc);
+ // psFree(tmpPeak);
+ psFree(imgData);
+ return(testStatus);
+}
+
+
+#define TST09_NUM_ROWS 70
+#define TST09_NUM_COLS 70
+#define TST09_SKY 5.0
+#define TST09_INNER_RADIUS 3
+#define TST09_OUTER_RADIUS 10
+#define LEVEL (TST09_SKY + 10.0)
+/******************************************************************************
+test09(): We first test pmSourceContour() with various NULL and unallowable
+input parameters.
+
+XXX: We don't verify the numbers.
+ *****************************************************************************/
+int test09( void )
+{
+ bool testStatus = true;
+ psImage *imgData = psImageAlloc(TST09_NUM_COLS, TST09_NUM_ROWS, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < imgData->numRows; i++) {
+ for (psS32 j = 0 ; j < imgData->numCols; j++) {
+ imgData->data.F32[i][j] = TST09_SKY;
+ }
+ }
+ psSource *mySrc = NULL;
+ psArray *rc = NULL;
+
+ psPeak *tmpPeak = pmPeakAlloc((psF32) (TST09_NUM_ROWS / 2),
+ (psF32) (TST09_NUM_COLS / 2),
+ 200.0,
+ PM_PEAK_LONE);
+
+ printf("Calling pmSourceLocalSky with valid data.\n");
+ tmpPeak->x = (psF32) (TST09_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST09_NUM_COLS / 2);
+ mySrc = pmSourceLocalSky(imgData,
+ tmpPeak,
+ PS_STAT_SAMPLE_MEAN,
+ (psF32) TST09_INNER_RADIUS,
+ (psF32) TST09_OUTER_RADIUS);
+
+ if (mySrc == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource.\n");
+ testStatus = false;
+ }
+
+ bool rcBool = pmSourceModelGuess(mySrc, imgData, PS_MODEL_GAUSS);
+ if (rcBool != true) {
+ printf("TEST ERROR: pmSourceModelGuess() returned FALSE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceContour with NULL psSource . Should generate error, return NULL.\n");
+ rc = pmSourceContour(NULL, imgData, LEVEL, PS_CONTOUR_CRUDE);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceContour() returned non-NULL.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceContour with NULL psImage . Should generate error, return NULL.\n");
+ rc = pmSourceContour(mySrc, NULL, LEVEL, PS_CONTOUR_CRUDE);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceContour() returned non-NULL.\n");
+ testStatus = false;
+ }
+
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceContour with acceptable data.\n");
+ printf("NOTE: must figure out the parameters for this test to be meaningful.\n");
+ mySrc->models->params[0] = TST09_SKY;
+ mySrc->models->params[1] = 15.0;
+ mySrc->models->params[2] = (psF32) (TST09_NUM_ROWS / 2);
+ mySrc->models->params[3] = (psF32) (TST09_NUM_COLS / 2);
+ mySrc->models->params[4] = 2.0;
+ mySrc->models->params[5] = 2.0;
+ mySrc->models->params[6] = 2.0;
+ rc = pmSourceContour(mySrc, imgData, LEVEL, PS_CONTOUR_CRUDE);
+ if (rc != NULL) {
+ printf("TEST ERROR: pmSourceContour() returned NULL.\n");
+ testStatus = false;
+ }
+
+ psFree(rc);
+ psFree(mySrc);
+ // psFree(tmpPeak);
+ psFree(imgData);
+ return(testStatus);
+}
+
+#define TST15_NUM_ROWS 100
+#define TST15_NUM_COLS 100
+#define TST15_SKY 10.0
+#define TST15_INNER_RADIUS 3
+#define TST15_OUTER_RADIUS 5
+/******************************************************************************
+test15(): We first test pmSourceAddModel() with various NULL and unallowable
+input parameters.
+
+XXX: We don't verify the numbers.
+ *****************************************************************************/
+int test15( void )
+{
+ bool testStatus = true;
+ psImage *imgData = psImageAlloc(TST15_NUM_COLS, TST15_NUM_ROWS, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < imgData->numRows; i++) {
+ for (psS32 j = 0 ; j < imgData->numCols; j++) {
+ imgData->data.F32[i][j] = TST15_SKY;
+ }
+ }
+ psSource *mySrc = NULL;
+ psBool rc = false;
+
+ psPeak *tmpPeak = pmPeakAlloc((psF32) (TST15_NUM_ROWS / 2),
+ (psF32) (TST15_NUM_COLS / 2),
+ 200.0,
+ PM_PEAK_LONE);
+
+ printf("Calling pmSourceLocalSky with valid data.\n");
+ tmpPeak->x = (psF32) (TST15_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST15_NUM_COLS / 2);
+ mySrc = pmSourceLocalSky(imgData,
+ tmpPeak,
+ PS_STAT_SAMPLE_MEAN,
+ (psF32) TST15_INNER_RADIUS,
+ (psF32) TST15_OUTER_RADIUS);
+
+ if (mySrc == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource.\n");
+ testStatus = false;
+ }
+
+ mySrc->models = pmModelAlloc(PS_MODEL_GAUSS);
+ mySrc->models->params[0] = 5.0;
+ mySrc->models->params[1] = 70.0;
+ mySrc->models->params[2] = (psF32) (TST15_NUM_ROWS / 2);
+ mySrc->models->params[3] = (psF32) (TST15_NUM_COLS / 2);
+ mySrc->models->params[4] = 1.0;
+ mySrc->models->params[5] = 1.0;
+ mySrc->models->params[6] = 2.0;
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceAddModel with NULL psImage. Should generate error, return FALSE.\n");
+ rc = pmSourceAddModel(NULL, mySrc, true);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceAddModel() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceAddModel with NULL psSrc. Should generate error, return FALSE.\n");
+ rc = pmSourceAddModel(imgData, NULL, true);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceAddModel() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceAddModel with acceptable data.\n");
+ rc = pmSourceAddModel(imgData, mySrc, true);
+ if (rc != true) {
+ printf("TEST ERROR: pmSourceAddModel() returned FALSE.\n");
+ testStatus = false;
+ }
+
+ psFree(mySrc);
+ psFree(imgData);
+ return(testStatus);
+}
+
+#define TST16_NUM_ROWS 100
+#define TST16_NUM_COLS 100
+#define TST16_SKY 10.0
+#define TST16_INNER_RADIUS 3
+#define TST16_OUTER_RADIUS 5
+/******************************************************************************
+test16(): We first test pmSourceSubModel() with various NULL and unallowable
+input parameters.
+
+XXX: We don't verify the numbers.
+ *****************************************************************************/
+int test16( void )
+{
+ bool testStatus = true;
+ psImage *imgData = psImageAlloc(TST16_NUM_COLS, TST16_NUM_ROWS, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < imgData->numRows; i++) {
+ for (psS32 j = 0 ; j < imgData->numCols; j++) {
+ imgData->data.F32[i][j] = TST16_SKY;
+ }
+ }
+ psSource *mySrc = NULL;
+ psBool rc = false;
+
+ psPeak *tmpPeak = pmPeakAlloc((psF32) (TST16_NUM_ROWS / 2),
+ (psF32) (TST16_NUM_COLS / 2),
+ 200.0,
+ PM_PEAK_LONE);
+
+ printf("Calling pmSourceLocalSky with valid data.\n");
+ tmpPeak->x = (psF32) (TST16_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST16_NUM_COLS / 2);
+ mySrc = pmSourceLocalSky(imgData,
+ tmpPeak,
+ PS_STAT_SAMPLE_MEAN,
+ (psF32) TST16_INNER_RADIUS,
+ (psF32) TST16_OUTER_RADIUS);
+
+ if (mySrc == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource.\n");
+ testStatus = false;
+ }
+
+ mySrc->models = pmModelAlloc(PS_MODEL_GAUSS);
+ mySrc->models->params[0] = 5.0;
+ mySrc->models->params[1] = 70.0;
+ mySrc->models->params[2] = (psF32) (TST16_NUM_ROWS / 2);
+ mySrc->models->params[3] = (psF32) (TST16_NUM_COLS / 2);
+ mySrc->models->params[4] = 1.0;
+ mySrc->models->params[5] = 1.0;
+ mySrc->models->params[6] = 2.0;
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSubModel with NULL psImage. Should generate error, return FALSE.\n");
+ rc = pmSourceSubModel(NULL, mySrc, true);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSubModel() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSubModel with NULL psSrc. Should generate error, return FALSE.\n");
+ rc = pmSourceSubModel(imgData, NULL, true);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceSubModel() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceSubModel with acceptable data.\n");
+ rc = pmSourceSubModel(imgData, mySrc, true);
+ if (rc != true) {
+ printf("TEST ERROR: pmSourceSubModel() returned FALSE.\n");
+ testStatus = false;
+ }
+
+ psFree(mySrc);
+ psFree(imgData);
+ return(testStatus);
+}
+
+#define TST20_NUM_ROWS 100
+#define TST20_NUM_COLS 100
+#define TST20_SKY 10.0
+#define TST20_INNER_RADIUS 3
+#define TST20_OUTER_RADIUS 5
+/******************************************************************************
+test20(): We first test pmSourceSubModel() with various NULL and unallowable
+input parameters.
+
+XXX: We don't verify the numbers.
+ *****************************************************************************/
+int test20( void )
+{
+ bool testStatus = true;
+ psImage *imgData = psImageAlloc(TST20_NUM_COLS, TST20_NUM_ROWS, PS_TYPE_F32);
+ for (psS32 i = 0 ; i < imgData->numRows; i++) {
+ for (psS32 j = 0 ; j < imgData->numCols; j++) {
+ imgData->data.F32[i][j] = TST20_SKY;
+ }
+ }
+ psSource *mySrc = NULL;
+ psBool rc = false;
+
+ psPeak *tmpPeak = pmPeakAlloc((psF32) (TST20_NUM_ROWS / 2),
+ (psF32) (TST20_NUM_COLS / 2),
+ 200.0,
+ PM_PEAK_LONE);
+
+ printf("Calling pmSourceLocalSky with valid data.\n");
+ tmpPeak->x = (psF32) (TST20_NUM_ROWS / 2);
+ tmpPeak->y = (psF32) (TST20_NUM_COLS / 2);
+ mySrc = pmSourceLocalSky(imgData,
+ tmpPeak,
+ PS_STAT_SAMPLE_MEAN,
+ (psF32) TST20_INNER_RADIUS,
+ (psF32) TST20_OUTER_RADIUS);
+
+ if (mySrc == NULL) {
+ printf("TEST ERROR: pmSourceLocalSky() returned a NULL psSource.\n");
+ testStatus = false;
+ }
+
+ mySrc->models = pmModelAlloc(PS_MODEL_GAUSS);
+
+
+ mySrc->models->params[0] = 5.0;
+ mySrc->models->params[1] = 70.0;
+ mySrc->models->params[2] = (psF32) (TST20_NUM_ROWS / 2);
+ mySrc->models->params[3] = (psF32) (TST20_NUM_COLS / 2);
+ mySrc->models->params[4] = 1.0;
+ mySrc->models->params[5] = 1.0;
+ mySrc->models->params[6] = 2.0;
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceFitModel with NULL psImage. Should generate error, return FALSE.\n");
+ rc = pmSourceFitModel(mySrc, NULL);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceFitModel() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceFitModel with NULL psSource. Should generate error, return FALSE.\n");
+ rc = pmSourceFitModel(NULL, imgData);
+ if (rc == true) {
+ printf("TEST ERROR: pmSourceFitModel() returned TRUE.\n");
+ testStatus = false;
+ }
+
+ printf("----------------------------------------------------------------------------------\n");
+ printf("Calling pmSourceFitModel with acceptable data.\n");
+ rc = pmSourceFitModel(mySrc, imgData);
+ printf("pmSourceFitModel returned %d\n", rc);
+
+ // XXX: Memory leaks are not being tested
+ psVector *junk = psVectorAlloc(10, PS_TYPE_F32);
+ junk->data.F32[0] = 0.0;
+
+ psFree(mySrc);
+ psFree(imgData);
+ return(testStatus);
+}
+
+
+// this code will
+
+
+
Index: /branches/rel5/psModules/test/verified/tst_pmObjects01.stderr
===================================================================
--- /branches/rel5/psModules/test/verified/tst_pmObjects01.stderr (revision 3634)
+++ /branches/rel5/psModules/test/verified/tst_pmObjects01.stderr (revision 3634)
@@ -0,0 +1,237 @@
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: structure allocators and deallocators} *
+* TestType: Positive *
+\**********************************************************************************/
+
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: structure allocators and deallocators} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: psFindVectorPeaks()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|pmFindVectorPeaks (FILE:LINENO)
+ Unallowable operation: psVector vector or its data is NULL.
+||E|pmFindVectorPeaks (FILE:LINENO)
+ Unallowable operation: psVector vector has no elements.
+||E|pmFindVectorPeaks (FILE:LINENO)
+ Unallowable operation: psVector vector has incorrect type.
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: psFindVectorPeaks()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: psFindImagePeaks()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|psImageAlloc (FILE:LINENO)
+ Specified number of rows (0) or columns (0) is invalid.
+||E|pmFindImagePeaks (FILE:LINENO)
+ Unallowable operation: psImage image or its data is NULL.
+||E|pmFindImagePeaks (FILE:LINENO)
+ Unallowable operation: psImage image or its data is NULL.
+||E|pmFindImagePeaks (FILE:LINENO)
+ Unallowable operation: psImage image has incorrect type.
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: psFindImagePeaks()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmCullPeaks()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|pmCullPeaks (FILE:LINENO)
+ Unallowable operation: peaks is NULL.
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmCullPeaks()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmSourceLocalSky()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|pmSourceLocalSky (FILE:LINENO)
+ Unallowable operation: psImage image or its data is NULL.
+||E|pmSourceLocalSky (FILE:LINENO)
+ Unallowable operation: psImage image has incorrect type.
+||E|pmSourceLocalSky (FILE:LINENO)
+ Unallowable operation: peak is NULL.
+||E|pmSourceLocalSky (FILE:LINENO)
+ Error: (0.0 > innerRadius) (0.000000 -10.000000)
+||E|pmSourceLocalSky (FILE:LINENO)
+ Error: (innerRadius > outerRadius) (10.000000 5.000000)
+||E|pmSourceLocalSky (FILE:LINENO)
+ Sub image startRow is outside image boundaries (-19).
+||E|pmSourceLocalSky (FILE:LINENO)
+ Sub image endRow is outside image boundaries (120).
+||E|pmSourceLocalSky (FILE:LINENO)
+ Sub image startCol is outside image boundaries (-19).
+||E|pmSourceLocalSky (FILE:LINENO)
+ Sub image endCol is outside image boundaries (105).
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmSourceLocalSky()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmSourceSetPixelCircle()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|pmSourceSetPixelCircle (FILE:LINENO)
+ Unallowable operation: source is NULL.
+||E|pmSourceSetPixelCircle (FILE:LINENO)
+ Unallowable operation: psImage image or its data is NULL.
+||E|pmSourceSetPixelCircle (FILE:LINENO)
+ Unallowable operation: psImage image has incorrect type.
+||E|pmSourceSetPixelCircle (FILE:LINENO)
+ Error: (0.0 > radius) (0.000000 -10.000000)
+||E|pmSourceSetPixelCircle (FILE:LINENO)
+ Sub image startCol is outside image boundaries (-9).
+||E|pmSourceSetPixelCircle (FILE:LINENO)
+ Sub image endCol is outside image boundaries (110).
+||E|pmSourceSetPixelCircle (FILE:LINENO)
+ Sub image startRow is outside image boundaries (-9).
+||E|pmSourceSetPixelCircle (FILE:LINENO)
+ Sub image endRow is outside image boundaries (110).
+||W|pmSourceSetPixelCircle
+ WARNING: pmSourceSetPixelCircle(): image->pixels not NULL. Freeing and reallocating.
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmSourceSetPixelCircle()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmSourceMoments()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|pmSourceMoments (FILE:LINENO)
+ Unallowable operation: source is NULL.
+||E|pmSourceMoments (FILE:LINENO)
+ Error: (0.0 > radius) (0.000000 -10.000000)
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmSourceMoments()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmMin()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|pmMinLM_Gauss2D (FILE:LINENO)
+ Unallowable operation: psVector params or its data is NULL.
+||E|pmMinLM_Gauss2D (FILE:LINENO)
+ Unallowable operation: psVector x or its data is NULL.
+||E|pmMinLM_PsuedoGauss2D (FILE:LINENO)
+ Unallowable operation: psVector params or its data is NULL.
+||E|pmMinLM_PsuedoGauss2D (FILE:LINENO)
+ Unallowable operation: psVector x or its data is NULL.
+||E|pmMinLM_Wauss2D (FILE:LINENO)
+ Unallowable operation: psVector params or its data is NULL.
+||E|pmMinLM_Wauss2D (FILE:LINENO)
+ Unallowable operation: psVector x or its data is NULL.
+||E|pmMinLM_TwistGauss2D (FILE:LINENO)
+ Unallowable operation: psVector params or its data is NULL.
+||E|pmMinLM_TwistGauss2D (FILE:LINENO)
+ Unallowable operation: psVector x or its data is NULL.
+||E|pmMinLM_Sersic (FILE:LINENO)
+ This function is not implemented yet.
+||E|pmMinLM_Sersic (FILE:LINENO)
+ Unallowable operation: psVector params or its data is NULL.
+||E|pmMinLM_Sersic (FILE:LINENO)
+ Unallowable operation: psVector x or its data is NULL.
+||E|pmMinLM_SersicCore (FILE:LINENO)
+ This function is not implemented yet.
+||E|pmMinLM_SersicCore (FILE:LINENO)
+ Unallowable operation: psVector params or its data is NULL.
+||E|pmMinLM_SersicCore (FILE:LINENO)
+ Unallowable operation: psVector x or its data is NULL.
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmMin()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmSourceModelGuess()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|pmSourceModelGuess (FILE:LINENO)
+ Unallowable operation: source is NULL.
+||E|pmSourceModelGuess (FILE:LINENO)
+ Unallowable operation: psImage image or its data is NULL.
+||E|pmModelAlloc (FILE:LINENO)
+ Undefined psModelType
+||E|pmSourceModelGuess (FILE:LINENO)
+ Undefined psModelType
+||W|pmSourceModelGuess
+ WARNING: source->models was non-NULL; calling psFree(source->models).
+||W|pmSourceModelGuess
+ WARNING: source->models was non-NULL; calling psFree(source->models).
+||W|pmSourceModelGuess
+ WARNING: source->models was non-NULL; calling psFree(source->models).
+||W|pmSourceModelGuess
+ WARNING: source->models was non-NULL; calling psFree(source->models).
+||W|pmSourceModelGuess
+ WARNING: source->models was non-NULL; calling psFree(source->models).
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmSourceModelGuess()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmSourceContour()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|pmSourceContour (FILE:LINENO)
+ Unallowable operation: source is NULL.
+||E|pmSourceContour (FILE:LINENO)
+ Unallowable operation: image is NULL.
+||E|pmSourceContour (FILE:LINENO)
+ Could not find contour edge (NAN)
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmSourceContour()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmSourceAddModel()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|p_pmSourceAddOrSubModel (FILE:LINENO)
+ Unallowable operation: psImage image or its data is NULL.
+||E|p_pmSourceAddOrSubModel (FILE:LINENO)
+ Unallowable operation: src is NULL.
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmSourceAddModel()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmSourceSubModel()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|p_pmSourceAddOrSubModel (FILE:LINENO)
+ Unallowable operation: psImage image or its data is NULL.
+||E|p_pmSourceAddOrSubModel (FILE:LINENO)
+ Unallowable operation: src is NULL.
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmSourceSubModel()} | tst_pmObjects01.c)
+
+/***************************** TESTPOINT ******************************************\
+* TestFile: tst_pmObjects01.c *
+* TestPoint: Test Point Driver{pmObjects: pmSourceSubModel()} *
+* TestType: Positive *
+\**********************************************************************************/
+
+||E|pmSourceFitModel (FILE:LINENO)
+ Unallowable operation: psImage image or its data is NULL.
+||E|pmSourceFitModel (FILE:LINENO)
+ Unallowable operation: source is NULL.
+
+---> TESTPOINT PASSED (Test Point Driver{pmObjects: pmSourceSubModel()} | tst_pmObjects01.c)
+
Index: /branches/rel5/psModules/test/verified/tst_pmObjects01.stdout
===================================================================
--- /branches/rel5/psModules/test/verified/tst_pmObjects01.stdout (revision 3634)
+++ /branches/rel5/psModules/test/verified/tst_pmObjects01.stdout (revision 3634)
@@ -0,0 +1,189 @@
+----------------------------------------------------------------------------------
+Calling pmFindVectorPeaks with NULL psVector. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmFindVectorPeaks with empty psVector. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmFindVectorPeaks with PS_TYPE_F64 psVector. Should generate error and return NULL.
+-------------- Calling test_pmFindVectorPeaks on an 1 size vector. --------------
+Test pmFindVectorPeaks() with a first-element peak.
+Test pmFindVectorPeaks() with a first-element peak, large threshold.
+-------------- Calling test_pmFindVectorPeaks on an 10 size vector. --------------
+Test pmFindVectorPeaks() with a first-element peak.
+Test pmFindVectorPeaks() with a first-element peak, large threshold.
+Test pmFindVectorPeaks() with a last-element peak.
+Test pmFindVectorPeaks() with a last-element peak, large threshold.
+Test pmFindVectorPeaks() with all even-numbered elements peak.
+Test pmFindVectorPeaks() with all even-numbered elements peak, large threshold.
+----------------------------------------------------------------------------------
+Calling pmFindImagePeaks with NULL psImage. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmFindImagePeaks with empty psImage. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmFindImagePeaks with PS_TYPE_F64 psImage. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+-------------- Calling test_pmFindVectorPeaks on an 10-by-10 image. --------------
+(200.0) (41.0) (34.0) (29.0) (26.0) (25.0) (26.0) (29.0) (34.0) (200.0)
+(41.0) (32.0) (25.0) (20.0) (17.0) (16.0) (17.0) (20.0) (25.0) (32.0)
+(34.0) (25.0) (18.0) (13.0) (10.0) (9.0) (10.0) (13.0) (18.0) (25.0)
+(29.0) (20.0) (13.0) (8.0) (5.0) (4.0) (5.0) (8.0) (13.0) (20.0)
+(26.0) (17.0) (10.0) (5.0) (2.0) (1.0) (2.0) (5.0) (10.0) (17.0)
+(25.0) (16.0) (9.0) (4.0) (1.0) (0.0) (1.0) (4.0) (9.0) (16.0)
+(26.0) (17.0) (10.0) (5.0) (2.0) (1.0) (2.0) (5.0) (10.0) (17.0)
+(29.0) (20.0) (13.0) (8.0) (5.0) (4.0) (5.0) (8.0) (13.0) (20.0)
+(34.0) (25.0) (18.0) (13.0) (10.0) (9.0) (10.0) (13.0) (18.0) (25.0)
+(200.0) (32.0) (25.0) (20.0) (17.0) (16.0) (17.0) (20.0) (25.0) (200.0)
+----------------------------------------------------------------------------------
+Calling pmCullPeaks with NULL psList. Should generate error and return NULL.
+(10.0) (0.0) (10.0) (0.0) (10.0) (0.0) (10.0) (0.0) (10.0) (0.0)
+(0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0)
+(12.0) (0.0) (12.0) (0.0) (12.0) (0.0) (12.0) (0.0) (12.0) (0.0)
+(0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0)
+(14.0) (0.0) (14.0) (0.0) (14.0) (0.0) (14.0) (0.0) (14.0) (0.0)
+(0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0)
+(16.0) (0.0) (16.0) (0.0) (16.0) (0.0) (16.0) (0.0) (16.0) (0.0)
+(0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0)
+(18.0) (0.0) (18.0) (0.0) (18.0) (0.0) (18.0) (0.0) (18.0) (0.0)
+(0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0) (0.0)
+Set 25 peaks
+----------------------------------------------------------------------------------
+Calling pmCullPeaks with large maxValue and NULL psRegion.
+----------------------------------------------------------------------------------
+Calling pmCullPeaks with tiny maxValue and NULL psRegion.
+pmFindImagePeaks found 25 peaks
+----------------------------------------------------------------------------------
+Calling pmCullPeaks with large maxValue and disjoint psRegion.
+pmFindImagePeaks found 25 peaks
+----------------------------------------------------------------------------------
+Calling pmCullPeaks with large maxValue and non-disjoint psRegion.
+pmFindImagePeaks found 25 peaks
+----------------------------------------------------------------------------------
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with NULL psImage. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with wrong-type psImage. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with NULL psPeak. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with innerRadius<0.0. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with innerRadius>outerRadius. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with subImage startRow < 0. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with subImage endRow > numRows. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with subImage startCol < 0. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with subImage endCol > numCols. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceLocalSky with valid data.
+----------------------------------------------------------------------------------
+----------------------------------------------------------------------------------
+Calling pmSourceSetPixelCircle with NULL psSource. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceSetPixelCircle with NULL psImage. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceSetPixelCircle with wrong type psImage. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceSetPixelCircle with radius < 0.0. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceSetPixelCircle with subImage startCol < 0. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceSetPixelCircle with subImage endCol > numCols. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceSetPixelCircle with subImage startRow < 0. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceSetPixelCircle with subImage endRow > numRows. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceSetPixelCircle with valid data.
+----------------------------------------------------------------------------------
+----------------------------------------------------------------------------------
+Calling pmSourceMoments with NULL psSource. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceMoments with radius < 0.0. Should generate error and return NULL.
+----------------------------------------------------------------------------------
+----------------------------------------------------------------------------------
+Calling pmMinLM_Gauss2D with NULL deriv vector. Should not generate error.
+----------------------------------------------------------------------------------
+Calling pmMinLM_Gauss2D with NULL params vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_Gauss2D with NULL x vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_PsuedoGauss2D with NULL deriv vector. Should not generate error.
+----------------------------------------------------------------------------------
+Calling pmMinLM_PsuedoGauss2D with NULL params vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_PsuedoGauss2D with NULL x vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_Wauss2D with NULL deriv vector. Should not generate error.
+----------------------------------------------------------------------------------
+Calling pmMinLM_Wauss2D with NULL params vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_Wauss2D with NULL x vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_TwistGauss2D with NULL deriv vector. Should not generate error.
+----------------------------------------------------------------------------------
+Calling pmMinLM_TwistGauss2D with NULL params vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_TwistGauss2D with NULL x vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_Sersic with NULL deriv vector. Should not generate error.
+----------------------------------------------------------------------------------
+Calling pmMinLM_Sersic with NULL params vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_Sersic with NULL x vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_SersicCore with NULL deriv vector. Should not generate error.
+----------------------------------------------------------------------------------
+Calling pmMinLM_SersicCore with NULL params vector. Should generate error and return NAN.
+----------------------------------------------------------------------------------
+Calling pmMinLM_SersicCore with NULL x vector. Should generate error and return NAN.
+Calling pmSourceLocalSky with valid data.
+----------------------------------------------------------------------------------
+Calling pmSourceModelGuess with NULL psSource. Should generate error, return FALSE.
+----------------------------------------------------------------------------------
+Calling pmSourceModelGuess with NULL psImage. Should generate error, return FALSE.
+----------------------------------------------------------------------------------
+Calling pmSourceModelGuess with bad model type. Should generate error, return FALSE.
+----------------------------------------------------------------------------------
+Calling pmSourceModelGuess with PS_MODEL_GAUSS
+----------------------------------------------------------------------------------
+Calling pmSourceModelGuess with PS_MODEL_PGAUSS
+----------------------------------------------------------------------------------
+Calling pmSourceModelGuess with PS_MODEL_TWIST_GAUSS
+----------------------------------------------------------------------------------
+Calling pmSourceModelGuess with PS_MODEL_WAUSS
+----------------------------------------------------------------------------------
+Calling pmSourceModelGuess with PS_MODEL_SERSIC
+----------------------------------------------------------------------------------
+Calling pmSourceModelGuess with PS_MODEL_SERSIC_CORE
+Calling pmSourceLocalSky with valid data.
+----------------------------------------------------------------------------------
+Calling pmSourceContour with NULL psSource . Should generate error, return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceContour with NULL psImage . Should generate error, return NULL.
+----------------------------------------------------------------------------------
+Calling pmSourceContour with acceptable data.
+NOTE: must figure out the parameters for this test to be meaningful.
+Calling pmSourceLocalSky with valid data.
+----------------------------------------------------------------------------------
+Calling pmSourceAddModel with NULL psImage. Should generate error, return FALSE.
+----------------------------------------------------------------------------------
+Calling pmSourceAddModel with NULL psSrc. Should generate error, return FALSE.
+----------------------------------------------------------------------------------
+Calling pmSourceAddModel with acceptable data.
+Calling pmSourceLocalSky with valid data.
+----------------------------------------------------------------------------------
+Calling pmSourceSubModel with NULL psImage. Should generate error, return FALSE.
+----------------------------------------------------------------------------------
+Calling pmSourceSubModel with NULL psSrc. Should generate error, return FALSE.
+----------------------------------------------------------------------------------
+Calling pmSourceSubModel with acceptable data.
+Calling pmSourceLocalSky with valid data.
+----------------------------------------------------------------------------------
+Calling pmSourceFitModel with NULL psImage. Should generate error, return FALSE.
+----------------------------------------------------------------------------------
+Calling pmSourceFitModel with NULL psSource. Should generate error, return FALSE.
+----------------------------------------------------------------------------------
+Calling pmSourceFitModel with acceptable data.
+pmSourceFitModel returned 0