Changeset 4029 for trunk/psLib/src/astronomy
- Timestamp:
- May 25, 2005, 10:26:55 AM (21 years ago)
- Location:
- trunk/psLib/src/astronomy
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/astronomy/psAstrometry.c
r3708 r4029 8 8 * @author GLG, MHPCC 9 9 * 10 * @version $Revision: 1.6 5$ $Name: not supported by cvs2svn $11 * @date $Date: 2005-0 4-19 00:17:05 $10 * @version $Revision: 1.66 $ $Name: not supported by cvs2svn $ 11 * @date $Date: 2005-05-25 20:26:55 $ 12 12 * 13 13 * Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii … … 41 41 psImage* tmpImage) 42 42 { 43 PS_ IMAGE_CHECK_NULL(tmpImage, 0);43 PS_ASSERT_IMAGE_NON_NULL(tmpImage, 0); 44 44 45 45 if ((x < 0.0) || (x > (double)tmpImage->numCols) || … … 154 154 psS32 j; 155 155 156 PS_ IMAGE_CHECK_NULL(x, NULL);157 PS_ IMAGE_CHECK_NULL(y, NULL);158 PS_ IMAGE_CHECK_TYPE(x, PS_TYPE_F64, NULL);159 PS_ IMAGE_CHECK_TYPE(y, PS_TYPE_F64, NULL);156 PS_ASSERT_IMAGE_NON_NULL(x, NULL); 157 PS_ASSERT_IMAGE_NON_NULL(y, NULL); 158 PS_ASSERT_IMAGE_TYPE(x, PS_TYPE_F64, NULL); 159 PS_ASSERT_IMAGE_TYPE(y, PS_TYPE_F64, NULL); 160 160 161 161 tmp = (psFixedPattern *) psAlloc(sizeof(psFixedPattern)); … … 211 211 const psObservatory* observatory) 212 212 { 213 PS_ PTR_CHECK_NULL(observatory, NULL);213 PS_ASSERT_PTR_NON_NULL(observatory, NULL); 214 214 215 215 psExposure* exp = psAlloc(sizeof(psExposure)); … … 271 271 const psExposure* exp) 272 272 { 273 PS_ INT_CHECK_NON_NEGATIVE(nChips, NULL);273 PS_ASSERT_INT_NONNEGATIVE(nChips, NULL); 274 274 275 275 psFPA* newFPA = psAlloc(sizeof(psFPA)); … … 315 315 psFPA *parentFPA) 316 316 { 317 PS_ INT_CHECK_NON_NEGATIVE(nCells, NULL);317 PS_ASSERT_INT_NONNEGATIVE(nCells, NULL); 318 318 319 319 psChip* chip = psAlloc(sizeof(psChip)); … … 350 350 psChip* parentChip) 351 351 { 352 PS_ INT_CHECK_NON_NEGATIVE(nReadouts, NULL);352 PS_ASSERT_INT_NONNEGATIVE(nReadouts, NULL); 353 353 354 354 psCell* cell = psAlloc(sizeof(psCell)); … … 406 406 psGrommit* psGrommitAlloc(const psExposure* exp) 407 407 { 408 PS_ PTR_CHECK_NULL(exp, NULL);408 PS_ASSERT_PTR_NON_NULL(exp, NULL); 409 409 410 410 double date = psTimeToMJD(exp->time); … … 436 436 const psFPA* FPA) 437 437 { 438 PS_ PTR_CHECK_NULL(fpaCoord, NULL);439 PS_ PTR_CHECK_NULL(FPA, NULL);438 PS_ASSERT_PTR_NON_NULL(fpaCoord, NULL); 439 PS_ASSERT_PTR_NON_NULL(FPA, NULL); 440 440 441 441 psChip* tmpChip = NULL; … … 461 461 const psFPA* FPA) 462 462 { 463 PS_ PTR_CHECK_NULL(fpaCoord, NULL);464 PS_ PTR_CHECK_NULL(FPA, NULL);465 PS_ PTR_CHECK_NULL(FPA->chips, NULL);463 PS_ASSERT_PTR_NON_NULL(fpaCoord, NULL); 464 PS_ASSERT_PTR_NON_NULL(FPA, NULL); 465 PS_ASSERT_PTR_NON_NULL(FPA->chips, NULL); 466 466 467 467 psArray* chips = FPA->chips; … … 476 476 for (psS32 i = 0; i < nChips; i++) { 477 477 psChip* tmpChip = chips->data[i]; 478 PS_ PTR_CHECK_NULL(tmpChip, NULL);479 PS_ PTR_CHECK_NULL(tmpChip->fromFPA, NULL);478 PS_ASSERT_PTR_NON_NULL(tmpChip, NULL); 479 PS_ASSERT_PTR_NON_NULL(tmpChip->fromFPA, NULL); 480 480 481 481 psPlaneTransformApply(&chipCoord, tmpChip->fromFPA, fpaCoord); … … 494 494 const psChip* chip) 495 495 { 496 PS_ PTR_CHECK_NULL(chipCoord, NULL);497 PS_ PTR_CHECK_NULL(chip, NULL);496 PS_ASSERT_PTR_NON_NULL(chipCoord, NULL); 497 PS_ASSERT_PTR_NON_NULL(chip, NULL); 498 498 499 499 psPlane cellCoord; … … 511 511 for (psS32 i = 0; i < cells->n; i++) { 512 512 psCell* tmpCell = (psCell* ) cells->data[i]; 513 PS_ PTR_CHECK_NULL(tmpCell, NULL);514 PS_ PTR_CHECK_NULL(tmpCell->fromChip, NULL);513 PS_ASSERT_PTR_NON_NULL(tmpCell, NULL); 514 PS_ASSERT_PTR_NON_NULL(tmpCell->fromChip, NULL); 515 515 psArray* readouts = tmpCell->readouts; 516 516 … … 518 518 for (psS32 j = 0; j < readouts->n; j++) { 519 519 psReadout* tmpReadout = readouts->data[j]; 520 PS_ READOUT_CHECK_NULL(tmpReadout, NULL);520 PS_ASSERT_READOUT_NON_NULL(tmpReadout, NULL); 521 521 522 522 psPlaneTransformApply(&cellCoord, … … 540 540 const psCell* cell) 541 541 { 542 PS_ PTR_CHECK_NULL(inCoord, NULL);543 PS_ PTR_CHECK_NULL(cell, NULL);542 PS_ASSERT_PTR_NON_NULL(inCoord, NULL); 543 PS_ASSERT_PTR_NON_NULL(cell, NULL); 544 544 545 545 return (psPlaneTransformApply(outCoord, cell->toChip, inCoord)); … … 550 550 const psChip* chip) 551 551 { 552 PS_ PTR_CHECK_NULL(inCoord, NULL);553 PS_ PTR_CHECK_NULL(chip, NULL);552 PS_ASSERT_PTR_NON_NULL(inCoord, NULL); 553 PS_ASSERT_PTR_NON_NULL(chip, NULL); 554 554 555 555 return (psPlaneTransformApply(outCoord, chip->toFPA, inCoord)); … … 562 562 const psFPA* fpa) 563 563 { 564 PS_ PTR_CHECK_NULL(inCoord, NULL);565 PS_ PTR_CHECK_NULL(fpa, NULL);564 PS_ASSERT_PTR_NON_NULL(inCoord, NULL); 565 PS_ASSERT_PTR_NON_NULL(fpa, NULL); 566 566 567 567 return(psPlaneDistortApply(outCoord, fpa->toTangentPlane, inCoord, … … 576 576 const psGrommit* grommit) 577 577 { 578 PS_ PTR_CHECK_NULL(tpCoord, NULL);579 PS_ PTR_CHECK_NULL(grommit, NULL);578 PS_ASSERT_PTR_NON_NULL(tpCoord, NULL); 579 PS_ASSERT_PTR_NON_NULL(grommit, NULL); 580 580 581 581 double AOB = 0.0; … … 596 596 const psCell* cell) 597 597 { 598 PS_ PTR_CHECK_NULL(cellCoord, NULL);599 PS_ PTR_CHECK_NULL(cell, NULL);598 PS_ASSERT_PTR_NON_NULL(cellCoord, NULL); 599 PS_ASSERT_PTR_NON_NULL(cell, NULL); 600 600 601 601 return (psPlaneTransformApply(fpaCoord, cell->toFPA, cellCoord)); … … 608 608 const psCell* cell) 609 609 { 610 PS_ PTR_CHECK_NULL(cellCoord, NULL);611 PS_ PTR_CHECK_NULL(cell, NULL);612 PS_ PTR_CHECK_NULL(cell->toFPA, NULL);613 PS_ PTR_CHECK_NULL(cell->parent, NULL);614 PS_ PTR_CHECK_NULL(cell->parent->parent, NULL);615 PS_ PTR_CHECK_NULL(cell->parent->parent->toTangentPlane, NULL);616 PS_ PTR_CHECK_NULL(cell->parent->parent->exposure, NULL);610 PS_ASSERT_PTR_NON_NULL(cellCoord, NULL); 611 PS_ASSERT_PTR_NON_NULL(cell, NULL); 612 PS_ASSERT_PTR_NON_NULL(cell->toFPA, NULL); 613 PS_ASSERT_PTR_NON_NULL(cell->parent, NULL); 614 PS_ASSERT_PTR_NON_NULL(cell->parent->parent, NULL); 615 PS_ASSERT_PTR_NON_NULL(cell->parent->parent->toTangentPlane, NULL); 616 PS_ASSERT_PTR_NON_NULL(cell->parent->parent->exposure, NULL); 617 617 618 618 psPlane* fpaCoord = NULL; … … 645 645 const psCell* cell) 646 646 { 647 PS_ PTR_CHECK_NULL(cellCoord, NULL);648 PS_ PTR_CHECK_NULL(cell, NULL);649 PS_ PTR_CHECK_NULL(cell->toSky, NULL);650 PS_ PTR_CHECK_NULL(cell->parent, NULL);651 PS_ PTR_CHECK_NULL(cell->parent->parent, NULL);652 PS_ PTR_CHECK_NULL(cell->parent->parent->projection, NULL);647 PS_ASSERT_PTR_NON_NULL(cellCoord, NULL); 648 PS_ASSERT_PTR_NON_NULL(cell, NULL); 649 PS_ASSERT_PTR_NON_NULL(cell->toSky, NULL); 650 PS_ASSERT_PTR_NON_NULL(cell->parent, NULL); 651 PS_ASSERT_PTR_NON_NULL(cell->parent->parent, NULL); 652 PS_ASSERT_PTR_NON_NULL(cell->parent->parent->projection, NULL); 653 653 if (cell->toSky) { 654 654 // XXX: Should we use toTP or toSky? … … 690 690 const psGrommit* grommit) 691 691 { 692 PS_ PTR_CHECK_NULL(in, NULL);693 PS_ PTR_CHECK_NULL(grommit, NULL);692 PS_ASSERT_PTR_NON_NULL(in, NULL); 693 PS_ASSERT_PTR_NON_NULL(grommit, NULL); 694 694 695 695 char* type = "RA"; … … 711 711 const psFPA* fpa) 712 712 { 713 PS_ PTR_CHECK_NULL(tpCoord, NULL);714 PS_ PTR_CHECK_NULL(fpa, NULL);715 PS_ PTR_CHECK_NULL(fpa->fromTangentPlane, NULL);713 PS_ASSERT_PTR_NON_NULL(tpCoord, NULL); 714 PS_ASSERT_PTR_NON_NULL(fpa, NULL); 715 PS_ASSERT_PTR_NON_NULL(fpa->fromTangentPlane, NULL); 716 716 717 717 return (psPlaneDistortApply(fpaCoord, fpa->fromTangentPlane, … … 723 723 const psChip* chip) 724 724 { 725 PS_ PTR_CHECK_NULL(fpaCoord, NULL);726 PS_ PTR_CHECK_NULL(chip, NULL);727 PS_ PTR_CHECK_NULL(chip->fromFPA, NULL);725 PS_ASSERT_PTR_NON_NULL(fpaCoord, NULL); 726 PS_ASSERT_PTR_NON_NULL(chip, NULL); 727 PS_ASSERT_PTR_NON_NULL(chip->fromFPA, NULL); 728 728 729 729 chipCoord = psPlaneTransformApply(chipCoord, chip->fromFPA, fpaCoord); … … 735 735 const psCell* cell) 736 736 { 737 PS_ PTR_CHECK_NULL(chipCoord, NULL);738 PS_ PTR_CHECK_NULL(cell, NULL);739 PS_ PTR_CHECK_NULL(cell->fromChip, NULL);737 PS_ASSERT_PTR_NON_NULL(chipCoord, NULL); 738 PS_ASSERT_PTR_NON_NULL(cell, NULL); 739 PS_ASSERT_PTR_NON_NULL(cell->fromChip, NULL); 740 740 741 741 cellCoord = psPlaneTransformApply(cellCoord, cell->fromChip, chipCoord); … … 749 749 const psCell* cell) 750 750 { 751 PS_ PTR_CHECK_NULL(skyCoord, NULL);752 PS_ PTR_CHECK_NULL(cell, NULL);753 PS_ PTR_CHECK_NULL(cell->parent, NULL);754 PS_ PTR_CHECK_NULL(cell->parent->parent, NULL);755 PS_ PTR_CHECK_NULL(cell->parent->parent->grommit, NULL);751 PS_ASSERT_PTR_NON_NULL(skyCoord, NULL); 752 PS_ASSERT_PTR_NON_NULL(cell, NULL); 753 PS_ASSERT_PTR_NON_NULL(cell->parent, NULL); 754 PS_ASSERT_PTR_NON_NULL(cell->parent->parent, NULL); 755 PS_ASSERT_PTR_NON_NULL(cell->parent->parent->grommit, NULL); 756 756 757 757 psChip *parChip = cell->parent; … … 783 783 const psCell* cell) 784 784 { 785 PS_ PTR_CHECK_NULL(skyCoord, NULL);786 PS_ PTR_CHECK_NULL(cell, NULL);787 PS_ PTR_CHECK_NULL(cell->parent, NULL);788 PS_ PTR_CHECK_NULL(cell->parent->parent, NULL);789 PS_ PTR_CHECK_NULL(cell->parent->parent->projection, NULL);790 PS_ PTR_CHECK_NULL(cell->toTP, NULL);785 PS_ASSERT_PTR_NON_NULL(skyCoord, NULL); 786 PS_ASSERT_PTR_NON_NULL(cell, NULL); 787 PS_ASSERT_PTR_NON_NULL(cell->parent, NULL); 788 PS_ASSERT_PTR_NON_NULL(cell->parent->parent, NULL); 789 PS_ASSERT_PTR_NON_NULL(cell->parent->parent->projection, NULL); 790 PS_ASSERT_PTR_NON_NULL(cell->toTP, NULL); 791 791 if (cell->toSky) { 792 792 // XXX: Should we use toTP or toSky? -
trunk/psLib/src/astronomy/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/astronomy/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.
