20220501_DX_Day DX AdSec Calculate Influence Functions, Day 12

08:00 Wrote IDL routine to check if shape removal was included in the 4D processing of the phase maps
  • /home/aoeng/doug/i4d/i4d_check_remove.pro

09:00 Confirm
  • Set_1 from 20220424 Shape Removal
  • Set_1 from 20220424 Shape Removal
  • Set_3 from 20220425 Shape Removal
  • Set_4 from 20220425 No Shape Removal!

09:40 Process Set_4 with 505 modes
IDL> produce =  ['20220425_182443','20220425_182601','20220425_182623','20220425_182654','20220425_182716','20220425_182739','20220425_182802','20220425_182828','20220425_182855','20220425_182922','20220425_182949','20220425_183015','20220425_183042','20220425_183109','20220425_183135','20220425_183202','20220425_183229','20220425_183255','20220425_183322','20220425_183349','20220425_183415','20220425_183506','20220425_183533','20220425_183600','20220425_183627','20220425_183653','20220425_183720','20220425_183746','20220425_183813','20220425_183840','20220425_183906','20220425_183933','20220425_184000','20220425_184026','20220425_184053','20220425_184120','20220425_184146','20220425_184213','20220425_184240','20220425_184306','20220425_184333','20220425_184433','20220425_184459','20220425_184526','20220425_184553','20220425_184619','20220425_184646','20220425_184712','20220425_184739','20220425_184806','20220425_184833','20220425_184900','20220425_184927','20220425_184953','20220425_185020','20220425_185047','20220425_185113','20220425_185140','20220425_185207','20220425_185234','20220425_185300','20220425_185431']
  • IDL> mm=obj_new('mirmodes', produce, binning=4)
  • IDL> mm->summary_opds
