Changeset 37706 for branches/eam_branches/ipp-20140904
- Timestamp:
- Dec 2, 2014, 10:05:38 AM (12 years ago)
- Location:
- branches/eam_branches/ipp-20140904/Ohana/src/relastro
- Files:
-
- 1 added
- 4 edited
-
Makefile (modified) (1 diff)
-
include/relastro.h (modified) (1 diff)
-
src/BrightCatalog.c (modified) (2 diffs)
-
src/client_logger.c (added)
-
src/relastro_parallel_images.c (modified) (7 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20140904/Ohana/src/relastro/Makefile
r37684 r37706 89 89 $(SRC)/indexCatalogs.$(ARCH).o \ 90 90 $(SRC)/syncfile.$(ARCH).o \ 91 $(SRC)/client_logger.$(ARCH).o \ 91 92 $(SRC)/relastroVisual.$(ARCH).o 92 93 -
branches/eam_branches/ipp-20140904/Ohana/src/relastro/include/relastro.h
r37699 r37706 653 653 void lockUpdateChips (); 654 654 void unlockUpdateChips (); 655 656 int client_logger_init (char *dirname); 657 int client_logger_message (char *format,...); -
branches/eam_branches/ipp-20140904/Ohana/src/relastro/src/BrightCatalog.c
r37655 r37706 43 43 if (!gfits_fread_header (f, &header)) { 44 44 if (VERBOSE) fprintf (stderr, "can't read image subset header\n"); 45 fclose (f);46 45 goto escape; 47 46 } 48 47 if (!gfits_fread_matrix (f, &matrix, &header)) { 49 48 if (VERBOSE) fprintf (stderr, "can't read image subset matrix\n"); 50 gfits_free_header (&header);51 fclose (f);52 49 goto escape; 53 50 } … … 316 313 if (catalog) free (catalog); 317 314 318 fclose (f);315 if (f) fclose (f); 319 316 return NULL; 320 317 } -
branches/eam_branches/ipp-20140904/Ohana/src/relastro/src/relastro_parallel_images.c
r37703 r37706 17 17 INITTIME; 18 18 19 client_logger_init (CATDIR); 20 19 21 // load the RegionTable (UserRegion should not be used at this level) 20 22 RegionHostTable *regionHosts = RegionHostTableLoad (CATDIR, REGION_FILE); 21 23 int myHost = regionHosts->index[REGION_HOST_ID]; 22 24 RegionHostFindNeighbors (regionHosts, myHost); 25 client_logger_message ("started parallel images on %s\n", regionHosts->hosts[myHost].hostname); 23 26 24 27 // load the subset images belonging to this host … … 44 47 } 45 48 put_astrom_table (table); 49 client_logger_message ("loaded images on %s\n", regionHosts->hosts[myHost].hostname); 46 50 47 51 // once we have read this table, we should remove it for repeat runs … … 63 67 free (syncfile); 64 68 69 client_logger_message ("loaded catalog data %s\n", regionHosts->hosts[myHost].hostname); 70 65 71 // find ICRF QSOs for reference downstream (only if USE_ICRF_CORRECT) 66 72 select_catalog_ICRF (catalog, Ncatalog); … … 84 90 createStarMap (catalog, Ncatalog); 85 91 92 client_logger_message ("set up image indexes on %s\n", regionHosts->hosts[myHost].hostname); 93 86 94 markObjects (catalog, Ncatalog); 87 95 88 96 SAVEPLOT = FALSE; 97 98 client_logger_message ("starting the loops: %s\n", regionHosts->hosts[myHost].hostname); 89 99 90 100 /* major modes */ … … 105 115 UpdateMeasures (catalog, Ncatalog); 106 116 MARKTIME("reset measures, host %d: %f sec\n", REGION_HOST_ID, dtime); 117 client_logger_message ("reset measures: %s\n", regionHosts->hosts[myHost].hostname); 107 118 } 108 119 for (i = 0; i < NLOOP; i++) { 109 120 UpdateObjects (catalog, Ncatalog, (i > 0)); 110 121 MARKTIME("update objects loop %d, host %d: %f sec\n", i, REGION_HOST_ID, dtime); 122 client_logger_message ("update objects loop %d: %s\n", i, regionHosts->hosts[myHost].hostname); 111 123 if ((i > 1) || !USE_GALAXY_MODEL) { 112 124 // if GALAXY_MODEL is selected, we want to delay the frame correction until we have … … 114 126 FrameCorrectionParallelSlave (catalog, Ncatalog, regionHosts, i); 115 127 MARKTIME("frame correction loop %d, host %d: %f sec\n", i, REGION_HOST_ID, dtime); 128 client_logger_message ("frame correction loop %d: %s\n", i, regionHosts->hosts[myHost].hostname); 116 129 } 117 130 share_mean_pos (catalog, Ncatalog, regionHosts, i); 118 131 slurp_mean_pos (catalog, Ncatalog, regionHosts, i); 132 client_logger_message ("slurp mags loop %d: %s\n", i, regionHosts->hosts[myHost].hostname); 119 133 120 134 UpdateChips (catalog, Ncatalog); 121 135 MARKTIME("update chips loop %d, host %d: %f sec\n", i, REGION_HOST_ID, dtime); 136 client_logger_message ("update chips loop %d: %s\n", i, regionHosts->hosts[myHost].hostname); 137 122 138 share_image_pos (regionHosts, i); 123 139 slurp_image_pos (catalog, Ncatalog, regionHosts, i); 140 client_logger_message ("slurp image pos loop %d: %s\n", i, regionHosts->hosts[myHost].hostname); 141 124 142 share_meas_pos (catalog, Ncatalog, regionHosts, i); 125 143 slurp_meas_pos (catalog, Ncatalog, regionHosts, i); 126 144 MARKTIME("done exchange loop %d, host %d: %f sec\n", i, REGION_HOST_ID, dtime); 145 client_logger_message ("slurp meas loop %d: %s\n", i, regionHosts->hosts[myHost].hostname); 127 146 } 128 147 // create summary plots of the process … … 159 178 160 179 share_image_pos (regionHosts, -1); 180 client_logger_message ("share image pos loop %d: %s\n", -1, regionHosts->hosts[myHost].hostname); 161 181 162 182 exit (0);
Note:
See TracChangeset
for help on using the changeset viewer.
