Changeset 4029 for trunk/psLib/src/astronomy/psAstrometry.c
- Timestamp:
- May 25, 2005, 10:26:55 AM (21 years ago)
- File:
-
- 1 edited
-
trunk/psLib/src/astronomy/psAstrometry.c (modified) (27 diffs)
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?
Note:
See TracChangeset
for help on using the changeset viewer.
