DX Calculation of Influence Functions for TS3

18:02 Working on offload matrix calculations

18:02 MMmatrix.sav contents from DX on 20160216

[AOeng@adsecdx optical-projection]$ cd 20160206_154227
[AOeng@adsecdx 20160206_154227]$ ls
MMmatrix.sav
[AOeng@adsecdx 20160206_154227]$ idl
IDL Version 7.1 (linux x86 m32). (c) 2009, ITT Visual Information Solutions
Installation number: 223182.
Licensed for use by: LBTO-University of Arizona

IDL> restore, 'MMmatrix.sav', /v
% RESTORE: Portable (XDR) SAVE/RESTORE file.
% RESTORE: Save file written by AOeng@adsecdx.mountain.lbto.org, Wed Feb 17 09:54:51 2016.
% RESTORE: IDL version 7.1 (linux, x86).
% 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.
% RESTORE: Restored variable: ACT_W_CL.
IDL> help, mmmatrix, dpix, xx, yy, mm2c, idx_mask, act_w_cl
MMMATRIX        FLOAT     = Array[641, 35969]
DPIX            LONG      =          218
XX              FLOAT     = Array[35969]
YY              FLOAT     = Array[35969]
MM2C            DOUBLE    = Array[641, 672]
IDX_MASK        LONG      = Array[35969]
ACT_W_CL        INT       = Array[641]

19:00 Modify im2mmmatrix_new
  • add
    • xx=xx[idx_mask]
    • yy=yy[idx_mask]

  • IDL> im2mmmatrix_new, '20180209_190244', rebfactor=4
  • IDL> help, MMMATRIX, DPIX, XX, YY, MM2C, IDX_MASK
    • MMMATRIX FLOAT = Array[201, 36653]
    • DPIX LONG = 216
    • XX FLOAT = Array[36653]
    • YY FLOAT = Array[36653]
    • MM2C DOUBLE = Array[201, 672]
    • IDX_MASK LONG = Array[36653]
  • IDL> print, name
  • IDL> print, file_dirname(fname)
    • $ADOPT_MEAS/adsec_calib/optical-projection/20180209_190244/MMmatrix.sav
    • $ADOPT_MEAS/adsec_calib/optical-projection/20180209_190244

  • IDL> make_offload_matrix_proj_new, '20180209_190244'
threshold used for cutting mirror modes:  0.000100000
number of mirror modes used to compiute the zernikes:
      0.00000      200.000      200.000      124.000      194.000      194.000      200.000      199.000      200.000      200.000
      200.000      198.000      184.000      194.000      192.000      82.0000      200.000      200.000      200.000      198.000
      198.000      192.000
pseudo-inverting z2c...
min/max Singular value: Max=0.0000000 (Max singular value=45.709410)

19:35 Try wih 643 modes
  • cd to 20190209_mirmodes
  • cp ../20180209_100244/adsec_save.sav .
  • IDL> im2mmmatrix_new, '20180209', rebfactor=4
    • MMMATRIX FLOAT = Array[643, 35850]
    • DPIX LONG = 216
    • XX FLOAT = Array[35850]
    • YY FLOAT = Array[35850]
    • MM2C DOUBLE = Array[643, 672]
    • IDX_MASK LONG = Array[35850]
    • $ADOPT_MEAS/adsec_calib/optical-projection/20180209/MMmatrix.sav
    • $ADOPT_MEAS/adsec_calib/optical-projection/20180209

  • IDL> make_offload_matrix_proj_new, '20180209'
threshold used for cutting mirror modes:  0.000100000
number of mirror modes used to compiute the zernikes:
      0.00000      366.000      578.000      569.000      177.000      178.000      636.000      627.000      197.000      202.000
      299.000      323.000      231.000      355.000      379.000      296.000      107.000      52.0000      355.000      360.000
      408.000      299.000
pseudo-inverting z2c...
min/max Singular value: Max=0.0000000 (Max singular value=11.143621)

21:01 Copy offload matrix files to data directory
  • [AOeng@dxadsec 20180209]$ pwd
    • /local/aomeas/adsec_calib/optical-projection/20180209
  • [AOeng@dxadsec 20180209]$ cp pmz.fits /home/aoeng/confcalib/calib/adsec/672a/data/.
  • [AOeng@dxadsec 20180209]$ cp pos2modes.fits /home/aoeng/confcalib/calib/adsec/672a/data/.

  • [AOeng@dxadsec data]$ pwd
    • /home/aoeng/confcalib/calib/adsec/672a/data
  • [AOeng@dxadsec data]$ ls -l
