Updating post-processing script to compute dT/dkappa and dg/dkappa.

ChristianZhao · May 27, 2020 · 1bcde73 · 1bcde73
1 parent b2f6eec
commit 1bcde73
Showing 1 changed file with 27 additions and 19 deletions.
diff --git a/tests/small/Thermal1D/compute_errs_source_no_movie2avi.m b/tests/small/Thermal1D/compute_errs_source_no_movie2avi.m
@@ -13,34 +13,40 @@
 x = [0:dx:1]; 
 T = ncread('thermal1D_with_source_out.exo', 'vals_nod_var1'); 
 %size(T)
-Texact = 0*T; 
+Texact = 0*T;
+dTdkappa = 0*T;  
 for i=1:length(x)
   Texact(i,:) = a*x(i)*(1.0-x(i))*cos(2.0*pi*kappa*t/rho/c); 
+  dTdkappaExact(i, :) = -2.0*pi/rho/c*a*x(i)*(1.0-x(i))*t.*sin(2.0*pi*kappa*t/rho/c); 
 end  
 %size(Texact) 
 
 K = length(t); 
-fig1=figure(1);
-winsize = get(fig1,'Position');
-winsize(1:2) = [0 0];
-Movie=moviein(K,fig1,winsize);
-set(fig1,'NextPlot','replacechildren')
+%fig1=figure(1);
+%winsize = get(fig1,'Position');
+%winsize(1:2) = [0 0];
+%Movie=moviein(K,fig1,winsize);
+%set(fig1,'NextPlot','replacechildren')
 
 rel_err = []; 
-
+soln_avg_rel_err = [];  
 for i=1:K
-  set(gca,'NextPlot','replacechildren') ; 
-  plot(x, T(:,i), 'b'); 
-  hold on; 
-  plot(x, Texact(:,i), '--r'); 
+  %set(gca,'NextPlot','replacechildren') ; 
+  %plot(x, T(:,i), 'b');
+  soln_avg_computed = mean(T(:,i));  
+  %hold on; 
+  %plot(x, Texact(:,i), '--r'); 
+  soln_avg_exact = mean(Texact(:,i));
+  dgdp_exact = mean(dTdkappaExact(:,i));  
   err = norm(T(:,i)-Texact(:,i))/norm(Texact(:,i)); 
   rel_err = [rel_err; err]; 
-  xlabel('x'); 
-  ylabel('Temp'); 
-  legend('Computed','Exact', 'Location','NorthWest'); 
-  title(['Temperature solution at t = ', num2str(t(i))]); 
-  axis([0 1 0 4])  
-  pause(0.1)
+  soln_avg_rel_err = [soln_avg_rel_err; norm(soln_avg_computed-soln_avg_exact)/norm(soln_avg_exact)]; 
+  %xlabel('x'); 
+  %ylabel('Temp'); 
+  %legend('Computed','Exact', 'Location','NorthWest'); 
+  %title(['Temperature solution at t = ', num2str(t(i))]); 
+  %axis([0 1 0 4])  
+  %pause(0.1)
   %Movie(:,i)=getframe(fig1,winsize);
   %mov(i) = getframe(gcf);
   %n = strcat("thermal_", num2str(i)); 
@@ -51,7 +57,10 @@
 %To convert the pngs to mp4, run:
 %  ffmpeg -framerate 10 -start_number 1 -i thermal_%1d.png -r 5 -vf scale=-2:1080,setsar=1 -pix_fmt yuv420p thermal_MMS.mp4
 
-X = ['Average relative error over time = ', num2str(mean(rel_err))]; 
+X = ['Average solution relative error over time = ', num2str(mean(rel_err))]; 
+disp(X) 
+disp([' ']); 
+X = ['Average response relative error over time = ', num2str(mean(soln_avg_rel_err))]; 
 disp(X) 
 disp([' ']); 
 figure(); 
@@ -60,4 +69,3 @@
 ylabel('Relative Error'); 
 title('Relative Errors w.r.t. Exact Solution');  
 
-