tracknum 		 nn 	 amp [nm] 	 first 	  last 	 missing    sampling
20220425_182443 	 125 	  500.0 	 0 	 2 	 0 	 892 x 892
20220425_182601 	 125 	  500.0 	 3 	 4 	 0 	 892 x 892
20220425_182623 	 125 	  350.0 	 5 	 7 	 0 	 892 x 892
20220425_182654 	 125 	  350.0 	 8 	 9 	 0 	 892 x 892
20220425_182716 	 51 	  250.0 	 10 	 14 	 0 	 892 x 892
20220425_182739 	 51 	  250.0 	 15 	 19 	 0 	 892 x 892
20220425_182802 	 35 	  150.0 	 20 	 28 	 0 	 892 x 892
20220425_182828 	 35 	  150.0 	 29 	 37 	 0 	 892 x 892
20220425_182855 	 35 	  150.0 	 38 	 46 	 0 	 892 x 892
20220425_182922 	 35 	  150.0 	 47 	 55 	 0 	 892 x 892
20220425_182949 	 35 	  150.0 	 56 	 64 	 0 	 892 x 892
20220425_183015 	 35 	  150.0 	 65 	 73 	 0 	 892 x 892
20220425_183042 	 35 	  150.0 	 74 	 82 	 0 	 892 x 892
20220425_183109 	 35 	  150.0 	 83 	 91 	 0 	 892 x 892
20220425_183135 	 35 	  100.0 	 92 	 100 	 0 	 892 x 892
20220425_183202 	 35 	  100.0 	 101 	 109 	 0 	 892 x 892
20220425_183229 	 35 	  100.0 	 110 	 118 	 0 	 892 x 892
20220425_183255 	 35 	  100.0 	 119 	 127 	 0 	 892 x 892
20220425_183322 	 35 	   80.0 	 128 	 136 	 0 	 892 x 892
20220425_183349 	 35 	   80.0 	 137 	 145 	 0 	 892 x 892
20220425_183415 	 35 	   80.0 	 146 	 154 	 0 	 892 x 892
20220425_183506 	 35 	   80.0 	 155 	 163 	 0 	 892 x 892
20220425_183533 	 35 	   80.0 	 164 	 172 	 0 	 892 x 892
20220425_183600 	 35 	   80.0 	 173 	 181 	 0 	 892 x 892
20220425_183627 	 35 	   80.0 	 182 	 190 	 0 	 892 x 892
20220425_183653 	 35 	   80.0 	 191 	 199 	 0 	 892 x 892
20220425_183720 	 35 	   70.0 	 200 	 208 	 0 	 892 x 892
20220425_183746 	 35 	   70.0 	 209 	 217 	 0 	 892 x 892
20220425_183813 	 35 	   50.0 	 218 	 226 	 0 	 892 x 892
20220425_183840 	 35 	   50.0 	 227 	 235 	 0 	 892 x 892
20220425_183906 	 35 	   50.0 	 236 	 244 	 0 	 892 x 892
20220425_183933 	 35 	   50.0 	 245 	 253 	 0 	 892 x 892
20220425_184000 	 35 	   50.0 	 254 	 262 	 0 	 892 x 892
20220425_184026 	 35 	   50.0 	 263 	 271 	 0 	 892 x 892
20220425_184053 	 35 	   50.0 	 272 	 280 	 0 	 892 x 892
20220425_184120 	 35 	   50.0 	 281 	 289 	 0 	 892 x 892
20220425_184146 	 35 	   40.0 	 290 	 298 	 0 	 892 x 892
20220425_184213 	 35 	   40.0 	 299 	 307 	 0 	 892 x 892
20220425_184240 	 35 	   40.0 	 308 	 316 	 0 	 892 x 892
20220425_184306 	 35 	   40.0 	 317 	 325 	 0 	 892 x 892
20220425_184333 	 35 	   40.0 	 326 	 334 	 0 	 892 x 892
20220425_184433 	 35 	   40.0 	 335 	 343 	 0 	 892 x 892
20220425_184459 	 35 	   40.0 	 344 	 352 	 0 	 892 x 892
20220425_184526 	 35 	   40.0 	 353 	 361 	 0 	 892 x 892
20220425_184553 	 35 	   40.0 	 362 	 370 	 0 	 892 x 892
20220425_184619 	 35 	   40.0 	 371 	 379 	 0 	 892 x 892
20220425_184646 	 35 	   40.0 	 380 	 388 	 0 	 892 x 892
20220425_184712 	 35 	   40.0 	 389 	 397 	 0 	 892 x 892
20220425_184739 	 35 	   40.0 	 398 	 406 	 0 	 892 x 892
20220425_184806 	 35 	   40.0 	 407 	 415 	 0 	 892 x 892
20220425_184833 	 35 	   40.0 	 416 	 424 	 0 	 892 x 892
20220425_184900 	 35 	   40.0 	 425 	 433 	 0 	 892 x 892
20220425_184927 	 35 	   40.0 	 434 	 442 	 0 	 892 x 892
20220425_184953 	 35 	   40.0 	 443 	 451 	 0 	 892 x 892
20220425_185020 	 35 	   40.0 	 452 	 460 	 0 	 892 x 892
20220425_185047 	 35 	   40.0 	 461 	 469 	 0 	 892 x 892
20220425_185113 	 35 	   40.0 	 470 	 478 	 0 	 892 x 892
20220425_185140 	 35 	   40.0 	 479 	 487 	 0 	 892 x 892
20220425_185207 	 35 	   40.0 	 488 	 496 	 0 	 892 x 892
20220425_185234 	 35 	   40.0 	 497 	 505 	 0 	 892 x 892
20220425_185300 	 35 	   40.0 	 506 	 514 	 0 	 892 x 892
20220425_185431 	 35 	   40.0 	 515 	 523 	 0 	 892 x 892
  • IDL> mm->collect_opds
    • Collecting OPDs....(go grab a coffee)
    • Collecting OPDs completed!
  • IDL> mm->mask_calibration, /SHOW
    • Plots look the same as previous analyses
  • IDL> mm->produce_ifmat
  • IDL> mm->check_noise
    • Plots look different. More spurious TT signal (???)
    • 2 plots
    • 505_surface_rms_s4.jpg:
      505_surface_rms_s4.jpg
    • 505_surface_s4.jpg:
      505_surface_s4.jpg
  • IDL> mm->inv_ifmat
    • last Eigenmode now looks flat (???)
    • Last eigenmode amplitudes are half the size of 2019 analysis
    • 2 plots
    • 505_last_eigenmode_s4.jpg:
      505_last_eigenmode_s4.jpg
    • 505_single_values.jpg:
      505_single_values.jpg
  • IDL> mm->interactive_tweak_mask
    • outer 0.29
    • inner 0.798
    • 2 plots
  • IDL> mm->tweak_ifmat, '.29.8'
    • 505_outer_s4.jpg:
      505_outer_s4.jpg
    • 505_inner_s4.jpg:
      505_inner_s4.jpg
  • IDL> mm->inv_ifmat, mask_id=".29.8"
    • 505_last_eigenmode.29.8_s4.jpg:
      505_last_eigenmode.29.8_s4.jpg
    • 505_single_values.29.8.jpg:
      505_single_values.29.8.jpg
  • IDL> mm->showme_the_modes, nrows=16, ncol=20, zoom=0.4
  • IDL> list=indgen(185)+320
  • IDL> mm->showme_the_modes, mode_num_idx=list, nrows=10, ncol=20, zoom=0.4
    • 2 plots
    • 505_mm_shapes_1.jpg:
      505_mm_shapes_1.jpg
    • 505_mm_shapes_2.jpg:
      505_mm_shapes_2.jpg
  • IDL> save, mm, file='/local/aomeas/adsec_calib/if_functions/20220425_mirmodes/mm_.29.8.sav'
  • IDL> klf = obj_new('klfit', mm, '.29.8')
  • IDL> klf->create_theo_kls
  • IDL> klf->create_ortho_fitted_kls
    • NOTE: Only 405 KL Shaped used (???)
    • Error
