Second DX Calculation of Influence Functions for TS3 20180302

23:31 Start completely over with data collected on 20180209
  • mv if_functions/20180209_mirmodes if_functions/20180209_mirmodes_orig
  • from TWiki 20180209
list_low_amp=['20180209_210838','20180209_211239','20180209_211609','20180209_215352','20180209_215709','20180209_220108','20180209_220305','20180209_220534','20180209_220724','20180209_220900','20180209_221054','20180209_221221','20180209_221358','20180209_221655','20180209_221909','20180209_222027','20180209_222154','20180209_222424','20180209_222627','20180209_222805','20180209_222931','20180209_223123','20180209_223306','20180209_223538','20180209_223640','20180209_223755','20180209_223916','20180209_224155','20180209_224344','20180209_224515','20180209_224716','20180209_224844','20180209_225045','20180209_225236','20180209_225411','20180209_225522','20180209_225654','20180209_225811','20180209_225942']

23:34
  • ssh AOeng@dxadsec
  • aoidl
 @startup

IDL> mm = obj_new('mirmodes', list_low_amp, data_dir = '/local/aomeas/adsec_calib/if_functions')

IDL> mm->summary_opds
% Compiled module: MINMAX.
tracknum 		 nn 	 amp [nm] 	 first 	  last 	 missing    sampling
20180209_210838 	 125 	  500.0 	 0 	 4 	 0 	 990 x 1002
20180209_211239 	 125 	  350.0 	 5 	 9 	 0 	 990 x 1002
20180209_211609 	 51 	  250.0 	 10 	 19 	 0 	 990 x 1002
20180209_215352 	 35 	  150.0 	 20 	 37 	 0 	 990 x 1002
20180209_215709 	 35 	  150.0 	 38 	 55 	 0 	 990 x 1002
20180209_220108 	 35 	  150.0 	 56 	 73 	 0 	 990 x 1002
20180209_220305 	 35 	  150.0 	 74 	 91 	 0 	 990 x 1002
20180209_220534 	 35 	  100.0 	 92 	 109 	 0 	 990 x 1002
20180209_220724 	 35 	  100.0 	 110 	 127 	 0 	 990 x 1002
20180209_220900 	 35 	   80.0 	 128 	 145 	 0 	 990 x 1002
20180209_221054 	 35 	   80.0 	 146 	 163 	 0 	 990 x 1002
20180209_221221 	 35 	   80.0 	 164 	 181 	 0 	 990 x 1002
20180209_221358 	 35 	   80.0 	 182 	 199 	 0 	 990 x 1002
20180209_221655 	 35 	   70.0 	 200 	 217 	 0 	 990 x 1002
20180209_221909 	 35 	   50.0 	 218 	 235 	 0 	 990 x 1002
20180209_222027 	 35 	   50.0 	 236 	 253 	 0 	 990 x 1002
20180209_222154 	 35 	   50.0 	 254 	 271 	 0 	 990 x 1002
20180209_222424 	 35 	   50.0 	 272 	 289 	 0 	 990 x 1002
20180209_222627 	 35 	   40.0 	 290 	 307 	 0 	 990 x 1002
20180209_222805 	 35 	   40.0 	 308 	 325 	 0 	 990 x 1002
20180209_222931 	 35 	   40.0 	 326 	 343 	 0 	 990 x 1002
20180209_223123 	 35 	   40.0 	 344 	 361 	 0 	 990 x 1002
20180209_223306 	 35 	   40.0 	 362 	 370 	 0 	 990 x 1002
20180209_223538 	 35 	   40.0 	 371 	 379 	 0 	 990 x 1002
20180209_223640 	 35 	   40.0 	 380 	 397 	 0 	 990 x 1002
20180209_223755 	 35 	   40.0 	 398 	 415 	 0 	 990 x 1002
20180209_223916 	 35 	   40.0 	 416 	 433 	 0 	 990 x 1002
20180209_224155 	 35 	   40.0 	 434 	 451 	 0 	 990 x 1002
20180209_224344 	 35 	   40.0 	 452 	 469 	 0 	 990 x 1002
20180209_224515 	 35 	   40.0 	 470 	 487 	 0 	 990 x 1002
20180209_224716 	 35 	   40.0 	 488 	 505 	 0 	 990 x 1002
20180209_224844 	 35 	   40.0 	 506 	 523 	 0 	 990 x 1002
20180209_225045 	 35 	   40.0 	 524 	 541 	 0 	 990 x 1002
20180209_225236 	 35 	   40.0 	 542 	 559 	 0 	 990 x 1002
20180209_225411 	 35 	   40.0 	 560 	 577 	 0 	 990 x 1002
20180209_225522 	 35 	   40.0 	 578 	 595 	 0 	 990 x 1002
20180209_225654 	 35 	   40.0 	 596 	 613 	 0 	 990 x 1002
20180209_225811 	 35 	   40.0 	 614 	 631 	 0 	 990 x 1002
20180209_225942 	 35 	   40.0 	 632 	 642 	 0 	 990 x 1002

