Changeset 9539
- Timestamp:
- Oct 13, 2006, 11:22:09 AM (20 years ago)
- Location:
- trunk
- Files:
-
- 24 edited
-
ippTools/src/chiptool.c (modified) (2 diffs)
-
ippTools/src/dettool.c (modified) (1 diff)
-
ippTools/src/pzgetexp.c (modified) (1 diff)
-
ippTools/src/pzgetimfiles.c (modified) (1 diff)
-
ippTools/src/pztool.c (modified) (1 diff)
-
ppFringe/src/ppFringeLoop.c (modified) (2 diffs)
-
ppFringe/src/ppFringeSetup.c (modified) (1 diff)
-
ppImage/src/ppImageLoop.c (modified) (1 diff)
-
ppImage/src/ppImageOptions.c (modified) (1 diff)
-
ppImage/src/ppTest.c (modified) (1 diff)
-
ppMerge/src/ppMergeOptions.c (modified) (1 diff)
-
ppStats/src/ppStatsLoop.c (modified) (2 diffs)
-
ppStats/src/ppStatsSetupFromRecipe.c (modified) (1 diff)
-
psModules/src/camera/pmFPAConstruct.c (modified) (7 diffs)
-
psModules/src/camera/pmFPAfileDefine.c (modified) (1 diff)
-
psModules/src/camera/pmHDUGenerate.c (modified) (1 diff)
-
psModules/src/concepts/pmConceptsRead.c (modified) (4 diffs)
-
psModules/src/concepts/pmConceptsStandard.c (modified) (7 diffs)
-
psModules/src/concepts/pmConceptsWrite.c (modified) (3 diffs)
-
psModules/src/config/pmConfig.c (modified) (11 diffs)
-
psModules/src/objects/pmPSF_IO.c (modified) (8 diffs)
-
psModules/src/objects/pmPSF_IO.h (modified) (2 diffs)
-
psModules/test/config/tst_pmConfig.c (modified) (2 diffs)
-
psphot/src/psphotModelTest.c (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ippTools/src/chiptool.c
r9392 r9539 113 113 if (strcmp(rawFrame->exposure->exp_tag, 114 114 pendingFrame->exposure->exp_tag) == 0) { 115 psArrayRemove (rawFrames, rawFrame);115 psArrayRemoveData(rawFrames, rawFrame); 116 116 // dec the counter as the array just got shorter 117 117 // and we don't want to skip elemnts … … 132 132 if (strcmp(rawFrame->exposure->exp_tag, 133 133 doneFrame->exposure->exp_tag) == 0) { 134 psArrayRemove (rawFrames, rawFrame);134 psArrayRemoveData(rawFrames, rawFrame); 135 135 // dec the counter as the array just got shorter 136 136 // and we don't want to skip elemnts -
trunk/ippTools/src/dettool.c
r9494 r9539 387 387 return false; 388 388 } 389 psMetadataRemove (md, 0, "det_id");389 psMetadataRemoveKey(md, "det_id"); 390 390 psString det_idStr = psDBIntToString((psU64)det_id); 391 391 psMetadataAddStr(mds->data[i], PS_LIST_HEAD, "det_id", 0, NULL, det_idStr); -
trunk/ippTools/src/pzgetexp.c
r9392 r9539 94 94 if (strcmp(newSummitExp->exp_id, 95 95 summitExp->exp_id) == 0) { 96 psArrayRemove (newSummitExps, newSummitExp);96 psArrayRemoveData(newSummitExps, newSummitExp); 97 97 // dec the counter as the array just got shorter 98 98 // and we don't want to skip elemnts -
trunk/ippTools/src/pzgetimfiles.c
r9392 r9539 83 83 if (strcmp(newImfile->exp_id, 84 84 pendingImfile->exp_id) == 0) { 85 psArrayRemove (newImfiles, newImfile);85 psArrayRemoveData(newImfiles, newImfile); 86 86 // dec the counter as the array just got shorter 87 87 // and we don't want to skip elemnts -
trunk/ippTools/src/pztool.c
r9392 r9539 97 97 newExpRow *newExp = new->data[j]; 98 98 if (strcmp(summitExp->exp_id, newExp->exp_id) == 0) { 99 psArrayRemove (summit, summitExp);99 psArrayRemoveData(summit, summitExp); 100 100 // dec the counter as the array just got shorter 101 101 //and we don't want to skip elemnts -
trunk/ppFringe/src/ppFringeLoop.c
r9069 r9539 81 81 if (contents && contents->list->n > 0) { 82 82 bool mdok; // Status of MD lookup 83 psMetadata *header = psMetadataLookupM D(&mdok, contents, extname); // Header83 psMetadata *header = psMetadataLookupMetadata(&mdok, contents, extname); // Header 84 84 if (mdok && header) { 85 85 regions = pmFringeRegionsReadFits(NULL, data->outFile, extname); … … 141 141 psString checkName = NULL; // Extension name to check 142 142 psStringAppend(&checkName, "%s:%d", extname, version); 143 header = psMetadataLookupM D(&mdok, contents, checkName);143 header = psMetadataLookupMetadata(&mdok, contents, checkName); 144 144 psFree(checkName); 145 145 } while (mdok && header); -
trunk/ppFringe/src/ppFringeSetup.c
r8785 r9539 77 77 // Now we have a camera format, we can read the recipe 78 78 bool mdok; // Status of MD lookup 79 psMetadata *recipe = psMetadataLookupM D(&mdok, config->recipes, RECIPE); // The recipe79 psMetadata *recipe = psMetadataLookupMetadata(&mdok, config->recipes, RECIPE); // The recipe 80 80 if (!recipe) { 81 81 psError(PS_ERR_IO, true, "Unable to find recipe %s.\n", RECIPE); -
trunk/ppImage/src/ppImageLoop.c
r8751 r9539 14 14 15 15 bool mdok; // Status of MD lookup 16 psMetadata *extraOpts = psMetadataLookupM D(&mdok, config->arguments, "PPIMAGE.OPTIONS"); // Extra options16 psMetadata *extraOpts = psMetadataLookupMetadata(&mdok, config->arguments, "PPIMAGE.OPTIONS"); // Extra options 17 17 const char *statsName = psMetadataLookupStr(&mdok, extraOpts, "STATS"); // Filename for statistics 18 18 psMetadata *stats = NULL; // Container for statistics -
trunk/ppImage/src/ppImageOptions.c
r9342 r9539 65 65 // select the recipe for this analysis 66 66 bool mdStatus = false; // Result of MD lookup 67 psMetadata *recipe = psMetadataLookupM D(&mdStatus, config->recipes, RECIPE_NAME);67 psMetadata *recipe = psMetadataLookupMetadata(&mdStatus, config->recipes, RECIPE_NAME); 68 68 if (! mdStatus || !recipe) { 69 69 psLogMsg("ppImage", PS_LOG_ERROR, "Can't find recipe %s in the RECIPES.\n", RECIPE_NAME); -
trunk/ppImage/src/ppTest.c
r8822 r9539 100 100 ////////////////////////////////////////////////////////////////////////////////////////////////////////////// 101 101 102 pmFPAPrint(stdout, fpa, false, true); 103 exit(0); 104 102 105 103 106 #if 0 -
trunk/ppMerge/src/ppMergeOptions.c
r9463 r9539 115 115 // Now we can read the recipe 116 116 bool mdok = true; // Status of MD lookup 117 psMetadata *recipe = psMetadataLookupM D(&mdok, config->recipes, PPMERGE_RECIPE); // Recipe information117 psMetadata *recipe = psMetadataLookupMetadata(&mdok, config->recipes, PPMERGE_RECIPE); // Recipe information 118 118 if (!mdok || !recipe) { 119 119 psError(PS_ERR_IO, true, "Unable to find recipe %s", PPMERGE_RECIPE); -
trunk/ppStats/src/ppStatsLoop.c
r8749 r9539 89 89 // Cell-level results 90 90 bool mdok; // Status of MD lookup 91 psMetadata *cellResults = psMemIncrRefCounter(psMetadataLookupM D(&mdok, chipResults, cellName));91 psMetadata *cellResults = psMemIncrRefCounter(psMetadataLookupMetadata(&mdok, chipResults, cellName)); 92 92 if (!mdok || !cellResults) { 93 93 cellResults = psMetadataAlloc(); … … 242 242 // Chip-level results 243 243 bool mdok; // Status of MD lookup 244 psMetadata *chipResults = psMemIncrRefCounter(psMetadataLookupM D(&mdok, fpaResults, chipName));244 psMetadata *chipResults = psMemIncrRefCounter(psMetadataLookupMetadata(&mdok, fpaResults, chipName)); 245 245 if (!mdok || !chipResults) { 246 246 chipResults = psMetadataAlloc(); -
trunk/ppStats/src/ppStatsSetupFromRecipe.c
r8749 r9539 69 69 // Determine recipe parameters 70 70 bool mdok; // Status of MD lookup 71 psMetadata *recipe = psMetadataLookupM D(&mdok, config->recipes, PPSTATS_RECIPE);71 psMetadata *recipe = psMetadataLookupMetadata(&mdok, config->recipes, PPSTATS_RECIPE); 72 72 if (!mdok || !recipe) { 73 73 psLogMsg(__func__, PS_LOG_WARN, "Unable to find recipe %s.\n", PPSTATS_RECIPE); -
trunk/psModules/src/camera/pmFPAConstruct.c
r9514 r9539 28 28 29 29 bool status = true; // Result of MD lookup 30 psMetadata *cells = psMetadataLookupM D(&status, format, "CELLS"); // The CELLS30 psMetadata *cells = psMetadataLookupMetadata(&status, format, "CELLS"); // The CELLS 31 31 if (!status || !cells) { 32 32 psError(PS_ERR_IO, true, "Unable to determine CELLS of camera.\n"); … … 34 34 } 35 35 36 psMetadata *cellData = psMetadataLookupM D(&status, cells, cellName); // The data for the particular cell36 psMetadata *cellData = psMetadataLookupMetadata(&status, cells, cellName); // The data for the particular cell 37 37 if (!status || !cellData) { 38 38 psLogMsg(__func__, PS_LOG_WARN, "Unable to find specs for cell %s: ignored\n", cellName); … … 389 389 390 390 bool mdok = true; // Status of MD lookup 391 psMetadata *file = psMetadataLookupM D(&mdok, format, "FILE"); // File information391 psMetadata *file = psMetadataLookupMetadata(&mdok, format, "FILE"); // File information 392 392 if (!mdok || !file) { 393 393 psError(PS_ERR_IO, true, "Unable to find FILE information in camera format configuration.\n"); … … 427 427 428 428 bool mdok = true; // Status from MD lookups 429 psMetadata *components = psMetadataLookupM D(&mdok, camera, "FPA"); // FPA components429 psMetadata *components = psMetadataLookupMetadata(&mdok, camera, "FPA"); // FPA components 430 430 if (!mdok || !components) { 431 431 psError(PS_ERR_IO, true, "Failed to lookup \"FPA\""); … … 475 475 476 476 bool mdok = true; // Status from metadata lookups 477 psMetadata *fileInfo = psMetadataLookupM D(&mdok, format, "FILE"); // The file information477 psMetadata *fileInfo = psMetadataLookupMetadata(&mdok, format, "FILE"); // The file information 478 478 if (!mdok || !fileInfo) { 479 479 psError(PS_ERR_IO, false, "Unable to find FILE in the camera format configuration.\n"); … … 553 553 // In cases 2 and 3, the CONTENTS is of type METADATA, and is either a menu (if EXTENSIONS=NONE), or a 554 554 // list of extensions otherwise. 555 psMetadata *contents = psMetadataLookupM D(&mdok, format, "CONTENTS"); // The contents of the FITS file555 psMetadata *contents = psMetadataLookupMetadata(&mdok, format, "CONTENTS"); // The contents of the FITS file 556 556 if (!mdok || !contents) { 557 557 if (mdok && !contents) { … … 775 775 776 776 bool mdok = true; // Status from metadata lookups 777 psMetadata *fileInfo = psMetadataLookupM D(&mdok, format, "FILE"); // The file information777 psMetadata *fileInfo = psMetadataLookupMetadata(&mdok, format, "FILE"); // The file information 778 778 if (!mdok || !fileInfo) { 779 779 psError(PS_ERR_IO, false, "Unable to find FILE in the camera format configuration.\n"); -
trunk/psModules/src/camera/pmFPAfileDefine.c
r9518 r9539 256 256 257 257 // select the format list from the selected camera 258 psMetadata *formats = psMetadataLookupM D(&status, file->camera, "FORMATS"); // List of formats258 psMetadata *formats = psMetadataLookupMetadata(&status, file->camera, "FORMATS"); // List of formats 259 259 260 260 // select the desired format -
trunk/psModules/src/camera/pmHDUGenerate.c
r8848 r9539 389 389 } 390 390 bool mdok = true; // Status of MD lookup 391 psMetadata *file = psMetadataLookupM D(&mdok, hdu->format, "FILE"); // File information for camera format391 psMetadata *file = psMetadataLookupMetadata(&mdok, hdu->format, "FILE"); // File information for camera format 392 392 if (!mdok || !file) { 393 393 psError(PS_ERR_UNEXPECTED_NULL, true, "Can't file FILE information for camera format " -
trunk/psModules/src/concepts/pmConceptsRead.c
r8848 r9539 191 191 psMetadata *cameraFormat = hdu->format; // The camera format 192 192 bool mdok = true; // Status of MD lookup 193 psMetadata *defaults = psMetadataLookupM D(&mdok, cameraFormat, "DEFAULTS"); // The DEFAULTS spec193 psMetadata *defaults = psMetadataLookupMetadata(&mdok, cameraFormat, "DEFAULTS"); // The DEFAULTS spec 194 194 if (!mdok || !defaults) { 195 195 psError(PS_ERR_IO, true, "Failed to find \"DEFAULTS\""); … … 279 279 psMetadata *cameraFormat = hduLow->format; // The camera format 280 280 bool mdok = true; // Status of MD lookup 281 psMetadata *transSpec = psMetadataLookupM D(&mdok, cameraFormat, "TRANSLATION"); // The TRANSLATION spec281 psMetadata *transSpec = psMetadataLookupMetadata(&mdok, cameraFormat, "TRANSLATION"); // The TRANSLATION spec 282 282 if (!mdok || !transSpec) { 283 283 psError(PS_ERR_IO, true, "Failed to find \"TRANSLATION\""); … … 381 381 psMetadata *cameraFormat = hdu->format; // The camera format 382 382 bool mdok = true; // Status of MD lookup 383 psMetadata *dbSpec = psMetadataLookupM D(&mdok, cameraFormat, "DATABASE"); // The DATABASE spec383 psMetadata *dbSpec = psMetadataLookupMetadata(&mdok, cameraFormat, "DATABASE"); // The DATABASE spec 384 384 if (!mdok || !dbSpec) { 385 385 psError(PS_ERR_IO, true, "Unable to find DATABASE in camera configuration.\n"); … … 394 394 psString name = specItem->name; // The concept name 395 395 396 psMetadata *dbLookup = psMetadataLookupM D(&mdok, dbSpec, name);396 psMetadata *dbLookup = psMetadataLookupMetadata(&mdok, dbSpec, name); 397 397 if (mdok && dbLookup) { 398 398 const char *tableName = psMetadataLookupStr(&mdok, dbLookup, "TABLE"); // Table name -
trunk/psModules/src/concepts/pmConceptsStandard.c
r9510 r9539 91 91 // How to interpret the coordinates 92 92 bool mdok = true; // Status of MD lookup 93 psMetadata *formats = psMetadataLookupM D(&mdok, cameraFormat, "FORMATS");93 psMetadata *formats = psMetadataLookupMetadata(&mdok, cameraFormat, "FORMATS"); 94 94 if (mdok && formats) { 95 95 psString format = psMetadataLookupStr(&mdok, formats, pattern->name); … … 128 128 // How to interpret the coordinates 129 129 bool mdok = true; // Status of MD lookup 130 psMetadata *formats = psMetadataLookupM D(&mdok,130 psMetadata *formats = psMetadataLookupMetadata(&mdok, 131 131 cameraFormat, 132 132 "FORMATS"); … … 350 350 351 351 // Get format 352 psMetadata *formats = psMetadataLookupM D(&mdok, cameraFormat, "FORMATS");352 psMetadata *formats = psMetadataLookupMetadata(&mdok, cameraFormat, "FORMATS"); 353 353 if (!mdok || !formats) { 354 354 psError(PS_ERR_UNKNOWN, true, "Unable to find FORMATS in camera configuration.\n"); … … 528 528 { 529 529 bool mdok = false; // Result of MD lookup 530 psMetadata *formats = psMetadataLookupM D(&mdok, cameraFormat, "FORMATS");530 psMetadata *formats = psMetadataLookupMetadata(&mdok, cameraFormat, "FORMATS"); 531 531 if (mdok && formats) { 532 532 psString format = psMetadataLookupStr(&mdok, formats, name); … … 622 622 assert(concept); 623 623 624 psMetadata *translation = psMetadataLookupM D(NULL, cameraFormat, "TRANSLATION");624 psMetadata *translation = psMetadataLookupMetadata(NULL, cameraFormat, "TRANSLATION"); 625 625 bool xBinOK = true, yBinOK = true; // Status of MD lookups 626 626 psString xKeyword = psMetadataLookupStr(&xBinOK, translation, "CELL.XBIN"); … … 649 649 assert(concept); 650 650 651 psMetadata *translation = psMetadataLookupM D(NULL, cameraFormat, "TRANSLATION");651 psMetadata *translation = psMetadataLookupMetadata(NULL, cameraFormat, "TRANSLATION"); 652 652 bool xBinOK = true, yBinOK = true; // Status of MD lookups 653 653 psString xKeyword = psMetadataLookupStr(&xBinOK, translation, "CELL.XBIN"); … … 710 710 711 711 bool mdok = true; // Status of MD lookup 712 psMetadata *formats = psMetadataLookupM D(&mdok, cameraFormat, "FORMATS"); // The formats712 psMetadata *formats = psMetadataLookupMetadata(&mdok, cameraFormat, "FORMATS"); // The formats 713 713 if (mdok && formats) { 714 714 psString format = psMetadataLookupStr(&mdok, formats, concept->name); // The formats for eg, CELL.TIME -
trunk/psModules/src/concepts/pmConceptsWrite.c
r8848 r9539 307 307 psMetadata *cameraFormat = hdu->format; // The camera format 308 308 bool mdok = true; // Status of MD lookup 309 psMetadata *defaults = psMetadataLookupM D(&mdok, cameraFormat, "DEFAULTS"); // The DEFAULTS spec309 psMetadata *defaults = psMetadataLookupMetadata(&mdok, cameraFormat, "DEFAULTS"); // The DEFAULTS spec 310 310 if (!mdok || !defaults) { 311 311 return false; … … 373 373 psMetadata *cameraFormat = hdu->format; // The camera format 374 374 bool mdok = true; // Status of MD lookup 375 psMetadata *translation = psMetadataLookupM D(&mdok, cameraFormat, "TRANSLATION"); // The TRANSLATION spec375 psMetadata *translation = psMetadataLookupMetadata(&mdok, cameraFormat, "TRANSLATION"); // The TRANSLATION spec 376 376 if (!mdok || !translation) { 377 377 return false; … … 430 430 psMetadata *cameraFormat = hdu->format; // The camera format 431 431 bool mdok = true; // Status of MD lookup 432 psMetadata *database = psMetadataLookupM D(&mdok, cameraFormat, "DATABASE"); // The DATABASE spec432 psMetadata *database = psMetadataLookupMetadata(&mdok, cameraFormat, "DATABASE"); // The DATABASE spec 433 433 if (!mdok || !database) { 434 434 return false; -
trunk/psModules/src/config/pmConfig.c
r9435 r9539 3 3 * @author PAP, IfA 4 4 * 5 * @version $Revision: 1.4 6$ $Name: not supported by cvs2svn $6 * @date $Date: 2006-10-1 0 01:03:13$5 * @version $Revision: 1.47 $ $Name: not supported by cvs2svn $ 6 * @date $Date: 2006-10-13 21:15:45 $ 7 7 * 8 8 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 277 277 // filename 278 278 bool mdok = true; // Status of MD lookup 279 psMetadata *cameras = psMetadataLookupM D(&mdok, config->site, "CAMERAS");279 psMetadata *cameras = psMetadataLookupMetadata(&mdok, config->site, "CAMERAS"); 280 280 if (mdok && cameras) { 281 281 char *newFile = psMetadataLookupStr(NULL, cameras, cameraFile); // The filename … … 377 377 // XXX: This is not spec'ed in the SDRS. 378 378 // 379 psMetadata *trace = psMetadataLookupM D(&mdok, config->site, "TRACE");379 psMetadata *trace = psMetadataLookupMetadata(&mdok, config->site, "TRACE"); 380 380 if (mdok && trace) { 381 381 psMetadataIterator *traceIter = psMetadataIteratorAlloc(trace, PS_LIST_HEAD, NULL); // Iterator … … 446 446 // Read the rule for that camera format 447 447 bool mdStatus = true; 448 psMetadata *rule = psMetadataLookupM D(&mdStatus, cameraFormat, "RULE");448 psMetadata *rule = psMetadataLookupMetadata(&mdStatus, cameraFormat, "RULE"); 449 449 if (! mdStatus || ! rule) { 450 450 psLogMsg(__func__, PS_LOG_WARN, "Unable to read rule for camera.\n"); … … 493 493 // Read the list of formats 494 494 bool mdok = true; // Status of MD lookup 495 psMetadata *formats = psMetadataLookupM D(&mdok, camera, "FORMATS"); // List of formats495 psMetadata *formats = psMetadataLookupMetadata(&mdok, camera, "FORMATS"); // List of formats 496 496 if (!mdok || !formats) { 497 497 psLogMsg(__func__, PS_LOG_WARN, "Unable to find list of FORMATS in camera %s --- ignored\n", … … 550 550 // If we don't know what sort of camera we have, we try all that we know 551 551 if (! config->camera) { 552 psMetadata *cameras = psMetadataLookupM D(&mdok, config->site, "CAMERAS");552 psMetadata *cameras = psMetadataLookupMetadata(&mdok, config->site, "CAMERAS"); 553 553 if (! mdok) { 554 554 psError(PS_ERR_IO, true, "Unable to find CAMERAS in the configuration."); … … 615 615 PS_ASSERT_PTR_NON_NULL(config, NULL); 616 616 617 psMetadata *cameras = psMetadataLookupM D(NULL, config->site, "CAMERAS");617 psMetadata *cameras = psMetadataLookupMetadata(NULL, config->site, "CAMERAS"); 618 618 if (!cameras) { 619 619 psError(PS_ERR_IO, true, "Unable to find CAMERAS in the configuration."); … … 659 659 660 660 bool mdok = true; // Status of MD lookup 661 psMetadata *recipes = psMetadataLookupM D(&mdok, source, "RECIPES"); // The list of recipes661 psMetadata *recipes = psMetadataLookupMetadata(&mdok, source, "RECIPES"); // The list of recipes 662 662 if (!mdok || !recipes) { 663 663 psLogMsg(__func__, PS_LOG_WARN, "RECIPES in the %s is not of type METADATA --- ignored\n", … … 720 720 if (sourceItem->data.S32 == P_PM_RECIPE_SOURCE_SYMBOLIC) { 721 721 const char *linkName = sourceItem->comment; // The name of the link 722 psMetadata *linkSource = psMetadataLookupM D(&mdok, config->recipes, linkName); // The source722 psMetadata *linkSource = psMetadataLookupMetadata(&mdok, config->recipes, linkName); // The source 723 723 if (!mdok || !linkSource) { 724 724 // Can't yet resolve it … … 794 794 795 795 // Is the source a symbolic reference? 796 psMetadata *extant = psMetadataLookupM D(&mdok, config->recipes, recipeSource); // Does it exist?796 psMetadata *extant = psMetadataLookupMetadata(&mdok, config->recipes, recipeSource); // Does it exist? 797 797 if (mdok && extant) { 798 798 psString comment = psStringCopy("Recipe added from command line as symbolic link to "); … … 889 889 890 890 bool mdok = true; // Status of MD lookup 891 psMetadata *rules = psMetadataLookupM D(&mdok, format, "RULE"); // How to identify this format891 psMetadata *rules = psMetadataLookupMetadata(&mdok, format, "RULE"); // How to identify this format 892 892 if (!mdok || !rules) { 893 893 psError(PS_ERR_IO, true, "Unable to find RULE in camera format.\n"); -
trunk/psModules/src/objects/pmPSF_IO.c
r8815 r9539 6 6 * @author EAM, IfA 7 7 * 8 * @version $Revision: 1. 5$ $Name: not supported by cvs2svn $9 * @date $Date: 2006- 09-15 09:49:01$8 * @version $Revision: 1.6 $ $Name: not supported by cvs2svn $ 9 * @date $Date: 2006-10-13 21:15:45 $ 10 10 * 11 11 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 37 37 #include "pmSourceIO.h" 38 38 39 psMetadata *pmPSFtoM D(psMetadata *metadata, pmPSF *psf)39 psMetadata *pmPSFtoMetadata (psMetadata *metadata, pmPSF *psf) 40 40 { 41 41 … … 54 54 for (int i = 0; i < nPar - 4; i++) { 55 55 psPolynomial2D *poly = psf->params->data[i]; 56 psPolynomial2DtoM D(metadata, poly, "PSF_PAR%02d", i);57 } 58 psPolynomial4DtoM D(metadata, psf->ApTrend, "APTREND");56 psPolynomial2DtoMetadata (metadata, poly, "PSF_PAR%02d", i); 57 } 58 psPolynomial4DtoMetadata (metadata, psf->ApTrend, "APTREND"); 59 59 60 60 psMetadataAdd (metadata, PS_LIST_TAIL, "PSF_AP_RESID", PS_DATA_F32, "aperture residual", psf->ApResid); … … 68 68 } 69 69 70 pmPSF *pmPSFfromM D(psMetadata *metadata)70 pmPSF *pmPSFfromMetadata (psMetadata *metadata) 71 71 { 72 72 … … 90 90 sprintf (keyword, "PSF_PAR%02d", i); 91 91 psMetadata *folder = psMetadataLookupPtr (&status, metadata, keyword); 92 psPolynomial2D *poly = psPolynomial2DfromM D(folder);92 psPolynomial2D *poly = psPolynomial2DfromMetadata (folder); 93 93 psFree (psf->params->data[i]); 94 94 psf->params->data[i] = poly; … … 96 96 sprintf (keyword, "APTREND"); 97 97 psMetadata *folder = psMetadataLookupPtr (&status, metadata, keyword); 98 psPolynomial4D *poly = psPolynomial4DfromM D(folder);98 psPolynomial4D *poly = psPolynomial4DfromMetadata (folder); 99 99 psFree (psf->ApTrend); 100 100 psf->ApTrend = poly; … … 197 197 case PM_FPA_FILE_PSF: 198 198 filename = pmFPAfileNameFromRule (file->filerule, file, view); 199 psMetadata *psfData = pmPSFtoM D(NULL, psf);199 psMetadata *psfData = pmPSFtoMetadata (NULL, psf); 200 200 psMetadataConfigWrite (psfData, filename); 201 201 psFree (psfData); … … 299 299 300 300 psMetadata *psfData = psMetadataConfigRead(NULL, &Nfail, filename, FALSE); 301 pmPSF *psf = pmPSFfromM D(psfData);301 pmPSF *psf = pmPSFfromMetadata (psfData); 302 302 psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PSPHOT.PSF", PS_DATA_UNKNOWN, "psphot psf", psf); 303 303 -
trunk/psModules/src/objects/pmPSF_IO.h
r6872 r9539 6 6 * @author EAM, IfA 7 7 * 8 * @version $Revision: 1. 2$ $Name: not supported by cvs2svn $9 * @date $Date: 2006- 04-17 18:01:05 $8 * @version $Revision: 1.3 $ $Name: not supported by cvs2svn $ 9 * @date $Date: 2006-10-13 21:15:45 $ 10 10 * 11 11 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 16 16 # define PM_PSF_IO_H 17 17 18 psMetadata *pmPSFtoM D(psMetadata *metadata, pmPSF *psf);19 pmPSF *pmPSFfromM D(psMetadata *metadata);18 psMetadata *pmPSFtoMetadata (psMetadata *metadata, pmPSF *psf); 19 pmPSF *pmPSFfromMetadata (psMetadata *metadata); 20 20 bool pmFPAviewWritePSFmodel (const pmFPAview *view, pmFPAfile *file); 21 21 bool pmFPAWritePSFmodel (pmFPA *fpa, const pmFPAview *view, pmFPAfile *file); -
trunk/psModules/test/config/tst_pmConfig.c
r6414 r9539 15 15 * XXXX: Must determine what to do with NULL arguments, then test it. 16 16 * 17 * @version $Revision: 1. 7$ $Name: not supported by cvs2svn $18 * @date $Date: 2006- 02-10 01:41:37$17 * @version $Revision: 1.8 $ $Name: not supported by cvs2svn $ 18 * @date $Date: 2006-10-13 21:15:45 $ 19 19 * 20 20 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 210 210 // Test the database camera metadata keywords. 211 211 // 212 psMetadata *tmpMeta = psMetadataLookupM D(&rc, site, "CAMERAS");212 psMetadata *tmpMeta = psMetadataLookupMetadata(&rc, site, "CAMERAS"); 213 213 if (rc == false) { 214 214 printf("TEST ERROR: failed to properly set metadata metadata CAMERAS.\n"); -
trunk/psphot/src/psphotModelTest.c
r9270 r9539 25 25 char *fitModeWord = psMetadataLookupStr (&status, recipe, "TEST_FIT_MODE"); 26 26 if (!status) { 27 fitModeWord = DEFAULT_MODE;27 fitModeWord = DEFAULT_MODE; 28 28 } 29 29 fitMode = PM_SOURCE_FIT_EXT; … … 32 32 // in fitMode, psf sets the model type 33 33 if (fitMode == PM_SOURCE_FIT_PSF) { 34 char *psfFile = psMetadataLookupStr (&status, recipe, "PSF_INPUT_FILE");35 if (!status) psAbort ("psphotModelTest", "PSF_INPUT_FILE not supplied");36 psMetadata *psfData = psMetadataConfigRead(NULL, &Nfail, psfFile, FALSE);37 psf = pmPSFfromMD(psfData);38 modelType = psf->type;34 char *psfFile = psMetadataLookupStr (&status, recipe, "PSF_INPUT_FILE"); 35 if (!status) psAbort ("psphotModelTest", "PSF_INPUT_FILE not supplied"); 36 psMetadata *psfData = psMetadataConfigRead(NULL, &Nfail, psfFile, FALSE); 37 psf = pmPSFfromMetadata (psfData); 38 modelType = psf->type; 39 39 } else { 40 // find the model: supplied by user or first in the PSF_MODEL list41 char *modelName = psMetadataLookupStr (&status, recipe, "TEST_FIT_MODEL");42 if (modelName == NULL) {43 // get the list pointers for the PSF_MODEL entries40 // find the model: supplied by user or first in the PSF_MODEL list 41 char *modelName = psMetadataLookupStr (&status, recipe, "TEST_FIT_MODEL"); 42 if (modelName == NULL) { 43 // get the list pointers for the PSF_MODEL entries 44 44 45 psList *list = NULL;46 psMetadataItem *mdi = psMetadataLookup (recipe, "PSF_MODEL");47 if (mdi == NULL) psAbort ("psphotModelTest", "missing PSF_MODEL selection");48 if (mdi->type == PS_DATA_STRING) {49 list = psListAlloc(NULL);50 psListAdd (list, PS_LIST_HEAD, mdi);51 } else {52 if (mdi->type != PS_DATA_METADATA_MULTI) psAbort ("psphotChoosePSF", "missing PSF_MODEL selection");53 list = psMemIncrRefCounter(mdi->data.list);54 }45 psList *list = NULL; 46 psMetadataItem *mdi = psMetadataLookup (recipe, "PSF_MODEL"); 47 if (mdi == NULL) psAbort ("psphotModelTest", "missing PSF_MODEL selection"); 48 if (mdi->type == PS_DATA_STRING) { 49 list = psListAlloc(NULL); 50 psListAdd (list, PS_LIST_HEAD, mdi); 51 } else { 52 if (mdi->type != PS_DATA_METADATA_MULTI) psAbort ("psphotChoosePSF", "missing PSF_MODEL selection"); 53 list = psMemIncrRefCounter(mdi->data.list); 54 } 55 55 56 // take the first list element57 item = psListGet (list, PS_LIST_HEAD);58 modelName = item->data.V;59 }60 modelType = pmModelSetType (modelName);61 if (modelType < 0) psAbort ("fitsource", "unknown model %s", modelName);56 // take the first list element 57 item = psListGet (list, PS_LIST_HEAD); 58 modelName = item->data.V; 59 } 60 modelType = pmModelSetType (modelName); 61 if (modelType < 0) psAbort ("fitsource", "unknown model %s", modelName); 62 62 } 63 63 … … 65 65 float INNER = psMetadataLookupF32 (&status, recipe, "TEST_FIT_INNER_RADIUS"); 66 66 if (!status) { 67 INNER = psMetadataLookupF32 (&status, recipe, "SKY_INNER_RADIUS");67 INNER = psMetadataLookupF32 (&status, recipe, "SKY_INNER_RADIUS"); 68 68 } 69 69 70 70 float OUTER = psMetadataLookupF32 (&status, recipe, "TEST_FIT_OUTER_RADIUS"); 71 71 if (!status) { 72 OUTER = psMetadataLookupF32 (&status, recipe, "SKY_OUTER_RADIUS");72 OUTER = psMetadataLookupF32 (&status, recipe, "SKY_OUTER_RADIUS"); 73 73 } 74 74 75 75 float RADIUS = psMetadataLookupF32 (&status, recipe, "TEST_FIT_RADIUS"); 76 76 if (!status) { 77 RADIUS = psMetadataLookupF32 (&status, recipe, "PSF_FIT_RADIUS");77 RADIUS = psMetadataLookupF32 (&status, recipe, "PSF_FIT_RADIUS"); 78 78 } 79 79 80 80 float mRADIUS = psMetadataLookupF32 (&status, recipe, "TEST_MOMENTS_RADIUS"); 81 81 if (!status) { 82 mRADIUS = psMetadataLookupF32 (&status, recipe, "PSF_MOMENTS_RADIUS");82 mRADIUS = psMetadataLookupF32 (&status, recipe, "PSF_MOMENTS_RADIUS"); 83 83 } 84 84 … … 117 117 psF32 *params = model->params->data.F32; 118 118 for (int i = 0; i < nParams; i++) { 119 if (i == 2) {120 params[i] = xObj;121 continue;122 }123 if (i == 3) {124 params[i] = yObj;125 continue;126 }127 sprintf (name, "TEST_FIT_PAR%d", i);128 value = psMetadataLookupF32 (&status, recipe, name);129 if (status) {130 params[i] = value;131 }119 if (i == 2) { 120 params[i] = xObj; 121 continue; 122 } 123 if (i == 3) { 124 params[i] = yObj; 125 continue; 126 } 127 sprintf (name, "TEST_FIT_PAR%d", i); 128 value = psMetadataLookupF32 (&status, recipe, name); 129 if (status) { 130 params[i] = value; 131 } 132 132 } 133 133 … … 136 136 137 137 if (fitMode == PM_SOURCE_FIT_PSF) { 138 pmModel *modelPSF = pmModelFromPSF (model, psf);139 psFree (model);140 model = modelPSF;141 params = model->params->data.F32;138 pmModel *modelPSF = pmModelFromPSF (model, psf); 139 psFree (model); 140 model = modelPSF; 141 params = model->params->data.F32; 142 142 } 143 143 … … 154 154 fprintf (stderr, "input parameters: \n"); 155 155 for (int i = 0; i < nParams; i++) { 156 fprintf (stderr, "%d : %f\n", i, params[i]);156 fprintf (stderr, "%d : %f\n", i, params[i]); 157 157 } 158 158 … … 162 162 char *fitset = psMetadataLookupStr (&status, recipe, "TEST_FIT_SET"); 163 163 if (status) { 164 status = psphotFitSet (source, model, fitset, fitMode);165 exit (0);164 status = psphotFitSet (source, model, fitset, fitMode); 165 exit (0); 166 166 } 167 167 … … 178 178 // subtract object, leave local sky 179 179 pmModelSub (source->pixels, source->mask, model, false, false); 180 180 181 181 fprintf (stderr, "output parameters: \n"); 182 182 for (int i = 0; i < nParams; i++) { 183 fprintf (stderr, "%d : %f\n", i, params[i]);183 fprintf (stderr, "%d : %f\n", i, params[i]); 184 184 } 185 185 186 // write out 186 // write out 187 187 psphotSaveImage (NULL, source->pixels, "resid.fits"); 188 188 psphotSaveImage (NULL, source->mask, "mask.fits");
Note:
See TracChangeset
for help on using the changeset viewer.