Restoring Mirror Modes Influence Function data...
% RESTORE: Portable (XDR) compressed SAVE/RESTORE file.
% RESTORE: Save file written by AOeng@dxadsec.mountain.lbto.org, Sun May  1 17:48:26 2022.
% RESTORE: IDL version 7.1 (linux, x86_64).
% RESTORE: Restored variable: DMCOUPLING.
% RESTORE: Restored variable: INV_IFMATRIX.
% RESTORE: Restored variable: IFMATRIX.
Restoring theoretical KL modes data...
Computing M2C...
Computing fitted KLs...
Computing fitting error...
Saving M2V data...
Re-orthonormalizing fitted KLs...
Number of final KLs:          405
Computing ortho-fitted KLs...
Computing ortho M2C...
Saving ortho M2V data...
% Object reference type required in this context: MM_OBJ.
% Execution halted at: KLFIT::CREATE_ORTHO_FITTED_KLS  249 /home/aoeng/soul/mixed_mode_control/aoroot/idl/adsec_lib/calibration/optical/klfit__define.pro
%                      $MAIN$          
  • IDL> klf->check_fitting
    • fitting_error plot (better ???)
    • Error

20220502

Xianyu fixed an error in klfit__define.pro

  • added _ to mm_obj (self._mm_obj->save_dir() now works)

10:30 re-run kl fit routine
  • IDL> klf->create_ortho_fitted_kls
Restoring Mirror Modes Influence Function data...
% RESTORE: Portable (XDR) compressed SAVE/RESTORE file.
% RESTORE: Save file written by AOeng@dxadsec.mountain.lbto.org, Sun May  1 17:48:26 2022.
% RESTORE: IDL version 7.1 (linux, x86_64).
% RESTORE: Restored variable: DMCOUPLING.
% RESTORE: Restored variable: INV_IFMATRIX.
% RESTORE: Restored variable: IFMATRIX.
Restoring theoretical KL modes data...
Computing M2C...
Computing fitted KLs...
Computing fitting error...
Saving M2V data...
Re-orthonormalizing fitted KLs...
Number of final KLs:          405
Computing ortho-fitted KLs...
Computing ortho M2C...
Saving ortho M2V data...
% Compiled module: MAKE_FITS_HEADER.
% Compiled module: FXADDPAR.
% Compiled module: FXPARPOS.
% Compiled module: DETABIFY.
% Compiled module: CHECK_FITS.
% Compiled module: FXPAR.
Computing Force requirements...
Finished!!!
Time required:       0.69599698 minutes
IDL> klf->check_fitting
  • 505_fitting_error_02.jpg:
    505_fitting_error_02.jpg
  • 505_fitting_error_02.jpg:
    505_fitting_error_02.jpg

Old version from 20220501

  • IDL> klf->check_fitting
% Compiled module: READFITS. % Compiled module: VALID_NUM. % READFITS: Now reading 702 element array % Compiled module: TEXTOIDL. % Compiled module: TEXTABLE. % Compiled module: STRTRANS. % Compiled module: TRANSLATE_SUB_SUPER. % Compiled module: NEXTTOK. % RESTORE: Error opening file. Unit: 100, File: /towerdata/adsec_calib/if_functions/20220425_mirmodes/KLdata_fit_.29.8_/force_data.sav No such file or directory

% Execution halted at: KLFIT::CHECK_FITTING 288 /home/aoeng/soul/mixed_mode_control/aoroot/idl/adsec_lib/calibration/optical/klfit__define.pro % $MAIN$
    • 505_fitting_error.jpg:
      505_fitting_error.jpg

