Changeset 4029 for trunk/psLib/src/astro
- Timestamp:
- May 25, 2005, 10:26:55 AM (21 years ago)
- Location:
- trunk/psLib/src/astro
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/astro/psCoord.c
r3991 r4029 10 10 * @author GLG, MHPCC 11 11 * 12 * @version $Revision: 1. 69$ $Name: not supported by cvs2svn $13 * @date $Date: 2005-05-2 0 00:15:50$12 * @version $Revision: 1.70 $ $Name: not supported by cvs2svn $ 13 * @date $Date: 2005-05-25 20:26:55 $ 14 14 * 15 15 * Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii … … 94 94 psPlaneTransform *p_psPlaneTransformLinearInvert(psPlaneTransform *transform) 95 95 { 96 PS_ PTR_CHECK_NULL(transform, 0);97 PS_ PTR_CHECK_NULL(transform->x, 0);98 PS_ PTR_CHECK_NULL(transform->y, 0);96 PS_ASSERT_PTR_NON_NULL(transform, 0); 97 PS_ASSERT_PTR_NON_NULL(transform->x, 0); 98 PS_ASSERT_PTR_NON_NULL(transform->y, 0); 99 99 100 100 psF64 A = 0.0; … … 160 160 psS32 p_psIsProjectionLinear(psPlaneTransform *transform) 161 161 { 162 PS_ PTR_CHECK_NULL(transform, 0);163 PS_ PTR_CHECK_NULL(transform->x, 0);164 PS_ PTR_CHECK_NULL(transform->y, 0);162 PS_ASSERT_PTR_NON_NULL(transform, 0); 163 PS_ASSERT_PTR_NON_NULL(transform->x, 0); 164 PS_ASSERT_PTR_NON_NULL(transform->y, 0); 165 165 166 166 for (psS32 i=0;i<(transform->x->nX);i++) { … … 223 223 psPlaneTransform* psPlaneTransformAlloc(psS32 n1, psS32 n2) 224 224 { 225 PS_ INT_CHECK_NON_NEGATIVE(n1, NULL);226 PS_ INT_CHECK_NON_NEGATIVE(n2, NULL);225 PS_ASSERT_INT_NONNEGATIVE(n1, NULL); 226 PS_ASSERT_INT_NONNEGATIVE(n2, NULL); 227 227 228 228 psPlaneTransform *pt = psAlloc(sizeof(psPlaneTransform)); … … 238 238 const psPlane* coords) 239 239 { 240 PS_ PTR_CHECK_NULL(transform, NULL);241 PS_ PTR_CHECK_NULL(transform->x, NULL);242 PS_ PTR_CHECK_NULL(transform->y, NULL);243 PS_ PTR_CHECK_NULL(coords, NULL);240 PS_ASSERT_PTR_NON_NULL(transform, NULL); 241 PS_ASSERT_PTR_NON_NULL(transform->x, NULL); 242 PS_ASSERT_PTR_NON_NULL(transform->y, NULL); 243 PS_ASSERT_PTR_NON_NULL(coords, NULL); 244 244 245 245 if (out == NULL) { … … 267 267 psPlaneDistort* psPlaneDistortAlloc(psS32 n1, psS32 n2, psS32 n3, psS32 n4) 268 268 { 269 PS_ INT_CHECK_NON_NEGATIVE(n1, NULL);270 PS_ INT_CHECK_NON_NEGATIVE(n2, NULL);271 PS_ INT_CHECK_NON_NEGATIVE(n3, NULL);272 PS_ INT_CHECK_NON_NEGATIVE(n4, NULL);269 PS_ASSERT_INT_NONNEGATIVE(n1, NULL); 270 PS_ASSERT_INT_NONNEGATIVE(n2, NULL); 271 PS_ASSERT_INT_NONNEGATIVE(n3, NULL); 272 PS_ASSERT_INT_NONNEGATIVE(n4, NULL); 273 273 274 274 psPlaneDistort *pt = psAlloc(sizeof(psPlaneDistort)); … … 290 290 float magnitude) 291 291 { 292 PS_ PTR_CHECK_NULL(transform, NULL);293 PS_ PTR_CHECK_NULL(transform->x, NULL);294 PS_ PTR_CHECK_NULL(transform->y, NULL);295 PS_ PTR_CHECK_NULL(coords, NULL);292 PS_ASSERT_PTR_NON_NULL(transform, NULL); 293 PS_ASSERT_PTR_NON_NULL(transform->x, NULL); 294 PS_ASSERT_PTR_NON_NULL(transform->y, NULL); 295 PS_ASSERT_PTR_NON_NULL(coords, NULL); 296 296 297 297 if (out == NULL) { … … 364 364 const psSphere* coord) 365 365 { 366 PS_ PTR_CHECK_NULL(transform, NULL);367 PS_ PTR_CHECK_NULL(coord, NULL);366 PS_ASSERT_PTR_NON_NULL(transform, NULL); 367 PS_ASSERT_PTR_NON_NULL(coord, NULL); 368 368 369 369 if (out == NULL) { … … 394 394 395 395 // Check for null parameter 396 PS_ PTR_CHECK_NULL(time, NULL);396 PS_ASSERT_PTR_NON_NULL(time, NULL); 397 397 398 398 // Convert psTime to MJD … … 427 427 428 428 // Check for null parameter 429 PS_ PTR_CHECK_NULL(time, NULL);429 PS_ASSERT_PTR_NON_NULL(time, NULL); 430 430 431 431 // Convert psTime to MJD … … 498 498 const psProjection* projection) 499 499 { 500 PS_ PTR_CHECK_NULL(coord, NULL);501 PS_ PTR_CHECK_NULL(projection, NULL);500 PS_ASSERT_PTR_NON_NULL(coord, NULL); 501 PS_ASSERT_PTR_NON_NULL(projection, NULL); 502 502 503 503 psF64 theta = 0.0; … … 550 550 const psProjection* projection) 551 551 { 552 PS_ PTR_CHECK_NULL(coord, NULL);553 PS_ PTR_CHECK_NULL(projection, NULL);552 PS_ASSERT_PTR_NON_NULL(coord, NULL); 553 PS_ASSERT_PTR_NON_NULL(projection, NULL); 554 554 555 555 psF64 theta = 0.0; … … 622 622 psSphereOffsetUnit unit) 623 623 { 624 PS_ PTR_CHECK_NULL(position1, NULL);625 PS_ PTR_CHECK_NULL(position2, NULL);624 PS_ASSERT_PTR_NON_NULL(position1, NULL); 625 PS_ASSERT_PTR_NON_NULL(position2, NULL); 626 626 627 627 // Check positions near 90 degree and issue warnings if necessary … … 721 721 psSphereOffsetUnit unit) 722 722 { 723 PS_ PTR_CHECK_NULL(position, NULL);724 PS_ PTR_CHECK_NULL(offset, NULL);723 PS_ASSERT_PTR_NON_NULL(position, NULL); 724 PS_ASSERT_PTR_NON_NULL(offset, NULL); 725 725 726 726 psSphere* tmp; … … 811 811 { 812 812 // Check input for NULL pointers 813 PS_ PTR_CHECK_NULL(coords, NULL);814 PS_ PTR_CHECK_NULL(fromTime, NULL);815 PS_ PTR_CHECK_NULL(toTime, NULL);813 PS_ASSERT_PTR_NON_NULL(coords, NULL); 814 PS_ASSERT_PTR_NON_NULL(fromTime, NULL); 815 PS_ASSERT_PTR_NON_NULL(toTime, NULL); 816 816 817 817 // Calculate Julian centuries … … 904 904 // XXX: This does not yet use region and nSamples: need to modify -rdd 905 905 906 PS_ PTR_CHECK_NULL(trans1, NULL);907 PS_ PTR_CHECK_NULL(trans2, NULL);906 PS_ASSERT_PTR_NON_NULL(trans1, NULL); 907 PS_ASSERT_PTR_NON_NULL(trans2, NULL); 908 908 //TRACE: printf("psPlaneTransformCombine(%d, %d, %d, %d: %d, %d, %d, %d)\n", trans1->x->nX, trans1->x->nY, trans1->y->nX, trans1->y->nY, trans2->x->nX, trans2->x->nY, trans2->y->nX, trans2->y->nY); 909 909 // … … 1047 1047 float sigmaClip) 1048 1048 { 1049 PS_ PTR_CHECK_NULL(trans, NULL);1050 PS_ PTR_CHECK_NULL(source, NULL);1051 PS_ PTR_CHECK_NULL(dest, NULL);1049 PS_ASSERT_PTR_NON_NULL(trans, NULL); 1050 PS_ASSERT_PTR_NON_NULL(source, NULL); 1051 PS_ASSERT_PTR_NON_NULL(dest, NULL); 1052 1052 1053 1053 psS32 numCoords = PS_MIN(source->n, dest->n); … … 1159 1159 int nSamples) 1160 1160 { 1161 PS_ PTR_CHECK_NULL(in, NULL);1161 PS_ASSERT_PTR_NON_NULL(in, NULL); 1162 1162 // 1163 1163 // If the transform is linear, then invert it exactly and return. -
trunk/psLib/src/astro/psTime.c
r3884 r4029 10 10 * @author Ross Harman, MHPCC 11 11 * 12 * @version $Revision: 1.6 0$ $Name: not supported by cvs2svn $13 * @date $Date: 2005-05- 11 22:02:15 $12 * @version $Revision: 1.61 $ $Name: not supported by cvs2svn $ 13 * @date $Date: 2005-05-25 20:26:55 $ 14 14 * 15 15 * Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii … … 178 178 } 179 179 table = (psLookupTable*)tableMetadataItem->data.V; 180 PS_ PTR_CHECK_NULL(table,0.0);180 PS_ASSERT_PTR_NON_NULL(table,0.0); 181 181 182 182 // Attempt to interpolate table … … 403 403 404 404 // Error checks 405 PS_ PTR_CHECK_NULL(time,NULL);406 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NULL);405 PS_ASSERT_PTR_NON_NULL(time,NULL); 406 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NULL); 407 407 408 408 if (time->type == type) { // time already right type. That was easy! … … 448 448 449 449 // Error checks 450 PS_ PTR_CHECK_NULL(time,NAN);451 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NAN);450 PS_ASSERT_PTR_NON_NULL(time,NAN); 451 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NAN); 452 452 453 453 // Calculate TAI or UTC time based on type of time user passes … … 519 519 520 520 // Error checks 521 PS_ PTR_CHECK_NULL(time,NAN);522 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NAN);521 PS_ASSERT_PTR_NON_NULL(time,NAN); 522 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NAN); 523 523 524 524 if(time->type != PS_TIME_TAI) { … … 560 560 } 561 561 dut = (psVector*)tableMetadataItem->data.V; 562 PS_ PTR_CHECK_NULL(dut,0.0);562 PS_ASSERT_PTR_NON_NULL(dut,0.0); 563 563 564 564 t = 2000.0 + (mjd - 51544.03)/365.2422; … … 592 592 593 593 // Error checks 594 PS_ PTR_CHECK_NULL(time,NULL);595 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NULL);594 PS_ASSERT_PTR_NON_NULL(time,NULL); 595 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NULL); 596 596 597 597 if(time->type != PS_TIME_TAI) … … 658 658 } 659 659 xp = (psVector*)tableMetadataItem->data.V; 660 PS_ PTR_CHECK_NULL(xp,NULL);660 PS_ASSERT_PTR_NON_NULL(xp,NULL); 661 661 662 662 // Get yp … … 667 667 } 668 668 yp = (psVector*)tableMetadataItem->data.V; 669 PS_ PTR_CHECK_NULL(yp,NULL);669 PS_ASSERT_PTR_NON_NULL(yp,NULL); 670 670 671 671 // Calculate "a" and "c" constants … … 715 715 716 716 // Error checks 717 PS_ PTR_CHECK_NULL(time,NAN);718 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NAN);717 PS_ASSERT_PTR_NON_NULL(time,NAN); 718 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NAN); 719 719 720 720 // Check time metadata … … 733 733 } 734 734 table = (psLookupTable*)tableMetadataItem->data.V; 735 PS_ PTR_CHECK_NULL(table,0);735 PS_ASSERT_PTR_NON_NULL(table,0); 736 736 737 737 // Determine Julian and modified Julian dates used in table lookup and time delta calculation … … 767 767 768 768 // Error checks 769 PS_ PTR_CHECK_NULL(time1,0);770 PS_ PTR_CHECK_NULL(time2,0);771 PS_ INT_CHECK_RANGE(time1->nsec,0,1e9-1,0);772 PS_ INT_CHECK_RANGE(time2->nsec,0,1e9-1,0);769 PS_ASSERT_PTR_NON_NULL(time1,0); 770 PS_ASSERT_PTR_NON_NULL(time2,0); 771 PS_ASSERT_INT_WITHIN_RANGE(time1->nsec,0,1e9-1,0); 772 PS_ASSERT_INT_WITHIN_RANGE(time2->nsec,0,1e9-1,0); 773 773 diff = abs((psS64)psTimeGetTAIDelta((psTime*)time1)-(psS64)psTimeGetTAIDelta((psTime*)time2)); 774 774 … … 781 781 782 782 // Error checks 783 PS_ PTR_CHECK_NULL(time,NAN);784 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NAN);783 PS_ASSERT_PTR_NON_NULL(time,NAN); 784 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NAN); 785 785 786 786 // Julian date conversion … … 800 800 801 801 // Error checks 802 PS_ PTR_CHECK_NULL(time,NAN);803 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NAN);802 PS_ASSERT_PTR_NON_NULL(time,NAN); 803 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NAN); 804 804 805 805 // Modified Julian date conversion … … 823 823 824 824 // Error checks 825 PS_ PTR_CHECK_NULL(time,NULL);826 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NULL);825 PS_ASSERT_PTR_NON_NULL(time,NULL); 826 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NULL); 827 827 828 828 tempString = psAlloc(MAX_TIME_STRING_LENGTH); … … 856 856 857 857 // Error checks 858 PS_ PTR_CHECK_NULL(time,timevalTime);859 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,timevalTime);858 PS_ASSERT_PTR_NON_NULL(time,timevalTime); 859 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,timevalTime); 860 860 861 861 timevalTime.tv_sec = time->sec; … … 879 879 880 880 // Error checks 881 PS_ PTR_CHECK_NULL(time,NULL);882 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NULL);881 PS_ASSERT_PTR_NON_NULL(time,NULL); 882 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NULL); 883 883 884 884 seconds = time->sec%60; … … 949 949 950 950 // Error check 951 PS_ INT_CHECK_RANGE(outTime->nsec,0,1e9-1,outTime);951 PS_ASSERT_INT_WITHIN_RANGE(outTime->nsec,0,1e9-1,outTime); 952 952 953 953 return outTime; … … 975 975 976 976 // Error check 977 PS_ INT_CHECK_RANGE(outTime->nsec,0,1e9-1,outTime);977 PS_ASSERT_INT_WITHIN_RANGE(outTime->nsec,0,1e9-1,outTime); 978 978 979 979 return outTime; … … 993 993 } 994 994 995 PS_ INT_CHECK_NON_NEGATIVE(tmTime.tm_year, outTime);996 PS_ INT_CHECK_RANGE(tmTime.tm_mon,1,12,outTime);997 PS_ INT_CHECK_RANGE(tmTime.tm_mday,1,31,outTime);998 PS_ INT_CHECK_RANGE(tmTime.tm_hour,0,23,outTime);999 PS_ INT_CHECK_RANGE(tmTime.tm_min,0,59,outTime);1000 PS_ INT_CHECK_RANGE(tmTime.tm_sec,0,59,outTime);1001 PS_ INT_CHECK_RANGE(millisecond,0,999,outTime);995 PS_ASSERT_INT_NONNEGATIVE(tmTime.tm_year, outTime); 996 PS_ASSERT_INT_WITHIN_RANGE(tmTime.tm_mon,1,12,outTime); 997 PS_ASSERT_INT_WITHIN_RANGE(tmTime.tm_mday,1,31,outTime); 998 PS_ASSERT_INT_WITHIN_RANGE(tmTime.tm_hour,0,23,outTime); 999 PS_ASSERT_INT_WITHIN_RANGE(tmTime.tm_min,0,59,outTime); 1000 PS_ASSERT_INT_WITHIN_RANGE(tmTime.tm_sec,0,59,outTime); 1001 PS_ASSERT_INT_WITHIN_RANGE(millisecond,0,999,outTime); 1002 1002 1003 1003 tmTime.tm_year -= 1900; … … 1018 1018 1019 1019 // Error check 1020 PS_ PTR_CHECK_NULL(time,NULL);1020 PS_ASSERT_PTR_NON_NULL(time,NULL); 1021 1021 1022 1022 // Allocate psTime struct … … 1028 1028 1029 1029 // Error check 1030 PS_ INT_CHECK_RANGE(outTime->nsec,0,1e9-1,outTime);1030 PS_ASSERT_INT_WITHIN_RANGE(outTime->nsec,0,1e9-1,outTime); 1031 1031 1032 1032 return outTime; … … 1046 1046 1047 1047 // Error check 1048 PS_ PTR_CHECK_NULL(time,NULL);1048 PS_ASSERT_PTR_NON_NULL(time,NULL); 1049 1049 1050 1050 // Allocate psTime struct … … 1087 1087 1088 1088 // Error check 1089 PS_ INT_CHECK_RANGE(outTime->nsec,0,1e9-1,outTime);1089 PS_ASSERT_INT_WITHIN_RANGE(outTime->nsec,0,1e9-1,outTime); 1090 1090 1091 1091 return outTime; … … 1100 1100 1101 1101 // Error checks 1102 PS_ PTR_CHECK_NULL(time,NULL);1103 PS_ INT_CHECK_RANGE(time->nsec,0,1e9-1,NULL);1102 PS_ASSERT_PTR_NON_NULL(time,NULL); 1103 PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,1e9-1,NULL); 1104 1104 1105 1105 // Convert time to TAI if necessary, but without changing input arguments … … 1120 1120 1121 1121 // Error check 1122 PS_ INT_CHECK_RANGE(outTime->nsec,0,1e9-1,outTime);1122 PS_ASSERT_INT_WITHIN_RANGE(outTime->nsec,0,1e9-1,outTime); 1123 1123 1124 1124 // Convert result to same time type as input … … 1142 1142 1143 1143 // Error checks 1144 PS_ PTR_CHECK_NULL(time1,0.0);1145 PS_ INT_CHECK_RANGE(time1->nsec,0,1e9-1,0.0);1146 PS_ PTR_CHECK_NULL(time2,0.0);1147 PS_ INT_CHECK_RANGE(time2->nsec,0,1e9-1,0.0);1144 PS_ASSERT_PTR_NON_NULL(time1,0.0); 1145 PS_ASSERT_INT_WITHIN_RANGE(time1->nsec,0,1e9-1,0.0); 1146 PS_ASSERT_PTR_NON_NULL(time2,0.0); 1147 PS_ASSERT_INT_WITHIN_RANGE(time2->nsec,0,1e9-1,0.0); 1148 1148 1149 1149 // Convert time to TAI if necessary, but without changing input arguments
Note:
See TracChangeset
for help on using the changeset viewer.