IDL> mm->collect_opds
Collecting OPDs....(go grab a coffee)
Collecting OPDs completed!

IDL> mm->mask_calibration, /SHOW

  • Mask Images: valid_mask_points
    valid_mask_points.jpg
  • Mask Images: intersection_mask
    intersection_mask.jpg
  • Mask Images: union_mask
    union_mask.jpg

00:32 Continue
IDL> mm->check_noise
  • Noise Plots: surface_nm
    surface_nm.jpg
  • Noise Plots: surface_nm_rms
    surface_nm_rms.jpg

01:39 Continue
IDL> mm->inv_ifmat  
Inverting IFmatrix...
% Compiled module: PSEUDO_INVERT.
% Compiled module: DO_SVD.
% Compiled module: SYSV_SET.
% Loaded DLM: LAPACK.
% Program caused arithmetic error: Floating underflow
% Detected at DO_SVD             26 /usr/local/exelisvis/idl71/lib/oaa_lib/utilities/do_svd.pro
min/max Singular value: Max=0.0010908856 (Max singular value=32.648881)
% Compiled module: REVERSE.
% Compiled module: DO_TRANSPOSE.
% Compiled module: MATMAT_MULTIPLY.
IFmatrix inverted.
Computing DMcoupling...
  • Matrix Invert Plots: single_values
    single_values.jpg
  • Matrix Invert Plots: eigenmodes
    eigenmodes.jpg

IDL> mm->interactive_tweak_mask
...Restoring modal matrices and mask data...
...Computing subset of Zernikes for fitting test...
% Compiled module: ZERN.
=================================
MASK OUTER BORDER TUNING:
press "D" to Reduce mask
press "I" to Increase the mask
press "A" to Accept current mask
=================================
Great! Mask Size Tuned!
=================================
CENTRAL OCCULTATION TUNING:
press "D" to Decrease mask
press "I" to Increase the mask
press "A" to Accept current mask
=================================
Great! Mask Central Obscuration Tuned!
Ellapsed time: 1.7551239
Now let's save the mask information!
Type a tag for the mask identification (only valid characters please):second
  • Mask Images: outer
    mask_outer.jpg
  • Mask Images: inner
    mask_inner.jpg

02:51 Continue
IDL> mm->showme_the_modes, nrows=20, ncol=30, zoom=0.11
  • MM Mode Plots:
    mm_modes.jpg

03:00 Continue
IDL> mm->tweak_ifmat, "second"

IDL> mm->inv_ifmat, mask_id="second"
  • Inversion Plot with Second pupil: single_values
    single_values_second.jpg
  • Inversion Plot with Second pupil: eigenmodes
    eigenmodes_seconds.jpg

03:38 Save the mm structure
  • IDL> save, mm, FILE='/local/aomeas/adsec_calib/if_functions/20180209_mirmodes/mm_struct.sav'

04:03 Create the klf object
IDL> klf = obj_new('klfit', mm, "second")

IDL> klf->create_theo_kls

IDL> klf->create_ortho_fitted_kls, klmax=500

IDL> klf->check_fitting
  • KL Fitting Plots: fitting_error
    fitting_error.jpg
  • KL Fitting Plots: peak_force
    peak_force.jpg

04:59 Continue
IDL> print, klf->kl_theo_dir()
/local/aomeas/adsec_calib/if_functions/20180209_mirmodes/KLdata_theo_second/
IDL> print, klf->kl_fit_dir()
/local/aomeas/adsec_calib/if_functions/20180209_mirmodes/KLdata_fit_second_/

