gamma=3) bf3Full = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='true level', gamma=3) bf4Full = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='true location', gamma=3) Lbf1Rem = L_p(bf1Rem.synthetic(4000, 1)) Lbf2Rem = L_p(bf2Rem.synthetic(4000, 1)) Lbf3Rem = L_p(bf3Rem.synthetic(4000, 1)) Lbf4Rem = L_p(bf4Rem.synthetic(4000, 1)) Lbf1Full = L_p(bf1Full.synthetic(4000, 1)) Lbf2Full = L_p(bf2Full.synthetic(4000, 1)) Lbf3Full = L_p(bf3Full.synthetic(4000, 1)) Lbf4Full = L_p(bf4Full.synthetic(4000, 1)) bca1Full = BeamformerCapon(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='classic') bca2Full = BeamformerCapon(freq_data=f,
eva32 = f32.eva[:] psf32 = PointSpreadFunction(grid=g, mpos=m, c=346.04, precision='float32') psf32Res = psf32.psf[:] psf64 = PointSpreadFunction(grid=g, mpos=m, c=346.04, precision='float64') psf64Res = psf64.psf[:] bb32 = BeamformerBase(freq_data=f, grid=g, mpos=m, r_diag=True, c=346.04, precision='float32') bb32Res = bb32.synthetic(cfreq,1) bb64 = BeamformerBase(freq_data=f, grid=g, mpos=m, r_diag=True, c=346.04, precision='float64') bb64Res = bb64.synthetic(cfreq,1) bf = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, gamma = 60, precision='float32') bfRes = bf.synthetic(cfreq,1) # 32 Bit PSF precision bd3232 = BeamformerDamas(beamformer=bb32, n_iter=100, psf_precision='float32') bd3232Res = bd3232.synthetic(cfreq,1) bc3232 = BeamformerClean(beamformer=bb32, psf_precision='float32') bc3232Res = bc3232.synthetic(cfreq,1) bdp3232 = BeamformerDamasPlus(beamformer=bb32, n_iter=100, psf_precision='float32') bdp3232Res = bdp3232.synthetic(cfreq,1) #64 Bit bd3264 = BeamformerDamas(beamformer=bb32, n_iter=100, psf_precision='float64') bd3264Res = bd3264.synthetic(cfreq,1) bc3264 = BeamformerClean(beamformer=bb32, psf_precision='float64') bc3264Res = bc3264.synthetic(cfreq,1) bdp3264 = BeamformerDamasPlus(beamformer=bb32, n_iter=100, psf_precision='float64')
Lbb3Rem = L_p(bb3Rem.synthetic(4000,1)) Lbb4Rem = L_p(bb4Rem.synthetic(4000,1)) Lbb1Full = L_p(bb1Full.synthetic(4000,1)) Lbb2Full = L_p(bb2Full.synthetic(4000,1)) Lbb3Full = L_p(bb3Full.synthetic(4000,1)) Lbb4Full = L_p(bb4Full.synthetic(4000,1)) bf1Rem = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=True, c=346.04, steer='classic', gamma=3) bf2Rem = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=True, c=346.04, steer='inverse', gamma=3) bf3Rem = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=True, c=346.04, steer='true level', gamma=3) bf4Rem = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=True, c=346.04, steer='true location', gamma=3) bf1Full = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='classic', gamma=3) bf2Full = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='inverse', gamma=3) bf3Full = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='true level', gamma=3) bf4Full = BeamformerFunctional(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='true location', gamma=3) Lbf1Rem = L_p(bf1Rem.synthetic(4000,1)) Lbf2Rem = L_p(bf2Rem.synthetic(4000,1)) Lbf3Rem = L_p(bf3Rem.synthetic(4000,1)) Lbf4Rem = L_p(bf4Rem.synthetic(4000,1)) Lbf1Full = L_p(bf1Full.synthetic(4000,1)) Lbf2Full = L_p(bf2Full.synthetic(4000,1)) Lbf3Full = L_p(bf3Full.synthetic(4000,1)) Lbf4Full = L_p(bf4Full.synthetic(4000,1)) bca1Full = BeamformerCapon(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='classic') bca2Full = BeamformerCapon(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='inverse') bca3Full = BeamformerCapon(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='true level') bca4Full = BeamformerCapon(freq_data=f, grid=g, mpos=m, r_diag=False, c=346.04, steer='true location') Lbca1Full = L_p(bca1Full.synthetic(4000,1)) Lbca2Full = L_p(bca2Full.synthetic(4000,1)) Lbca3Full = L_p(bca3Full.synthetic(4000,1))