From 19dec53f0319b92f4aff0ec2f1be3dfd67d12824 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Fri, 13 Jul 2012 15:29:05 +0200
Subject: [PATCH] Add a helper for velocity stepping plots
---
src/commands.m | 22 ++++++++++++++++++++++
1 file changed, 22 insertions(+)
create mode 100644 src/commands.m
diff --git a/src/commands.m b/src/commands.m
new file mode 100644
index 00000000..da7457d9
--- /dev/null
+++ b/src/commands.m
@@ -0,0 +1,22 @@
+## State evolution
+load ../data;
+graphics_toolkit fltk;
+path(path, '../../../');
+plot_state
+
+## Comparison of steady-state convergence rates, to and fro
+# 600 time steps
+load ../data;
+graphics_toolkit fltk;
+path(path, '../../../');
+plot(1:144,A(121:(120+144),1)-A(121,1),1:144,-interp1(1:240, A(361:600,1), linspace(1,240,144))+A(361,1));
+axis tight;
+legend('decline (first)','incline (second)');
+
+# 6000 time steps
+load ../data;
+graphics_toolkit fltk;
+path(path, '../../../');
+plot(1:1440,A(1201:(1200+1440),1)-A(1201,1),1:1440,-interp1(1:2400, A(3601:6000,1), linspace(1,2400,1440))+A(3601,1));
+axis tight;
+legend('decline (first)','incline (second)');
--
GitLab