IDL> save, klf, FILE='/local/aomeas/adsec_calib/if_functions/20180209_mirmodes/klf_object.sav'

06:26 Create KL_v23 with 500 mode KL basis
[AOeng@dxadsec M2C]$ make_m2c.py KL_v23
Configuration directory: /home/aoeng/aoroot/conf/adsec/current/processConf

Copying data...

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

M2C file should be copied manually with this command:

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

06:30 aoidl killed by adsc_stop
$ aoidl
IDL> @startup

IDL> restore, '/local/aomeas/adsec_calib/if_functions/20180209_mirmodes/klf_object.sav', /v
% RESTORE: Portable (XDR) SAVE/RESTORE file.
% RESTORE: Save file written by AOeng@dxadsec.mountain.lbto.org, Fri Mar  2 05:24:05 2018.
% RESTORE: IDL version 7.1 (linux, x86_64).
% RESTORE: Restored variable: KLF.
IDL> restore, '/local/aomeas/adsec_calib/if_functions/20180209_mirmodes/mm_struct.sav', /v 
% RESTORE: Portable (XDR) SAVE/RESTORE file.
% RESTORE: Save file written by AOeng@dxadsec.mountain.lbto.org, Fri Mar  2 03:36:53 2018.
% RESTORE: IDL version 7.1 (linux, x86_64).
% RESTORE: Restored variable: MM.

IDL> mm->save_MMmatrix, "FLAO1", mask_id="second"
Loading IF matrix...
Loading MM2C....
Saving MMmatrix file: /local/aomeas/adsec_calib/if_functions/20180209_mirmodes/MMmatrix_FLAO1_20180209.sav

IDL> klf->save_KLmatrix, "KL_v23"
Restoring KL data...
% READFITS: Now reading 672 by 672 array
Saving KLmatrix file: /local/aomeas/adsec_calib/if_functions/20180209_mirmodes/KLdata_fit_second_/KLmatrix_KL_v23.sav

-- %USERSIG{DougMiller - 2018-03-01}%
I Attachment Action Size Date Who Comment
eigenmodes.jpgjpg eigenmodes.jpg manage 52 K 02 Mar 2018 - 02:04 DougMiller Matrix Invert Plots
eigenmodes_seconds.jpgjpg eigenmodes_seconds.jpg manage 50 K 02 Mar 2018 - 03:28 DougMiller Inversion Plot with Second pupil
fitting_error.jpgjpg fitting_error.jpg manage 81 K 02 Mar 2018 - 05:06 DougMiller KL Fitting Plots
intersection_mask.jpgjpg intersection_mask.jpg manage 66 K 02 Mar 2018 - 00:30 DougMiller Mask Images
mask_inner.jpgjpg mask_inner.jpg manage 92 K 02 Mar 2018 - 02:22 DougMiller Mask Images:
mask_outer.jpgjpg mask_outer.jpg manage 91 K 02 Mar 2018 - 02:22 DougMiller Mask Images:
mm_modes.jpgjpg mm_modes.jpg manage 1 MB 02 Mar 2018 - 02:53 DougMiller MM Mode Plots
peak_force.jpgjpg peak_force.jpg manage 81 K 02 Mar 2018 - 05:06 DougMiller KL Fitting Plots
single_values.jpgjpg single_values.jpg manage 37 K 02 Mar 2018 - 02:04 DougMiller Matrix Invert Plots
single_values_second.jpgjpg single_values_second.jpg manage 36 K 02 Mar 2018 - 03:27 DougMiller Inversion Plot with Second pupil
surface_nm.jpgjpg surface_nm.jpg manage 78 K 02 Mar 2018 - 01:13 DougMiller Noise Plots
surface_nm_rms.jpgjpg surface_nm_rms.jpg manage 95 K 02 Mar 2018 - 01:13 DougMiller Noise Plots
union_mask.jpgjpg union_mask.jpg manage 70 K 02 Mar 2018 - 00:30 DougMiller Mask Images
valid_mask_points.jpgjpg valid_mask_points.jpg manage 133 K 02 Mar 2018 - 00:30 DougMiller Mask Images
Topic revision: r3 - 02 Mar 2018, 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