log-log plots

fixed_domain/diffusion_fixed_domain_solution.py

@@ -7,7 +7,7 @@ import progressbar
 df.set_log_level(df.LogLevel.ERROR)
 df.set_log_level(df.LogLevel.ERROR)
 df.parameters['form_compiler']['optimize'] = True
 df.parameters['form_compiler']['optimize'] = True
 
 
-Nx, L, D, tmax, dt = 32, 1, 0.1, 0.25, 0.25/50
+Nx, L, D, tmax, dt = 32, 1, 0.1, 0.25, 0.25/400
 Nt = int(tmax/dt)
 Nt = int(tmax/dt)
 
 
 mesh = df.IntervalMesh(Nx, 0, L)
 mesh = df.IntervalMesh(Nx, 0, L)
@@ -45,7 +45,7 @@ ax.set_ylabel(r'$c(x,t)$')
 ax.legend(loc=1)
 ax.legend(loc=1)
 
 
 error = np.abs(c_array - c_exact)
 error = np.abs(c_array - c_exact)
-print('Error at t=0.25 is',np.max(error[-1]))
+print('Max error at t=0.25 is',np.max(error[-1]))
 
 
 err_plot, = ax1.plot(x, error[1], 'bo', mfc='none', ms=6)
 err_plot, = ax1.plot(x, error[1], 'bo', mfc='none', ms=6)
 ax1.set_ylim(np.min(error), np.max(error))
 ax1.set_ylim(np.min(error), np.max(error))