-rw-r--r--. 1 AOeng aoeng   3617280 Dec 11 18:52 a_delay_matrix.fits
-rw-r--r--. 1 AOeng aoeng      2412 Feb  4  2017 astig_lut.sav
-rw-r--r--. 1 AOeng aoeng   4305600 Dec 11 18:51 b0_a_matrix.fits
-rw-r--r--. 1 AOeng aoeng   4305600 Dec 11 18:52 b0_b_matrix.fits
-rw-r--r--. 1 AOeng aoeng 118845056 Feb  5 20:34 ff_+5.99e-05_3.40e-06m.sav
-rw-r--r--. 1 AOeng aoeng   3343680 Feb  5 20:34 ff_f_svec.fits
-rw-r--r--. 1 AOeng aoeng      5760 Feb  5 20:34 ff_list.fits
-rw-r--r--. 1 AOeng aoeng   3343680 Feb  5 20:34 ff_matrix.fits
-rw-r--r--. 1 AOeng aoeng   3343680 Feb  5 20:34 ff_p_svec.fits
-rw-r--r--. 1 AOeng aoeng      8640 Feb  5 20:34 ff_sval.fits
-rw-r--r--. 1 AOeng aoeng   3617280 Jan 20  2017 null_pmhofor.fits
-rw-r--r--. 1 AOeng aoeng     14400 Jan 20  2017 pbgainmatrix.fits
-rw-r--r--. 1 AOeng aoeng   3617280 Jan 20  2017 pmhofor.fits
-rw-r--r--. 1 AOeng aoeng    123840 Feb 13 21:01 pmz.fits
-rw-r--r--. 1 AOeng aoeng    123840 Feb 13 21:02 pos2modes.fits

20:56 Make all the needed KL directories
  • $ make_m2c.py KL_v22
  • /local/aomeas/adsec_calib/M2C/KL_v22
    • $ ls
      • disturb filtering gain intmatAcq modesAmp RECs
      • All directories empty except for gain

21:37 more auxiliary files
  • IDL> mm->save_MMmatrix, "FLAO1", mask_id="best"
    • Loading IF matrix...
    • Loading MM2C....
    • Saving MMmatrix file: /local/aomeas/adsec_calib/if_functions/20180209_mirmodes/MMmatrix_FLAO1_20180209.sav

21:40 more files
  • IDL> klf->save_KLmatrix, "KL_v22"
    • Restoring KL data...
    • % READFITS: Now reading 672 by 672 array
    • Saving KLmatrix file: /local/aomeas/adsec_calib/if_functions/20180209_mirmodes/KLdata_fit_best_/KLmatrix_KL_v22.sav

Create Seeing disturbance files

09:00 Calculate disturbance files for a range of seeing
  • Edit file: /home/aoeng/aoroot/idl/elab_lib/batch_generate_disturb.pro
    • disturb_dir = '/local/aomeas/adsec_calib/CMD/disturb'
    • seeing = [0.4, 0.6, 0.8, 1.0, 1.2, 1.5, 2.0, 2.4, 3.0]
    • hz = [1000.]
    • mirmodes_file = '/local/phase_maps/MMmatrix_TS3_20180209.sav'
      • changed name from FLAO1 to TS3, not to insult LBTI
    • dpix = 216
      • From IDL> im2mmmatrix_new, '20180209_190244', rebfactor=4 (see above)

13:00 Create 1.0" seeing at 1000 Hz
  • $ aoidl
  • IDL> @startup
  • IDL> @batch_generate_disturb
  • files:
[AOeng@dxadsec disturb]$ ls -l
total 360180
-rw-r--r--. 1 AOeng aoeng  10396800 Feb 14 12:45 atm_np1140_sz00042.41_L0040_seed3892_nosha_nl2.dat
-rw-r--r--. 1 AOeng aoeng      2540 Feb 14 12:45 atm_np1140_sz00042.41_L0040_seed3892_nosha_nl2.dat.sav
-rw-r--r--. 1 AOeng aoeng  10756800 Feb 14 13:08 dist_flao1_ts4_atm_s1.0_L040.0_v15.0_ovfreq1000.00_sd3892.fits
-rw-r--r--. 1 AOeng aoeng 347652814 Feb 14 12:45 inv_MMmatrix_TS3_20180209.sav
  • Need to change the name structure to seeing and Hz