Old version from 20220501

  • IDL> klf->showme_the_modes, 'ORTHOFIT', nrows=16, ncol=20, zoom=0.4
  • IDL> list=indgen(185)+320 (too many. This is for 505 but only 405 KL fits to display)
  • IDL> klf->showme_the_modes, 'ORTHOFIT', nrows=5, ncol=20, zoom=0.4, mode_num_idx=list
    • 2 plots
    • 505_kl_shapes_1.jpg:
      505_kl_shapes_1.jpg
    • 505_kl_shapes_2.jpg:
      505_kl_shapes_2.jpg
  • IDL> print, klf->kl_theo_dir()
    • /towerdata/adsec_calib/if_functions/20220425_mirmodes/KLdata_theo_.29.8/

  • [AOeng@dxadsec M2C]$ pwd
    • /local/aomeas/adsec_calib/M2C
  • [AOeng@dxadsec M2C]$ make_m2c.py KL_v28
Copying data...

Created directory: /local/aomeas/adsec_calib/M2C/KL_v28

M2C file should be copied manually with this command:

cp m2c.fits /local/aomeas/adsec_calib/M2C/KL_v28

  • IDL> print, klf->kl_fit_dir()
    • /towerdata/adsec_calib/if_functions/20220425_mirmodes/KLdata_fit_.29.8_/

  • [AOeng@dxadsec M2C]$ cd /towerdata/adsec_calib/if_functions/20220425_mirmodes/KLdata_fit_.29.8_/
  • [AOeng@dxadsec KLdata_fit_.29.8_]$ ls
    • fitting_error.fits M2V_data.sav orthofit_M2V_data.sav
    • missing files: force_data.sav, KLmatrix_KL_v27.sav, m2c.fits. orthofit_M2V_data.sav
    • must be due to error in klf->create_ortho_fitted_kls (see above) when "Saving ortho M2V data..."
  • [AOeng@dxadsec KLdata_fit_.29.8_]$ ls /local/aomeas/adsec_calib/tt_functions/20190907_mirmodes/KLdata_fit_.29.29_
    • fitting_error.fits force_data.sav KLmatrix_KL_v27.sav m2c.fits M2V_data.sav orthofit_M2V_data.sav

12:17 Next analysis tests
  • Use ADOPT_ROOT with int4d, used in 2019.
    • did routines used for KL calculation get modified by MMC upgrade?
  • Use 600 modes instead for 505 or 632.

