Changeset 21183 for trunk/psLib/src/math/psPolynomialUtils.c
- Timestamp:
- Jan 26, 2009, 8:40:07 PM (17 years ago)
- File:
-
- 1 edited
-
trunk/psLib/src/math/psPolynomialUtils.c (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/math/psPolynomialUtils.c
r19844 r21183 21 21 psStats *stats, 22 22 const psVector *mask, 23 ps MaskType maskValue,23 psVectorMaskType maskValue, 24 24 const psVector *f, 25 25 const psVector *fErr, … … 36 36 if (mask != NULL) { 37 37 PS_ASSERT_VECTORS_SIZE_EQUAL(f, mask, NULL); 38 PS_ASSERT_VECTOR_TYPE(mask, PS_TYPE_ U8, NULL);38 PS_ASSERT_VECTOR_TYPE(mask, PS_TYPE_VECTOR_MASK, NULL); 39 39 } 40 40 PS_ASSERT_VECTOR_NON_NULL(x, NULL); … … 106 106 // recovery is not allowed with this scheme 107 107 for (int i = 0; i < resid->n; i++) { 108 if ((mask != NULL) && (mask->data. U8[i] & maskValue)) {108 if ((mask != NULL) && (mask->data.PS_TYPE_VECTOR_MASK_DATA[i] & maskValue)) { 109 109 continue; 110 110 } … … 112 112 if (resid->data.F64[i] - stats->sampleMedian > sigma*maxClipSigma) { 113 113 if (mask != NULL) { 114 mask->data. U8[i] |= 0x01;114 mask->data.PS_TYPE_VECTOR_MASK_DATA[i] |= 0x01; 115 115 } 116 116 continue; … … 118 118 if (resid->data.F64[i] - stats->sampleMedian < sigma*minClipSigma) { 119 119 if (mask != NULL) { 120 mask->data. U8[i] |= 0x01;120 mask->data.PS_TYPE_VECTOR_MASK_DATA[i] |= 0x01; 121 121 } 122 122 continue;
Note:
See TracChangeset
for help on using the changeset viewer.
