Mercurial > repos > public > sbplib

changeset 64:7067bf8adbfa

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed some typing and added function to plot eigen values of a matrix.
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 18 Nov 2015 17:14:19 +0100
parents 21b18edb9667
children 33f0654a2413
files +noname/Discretization.m +util/get_grid.m plotEig.m
diffstat 3 files changed, 10 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/+noname/Discretization.m	Sat Nov 14 19:48:23 2015 -0800
+++ b/+noname/Discretization.m	Wed Nov 18 17:14:19 2015 +0100
@@ -3,7 +3,7 @@
         name         %Short description
         description  %Longer description
         order        %Order of accuracy
-        h            % scalar desciribing the grid spacing.. (IS THIS THE RIGHT PLACE FOR THIS?)
+        % h            % scalar desciribing the grid spacing.. (IS THIS THE RIGHT PLACE FOR THIS?)
     end
 
     methods (Abstract)
@@ -19,7 +19,7 @@
         %          the appropriate timestepper. It should also provide a default value.
         %     k is a desired timestep
         %     cfl is a choses cfl constant used to set the timestep. ignored if k is set.
-        ts = getTimestepper(obj,method,k,cfl)
+        ts = getTimestepper(obj, method, k, cfl)
 
         % Calculates a timestep for the discretization and a given timestepping method.
         % Can take order, differnt types of scaling in h, or other parameters in Discr into
@@ -31,7 +31,7 @@
 
         % getTimeSnapshot returns a struct which represents the solution in ts at current time.
         % if ts is empty or 0 a representation of the initial conditions be returned.
-        repr = getTimeSnapshot(obj,ts)
+        repr = getTimeSnapshot(obj, ts)
 
 
         % Sets up movie recording to a given file.
--- a/+util/get_grid.m	Sat Nov 14 19:48:23 2015 -0800
+++ b/+util/get_grid.m	Wed Nov 18 17:14:19 2015 +0100
@@ -1,3 +1,4 @@
+% [x,h] = get_grid(a,b,m)
 function [x,h] = get_grid(a,b,m)
     % TODO: allow the interval to be a vector
     x = linspace(a,b,m)';
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/plotEig.m	Wed Nov 18 17:14:19 2015 +0100
@@ -0,0 +1,6 @@
+% Plots all the eigen values for a sparse matrix.
+function plotEig(A)
+    e = eig(full(A));
+    e = sort(e);
+    plot(e,'.')
+end
\ No newline at end of file