IDL> help, /recall_commands       
Recall buffer length: 500
1	help, /recall_commands
2	restore, '/towerdata/adsec_calib/if_functions/20220425_mirmodes/adsec_shell.sav', /v
3	help, adsec.act_w_cl
4	ADSEC=ADSEC_SAVE
5	restore, '/towerdata/adsec_calib/if_functions/20220425_mirmodes/adsec_save.sav', /v
6	help, adsec.act_w_cl
7	restore, '/towerdata/adsec_calib/if_functions/20220425_mirmodes/adsec_shell.sav', /v
8	restore, '/towerdata/adsec_calib/if_functions/20220425_182443/adsec_shell.sav', /v
9	retall
10	restore, '/towerdata/adsec_calib/if_functions/20220425_182443/adsec_shell.sav', /v
11	make_offload_matrix_proj_new,'20220425_mirmodes', /save, tracknum='20220425_mirmodes'
12	make_offload_matrix_proj_new,'20220425_mirmodes', /do_save, tracknum='20220425_mirmodes'
13	ADSEC=ADSEC_SAVE
14	ADSEC=ADSEC_SAV
15	restore, '/towerdata/adsec_calib/if_functions/20220425_182443/adsec_save.sav', /v
16	make_offload_matrix_proj_new,'20220425_mirmodes', /save, tracknum='20220425_mirmodes'
17	.comp /home/aoeng/soul/mixed_mode_control/aoroot/idl/adsec_lib/mmc_removal/make_offload_matrix_proj_new.pro
18	make_offload_matrix_proj_new,'20220425_mirmodes', /do_save, tracknum='20220425_mirmodes'
19	retall
20	make_offload_matrix_proj_new,'20220425_mirmodes', /do_save, tracknum='20220425_mirmodes'
21	make_offload_matrix_proj_new,'20220425_mirmodes', /save, tracknum='20220425_mirmodes'
22	im2mmmatrix_new, '20220425_mirmodes', rebfactor=4
23	$pwd
24	im2mmmatrix_new, '20220425', rebfactor=4
25	im2mmmatrix_new, '../if_functions/20220425_mirmodes', rebfactor=4
26	.compile /home/aoeng/aoroot/idl/adsec_lib/test/im2mmmatrix_new.pro
27	im2mmmatrix_new, '../tt_functions/20190907_mirmodes', rebfactor=4
28	.compile /home/aoeng/aoroot/idl/adsec_lib/test/im2mmmatrix_new.pro
29	print, dpix
30	restore, '/local/aomeas/adsec_calib/if_functions/20220425_mirmodes/KLdata_fit_.29.8_/KLmatrix_KL_v28.sav', /v
31	print, dpix
32	restore, '/local/aomeas/adsec_calib/tt_functions/20190907_mirmodes/MMmatrix_FLAO1_20190907.sav', /v
33	klf->save_KLmatrix, "KL_v28"
34	mm->save_MMmatrix, 'FLAO1', mask_id='.29.8'
35	klf->showme_the_modes, 'ORTHOFIT', nrows=16, ncol=20, zoom=0.4
36	klf->check_fitting
37	klf->create_ortho_fitted_kls
38	.compile klfit__define
39	retall
40	print, mm->save_dir()
41	retall
42	print, klf->kl_fit_dir()
43	print, klf->kl_theo_dir()
44	retall
45	klf->showme_the_modes, 'ORTHOFIT', nrows=5, ncol=20, zoom=0.4, mode_num_idx=list
46	list=indgen(185)+320
47	klf->showme_the_modes, 'ORTHOFIT', nrows=16, ncol=20, zoom=0.4
48	retall
49	klf->showme_the_modes, 'ORTHOFIT', nrows=16, ncol=20, zoom=0.4
50	klf->check_fitting
51	retall
52	klf->create_ortho_fitted_kls
53	retall
54	klf->check_fitting
55	klf->create_ortho_fitted_kls
56	IDL> klf->create_ortho_fitted_kls
57	klf->create_theo_kls
58	klf = obj_new('klfit', mm, '.29.8')
59	save, mm, file='/local/aomeas/adsec_calib/if_functions/20220425_mirmodes/mm_.29.8.sav'
60	mm->showme_the_modes, mode_num_idx=list, nrows=10, ncol=20, zoom=0.4
61	mm->showme_the_modes, mode_num_idx=list, nrows=16, ncol=20, zoom=0.4
62	list=indgen(185)+320
63	mm->showme_the_modes, nrows=16, ncol=20, zoom=0.4
64	mm->inv_ifmat, mask_id=".29.8"
65	mm->tweak_ifmat, '.29.8'
66	mm->interactive_tweak_mask
67	mm->inv_ifmat
68	mm->check_noise
69	 mm->produce_ifmat
70	mm->mask_calibration, /SHOW
71	mm->collect_opds
72	mm->summary_opds
73	mm=obj_new('mirmodes', produce, binning=4)
74	mm->summary_opds
75	retall
76	mm->summary_opds
77	IDL> mm->summary_opds
78	IDLmm=obj_new('mirmodes', produce, binning=4)
79	IDL> mm=obj_new('mirmodes', produce, binning=4)
80	produce = ['20220425_182443','20220425_182601','20220425_182623','20220425_182654','20220425_182716','20220425_182739','20220425_182802','20220425_182828','20220425_182855','20220425_182922','20220425_182949','20220425_183015','20220425_183042','20220425_183109','20220425_183135','20220425_183202','20220425_183229','20220425_183255','20220425_183322','20220425_183349','20220425_183415','20220425_183506','20220425_183533','20220425_183600','20220425_183627','20220425_183653','20220425_183720','20220425_183746','20220425_183813','20220425_183840','20220425_183906','20220425_183933','20220425_184000','20220425_184026','20220425_184053','20220425_184120','20220425_184146','20220425_184213','20220425_184240','20220425_184306','20220425_184333','20220425_184433','20220425_184459','20220425_184526','20220425_184553','20220425_184619','20220425_184646','20220425_184712','20220425_184739','20220425_184806','20220425_184833','20220425_184900','20220425_184927','20220425_184953','20220425_185020','20220425_185047','20220425_185113','20220425_185140','20220425_185207','20220425_185234','20220425_185300','20220425_185431']

IDL> make_offload_matrix_proj_new,'20220425_mirmodes', /save, tracknum='20220425_mirmodes'
% RESTORE: Portable (XDR) SAVE/RESTORE file.
% RESTORE: Save file written by AOeng@dxadsec.mountain.lbto.org, Mon May  2 19:11:39 2022.
% RESTORE: IDL version 7.1 (linux, x86_64).
% RESTORE: Restored variable: MMMATRIX.
% RESTORE: Restored variable: DPIX.
% RESTORE: Restored variable: XX.
% RESTORE: Restored variable: YY.
% RESTORE: Restored variable: MM2C.
% RESTORE: Restored variable: IDX_MASK.
pseudo-inverting MMmatrix...
MMMATRIX        FLOAT     = Array[632, 29944]
min/max Singular value: Max=1.8029689e-05 (Max singular value=380.51190)
  (total(max_ifunc[act_w_cl] gt 0)) =       0.00000
  n_elements(act_w_cl) =          632
  n_neg =       632.000
