Mercurial > repos > public > sbplib

changeset 264:8a625c5a3633 operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Changed input parameter L (domain length) to lim (cell with domain boundaries)
author Martin Almquist <martin.almquist@it.uu.se>
date Fri, 09 Sep 2016 10:39:12 +0200
parents 21a180acbd49
children 3e0432b3cb68
files +sbp/D1Nonequidistant.m +sbp/D1Upwind.m +sbp/D2BlockNorm.m +sbp/D2Standard.m +sbp/D2Variable.m +sbp/D4Compatible.m +sbp/D4CompatibleVariable.m +sbp/D4Standard.m
diffstat 8 files changed, 43 insertions(+), 15 deletions(-) [+]
line wrap: on
line diff
--- a/+sbp/D1Nonequidistant.m	Thu Sep 08 17:50:30 2016 +0200
+++ b/+sbp/D1Nonequidistant.m	Fri Sep 09 10:39:12 2016 +0200
@@ -13,13 +13,17 @@
     end
     
     methods
-        function obj = D1Nonequidistant(m,L,order,option)
+        function obj = D1Nonequidistant(m,lim,order,option)
             
             default_arg('option','Accurate');
             % 'Accurate' operators are optimized for accuracy
             % 'Minimal' operators have the smallest possible boundary
             %  closure
             
+            x_l = lim{1};
+            x_r = lim{2};
+            L = x_r-x_l;
+            
             switch option
                 
                 case {'Accurate','accurate','A'}
@@ -66,6 +70,8 @@
                     
             end
             
+            obj.x = obj.x + x_l;
+            
             obj.e_1 = sparse(m,1);
             obj.e_m = sparse(m,1);
             obj.e_1(1) = 1;
--- a/+sbp/D1Upwind.m	Thu Sep 08 17:50:30 2016 +0200
+++ b/+sbp/D1Upwind.m	Fri Sep 09 10:39:12 2016 +0200
@@ -13,10 +13,13 @@
     end
 
     methods
-        function obj = D1Upwind(m,L,order)
+        function obj = D1Upwind(m,lim,order)
             
+            x_l = lim{1};
+            x_r = lim{2};
+            L = x_r-x_l;
             obj.h = L/(m-1);
-            obj.x = linspace(0,L,m)';
+            obj.x = linspace(x_l,x_r,m)';
 
             switch order
                 case 2
--- a/+sbp/D2BlockNorm.m	Thu Sep 08 17:50:30 2016 +0200
+++ b/+sbp/D2BlockNorm.m	Fri Sep 09 10:39:12 2016 +0200
@@ -19,10 +19,13 @@
 
 
     methods
-        function obj = D2BlockNorm(m,L,order)
+        function obj = D2BlockNorm(m,lim,order)
 
+            x_l = lim{1};
+            x_r = lim{2};
+            L = x_r-x_l;
             obj.h = L/(m-1);
-            obj.x = linspace(0,L,m)';
+            obj.x = linspace(x_l,x_r,m)';
 
             if order == 4
                 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,...
--- a/+sbp/D2Standard.m	Thu Sep 08 17:50:30 2016 +0200
+++ b/+sbp/D2Standard.m	Fri Sep 09 10:39:12 2016 +0200
@@ -18,10 +18,13 @@
     end
 
     methods
-        function obj = D2Standard(m,L,order)
+        function obj = D2Standard(m,lim,order)
 
+            x_l = lim{1};
+            x_r = lim{2};
+            L = x_r-x_l;
             obj.h = L/(m-1);
-            obj.x = linspace(0,L,m)';
+            obj.x = linspace(x_l,x_r,m)';
 
             if order == 2
                 [obj.H, obj.HI, obj.D1, obj.D2, obj.e_1,...
--- a/+sbp/D2Variable.m	Thu Sep 08 17:50:30 2016 +0200
+++ b/+sbp/D2Variable.m	Fri Sep 09 10:39:12 2016 +0200
@@ -17,10 +17,13 @@
     end
 
     methods
-        function obj = D2Variable(m,L,order)
+        function obj = D2Variable(m,lim,order)
             
+            x_l = lim{1};
+            x_r = lim{2};
+            L = x_r-x_l;
             obj.h = L/(m-1);
-            obj.x = linspace(0,L,m)';
+            obj.x = linspace(x_l,x_r,m)';
 
             switch order
                 case 4
--- a/+sbp/D4Compatible.m	Thu Sep 08 17:50:30 2016 +0200
+++ b/+sbp/D4Compatible.m	Fri Sep 09 10:39:12 2016 +0200
@@ -27,10 +27,14 @@
 
 
     methods
-        function obj = D4Compatible(m,L,order)
+        function obj = D4Compatible(m,lim,order)
+            
             
+            x_l = lim{1};
+            x_r = lim{2};
+            L = x_r-x_l;
             obj.h = L/(m-1);
-            obj.x = linspace(0,L,m)';
+            obj.x = linspace(x_l,x_r,m)';
 
             if order == 2
                 [obj.H, obj.HI, obj.D1, obj.D4, obj.e_1, obj.e_m, obj.M4,...
--- a/+sbp/D4CompatibleVariable.m	Thu Sep 08 17:50:30 2016 +0200
+++ b/+sbp/D4CompatibleVariable.m	Fri Sep 09 10:39:12 2016 +0200
@@ -27,10 +27,13 @@
 
 
     methods
-        function obj = D4CompatibleVariable(m,L,order)
+        function obj = D4CompatibleVariable(m,lim,order)
             
+            x_l = lim{1};
+            x_r = lim{2};
+            L = x_r-x_l;
             obj.h = L/(m-1);
-            obj.x = linspace(0,L,m)';
+            obj.x = linspace(x_l,x_r,m)';
 
             if order == 2
                 [obj.H, obj.HI, ~, obj.D2, ~, obj.D4, obj.e_1, obj.e_m,...
--- a/+sbp/D4Standard.m	Thu Sep 08 17:50:30 2016 +0200
+++ b/+sbp/D4Standard.m	Fri Sep 09 10:39:12 2016 +0200
@@ -27,10 +27,13 @@
 
 
     methods
-        function obj = D4Standard(m,L,order)
+        function obj = D4Standard(m,lim,order)
             
+            x_l = lim{1};
+            x_r = lim{2};
+            L = x_r-x_l;
             obj.h = L/(m-1);
-            obj.x = linspace(0,L,m)';
+            obj.x = linspace(x_l,x_r,m)';
 
             if order == 4
                 [obj.H, obj.HI, obj.D1, obj.D2, obj.D3, obj.D4,...