From 5bf27846558d3d760c186b716e47980dac605dad Mon Sep 17 00:00:00 2001 From: Lalit Rajendran Date: Fri, 24 Jan 2020 11:33:05 -0500 Subject: [PATCH] added lines in the sample script to save results to a .mat file --- sample-result.mat | Bin 0 -> 8560 bytes sample_script.py | 7 ++++++- 2 files changed, 6 insertions(+), 1 deletion(-) create mode 100644 sample-result.mat diff --git a/sample-result.mat b/sample-result.mat new file mode 100644 index 0000000000000000000000000000000000000000..0f1f40ba72a7ae2c387aae318c7198d77a42c28e GIT binary patch literal 8560 zcmeHLdss~C9xf>rl}mC9ks`WGA<5;X8>P~ANix%1(o#ugMg}Q_q8&*ag-yz(+)HjH z6+$jWN<|@;okE1QW-8)*&Th^hI{Rs#^E^lA4>Qlq`>ge>?|Z-B?Ong`bD2ELd9sb( zXj7J+%j8)mlwW|4o?8Gvh*~UOWUVJ$yu@$0v7W8ihacoKUvF`swVs36PtS=TsAn-s z&)nR4q=mJGxt;~ff~D|Zd=wpBG|4cQs)B+FK?488Z!;;kTz4h-?Wah-a|x=E-?s<* za0YtBhKh3>UzJ{2q$13D|5le9J78hXFKU}C9a&UPn*JvZ1si_OppcTAi;_KZF7bUn zp55+}6Sl)PX-LLMIilr0TAoA8IJC@zmU+>_gBD&n{^bAZ|DXByv;KbS&rkjR`TqHV z@8AEVZ}9jxeS^Gz_|YiuSK15ZeNua&ygzF%l=p@0h4Oy1y-?la#!l#G)U6pFh)wy7KqgoHz6UCIY}#= z@F`GaOyVV$v9`gw@}@AYSRJxZjIzWwCc(t)Bj?Wg8j?6EKsfvx&%ml_#& zGbB%3qa#4ls(Bbw#>K@8%^BcjQ3tSO4r!6c}eX?HvdM?rtw=yUa#n>9(teT(+s<*FV#Nt(Wr1wIdlOR<(D*T#hM&s?NLw9 z(_rImXs%M66B`x^hvsc{=0GhfIeuIm7nhISzNJC(!ZIt#JnRw$`q%+HC-chd?62}a zqf|41Z6Lzd*c+?~$0?k8X4^@BfB<{8UFaRwiHGS=CBMyT;=ub|;>IZq8?vj_#s2*{ z2&owUSX#)zfWnQ#GE{gNGB)|i>`}!2OUA#QKTYyDJTa_lMBDt``}=q}=#cm?{bNn+ zZxl9F_H;ewEI=r_2CYlufywo)aJbGz;B?XO=s+&gPS4r4Duj!io;u4O_2ps1Zrz=| zuJSO$YoSKeDgnwI7r$&uq2PXb$gHD8Kiaat?|f~yBT}L_Qn!M{0t%Jy+(Szic%zob z)_v+JfX+>u!5>LJ{t_~AZ(@8MV55YpCE$G z;fY!X<{~IHtr;6Z>Rfxn=*LSKA$D6{%JHuf;-#us;rOk#^;hnn|22M}|Lt}C_s<`c z6uxy2|8IRm&xKY4avpxf#t-4M=W5arE87a6gGPN=OkQpTjl4AdT*)aK`jNGhhTWkd z6KLN|tfcXvpWEv(Pie$PA1qHSrC~Zz%c18v8Xp3UZup33C{Eekq1lSY&PVI7Ol?ZS_i?&&Vkgy0n?~h=&OwILX@r@+O25B@hT5%;ai3oh-W#(ePlht+`Yd*} zt0#k*yLGF^NEmdn4l}7;#^AuPVb`~MF?g0KG1MK-K;0oNhIvfmd69~4%vKtHT+8QW6-_3QdJK!F53R$&L%Ps$E7Nx zb|;PNn}bTW%xulysQDxx5*iM?z3al^Qe2?sA6-b+@jT>Uv+H>V^*K&P z4I>z&zSSL8MDpYN_$y1|wH{ofp{^b`u{*JQucKkAZUcki;S~nviBfomo4bi@WO$M) zDJsd8VYE&yudHJO5~k-*YwXZ~erXC-vTPah#KE6WkaZmE)alxG;*W$mJ{N=w86ZWt?=)lQVkxSQu3Wb{R0h9a zCD9Iw4R~EY(9_nb0Y}#+9C4b~fWE5B@B5N<1fOo)*|bK6&nJCfE+_FxI(>TZAd(-v zOm0lXat1M%5|o-re80E<)74`3hX*v?3`!_&2qboQyWyeQUkdMx^$&DPq;Qi_aa&f& zFl@g;Sdl^l%x0#SluT`aqD5Snp=2HJS89AMB>Z|GH=47DgSv@iqQ~+f5Z3lF5Rscsm+KlW(o?kEStuw$|vgCu!&mI@$Z? z6&m9<(D&(jqNi~YYeVXY&cF9}Tut(-gQGPyGl|BhBPMYZgUNfVEWc)`CyfKKgZIe# V5&x$Qn%eNW_50`hdi+t>{{Zn2NVWg~ literal 0 HcmV?d00001 diff --git a/sample_script.py b/sample_script.py index fbc9b85..6526ab9 100755 --- a/sample_script.py +++ b/sample_script.py @@ -17,6 +17,7 @@ # Load the data: data = loadmat('sample-data.mat', squeeze_me=True) +# set the density uncertainty at the boundary points sigma_rho_dirichlet = 0.01 # calculate density and uncertainty @@ -27,6 +28,9 @@ sigma_grad_x=data['sigma_rho_x'], sigma_grad_y=data['sigma_rho_y'], sigma_dirichlet=sigma_rho_dirichlet) +# save the results to file +savemat(file_name='sample-result.mat', mdict={'X': data['X'], 'Y': data['Y'], 'rho': rho, 'sigma_rho': sigma_rho}, long_field_names=True) + # Plot the results fig1 = plt.figure(1, figsize=(12,8)) plt.figure(1) @@ -60,7 +64,8 @@ plt.title('sigma rho') plt.tight_layout() -# save results to file + +# save plot to file plt.savefig('sample-result.png') plt.close()