psf3 = PointSpreadFunction(grid=g, mpos=m, c=346.04, steer='true level') psf4 = PointSpreadFunction(grid=g, mpos=m, c=346.04, steer='true location') Lpsf1 = L_p(psf1.psf[:]) Lpsf2 = L_p(psf2.psf[:]) Lpsf3 = L_p(psf3.psf[:]) Lpsf4 = L_p(psf4.psf[:]) bcpsf1Rem = BeamformerClean(beamformer=bb1Rem) bcpsf2Rem = BeamformerClean(beamformer=bb2Rem) bcpsf3Rem = BeamformerClean(beamformer=bb3Rem) bcpsf4Rem = BeamformerClean(beamformer=bb4Rem) bcpsf1Full = BeamformerClean(beamformer=bb1Full) bcpsf2Full = BeamformerClean(beamformer=bb2Full) bcpsf3Full = BeamformerClean(beamformer=bb3Full) bcpsf4Full = BeamformerClean(beamformer=bb4Full) Lbcpsf1Rem = L_p(bcpsf1Rem.synthetic(4000, 1)) Lbcpsf2Rem = L_p(bcpsf2Rem.synthetic(4000, 1)) Lbcpsf3Rem = L_p(bcpsf3Rem.synthetic(4000, 1)) Lbcpsf4Rem = L_p(bcpsf4Rem.synthetic(4000, 1)) Lbcpsf1Full = L_p(bcpsf1Full.synthetic(4000, 1)) Lbcpsf2Full = L_p(bcpsf2Full.synthetic(4000, 1)) Lbcpsf3Full = L_p(bcpsf3Full.synthetic(4000, 1)) Lbcpsf4Full = L_p(bcpsf4Full.synthetic(4000, 1)) bd1Rem = BeamformerDamas(beamformer=bb1Rem, n_iter=100) bd2Rem = BeamformerDamas(beamformer=bb2Rem, n_iter=100) bd3Rem = BeamformerDamas(beamformer=bb3Rem, n_iter=100) bd4Rem = BeamformerDamas(beamformer=bb4Rem, n_iter=100) bd1Full = BeamformerDamas(beamformer=bb1Full, n_iter=100) bd2Full = BeamformerDamas(beamformer=bb2Full, n_iter=100) bd3Full = BeamformerDamas(beamformer=bb3Full, n_iter=100)
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') bdp3264Res = bdp3264.synthetic(cfreq,1) # 32 Bit PSF precision bd6432 = BeamformerDamas(beamformer=bb64, n_iter=100, psf_precision='float32') bd6432Res = bd6432.synthetic(cfreq,1) bc6432 = BeamformerClean(beamformer=bb64, psf_precision='float32')
psf3 = PointSpreadFunction(grid=g, mpos=m, c=346.04, steer='true level') psf4 = PointSpreadFunction(grid=g, mpos=m, c=346.04, steer='true location') Lpsf1 = L_p(psf1.psf[:]) Lpsf2 = L_p(psf2.psf[:]) Lpsf3 = L_p(psf3.psf[:]) Lpsf4 = L_p(psf4.psf[:]) bcpsf1Rem = BeamformerClean(beamformer=bb1Rem) bcpsf2Rem = BeamformerClean(beamformer=bb2Rem) bcpsf3Rem = BeamformerClean(beamformer=bb3Rem) bcpsf4Rem = BeamformerClean(beamformer=bb4Rem) bcpsf1Full = BeamformerClean(beamformer=bb1Full) bcpsf2Full = BeamformerClean(beamformer=bb2Full) bcpsf3Full = BeamformerClean(beamformer=bb3Full) bcpsf4Full = BeamformerClean(beamformer=bb4Full) Lbcpsf1Rem = L_p(bcpsf1Rem.synthetic(4000,1)) Lbcpsf2Rem = L_p(bcpsf2Rem.synthetic(4000,1)) Lbcpsf3Rem = L_p(bcpsf3Rem.synthetic(4000,1)) Lbcpsf4Rem = L_p(bcpsf4Rem.synthetic(4000,1)) Lbcpsf1Full = L_p(bcpsf1Full.synthetic(4000,1)) Lbcpsf2Full = L_p(bcpsf2Full.synthetic(4000,1)) Lbcpsf3Full = L_p(bcpsf3Full.synthetic(4000,1)) Lbcpsf4Full = L_p(bcpsf4Full.synthetic(4000,1)) bd1Rem = BeamformerDamas(beamformer=bb1Rem, n_iter=100) bd2Rem = BeamformerDamas(beamformer=bb2Rem, n_iter=100) bd3Rem = BeamformerDamas(beamformer=bb3Rem, n_iter=100) bd4Rem = BeamformerDamas(beamformer=bb4Rem, n_iter=100) bd1Full = BeamformerDamas(beamformer=bb1Full, n_iter=100) bd2Full = BeamformerDamas(beamformer=bb2Full, n_iter=100) bd3Full = BeamformerDamas(beamformer=bb3Full, n_iter=100)