Angular spread [deg] wo x-mirroring:      103.773, w x-mirroring:     0.763476
rotation [deg]:      146.650(sigma:     0.763476)
scale_factor:   0.00226665(sigma:  5.26674e-05)
Zernike to analyze: 00
Zernike to analyze: 01
Zernike to analyze: 02
Zernike to analyze: 03
Zernike to analyze: 04
Zernike to analyze: 05
Zernike to analyze: 06
Zernike to analyze: 07
Zernike to analyze: 08
Zernike to analyze: 09
Zernike to analyze: 10
Zernike to analyze: 11
Zernike to analyze: 12
Zernike to analyze: 13
Zernike to analyze: 14
Zernike to analyze: 15
Zernike to analyze: 16
Zernike to analyze: 17
Zernike to analyze: 18
Zernike to analyze: 19
Zernike to analyze: 20
Zernike to analyze: 21
threshold used for cutting mirror modes:  0.000100000
number of mirror modes used to compiute the zernikes:
      0.00000      98.0000      444.000      573.000      11.0000      7.00000      345.000      152.000      7.00000      161.000      556.000      120.000      555.000      279.000      359.000      160.000      105.000
      314.000      286.000      348.000      363.000      331.000
pseudo-inverting z2c...
min/max Singular value: Max=0.0000000 (Max singular value=13.965402)

[AOeng@dxadsec 20220425_mirmodes]$ ls -ltr ../../optical-projection/20220425_mirmodes/
total 77844
-rw-r--r--. 1 AOeng aoeng 79457996 May  2 19:11 MMmatrix.sav
-rw-r--r--. 1 AOeng aoeng   123840 May  2 20:18 pos2modes.fits
-rw-r--r--. 1 AOeng aoeng   123840 May  2 20:18 pmz.fits

Needed files ???

From Guido:
  • AdOpt> tmp = comp_ho_offload_matrix(CUTMODE=cutmode, /SAV)
  • Data saved in $ADOPT_MEAS/adsec_calib/electric-projection/
  • “CUTMODE” parameter has to be adjusted.
  • Outputs: • pmhofor.fits (NOT USED) • curr2modes.fits (USED for OFFLOAD)

  • [AOeng@dxadsec 20220425_mirmodes]$ ls -l /home/aoeng/soul/mixed_mode_control/aoroot/calib/adsec/current/data.20220421.telescope/