13:10 Tried to load the file dist_flao1_ts4_atm_s1.0_L040.0_v15.0_ovfreq1000.00_sd3892.fits but get a error
  • IDL Invalid/Uncoherent input parameter
  • Says it is trying to read the correct file.
  • Will try manually tomorrow

Collect LBTI Bin 1 Reconstructors with KL_v22 basis set

10:00 Try LBTI data collection
  • in ~/.bashrc set umask 000

10:00 intmatAcquireGUI
  • IDL error file: permissions problem

10:09 intmatDisturbGUI
  • KL_v22, 20, 2, 10 Open loop
  • 030805

10:10 intmatAcquireGUI
  • Had to stop current disturbance by loading recon, set gain 0, close loop, set Disturbance "Disabled" on WFS Control GUI
  • 031247

10:17 intmatAnalyseGUI
  • 031613

10:20 Close loop
  • flat lbti_20180314_20mode_ol.sav
  • load it

10:27 Open Loop, 20 modes
  • dist 030805
  • acq 032650
  • rec 032818
  • closed loop
  • flat lbti_20180314_20modes_ol_2
  • apply it

10:54 Closed Loop, 20 modes
  • dist 030805
  • acq 035436 (gain 0.3)
  • rec 040038
  • closed loop
  • flat lbti_20180214_20modes_cl.sav (not much better)

11:10 66 modes
  • dist 041007 (amp file 672_0.3.fits)
  • acq 041308 (gain0.1.fits)
  • rec 041308
  • Closed Loop
  • flat lbti_20180214_66modes_cl (flatter, but still edges missing on each pupil)

11:25 153 modes
  • dist 042510 (optimized_amano_6l.fits)
  • acq 041601
  • rec 042810
  • closed loop
  • flat lbti_20180214_153modes.sav (much better, edges almost fixed)

10:38 400 modes
  • dist 043808 (optimized_amano_6l.fits)
  • acq 044018
  • rec 044203
  • closed loop
  • flat lbti_20140214_400modes.sav (a bit better)

Test TT forces in closed loop

In closed loop, turned off offloading to the hexapod so that any TT will remain on the Thin Shell. To force a TT, the bayside stage will be move a small amount. We will record the total force (min/max) from the Mirror GUI

12:05 closed loop, 400 modes, gain .2, 0.05, 0.05
  • cl_400m.mov
  • cl_400m.jpg

12:18 Start at X=-44.3
  • Offload on, delta cmd = -0.3 to +0.3 um
  • force = -0.2 to +0.26 N
  • Closed Loop, bin 1, 400 modes, 600 Hz:
    cl_400m.jpg
  • cl_400m.mov: Closed Loop, bin 1, 400 modes, 600 Hz

12:20 Move to X=46 (in 0.5 steps)
  • offload Off
  • Delta cmd limits set to -10 to +10 um
  • tip_10um.mov
  • tip_10um.jpg
  • tip_tilt_0um.jpg
  • Tip/Tilt 0 um:
    tip_tilt_0um.jpg
  • Tip 11 um:
    tip_10um.jpg
  • tip_10um.mov: Tip 11um

12:28 Start at Y=-36.3. Move to Y=-38.
  • offload off
  • tip_10um.mov
  • tip_10um.jpg
  • Tilt 11 um:
    tilt_10um.jpg
  • tilt_10um.mov: Tilt 11 um

-- %USERSIG{DougMiller - 2018-02-13}%
I Attachment Action Size Date Who Comment
cl_400m.jpgjpg cl_400m.jpg manage 300 K 14 Feb 2018 - 13:21 DougMiller Closed Loop, bin 1, 400 modes, 600 Hz
cl_400m.movmov cl_400m.mov manage 11 MB 14 Feb 2018 - 13:21 DougMiller Closed Loop, bin 1, 400 modes, 600 Hz
tilt_10um.jpgjpg tilt_10um.jpg manage 290 K 14 Feb 2018 - 13:23 DougMiller Tilt 13 um
tilt_10um.movmov tilt_10um.mov manage 18 MB 14 Feb 2018 - 13:23 DougMiller Tilt 13 um
tip_10um.jpgjpg tip_10um.jpg manage 293 K 14 Feb 2018 - 13:22 DougMiller Tip 13um
tip_10um.movmov tip_10um.mov manage 20 MB 14 Feb 2018 - 13:22 DougMiller Tip 13um
tip_tilt_0um.jpgjpg tip_tilt_0um.jpg manage 283 K 14 Feb 2018 - 13:22 DougMiller Tip 13um
Topic revision: r5 - 14 Feb 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