Changeset 1022 for trunk/psLib/src/dataManip
- Timestamp:
- Jun 14, 2004, 9:40:15 AM (22 years ago)
- Location:
- trunk/psLib/src/dataManip
- Files:
-
- 3 edited
-
psFunctions.c (modified) (5 diffs)
-
psMinimize.c (modified) (2 diffs)
-
psStats.c (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/dataManip/psFunctions.c
r1020 r1022 7 7 * polynomials. It also contains a Gaussian functions. 8 8 * 9 * @version $Revision: 1.1 2$ $Name: not supported by cvs2svn $10 * @date $Date: 2004-06-14 19: 32:42$9 * @version $Revision: 1.13 $ $Name: not supported by cvs2svn $ 10 * @date $Date: 2004-06-14 19:40:14 $ 11 11 * 12 12 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 73 73 Gaussian distribution as well. 74 74 75 GUS: There is no way to seed the random generator.75 NOTE: There is no way to seed the random generator. 76 76 *****************************************************************************/ 77 77 psVector *psGaussianDev(float mean, … … 94 94 } 95 95 96 // GUS: Should I free r as well?96 // NOTE: Should I free r as well? 97 97 return(gauss); 98 98 } … … 320 320 float xSum = 1.0; 321 321 322 // GUS: Do we want to flag this case?322 // NOTE: Do we want to flag this case? 323 323 if (myPoly->n == 0) { 324 324 return(1.0); … … 646 646 double xSum = 1.0; 647 647 648 // GUS: Do we want to flag this case?648 // NOTE: Do we want to flag this case? 649 649 if (myPoly->n == 0) { 650 650 return(1.0); -
trunk/psLib/src/dataManip/psMinimize.c
r845 r1022 15 15 #include "float.h" 16 16 #include <math.h> 17 // GUS: rewrite so there is no maximum order for the polynomials.17 // NOTE: rewrite so there is no maximum order for the polynomials. 18 18 #define MAX_POLY_ORDER 10 19 19 #define MAX_POLYNOMIAL_TERMS (((MAX_POLY_ORDER+1) * (MAX_POLY_ORDER + 2)) / 2) … … 243 243 } 244 244 col[j] = 1.0; 245 // GUS: substitue the LUD rotine245 // NOTE: substitue the LUD rotine 246 246 // lubksb(A, N, indx, col); 247 247 for(i=1;i<=N;i++) { -
trunk/psLib/src/dataManip/psStats.c
r1020 r1022 95 95 float binSize = 0.0; // The histogram bin size 96 96 97 // GUS: Verify that this is the correct action.97 // NOTE: Verify that this is the correct action. 98 98 if (n == 0) { 99 99 return(NULL); 100 100 } 101 101 102 // GUS: Verify that this is the correct action.102 // NOTE: Verify that this is the correct action. 103 103 if (lower > upper) { 104 104 return(NULL); … … 146 146 int i; // Loop index variable 147 147 148 // GUS: Verify that this is the correct action.148 // NOTE: Verify that this is the correct action. 149 149 if (bounds == NULL) { 150 150 // psAbort(__func__, "psHistogram requested with NULL bounds"); … … 152 152 } 153 153 154 // GUS: Verify that this is the correct action.154 // NOTE: Verify that this is the correct action. 155 155 if (bounds->n <= 1) { 156 156 // psAbort(__func__, "psHistogram requested with NULL bounds"); … … 215 215 int numBins = 0; // The total number of bins 216 216 217 // GUS: Verify that this is the correct action.217 // NOTE: Verify that this is the correct action. 218 218 if (out == NULL) { 219 219 return(NULL); 220 220 } 221 221 222 // GUS: Verify that this is the correct action.222 // NOTE: Verify that this is the correct action. 223 223 if (in == NULL) { 224 224 return(out); … … 240 240 } 241 241 } 242 // GUS: determine the correct action for a variety of other cases:242 // NOTE: determine the correct action for a variety of other cases: 243 243 // in vector has 0 elements, and histogram structure has zero bins. 244 244 … … 268 268 // bin number requires a bit more work. 269 269 } else { 270 // GUS: This is slow. Put a smarter algorithm here to270 // NOTE: This is slow. Put a smarter algorithm here to 271 271 // find the correct bin number (bin search, probably) 272 272 for (j=0;j<(out->bounds->n)-1;j++) { … … 841 841 psVectorFree(unsortedVector); 842 842 psVectorFree(sortedVector); 843 // GUS: This is the843 // NOTE: This is the 844 844 } 845 845 … … 950 950 951 951 // Fit a Gaussian to the bins in the range MODE-dL to Mode+dL 952 // GUS: This step is dependent on the functions in psMinimize.c being952 // NOTE: This step is dependent on the functions in psMinimize.c being 953 953 // implemented. Currently, they are not. 954 954 … … 1122 1122 1123 1123 // 1. Compute the sample median. 1124 // GUS: This seems odd. Verify with IfA that we want to calculate the1124 // NOTE: This seems odd. Verify with IfA that we want to calculate the 1125 1125 // median here, not the mean. 1126 1126 p_psVectorSampleMedian(myVector, maskVector, maskVal, stats); … … 1198 1198 unsigned int maskVal) 1199 1199 { 1200 // GUS: Verify that this is the correct action.1200 // NOTE: Verify that this is the correct action. 1201 1201 if (in == NULL) { 1202 1202 return(stats); … … 1237 1237 1238 1238 // ************************************************************************ 1239 // GUS: The Stdev calculation requires the mean. Should we assume the1239 // NOTE: The Stdev calculation requires the mean. Should we assume the 1240 1240 // mean has already been calculated? Or should we always calculate it? 1241 1241 if (stats->options & PS_STAT_SAMPLE_STDEV) {
Note:
See TracChangeset
for help on using the changeset viewer.