total 178364
-rw-r--r--. 1 AOeng aoeng   3617280 Apr 22 17:06 a_delay_matrix.fits
-rw-r--r--. 1 AOeng aoeng      2412 Apr 22 17:06 astig_lut.sav
-rw-r--r--. 1 AOeng aoeng   7660800 Apr 22 17:06 b0_a_matrix.fits
-rw-r--r--. 1 AOeng aoeng   7660800 Apr 22 17:06 b0_b_matrix.fits
-rw-r--r--. 1 AOeng aoeng 118845056 Apr 22 17:06 ff_+5.99e-05_3.40e-06m.sav
-rw-r--r--. 1 AOeng aoeng   3211200 Apr 22 17:06 ff_f_svec.fits
-rw-r--r--. 1 AOeng aoeng      5760 Apr 22 17:06 ff_list.fits
-rw-r--r--. 1 AOeng aoeng   3211200 Apr 22 17:06 ff_matrix.fits
-rw-r--r--. 1 AOeng aoeng   3211200 Apr 22 17:06 ff_p_svec.fits
-rw-r--r--. 1 AOeng aoeng      8640 Apr 22 17:06 ff_sval.fits
drwxr-xr-x. 2 AOeng aoeng      4096 Apr 22 17:06 flat/
-rw-r--r--. 1 AOeng aoeng  27673900 Apr 22 17:06 layer_phi_000000_time_hist_armando_672a.sav
-rw-r--r--. 1 AOeng aoeng   3617280 Apr 22 17:06 null_pmhofor.fits
-rw-r--r--. 1 AOeng aoeng     14400 Apr 22 17:06 pbgainmatrix.fits
-rw-r--r--. 1 AOeng aoeng   3617280 Apr 22 17:06 pmhofor.fits
-rw-r--r--. 1 AOeng aoeng    123840 Apr 22 17:06 pmz.fits
-rw-r--r--. 1 AOeng aoeng    123840 Apr 22 17:06 pos2modes.fits
[AOeng@dxadsec 20220425_mirmodes]$ more  /home/aoeng/soul/mixed_mode_control/aoroot/calib/adsec/current/data.20220421.telescope/pmhofor.fits 
SIMPLE  =                    T / Written by IDL:  Fri Aug 30 12:57:01 2013      
BITPIX  =                  -64 / Number of bits per data pixel                  
NAXIS   =                    2 / Number of data axes                            
NAXIS1  =                  672 /                                                
NAXIS2  =                  672 /                                                
EXTEND  =                    T / FITS data may contain extensions               
DATE    = '2013-08-30'         / Creation UTC (CCCC-MM-DD) date of FITS header  
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy  
COMMENT and Astrophysics', volume 376, page 359; bibcode 2001A&A...376..359H    
TRACKNUM= '20130830_125701'    /                                                
UNIT    = '672a    '           /                                                
SHELL   = 'TS4     '           /                                                
TYPE    = 'SL      '           /AO or SL mode                                   
M2C     = 'UNKNOW  '           /modes to command recontructor matrix            
END                                                                             
[AOeng@dxadsec 20220425_mirmodes]$ more  /home/aoeng/soul/mixed_mode_control/aoroot/calib/adsec/current/data.20220421.telescope/pmz.fits 
SIMPLE  =                    T / Written by IDL:  Sat Aug 31 09:37:23 2019      
BITPIX  =                  -64 / Number of bits per data pixel                  
NAXIS   =                    2 / Number of data axes                            
NAXIS1  =                   22 /                                                
NAXIS2  =                  672 /                                                
EXTEND  =                    T / FITS data may contain extensions               
DATE    = '2019-08-31'         / Creation UTC (CCCC-MM-DD) date of FITS header  
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy  
COMMENT and Astrophysics', volume 376, page 359; bibcode 2001A&A...376..359H    
TRACKNUM= '20190831_093723'    /                                                
UNIT    = '672a    '           /                                                
SHELL   = 'TS3     '           /                                                
MMMATRIX= '20190831_040458'    /filename of used matrix of FF mirror modes influ
THRESHOL=          0.000100000 /threshold on relative rms used to remove high or
END                                                                             
[AOeng@dxadsec 20220425_mirmodes]$ more  /home/aoeng/soul/mixed_mode_control/aoroot/calib/adsec/current/data.20220421.telescope/pos2modes.fits 
SIMPLE  =                    T / Written by IDL:  Sat Aug 31 09:37:23 2019      
BITPIX  =                  -64 / Number of bits per data pixel                  
NAXIS   =                    2 / Number of data axes                            
NAXIS1  =                  672 /                                                
NAXIS2  =                   22 /                                                
EXTEND  =                    T / FITS data may contain extensions               
DATE    = '2019-08-31'         / Creation UTC (CCCC-MM-DD) date of FITS header  
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy  
COMMENT and Astrophysics', volume 376, page 359; bibcode 2001A&A...376..359H    
TRACKNUM= '20190831_093723'    /                                                
UNIT    = '672a    '           /                                                
SHELL   = 'TS3     '           /                                                
MMMATRIX= '20190831_040458'    /filename of used matrix of FF mirror modes influ
THRESHOL=          0.000100000 /threshold on relative rms used to remove high or
END                                                                             

  • check pos2modes.fits
IDL> b=readfits("/local/aomeas/adsec_calib/optical-projection/20220425_mirmodes/pos2modes.fits", head2)
% READFITS: Now reading 672 by 22 array
IDL> blb= where(max(abs(b), dim=2) EQ 0)                                                               
IDL> print, n_elements(blb)                                                                            
          40
IDL> display, b[*,9], /sh, spot=1, map=512, /no_n                                                      
IDL> display, b[*,10], /sh, spot=1, map=512, /no_n
IDL> blb= where(max(abs(b), dim=2) NE 0)          
IDL> display, b[blb,10], blb, /sh, spot=1, map=512, /no_n
IDL> display, b[blb,10], blb, /sh, spot=1, map=512, /no_n, /sm

