Changeset 4029 for trunk/psLib/src/math/psStats.c
- Timestamp:
- May 25, 2005, 10:26:55 AM (21 years ago)
- File:
-
- 1 edited
-
trunk/psLib/src/math/psStats.c (modified) (12 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/math/psStats.c
r4028 r4029 9 9 * @author GLG, MHPCC 10 10 * 11 * @version $Revision: 1.12 8$ $Name: not supported by cvs2svn $12 * @date $Date: 2005-05-25 19:25:29$11 * @version $Revision: 1.129 $ $Name: not supported by cvs2svn $ 12 * @date $Date: 2005-05-25 20:26:55 $ 13 13 * 14 14 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 424 424 psStats* stats) 425 425 { 426 PS_ VECTOR_CHECK_NULL(myVector, -1);427 PS_ PTR_CHECK_NULL(stats, -1);426 PS_ASSERT_VECTOR_NON_NULL(myVector, -1); 427 PS_ASSERT_PTR_NON_NULL(stats, -1); 428 428 psS32 i = 0; // Loop index variable 429 429 … … 472 472 psStats* stats) 473 473 { 474 PS_ VECTOR_CHECK_NULL(myVector, -1);475 PS_ PTR_CHECK_NULL(stats, -1);474 PS_ASSERT_VECTOR_NON_NULL(myVector, -1); 475 PS_ASSERT_PTR_NON_NULL(stats, -1); 476 476 psS32 i = 0; // Loop index variable 477 477 psS32 numData = 0; // The number of data points … … 619 619 psF32 sigma) 620 620 { 621 PS_ PTR_CHECK_NULL(robustHistogram, NULL);622 PS_ PTR_CHECK_NULL(robustHistogram->bounds, NULL);621 PS_ASSERT_PTR_NON_NULL(robustHistogram, NULL); 622 PS_ASSERT_PTR_NON_NULL(robustHistogram->bounds, NULL); 623 623 624 624 psS32 i = 0; // Loop index variable … … 1029 1029 1030 1030 // Ensure that stats->clipIter is within the proper range. 1031 PS_ INT_CHECK_RANGE(stats->clipIter,1032 PS_CLIPPED_NUM_ITER_LB,1033 PS_CLIPPED_NUM_ITER_UB, -1);1031 PS_ASSERT_INT_WITHIN_RANGE(stats->clipIter, 1032 PS_CLIPPED_NUM_ITER_LB, 1033 PS_CLIPPED_NUM_ITER_UB, -1); 1034 1034 1035 1035 // Ensure that stats->clipSigma is within the proper range. 1036 PS_ INT_CHECK_RANGE(stats->clipSigma,1037 PS_CLIPPED_SIGMA_LB,1038 PS_CLIPPED_SIGMA_UB, -1);1036 PS_ASSERT_INT_WITHIN_RANGE(stats->clipSigma, 1037 PS_CLIPPED_SIGMA_LB, 1038 PS_CLIPPED_SIGMA_UB, -1); 1039 1039 1040 1040 // Allocate a psStats structure for calculating the mean, median, and … … 1263 1263 psF32 getThisValue) 1264 1264 { 1265 PS_ POLY_CHECK_NULL(myPoly, NAN);1265 PS_ASSERT_POLY_NON_NULL(myPoly, NAN); 1266 1266 PS_FLOAT_COMPARE(rangeLow, rangeHigh, NAN); 1267 1267 // We ensure that the requested f(y) value, which is getThisValue, is … … 1332 1332 psF32 yVal) 1333 1333 { 1334 PS_ VECTOR_CHECK_NULL(xVec, NAN);1335 PS_ VECTOR_CHECK_NULL(yVec, NAN);1336 PS_ VECTOR_CHECK_TYPE(xVec, PS_TYPE_F32, NAN);1337 PS_ VECTOR_CHECK_TYPE(yVec, PS_TYPE_F32, NAN);1338 // PS_ VECTOR_CHECK_SIZE_EQUAL(xVec, yVec, NAN);1339 PS_ INT_CHECK_RANGE(binNum, 0, (xVec->n - 1), NAN);1340 PS_ INT_CHECK_RANGE(binNum, 0, (yVec->n - 1), NAN);1334 PS_ASSERT_VECTOR_NON_NULL(xVec, NAN); 1335 PS_ASSERT_VECTOR_NON_NULL(yVec, NAN); 1336 PS_ASSERT_VECTOR_TYPE(xVec, PS_TYPE_F32, NAN); 1337 PS_ASSERT_VECTOR_TYPE(yVec, PS_TYPE_F32, NAN); 1338 // PS_ASSERT_VECTORS_SIZE_EQUAL(xVec, yVec, NAN); 1339 PS_ASSERT_INT_WITHIN_RANGE(binNum, 0, (xVec->n - 1), NAN); 1340 PS_ASSERT_INT_WITHIN_RANGE(binNum, 0, (yVec->n - 1), NAN); 1341 1341 1342 1342 // PS_VECTOR_DECLARE_ALLOC_STATIC(x, 3, PS_TYPE_F64); … … 1827 1827 psHistogram* psHistogramAlloc(psF32 lower, psF32 upper, psS32 n) 1828 1828 { 1829 PS_ INT_CHECK_POSITIVE(n, NULL);1829 PS_ASSERT_INT_POSITIVE(n, NULL); 1830 1830 PS_FLOAT_COMPARE(lower, upper, NULL); 1831 1831 … … 1875 1875 psHistogram* psHistogramAllocGeneric(const psVector* bounds) 1876 1876 { 1877 PS_ VECTOR_CHECK_NULL(bounds, NULL);1878 PS_ VECTOR_CHECK_TYPE(bounds, PS_TYPE_F32, NULL);1877 PS_ASSERT_VECTOR_NON_NULL(bounds, NULL); 1878 PS_ASSERT_VECTOR_TYPE(bounds, PS_TYPE_F32, NULL); 1879 1879 PS_INT_COMPARE(2, bounds->n, NULL); 1880 1880 … … 1930 1930 psF32 error) 1931 1931 { 1932 PS_ PTR_CHECK_NULL(out, -1);1933 PS_ PTR_CHECK_NULL(out->bounds, -1);1934 PS_ PTR_CHECK_NULL(out->nums, -1);1935 PS_ INT_CHECK_RANGE(binNum, 0, ((out->nums->n)-1), -2);1932 PS_ASSERT_PTR_NON_NULL(out, -1); 1933 PS_ASSERT_PTR_NON_NULL(out->bounds, -1); 1934 PS_ASSERT_PTR_NON_NULL(out->nums, -1); 1935 PS_ASSERT_INT_WITHIN_RANGE(binNum, 0, ((out->nums->n)-1), -2); 1936 1936 PS_FLOAT_COMPARE(0.0, error, -3); 1937 PS_ FLOAT_CHECK_RANGE(data, out->bounds->data.F32[0], out->bounds->data.F32[(out->bounds->n)-1], -4);1937 PS_ASSERT_FLOAT_WITHIN_RANGE(data, out->bounds->data.F32[0], out->bounds->data.F32[(out->bounds->n)-1], -4); 1938 1938 1939 1939 psF32 boxcarWidth = 2.35 * error; … … 2019 2019 psU32 maskVal) 2020 2020 { 2021 PS_ PTR_CHECK_NULL(out, NULL);2022 PS_ VECTOR_CHECK_NULL(out->bounds, NULL);2023 PS_ VECTOR_CHECK_TYPE(out->bounds, PS_TYPE_F32, NULL);2024 PS_ INT_CHECK_NON_NEGATIVE(out->bounds->n, NULL);2025 PS_ VECTOR_CHECK_NULL(out->nums, NULL);2026 PS_ VECTOR_CHECK_TYPE(out->nums, PS_TYPE_F32, NULL);2027 PS_ INT_CHECK_NON_NEGATIVE(out->nums->n, NULL);2028 PS_ VECTOR_CHECK_NULL(in, out);2021 PS_ASSERT_PTR_NON_NULL(out, NULL); 2022 PS_ASSERT_VECTOR_NON_NULL(out->bounds, NULL); 2023 PS_ASSERT_VECTOR_TYPE(out->bounds, PS_TYPE_F32, NULL); 2024 PS_ASSERT_INT_NONNEGATIVE(out->bounds->n, NULL); 2025 PS_ASSERT_VECTOR_NON_NULL(out->nums, NULL); 2026 PS_ASSERT_VECTOR_TYPE(out->nums, PS_TYPE_F32, NULL); 2027 PS_ASSERT_INT_NONNEGATIVE(out->nums->n, NULL); 2028 PS_ASSERT_VECTOR_NON_NULL(in, out); 2029 2029 if (mask != NULL) { 2030 PS_ VECTOR_CHECK_SIZE_EQUAL(in, mask, NULL);2031 PS_ VECTOR_CHECK_TYPE(mask, PS_TYPE_U8, NULL);2030 PS_ASSERT_VECTORS_SIZE_EQUAL(in, mask, NULL); 2031 PS_ASSERT_VECTOR_TYPE(mask, PS_TYPE_U8, NULL); 2032 2032 } 2033 2033 if (errors != NULL) { 2034 PS_ VECTOR_CHECK_SIZE_EQUAL(in, errors, NULL);2035 PS_ VECTOR_CHECK_TYPE(errors, in->type.type, NULL);2034 PS_ASSERT_VECTORS_SIZE_EQUAL(in, errors, NULL); 2035 PS_ASSERT_VECTOR_TYPE(errors, in->type.type, NULL); 2036 2036 } 2037 2037 … … 2216 2216 psU32 maskVal) 2217 2217 { 2218 PS_ PTR_CHECK_NULL(stats, NULL);2219 PS_ VECTOR_CHECK_NULL(in, stats);2218 PS_ASSERT_PTR_NON_NULL(stats, NULL); 2219 PS_ASSERT_VECTOR_NON_NULL(in, stats); 2220 2220 if (mask != NULL) { 2221 PS_ VECTOR_CHECK_SIZE_EQUAL(mask, in, stats);2222 PS_ VECTOR_CHECK_TYPE(mask, PS_TYPE_U8, stats);2221 PS_ASSERT_VECTORS_SIZE_EQUAL(mask, in, stats); 2222 PS_ASSERT_VECTOR_TYPE(mask, PS_TYPE_U8, stats); 2223 2223 } 2224 2224 if (errors != NULL) { 2225 PS_ VECTOR_CHECK_SIZE_EQUAL(errors, in, stats);2226 PS_ VECTOR_CHECK_TYPE(errors, in->type.type, stats);2225 PS_ASSERT_VECTORS_SIZE_EQUAL(errors, in, stats); 2226 PS_ASSERT_VECTOR_TYPE(errors, in->type.type, stats); 2227 2227 } 2228 2228
Note:
See TracChangeset
for help on using the changeset viewer.
