Mao_Chengfeng_lab12.ipynb Bai_Hongyang_lab12.ipynb: 0085 y_e = 38.297 Bai_Hongyang_lab12.ipynb: 0089 x_e = 142.372 Bai_Hongyang_lab12.ipynb: 0107 return J Li_Boyu_Lab12.ipynb: 0017 print cov_x feddersen_matthew_lab12.ipynb: 0101 J = [] Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0017 print cov_x snyder_bradley_Lab12_FA2014.ipynb: 0017 print cov_x Ling_Wang-Jung_Lab12.ipynb: 0037 for i in range(10): manuel_devasia_Lab12.ipynb: 0080 delta_t = 30 deVasconcelos_Thiago_Lab12.ipynb: 0055 for i in range(len(lats)): deVasconcelos_Thiago_Lab12.ipynb: 0080 delta_t = 30 Lei_TsungHan_Lab12.ipynb: 0037 for i in range(10): Lei_TsungHan_Lab12.ipynb: 0055 for i in range(len(lats)): Lei_TsungHan_Lab12.ipynb: 0107 return J vashistha_siddhant_lab12.ipynb: 0085 y_e = 38.297 survila2_Lab12_FA2014.ipynb: 0085 y_e = 38.297 survila2_Lab12_FA2014.ipynb: 0089 x_e = 142.372 survila2_Lab12_FA2014.ipynb: 0107 return J Wu_Yi-Luen_Lab12.ipynb: 0107 return J kabaghe_chuma_lab12.ipynb: 0012 x_e = mean(lons) kabaghe_chuma_lab12.ipynb: 0037 for i in range(10): kabaghe_chuma_lab12.ipynb: 0075 y_e = mean(lats) kabaghe_chuma_lab12.ipynb: 0085 y_e = 38.297 kabaghe_chuma_lab12.ipynb: 0089 x_e = 142.372 Bai_Hongyang_lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0086 y_e = 38.297 Mao_Chengfeng_lab12.ipynb: 0087 x_e = 142.372 Mao_Chengfeng_lab12.ipynb: 0117 return J Shih_ChienTung_Lab12.ipynb: 0088 print J Li_Boyu_Lab12.ipynb: 0088 print J tang_yuhan_lab12.ipynb: 0088 print J snyder_bradley_Lab12_FA2014.ipynb: 0088 print J Maier_Christopher_Lab12_FA2014.ipynb: 0017 y_e = 38.322 Maier_Christopher_Lab12_FA2014.ipynb: 0107 x_e = 142.369 Deng_Fei_lab12_SEENPZ.ipynb: 0088 print J Ling_Wang-Jung_Lab12.ipynb: 0088 print J McLaughlin_Sean_Lab12_FA2014.ipynb: 0119 figure(figsize=(10,10)) Lei_TsungHan_Lab12.ipynb: 0117 return J Wang_Hongru_Lab12.ipynb: 0088 print J vashistha_siddhant_lab12.ipynb: 0086 y_e = 38.297 survila2_Lab12_FA2014.ipynb: 0017 y_e = 38.322 survila2_Lab12_FA2014.ipynb: 0086 y_e = 38.297 survila2_Lab12_FA2014.ipynb: 0087 x_e = 142.372 survila2_Lab12_FA2014.ipynb: 0107 x_e = 142.369 survila2_Lab12_FA2014.ipynb: 0117 return J Wu_Yi-Luen_Lab12.ipynb: 0117 return J kabaghe_chuma_lab12.ipynb: 0086 y_e = 38.297 kabaghe_chuma_lab12.ipynb: 0087 x_e = 142.372 crawford_melissa_macrawf2_lab12.ipynb Shih_ChienTung_Lab12.ipynb: 0100 print "" Antonacci_Dominic_genetic_algorithm.ipynb: 0114 return 1 Antonacci_Dominic_lab12_SEENPZ.ipynb: 0114 return 1 Wu_Yi-Luen_Lab12.ipynb: 0100 print "" Lab12_CheatDetect.ipynb Shih_ChienTung_Lab12.ipynb Bai_Hongyang_lab12.ipynb: 0006 print J crawford_melissa_macrawf2_lab12.ipynb: 0037 print "" Li_Boyu_Lab12.ipynb: 0006 print J feddersen_matthew_lab12.ipynb: 0014 print sigma tang_yuhan_lab12.ipynb: 0006 print J snyder_bradley_Lab12_FA2014.ipynb: 0006 print J Deng_Fei_lab12_SEENPZ.ipynb: 0006 print J Ling_Wang-Jung_Lab12.ipynb: 0006 print J Ling_Wang-Jung_Lab12.ipynb: 0014 print sigma Ling_Wang-Jung_Lab12.ipynb: 0050 v = 1.39 Ling_Wang-Jung_Lab12.ipynb: 0064 for j in range(4): Ling_Wang-Jung_Lab12.ipynb: 0071 print cov deVasconcelos_Thiago_Lab12.ipynb: 0070 figure(figsize=(6,6)) Lei_TsungHan_Lab12.ipynb: 0070 figure(figsize=(6,6)) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0070 figure(figsize=(6,6)) Hsu_HengWei_Lab12_FA2014.ipynb: 0004 sigma = (3.3*3.3/12)**0.5 Hsu_HengWei_Lab12_FA2014.ipynb: 0005 xe = sum(lons)/len(lons) Hsu_HengWei_Lab12_FA2014.ipynb: 0008 elif j == 1: Hsu_HengWei_Lab12_FA2014.ipynb: 0011 J[i][j] = k Hsu_HengWei_Lab12_FA2014.ipynb: 0012 pv = [] Hsu_HengWei_Lab12_FA2014.ipynb: 0014 print sigma Hsu_HengWei_Lab12_FA2014.ipynb: 0015 if 143.55>i>141.45: Hsu_HengWei_Lab12_FA2014.ipynb: 0019 print "Uncertainty of time arrival:",sigma Hsu_HengWei_Lab12_FA2014.ipynb: 0025 sigma = (30*30/12)**0.5 Hsu_HengWei_Lab12_FA2014.ipynb: 0029 for k in px: Hsu_HengWei_Lab12_FA2014.ipynb: 0034 cov = [[0 for j in range(4)] for i in range(4)] Hsu_HengWei_Lab12_FA2014.ipynb: 0038 i += 1 Hsu_HengWei_Lab12_FA2014.ipynb: 0040 py = [] Hsu_HengWei_Lab12_FA2014.ipynb: 0042 elif j == 2: Hsu_HengWei_Lab12_FA2014.ipynb: 0043 ye = sum(lats)/len(lats) Hsu_HengWei_Lab12_FA2014.ipynb: 0047 sigma = (17*17/12)**0.5 Hsu_HengWei_Lab12_FA2014.ipynb: 0050 v = 1.39 Hsu_HengWei_Lab12_FA2014.ipynb: 0053 pt = [] Hsu_HengWei_Lab12_FA2014.ipynb: 0061 if j == 0: Hsu_HengWei_Lab12_FA2014.ipynb: 0064 for j in range(4): Hsu_HengWei_Lab12_FA2014.ipynb: 0065 px = [] Hsu_HengWei_Lab12_FA2014.ipynb: 0073 while i < len(lats): Hsu_HengWei_Lab12_FA2014.ipynb: 0076 plot( 142.372,38.297,'b.') Hsu_HengWei_Lab12_FA2014.ipynb: 0077 J = np.zeros((len(px),4)) Hsu_HengWei_Lab12_FA2014.ipynb: 0078 else : Wang_Hongru_Lab12.ipynb: 0006 print J vashistha_siddhant_lab12.ipynb: 0070 figure(figsize=(6,6)) Wu_Yi-Luen_Lab12.ipynb: 0037 print "" Wu_Yi-Luen_Lab12.ipynb: 0070 figure(figsize=(6,6)) Wu_Yi-Luen_Lab12.ipynb: 0071 print cov Li_Boyu_Lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0034 print cov_x Bai_Hongyang_lab12.ipynb: 0004 print J Shih_ChienTung_Lab12.ipynb: 0004 print J Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0011 print sigma_t Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0034 print cov_x tang_yuhan_lab12.ipynb: 0004 print J snyder_bradley_Lab12_FA2014.ipynb: 0004 print J snyder_bradley_Lab12_FA2014.ipynb: 0034 print cov_x Deng_Fei_lab12_SEENPZ.ipynb: 0004 print J Ling_Wang-Jung_Lab12.ipynb: 0004 print J deVasconcelos_Thiago_Lab12.ipynb: 0008 lats = [] deVasconcelos_Thiago_Lab12.ipynb: 0033 lons = [] Antonacci_Dominic_lab12_SEENPZ.ipynb: 0008 lats = [] Antonacci_Dominic_lab12_SEENPZ.ipynb: 0033 lons = [] Wang_Hongru_Lab12.ipynb: 0004 print J Wang_Xinao_Lab12_SEENPZ.ipynb: 0008 lats = [] Wang_Xinao_Lab12_SEENPZ.ipynb: 0033 lons = [] kabaghe_chuma_lab12.ipynb: 0007 print N feddersen_matthew_lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0085 J = [] Shih_ChienTung_Lab12.ipynb: 0049 print sigma Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0019 #print J tang_yuhan_lab12.ipynb: 0023 xe = 142.372 tang_yuhan_lab12.ipynb: 0068 ye = 38.297 Ling_Wang-Jung_Lab12.ipynb: 0019 #print J Ling_Wang-Jung_Lab12.ipynb: 0049 print sigma Hsu_HengWei_Lab12_FA2014.ipynb: 0023 xe = 142.372 Hsu_HengWei_Lab12_FA2014.ipynb: 0049 print sigma Hsu_HengWei_Lab12_FA2014.ipynb: 0068 ye = 38.297 Wang_Hongru_Lab12.ipynb: 0023 xe = 142.372 Wang_Hongru_Lab12.ipynb: 0068 ye = 38.297 Wang_Xinao_Lab12_SEENPZ.ipynb: 0023 xe = 142.372 Wang_Xinao_Lab12_SEENPZ.ipynb: 0068 ye = 38.297 Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb Mao_Chengfeng_lab12.ipynb: 0018 print cov_x Li_Boyu_Lab12.ipynb: 0018 print cov_x Li_Boyu_Lab12.ipynb: 0037 print sigma_t feddersen_matthew_lab12.ipynb: 0036 #print J snyder_bradley_Lab12_FA2014.ipynb: 0018 print cov_x snyder_bradley_Lab12_FA2014.ipynb: 0093 N = 100 snyder_bradley_Lab12_FA2014.ipynb: 0114 N = 10 Maier_Christopher_Lab12_FA2014.ipynb: 0114 N = 10 Deng_Fei_lab12_SEENPZ.ipynb: 0026 x = [] Deng_Fei_lab12_SEENPZ.ipynb: 0051 y = [] Ling_Wang-Jung_Lab12.ipynb: 0036 #print J manuel_devasia_Lab12.ipynb: 0093 N = 100 McLaughlin_Sean_Lab12_FA2014.ipynb: 0114 N = 10 Vankayala_Anuraag_Lab12.ipynb: 0102 print sampling_rate Antonacci_Dominic_genetic_algorithm.ipynb crawford_melissa_macrawf2_lab12.ipynb: 0101 return 1 McLaughlin_Sean_Lab12_FA2014.ipynb: 0054 v = 4 #km/s Antonacci_Dominic_lab12_SEENPZ.ipynb: 0002 # Velocity is in km/s Antonacci_Dominic_lab12_SEENPZ.ipynb: 0003 # Jacobian result is in degrees of latitude and longitude Antonacci_Dominic_lab12_SEENPZ.ipynb: 0009 dtdv(lats[i], lons[i], lat_e, lon_e, v)]) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0010 # Iterate over all lats and lons Antonacci_Dominic_lab12_SEENPZ.ipynb: 0032 J = computeJacobian(lats, lons, lat_e, lon_e, v) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0038 #return (lon_i-lon_e)*111*cos(lat_e * pi / 180.) / (v * 6371 * distance_on_unit_sphere(lat_i, lon_i, lat_e, lon_e)) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0042 jacobian = [] Antonacci_Dominic_lab12_SEENPZ.ipynb: 0048 # Compute the uncertainty of the estimate Antonacci_Dominic_lab12_SEENPZ.ipynb: 0049 #return (lat_i-lat_e)*111 / (v * 6371 * distance_on_unit_sphere(lat_i, lon_i, lat_e, lon_e)) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0050 return -sqrt( (lat_i - lat_e)**2 + (lon_i - lon_e)**2) / (v/111.)**2 Antonacci_Dominic_lab12_SEENPZ.ipynb: 0054 v = 4 #km/s Antonacci_Dominic_lab12_SEENPZ.ipynb: 0057 for i in xrange(len(lats)): Antonacci_Dominic_lab12_SEENPZ.ipynb: 0059 return (lat_i - lat_e) / (v/111. * sqrt( (lat_i - lat_e)**2 + (lon_i - lon_e)**2)) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0065 def dtdt(lat_i, lon_i, lat_e, lon_e, v): Antonacci_Dominic_lab12_SEENPZ.ipynb: 0070 def dtdv(lat_i, lon_i, lat_e, lon_e, v): Antonacci_Dominic_lab12_SEENPZ.ipynb: 0073 unc = unc[ind] Antonacci_Dominic_lab12_SEENPZ.ipynb: 0074 # Computes the Jacobian given receivers at given latitudes and longitudes, for an estimated epicenter in latitude and longitude Antonacci_Dominic_lab12_SEENPZ.ipynb: 0076 # Define a new Jacobian, better matching what is being used to grade it Antonacci_Dominic_lab12_SEENPZ.ipynb: 0077 jacobian.append([ dtdx(lats[i], lons[i], lat_e, lon_e, v), Antonacci_Dominic_lab12_SEENPZ.ipynb: 0078 def computeJacobian(lats, lons, lat_e, lon_e, v): Antonacci_Dominic_lab12_SEENPZ.ipynb: 0083 lon_e = 142.372 #degrees longitude Antonacci_Dominic_lab12_SEENPZ.ipynb: 0094 dtdy(lats[i], lons[i], lat_e, lon_e, v), Antonacci_Dominic_lab12_SEENPZ.ipynb: 0100 lat_e = 38.297 #degrees latitude Antonacci_Dominic_lab12_SEENPZ.ipynb: 0101 return 1 Antonacci_Dominic_lab12_SEENPZ.ipynb: 0102 dtdt(lats[i], lons[i], lat_e, lon_e, v), Antonacci_Dominic_lab12_SEENPZ.ipynb: 0104 def dtdy(lat_i, lon_i, lat_e, lon_e, v): Antonacci_Dominic_lab12_SEENPZ.ipynb: 0108 return (lon_i - lon_e) / (v/111. * sqrt( (lat_i - lat_e)**2 + (lon_i - lon_e)**2)) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0116 return np.array(jacobian) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0126 def dtdx(lat_i, lon_i, lat_e, lon_e, v): Antonacci_Dominic_lab12_SEENPZ.ipynb: 0131 #return -6371 * distance_on_unit_sphere(lat_i, lon_i, lat_e, lon_e) / v**2 Wang_Hongru_Lab12.ipynb: 0042 jacobian = [] Vankayala_Anuraag_Lab12.ipynb: 0057 for i in xrange(len(lats)): tang_yuhan_lab12.ipynb Bai_Hongyang_lab12.ipynb: 0094 print J Shih_ChienTung_Lab12.ipynb: 0094 print J Li_Boyu_Lab12.ipynb: 0094 print J feddersen_matthew_lab12.ipynb: 0034 xe = 142.372 feddersen_matthew_lab12.ipynb: 0098 ye = 38.297 snyder_bradley_Lab12_FA2014.ipynb: 0094 print J Deng_Fei_lab12_SEENPZ.ipynb: 0004 center_x = np.mean(lons) Deng_Fei_lab12_SEENPZ.ipynb: 0083 while True: Deng_Fei_lab12_SEENPZ.ipynb: 0094 print J Deng_Fei_lab12_SEENPZ.ipynb: 0123 center_y = np.mean(lats) Ling_Wang-Jung_Lab12.ipynb: 0094 print J Ling_Wang-Jung_Lab12.ipynb: 0116 for i in range(N): LaBelle_Julia_lab12.ipynb: 0084 figure(figsize=(15,15)) Hsu_HengWei_Lab12_FA2014.ipynb: 0034 xe = 142.372 Hsu_HengWei_Lab12_FA2014.ipynb: 0098 ye = 38.297 Wang_Hongru_Lab12.ipynb: 0034 xe = 142.372 Wang_Hongru_Lab12.ipynb: 0094 print J Wang_Hongru_Lab12.ipynb: 0098 ye = 38.297 Wang_Xinao_Lab12_SEENPZ.ipynb: 0034 xe = 142.372 Wang_Xinao_Lab12_SEENPZ.ipynb: 0098 ye = 38.297 Wu_Yi-Luen_Lab12.ipynb: 0116 for i in range(N): kabaghe_chuma_lab12.ipynb: 0116 for i in range(N): snyder_bradley_Lab12_FA2014.ipynb Mao_Chengfeng_lab12.ipynb: 0015 print cov_x Bai_Hongyang_lab12.ipynb: 0087 print J Shih_ChienTung_Lab12.ipynb: 0087 print J Li_Boyu_Lab12.ipynb: 0015 print cov_x Li_Boyu_Lab12.ipynb: 0087 print J Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0008 N = 100 Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0015 print cov_x Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0073 N = 10 tang_yuhan_lab12.ipynb: 0087 print J Maier_Christopher_Lab12_FA2014.ipynb: 0073 N = 10 Deng_Fei_lab12_SEENPZ.ipynb: 0087 print J Ling_Wang-Jung_Lab12.ipynb: 0087 print J manuel_devasia_Lab12.ipynb: 0008 N = 100 McLaughlin_Sean_Lab12_FA2014.ipynb: 0073 N = 10 LaBelle_Julia_lab12.ipynb: 0020 xe = 142.369 LaBelle_Julia_lab12.ipynb: 0041 ye = 38.322 Wang_Hongru_Lab12.ipynb: 0087 print J vashistha_siddhant_lab12.ipynb: 0097 for i in arange(N): Maier_Christopher_Lab12_FA2014.ipynb Bai_Hongyang_lab12.ipynb: 0013 y_e = 38.322 Bai_Hongyang_lab12.ipynb: 0091 x_e = 142.369 Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0067 N = 10 snyder_bradley_Lab12_FA2014.ipynb: 0067 N = 10 McLaughlin_Sean_Lab12_FA2014.ipynb: 0067 N = 10 vashistha_siddhant_lab12.ipynb: 0026 legend() survila2_Lab12_FA2014.ipynb: 0013 y_e = 38.322 survila2_Lab12_FA2014.ipynb: 0091 x_e = 142.369 Deng_Fei_lab12_SEENPZ.ipynb Bai_Hongyang_lab12.ipynb: 0102 print J Shih_ChienTung_Lab12.ipynb: 0102 print J Li_Boyu_Lab12.ipynb: 0102 print J Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0048 x = [] Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0054 y = [] tang_yuhan_lab12.ipynb: 0092 while True: tang_yuhan_lab12.ipynb: 0102 print J tang_yuhan_lab12.ipynb: 0116 center_x = np.mean(lons) tang_yuhan_lab12.ipynb: 0129 center_y = np.mean(lats) snyder_bradley_Lab12_FA2014.ipynb: 0102 print J Ling_Wang-Jung_Lab12.ipynb: 0102 print J Wang_Hongru_Lab12.ipynb: 0088 for i in xrange(N): Wang_Hongru_Lab12.ipynb: 0102 print J Wu_Yi-Luen_Lab12.ipynb: 0035 N = len(x) Ling_Wang-Jung_Lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0034 for i in range(10): Bai_Hongyang_lab12.ipynb: 0076 print J Shih_ChienTung_Lab12.ipynb: 0061 print sigma Shih_ChienTung_Lab12.ipynb: 0065 v = 1.39 Shih_ChienTung_Lab12.ipynb: 0076 print J Shih_ChienTung_Lab12.ipynb: 0099 for j in range(4): Shih_ChienTung_Lab12.ipynb: 0100 print cov Li_Boyu_Lab12.ipynb: 0076 print J feddersen_matthew_lab12.ipynb: 0030 #print J feddersen_matthew_lab12.ipynb: 0061 print sigma Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0030 #print J tang_yuhan_lab12.ipynb: 0076 print J tang_yuhan_lab12.ipynb: 0095 for i in range(N): snyder_bradley_Lab12_FA2014.ipynb: 0076 print J Deng_Fei_lab12_SEENPZ.ipynb: 0076 print J manuel_devasia_Lab12.ipynb: 0055 plt.show() McLaughlin_Sean_Lab12_FA2014.ipynb: 0055 plt.show() Lei_TsungHan_Lab12.ipynb: 0034 for i in range(10): Hsu_HengWei_Lab12_FA2014.ipynb: 0061 print sigma Hsu_HengWei_Lab12_FA2014.ipynb: 0065 v = 1.39 Hsu_HengWei_Lab12_FA2014.ipynb: 0099 for j in range(4): Wang_Hongru_Lab12.ipynb: 0076 print J Wu_Yi-Luen_Lab12.ipynb: 0095 for i in range(N): Wu_Yi-Luen_Lab12.ipynb: 0100 print cov kabaghe_chuma_lab12.ipynb: 0034 for i in range(10): kabaghe_chuma_lab12.ipynb: 0095 for i in range(N): manuel_devasia_Lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0059 delta_t = 30 Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0061 N = 100 snyder_bradley_Lab12_FA2014.ipynb: 0061 N = 100 Ling_Wang-Jung_Lab12.ipynb: 0035 plt.show() McLaughlin_Sean_Lab12_FA2014.ipynb: 0026 %matplotlib inline McLaughlin_Sean_Lab12_FA2014.ipynb: 0035 plt.show() deVasconcelos_Thiago_Lab12.ipynb: 0059 delta_t = 30 Gleason_Cole_Lab12_ALSO_HAS_NPZ_IF_NEEDED.ipynb: 0026 %matplotlib inline Gleason_Cole_Lab12_ALSO_HAS_NPZ_IF_NEEDED.ipynb: 0028 lats = np.array(npzfile['lats']) Gleason_Cole_Lab12_ALSO_HAS_NPZ_IF_NEEDED.ipynb: 0047 npzfile = np.load('receiver_positions.npz') Gleason_Cole_Lab12_ALSO_HAS_NPZ_IF_NEEDED.ipynb: 0085 lons = np.array(npzfile['lons']) McLaughlin_Sean_Lab12_FA2014.ipynb Bai_Hongyang_lab12.ipynb: 0215 figure(figsize=(10,10)) Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0137 N = 10 Antonacci_Dominic_genetic_algorithm.ipynb: 0081 v = 4 #km/s snyder_bradley_Lab12_FA2014.ipynb: 0137 N = 10 Maier_Christopher_Lab12_FA2014.ipynb: 0137 N = 10 Ling_Wang-Jung_Lab12.ipynb: 0087 plt.show() manuel_devasia_Lab12.ipynb: 0067 %matplotlib inline manuel_devasia_Lab12.ipynb: 0087 plt.show() Antonacci_Dominic_lab12_SEENPZ.ipynb: 0021 cos = (math.sin(phi1)*math.sin(phi2)*math.cos(theta1 - theta2) + Antonacci_Dominic_lab12_SEENPZ.ipynb: 0028 # phi = 90 - latitude Antonacci_Dominic_lab12_SEENPZ.ipynb: 0030 # theta = longitude Antonacci_Dominic_lab12_SEENPZ.ipynb: 0035 arc = math.acos( cos ) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0038 phi2 = (90.0 - lat2)*degrees_to_radians Antonacci_Dominic_lab12_SEENPZ.ipynb: 0039 theta2 = long2*degrees_to_radians Antonacci_Dominic_lab12_SEENPZ.ipynb: 0040 # spherical coordinates in radians. Antonacci_Dominic_lab12_SEENPZ.ipynb: 0064 # Compute spherical distance from spherical coordinates. Antonacci_Dominic_lab12_SEENPZ.ipynb: 0068 theta1 = long1*degrees_to_radians Antonacci_Dominic_lab12_SEENPZ.ipynb: 0078 # sin phi sin phi' cos(theta-theta') + cos phi cos phi' Antonacci_Dominic_lab12_SEENPZ.ipynb: 0081 v = 4 #km/s Antonacci_Dominic_lab12_SEENPZ.ipynb: 0103 # (1, theta, phi) and (1, theta, phi) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0109 degrees_to_radians = math.pi/180.0 Antonacci_Dominic_lab12_SEENPZ.ipynb: 0118 phi1 = (90.0 - lat1)*degrees_to_radians Antonacci_Dominic_lab12_SEENPZ.ipynb: 0124 # For two locations in spherical coordinates Antonacci_Dominic_lab12_SEENPZ.ipynb: 0146 math.cos(phi1)*math.cos(phi2)) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0164 for i in xrange(10): Antonacci_Dominic_lab12_SEENPZ.ipynb: 0170 # Convert latitude and longitude to Antonacci_Dominic_lab12_SEENPZ.ipynb: 0181 # cosine( arc length ) = Antonacci_Dominic_lab12_SEENPZ.ipynb: 0200 # distance = rho * arc length Gleason_Cole_Lab12_ALSO_HAS_NPZ_IF_NEEDED.ipynb: 0067 %matplotlib inline deVasconcelos_Thiago_Lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0050 for i in range(len(lats)): Mao_Chengfeng_lab12.ipynb: 0079 delta_t = 30 Shih_ChienTung_Lab12.ipynb: 0052 figure(figsize=(6,6)) Li_Boyu_Lab12.ipynb: 0046 lons = [] Li_Boyu_Lab12.ipynb: 0060 lats = [] manuel_devasia_Lab12.ipynb: 0079 delta_t = 30 Lei_TsungHan_Lab12.ipynb: 0050 for i in range(len(lats)): Lei_TsungHan_Lab12.ipynb: 0052 figure(figsize=(6,6)) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0046 lons = [] Antonacci_Dominic_lab12_SEENPZ.ipynb: 0052 figure(figsize=(6,6)) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0060 lats = [] Wang_Xinao_Lab12_SEENPZ.ipynb: 0046 lons = [] Wang_Xinao_Lab12_SEENPZ.ipynb: 0060 lats = [] vashistha_siddhant_lab12.ipynb: 0052 figure(figsize=(6,6)) Wu_Yi-Luen_Lab12.ipynb: 0052 figure(figsize=(6,6)) Lei_TsungHan_Lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0030 for i in range(10): Mao_Chengfeng_lab12.ipynb: 0060 return J Mao_Chengfeng_lab12.ipynb: 0064 for i in range(len(lats)): Bai_Hongyang_lab12.ipynb: 0060 return J Shih_ChienTung_Lab12.ipynb: 0058 figure(figsize=(6,6)) Ling_Wang-Jung_Lab12.ipynb: 0030 for i in range(10): deVasconcelos_Thiago_Lab12.ipynb: 0058 figure(figsize=(6,6)) deVasconcelos_Thiago_Lab12.ipynb: 0064 for i in range(len(lats)): Antonacci_Dominic_lab12_SEENPZ.ipynb: 0058 figure(figsize=(6,6)) vashistha_siddhant_lab12.ipynb: 0058 figure(figsize=(6,6)) survila2_Lab12_FA2014.ipynb: 0060 return J Wu_Yi-Luen_Lab12.ipynb: 0058 figure(figsize=(6,6)) Wu_Yi-Luen_Lab12.ipynb: 0060 return J kabaghe_chuma_lab12.ipynb: 0030 for i in range(10): LaBelle_Julia_lab12.ipynb tang_yuhan_lab12.ipynb: 0014 figure(figsize=(15,15)) snyder_bradley_Lab12_FA2014.ipynb: 0019 xe = 142.369 snyder_bradley_Lab12_FA2014.ipynb: 0040 ye = 38.322 Antonacci_Dominic_lab12_SEENPZ.ipynb crawford_melissa_macrawf2_lab12.ipynb: 0114 return 1 Shih_ChienTung_Lab12.ipynb: 0067 figure(figsize=(6,6)) Li_Boyu_Lab12.ipynb: 0057 lons = [] Li_Boyu_Lab12.ipynb: 0074 lats = [] Antonacci_Dominic_genetic_algorithm.ipynb: 0001 # Velocity is in km/s Antonacci_Dominic_genetic_algorithm.ipynb: 0003 # Jacobian result is in degrees of latitude and longitude Antonacci_Dominic_genetic_algorithm.ipynb: 0007 dtdv(lats[i], lons[i], lat_e, lon_e, v)]) Antonacci_Dominic_genetic_algorithm.ipynb: 0008 # Iterate over all lats and lons Antonacci_Dominic_genetic_algorithm.ipynb: 0020 lat_e = 38.297 #degrees latitude Antonacci_Dominic_genetic_algorithm.ipynb: 0025 def dtdv(lat_i, lon_i, lat_e, lon_e, v): Antonacci_Dominic_genetic_algorithm.ipynb: 0032 J = computeJacobian(lats, lons, lat_e, lon_e, v) Antonacci_Dominic_genetic_algorithm.ipynb: 0035 #return (lon_i-lon_e)*111*cos(lat_e * pi / 180.) / (v * 6371 * distance_on_unit_sphere(lat_i, lon_i, lat_e, lon_e)) Antonacci_Dominic_genetic_algorithm.ipynb: 0037 jacobian = [] Antonacci_Dominic_genetic_algorithm.ipynb: 0044 # Compute the uncertainty of the estimate Antonacci_Dominic_genetic_algorithm.ipynb: 0045 #return (lat_i-lat_e)*111 / (v * 6371 * distance_on_unit_sphere(lat_i, lon_i, lat_e, lon_e)) Antonacci_Dominic_genetic_algorithm.ipynb: 0049 return -sqrt( (lat_i - lat_e)**2 + (lon_i - lon_e)**2) / (v/111.)**2 Antonacci_Dominic_genetic_algorithm.ipynb: 0058 v = 4 #km/s Antonacci_Dominic_genetic_algorithm.ipynb: 0061 for i in xrange(len(lats)): Antonacci_Dominic_genetic_algorithm.ipynb: 0065 jacobian.append([ dtdx(lats[i], lons[i], lat_e, lon_e, v), Antonacci_Dominic_genetic_algorithm.ipynb: 0066 return (lat_i - lat_e) / (v/111. * sqrt( (lat_i - lat_e)**2 + (lon_i - lon_e)**2)) Antonacci_Dominic_genetic_algorithm.ipynb: 0071 return np.array(jacobian) Antonacci_Dominic_genetic_algorithm.ipynb: 0072 def dtdt(lat_i, lon_i, lat_e, lon_e, v): Antonacci_Dominic_genetic_algorithm.ipynb: 0081 unc = unc[ind] Antonacci_Dominic_genetic_algorithm.ipynb: 0082 # Computes the Jacobian given receivers at given latitudes and longitudes, for an estimated epicenter in latitude and longitude Antonacci_Dominic_genetic_algorithm.ipynb: 0085 # Define a new Jacobian, better matching what is being used to grade it Antonacci_Dominic_genetic_algorithm.ipynb: 0089 def computeJacobian(lats, lons, lat_e, lon_e, v): Antonacci_Dominic_genetic_algorithm.ipynb: 0094 lon_e = 142.372 #degrees longitude Antonacci_Dominic_genetic_algorithm.ipynb: 0109 dtdy(lats[i], lons[i], lat_e, lon_e, v), Antonacci_Dominic_genetic_algorithm.ipynb: 0114 return 1 Antonacci_Dominic_genetic_algorithm.ipynb: 0115 dtdt(lats[i], lons[i], lat_e, lon_e, v), Antonacci_Dominic_genetic_algorithm.ipynb: 0116 def dtdy(lat_i, lon_i, lat_e, lon_e, v): Antonacci_Dominic_genetic_algorithm.ipynb: 0119 return (lon_i - lon_e) / (v/111. * sqrt( (lat_i - lat_e)**2 + (lon_i - lon_e)**2)) Antonacci_Dominic_genetic_algorithm.ipynb: 0130 def dtdx(lat_i, lon_i, lat_e, lon_e, v): Antonacci_Dominic_genetic_algorithm.ipynb: 0137 #return -6371 * distance_on_unit_sphere(lat_i, lon_i, lat_e, lon_e) / v**2 McLaughlin_Sean_Lab12_FA2014.ipynb: 0015 cos = (math.sin(phi1)*math.sin(phi2)*math.cos(theta1 - theta2) + McLaughlin_Sean_Lab12_FA2014.ipynb: 0016 # Convert latitude and longitude to McLaughlin_Sean_Lab12_FA2014.ipynb: 0019 # phi = 90 - latitude McLaughlin_Sean_Lab12_FA2014.ipynb: 0021 # theta = longitude McLaughlin_Sean_Lab12_FA2014.ipynb: 0024 arc = math.acos( cos ) McLaughlin_Sean_Lab12_FA2014.ipynb: 0027 phi2 = (90.0 - lat2)*degrees_to_radians McLaughlin_Sean_Lab12_FA2014.ipynb: 0028 theta2 = long2*degrees_to_radians McLaughlin_Sean_Lab12_FA2014.ipynb: 0030 # spherical coordinates in radians. McLaughlin_Sean_Lab12_FA2014.ipynb: 0039 # Compute spherical distance from spherical coordinates. McLaughlin_Sean_Lab12_FA2014.ipynb: 0047 theta1 = long1*degrees_to_radians McLaughlin_Sean_Lab12_FA2014.ipynb: 0055 # sin phi sin phi' cos(theta-theta') + cos phi cos phi' McLaughlin_Sean_Lab12_FA2014.ipynb: 0058 v = 4 #km/s McLaughlin_Sean_Lab12_FA2014.ipynb: 0063 # For two locations in spherical coordinates McLaughlin_Sean_Lab12_FA2014.ipynb: 0076 # (1, theta, phi) and (1, theta, phi) McLaughlin_Sean_Lab12_FA2014.ipynb: 0077 degrees_to_radians = math.pi/180.0 McLaughlin_Sean_Lab12_FA2014.ipynb: 0080 phi1 = (90.0 - lat1)*degrees_to_radians McLaughlin_Sean_Lab12_FA2014.ipynb: 0112 for i in xrange(10): McLaughlin_Sean_Lab12_FA2014.ipynb: 0124 # cosine( arc length ) = McLaughlin_Sean_Lab12_FA2014.ipynb: 0135 # distance = rho * arc length McLaughlin_Sean_Lab12_FA2014.ipynb: 0140 math.cos(phi1)*math.cos(phi2)) deVasconcelos_Thiago_Lab12.ipynb: 0057 lons = [] deVasconcelos_Thiago_Lab12.ipynb: 0067 figure(figsize=(6,6)) deVasconcelos_Thiago_Lab12.ipynb: 0074 lats = [] Lei_TsungHan_Lab12.ipynb: 0067 figure(figsize=(6,6)) Wang_Hongru_Lab12.ipynb: 0037 jacobian = [] Vankayala_Anuraag_Lab12.ipynb: 0061 for i in xrange(len(lats)): Wang_Xinao_Lab12_SEENPZ.ipynb: 0057 lons = [] Wang_Xinao_Lab12_SEENPZ.ipynb: 0074 lats = [] vashistha_siddhant_lab12.ipynb: 0067 figure(figsize=(6,6)) Wu_Yi-Luen_Lab12.ipynb: 0067 figure(figsize=(6,6)) Hsu_HengWei_Lab12_FA2014.ipynb Shih_ChienTung_Lab12.ipynb: 0002 sigma = (3.3*3.3/12)**0.5 Shih_ChienTung_Lab12.ipynb: 0003 xe = sum(lons)/len(lons) Shih_ChienTung_Lab12.ipynb: 0006 elif j == 1: Shih_ChienTung_Lab12.ipynb: 0010 J[i][j] = k Shih_ChienTung_Lab12.ipynb: 0012 pv = [] Shih_ChienTung_Lab12.ipynb: 0015 print sigma Shih_ChienTung_Lab12.ipynb: 0016 if 143.55>i>141.45: Shih_ChienTung_Lab12.ipynb: 0022 print "Uncertainty of time arrival:",sigma Shih_ChienTung_Lab12.ipynb: 0024 sigma = (30*30/12)**0.5 Shih_ChienTung_Lab12.ipynb: 0027 for k in px: Shih_ChienTung_Lab12.ipynb: 0032 cov = [[0 for j in range(4)] for i in range(4)] Shih_ChienTung_Lab12.ipynb: 0033 plot( 142.372,38.297,'b.') Shih_ChienTung_Lab12.ipynb: 0036 i += 1 Shih_ChienTung_Lab12.ipynb: 0038 py = [] Shih_ChienTung_Lab12.ipynb: 0039 elif j == 2: Shih_ChienTung_Lab12.ipynb: 0040 ye = sum(lats)/len(lats) Shih_ChienTung_Lab12.ipynb: 0042 sigma = (17*17/12)**0.5 Shih_ChienTung_Lab12.ipynb: 0045 v = 1.39 Shih_ChienTung_Lab12.ipynb: 0047 pt = [] Shih_ChienTung_Lab12.ipynb: 0054 if j == 0: Shih_ChienTung_Lab12.ipynb: 0057 for j in range(4): Shih_ChienTung_Lab12.ipynb: 0060 px = [] Shih_ChienTung_Lab12.ipynb: 0064 while i < len(lats): Shih_ChienTung_Lab12.ipynb: 0068 J = np.zeros((len(px),4)) Shih_ChienTung_Lab12.ipynb: 0070 else : feddersen_matthew_lab12.ipynb: 0009 ye = 38.297 feddersen_matthew_lab12.ipynb: 0013 xe = 142.372 feddersen_matthew_lab12.ipynb: 0015 print sigma tang_yuhan_lab12.ipynb: 0009 ye = 38.297 tang_yuhan_lab12.ipynb: 0013 xe = 142.372 Ling_Wang-Jung_Lab12.ipynb: 0015 print sigma Ling_Wang-Jung_Lab12.ipynb: 0045 v = 1.39 Ling_Wang-Jung_Lab12.ipynb: 0057 for j in range(4): Wang_Hongru_Lab12.ipynb: 0009 ye = 38.297 Wang_Hongru_Lab12.ipynb: 0013 xe = 142.372 Wang_Xinao_Lab12_SEENPZ.ipynb: 0009 ye = 38.297 Wang_Xinao_Lab12_SEENPZ.ipynb: 0013 xe = 142.372 Gleason_Cole_Lab12_ALSO_HAS_NPZ_IF_NEEDED.ipynb manuel_devasia_Lab12.ipynb: 0019 lats = np.array(npzfile['lats']) manuel_devasia_Lab12.ipynb: 0031 %matplotlib inline manuel_devasia_Lab12.ipynb: 0055 npzfile = np.load('receiver_positions.npz') manuel_devasia_Lab12.ipynb: 0095 lons = np.array(npzfile['lons']) McLaughlin_Sean_Lab12_FA2014.ipynb: 0031 %matplotlib inline Wu_Yi-Luen_Lab12.ipynb: 0064 legend(loc='best') Wang_Hongru_Lab12.ipynb Bai_Hongyang_lab12.ipynb: 0068 print J Shih_ChienTung_Lab12.ipynb: 0068 print J Li_Boyu_Lab12.ipynb: 0068 print J feddersen_matthew_lab12.ipynb: 0023 xe = 142.372 feddersen_matthew_lab12.ipynb: 0069 ye = 38.297 Antonacci_Dominic_genetic_algorithm.ipynb: 0022 jacobian = [] tang_yuhan_lab12.ipynb: 0023 xe = 142.372 tang_yuhan_lab12.ipynb: 0068 print J tang_yuhan_lab12.ipynb: 0069 ye = 38.297 snyder_bradley_Lab12_FA2014.ipynb: 0068 print J Deng_Fei_lab12_SEENPZ.ipynb: 0056 for i in xrange(N): Deng_Fei_lab12_SEENPZ.ipynb: 0068 print J Ling_Wang-Jung_Lab12.ipynb: 0068 print J Antonacci_Dominic_lab12_SEENPZ.ipynb: 0022 jacobian = [] Hsu_HengWei_Lab12_FA2014.ipynb: 0023 xe = 142.372 Hsu_HengWei_Lab12_FA2014.ipynb: 0069 ye = 38.297 Wang_Xinao_Lab12_SEENPZ.ipynb: 0023 xe = 142.372 Wang_Xinao_Lab12_SEENPZ.ipynb: 0069 ye = 38.297 Vankayala_Anuraag_Lab12.ipynb Dandala_Kris_ECE398BD_Lab12_FA2014.ipynb: 0087 print sampling_rate Antonacci_Dominic_genetic_algorithm.ipynb: 0027 for i in xrange(len(lats)): Antonacci_Dominic_lab12_SEENPZ.ipynb: 0027 for i in xrange(len(lats)): Wang_Xinao_Lab12_SEENPZ.ipynb Li_Boyu_Lab12.ipynb: 0030 lats = [] Li_Boyu_Lab12.ipynb: 0034 lons = [] feddersen_matthew_lab12.ipynb: 0009 xe = 142.372 feddersen_matthew_lab12.ipynb: 0024 ye = 38.297 tang_yuhan_lab12.ipynb: 0009 xe = 142.372 tang_yuhan_lab12.ipynb: 0024 ye = 38.297 deVasconcelos_Thiago_Lab12.ipynb: 0030 lats = [] deVasconcelos_Thiago_Lab12.ipynb: 0034 lons = [] Antonacci_Dominic_lab12_SEENPZ.ipynb: 0030 lats = [] Antonacci_Dominic_lab12_SEENPZ.ipynb: 0034 lons = [] Hsu_HengWei_Lab12_FA2014.ipynb: 0009 xe = 142.372 Hsu_HengWei_Lab12_FA2014.ipynb: 0024 ye = 38.297 Wang_Hongru_Lab12.ipynb: 0009 xe = 142.372 Wang_Hongru_Lab12.ipynb: 0024 ye = 38.297 survila2_Lab12_FA2014.ipynb: 0014 y_center = mean(lats) survila2_Lab12_FA2014.ipynb: 0061 x_center = mean(lons) vashistha_siddhant_lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0082 y_e = 38.297 Bai_Hongyang_lab12.ipynb: 0082 y_e = 38.297 Shih_ChienTung_Lab12.ipynb: 0048 figure(figsize=(6,6)) snyder_bradley_Lab12_FA2014.ipynb: 0103 for i in arange(N): Maier_Christopher_Lab12_FA2014.ipynb: 0070 legend() deVasconcelos_Thiago_Lab12.ipynb: 0048 figure(figsize=(6,6)) Lei_TsungHan_Lab12.ipynb: 0048 figure(figsize=(6,6)) Antonacci_Dominic_lab12_SEENPZ.ipynb: 0048 figure(figsize=(6,6)) survila2_Lab12_FA2014.ipynb: 0082 y_e = 38.297 Wu_Yi-Luen_Lab12.ipynb: 0048 figure(figsize=(6,6)) kabaghe_chuma_lab12.ipynb: 0082 y_e = 38.297 survila2_Lab12_FA2014.ipynb Mao_Chengfeng_lab12.ipynb: 0017 x_e = 142.372 Mao_Chengfeng_lab12.ipynb: 0044 y_e = 38.297 Mao_Chengfeng_lab12.ipynb: 0070 return J Bai_Hongyang_lab12.ipynb: 0017 x_e = 142.372 Bai_Hongyang_lab12.ipynb: 0044 y_e = 38.297 Bai_Hongyang_lab12.ipynb: 0045 y_e = 38.322 Bai_Hongyang_lab12.ipynb: 0046 x_e = 142.369 Bai_Hongyang_lab12.ipynb: 0070 return J Maier_Christopher_Lab12_FA2014.ipynb: 0045 y_e = 38.322 Maier_Christopher_Lab12_FA2014.ipynb: 0046 x_e = 142.369 Lei_TsungHan_Lab12.ipynb: 0070 return J Wang_Xinao_Lab12_SEENPZ.ipynb: 0013 y_center = mean(lats) Wang_Xinao_Lab12_SEENPZ.ipynb: 0081 x_center = mean(lons) vashistha_siddhant_lab12.ipynb: 0044 y_e = 38.297 Wu_Yi-Luen_Lab12.ipynb: 0070 return J kabaghe_chuma_lab12.ipynb: 0017 x_e = 142.372 kabaghe_chuma_lab12.ipynb: 0044 y_e = 38.297 Wu_Yi-Luen_Lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0052 return J Bai_Hongyang_lab12.ipynb: 0052 return J crawford_melissa_macrawf2_lab12.ipynb: 0022 print "" Shih_ChienTung_Lab12.ipynb: 0022 print "" Shih_ChienTung_Lab12.ipynb: 0041 figure(figsize=(6,6)) Shih_ChienTung_Lab12.ipynb: 0053 print cov tang_yuhan_lab12.ipynb: 0033 for i in range(N): Deng_Fei_lab12_SEENPZ.ipynb: 0058 N = len(x) Ling_Wang-Jung_Lab12.ipynb: 0033 for i in range(N): Ling_Wang-Jung_Lab12.ipynb: 0053 print cov deVasconcelos_Thiago_Lab12.ipynb: 0041 figure(figsize=(6,6)) Lei_TsungHan_Lab12.ipynb: 0041 figure(figsize=(6,6)) Lei_TsungHan_Lab12.ipynb: 0052 return J Antonacci_Dominic_lab12_SEENPZ.ipynb: 0041 figure(figsize=(6,6)) Gleason_Cole_Lab12_ALSO_HAS_NPZ_IF_NEEDED.ipynb: 0061 legend(loc='best') vashistha_siddhant_lab12.ipynb: 0041 figure(figsize=(6,6)) survila2_Lab12_FA2014.ipynb: 0052 return J kabaghe_chuma_lab12.ipynb: 0033 for i in range(N): kabaghe_chuma_lab12.ipynb Mao_Chengfeng_lab12.ipynb: 0013 x_e = mean(lons) Mao_Chengfeng_lab12.ipynb: 0038 for i in range(10): Mao_Chengfeng_lab12.ipynb: 0063 y_e = mean(lats) Mao_Chengfeng_lab12.ipynb: 0079 y_e = 38.297 Mao_Chengfeng_lab12.ipynb: 0084 x_e = 142.372 Bai_Hongyang_lab12.ipynb: 0079 y_e = 38.297 Bai_Hongyang_lab12.ipynb: 0084 x_e = 142.372 Li_Boyu_Lab12.ipynb: 0085 print N tang_yuhan_lab12.ipynb: 0098 for i in range(N): Ling_Wang-Jung_Lab12.ipynb: 0038 for i in range(10): Ling_Wang-Jung_Lab12.ipynb: 0098 for i in range(N): Lei_TsungHan_Lab12.ipynb: 0038 for i in range(10): vashistha_siddhant_lab12.ipynb: 0079 y_e = 38.297 survila2_Lab12_FA2014.ipynb: 0079 y_e = 38.297 survila2_Lab12_FA2014.ipynb: 0084 x_e = 142.372 Wu_Yi-Luen_Lab12.ipynb: 0098 for i in range(N):