[AOeng@dxadsec 20220425_mirmodes]$ cp pmz.fits /home/aoeng/soul/mixed_mode_control/aoroot/calib/adsec/current/data.20220421.telescope/.
[AOeng@dxadsec 20220425_mirmodes]$ cp pos2modes.fits /home/aoeng/soul/mixed_mode_control/aoroot/calib/adsec/current/data.20220421.telescope/.
[AOeng@dxadsec 20220425_mirmodes]$ pwd
/local/aomeas/adsec_calib/optical-projection/20220425_mirmodes

[AOeng@dxadsec phase_maps]$ mkdir /local/aomeas/adsec_calib/if_functions/20190907_mirmodes/phase_maps
[AOeng@dxadsec phase_maps]$ mv KLmatrix_KL_v27.sav /local/aomeas/adsec_calib/if_functions/20190907_mirmodes/phase_maps/.
[AOeng@dxadsec phase_maps]$ mv MMmatrix_FLAO1_20190907.sav /local/aomeas/adsec_calib/if_functions/20190907_mirmodes/phase_maps/.
[AOeng@dxadsec phase_maps]$ mv MMmatrix_FLAO1_20190907.sav /local/aomeas/adsec_calib/if_functions/202204_mirmodes/
Display all 357 possibilities? (y or n)
[AOeng@dxadsec phase_maps]$ ls /local/aomeas/adsec_calib/if_functions/20220425_mirmodes/MM*
/local/aomeas/adsec_calib/if_functions/20220425_mirmodes/MMmatrix_FLAO1_20220425.sav
/local/aomeas/adsec_calib/if_functions/20220425_mirmodes/MMmatrix.sav
[AOeng@dxadsec phase_maps]$ cp /local/aomeas/adsec_calib/if_functions/20220425_mirmodes/MMmatrix_FLAO1_20220425.sav .
[AOeng@dxadsec phase_maps]$ ls
MMmatrix_FLAO1_20220425.sav  phase_screens_test/  phase_screens_test_1000Hz/
[AOeng@dxadsec phase_maps]$ ls /local/aomeas/adsec_calib/if_functions/20220425_mirmodes/KLdata_fit_.29.8_/
fitting_error.fits  force_data.sav  KLmatrix_KL_v28.sav  m2c.fits  M2V_data.sav  orthofit_M2V_data.sav
[AOeng@dxadsec phase_maps]$ cp  /local/aomeas/adsec_calib/if_functions/20220425_mirmodes/KLdata_fit_.29.8_/KLmatrix_KL_v28.sav  .

-- DougMiller - 01 May 2022
I Attachment Action Size Date Who Comment
505_fitting_error.jpgjpg 505_fitting_error.jpg manage 79 K 01 May 2022 - 19:03 DougMiller  
505_fitting_error_02.jpgjpg 505_fitting_error_02.jpg manage 71 K 02 May 2022 - 18:02 DougMiller  
505_inner_s4.jpgjpg 505_inner_s4.jpg manage 101 K 01 May 2022 - 19:03 DougMiller  
505_kl_shapes_1.jpgjpg 505_kl_shapes_1.jpg manage 688 K 01 May 2022 - 19:03 DougMiller  
505_kl_shapes_2.jpgjpg 505_kl_shapes_2.jpg manage 209 K 01 May 2022 - 19:03 DougMiller  
505_last_eigenmode.29.8_s4.jpgjpg 505_last_eigenmode.29.8_s4.jpg manage 60 K 01 May 2022 - 19:04 DougMiller  
505_last_eigenmode_s4.jpgjpg 505_last_eigenmode_s4.jpg manage 56 K 01 May 2022 - 19:04 DougMiller  
505_mm_shapes_1.jpgjpg 505_mm_shapes_1.jpg manage 672 K 01 May 2022 - 19:04 DougMiller  
505_mm_shapes_2.jpgjpg 505_mm_shapes_2.jpg manage 461 K 01 May 2022 - 19:04 DougMiller  
505_outer_s4.jpgjpg 505_outer_s4.jpg manage 86 K 01 May 2022 - 19:04 DougMiller  
505_single_values.29.8.jpgjpg 505_single_values.29.8.jpg manage 37 K 01 May 2022 - 19:05 DougMiller  
505_single_values.jpgjpg 505_single_values.jpg manage 39 K 01 May 2022 - 19:05 DougMiller  
505_surface_rms_s4.jpgjpg 505_surface_rms_s4.jpg manage 97 K 01 May 2022 - 19:05 DougMiller  
505_surface_s4.jpgjpg 505_surface_s4.jpg manage 78 K 01 May 2022 - 19:05 DougMiller  
Topic revision: r4 - 02 May 2022, DougMiller
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding Foswiki? Send feedback