Mercurial > repos > public > sbplib

changeset 55:a8ed986fcf57

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Minor renaming and clean up in 2d wave schemes.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 12 Nov 2015 09:47:04 -0800
parents 2194cd385419
children ce90abc350c5
files +scheme/Wave2d.m +scheme/Wave2dCurve.m
diffstat 2 files changed, 8 insertions(+), 9 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Wave2d.m	Wed Nov 11 19:21:38 2015 -0800
+++ b/+scheme/Wave2d.m	Thu Nov 12 09:47:04 2015 -0800
@@ -1,4 +1,4 @@
-classdef SchmWave2d < noname.Scheme
+classdef Wave2d < scheme.Scheme
     properties
         m % Number of points in each direction, possibly a vector
         h % Grid spacing
@@ -22,8 +22,11 @@
     end
 
     methods
-        function obj = SchmWave2d(m,xlim,ylim,order,alpha)
-            default_arg('a',1);
+        function obj = Wave2d(m,lim,order,alpha)
+            default_arg('alpha',1);
+
+            xlim = lim{1};
+            ylim = lim{2};
 
             if length(m) == 1
                 m = [m m];
--- a/+scheme/Wave2dCurve.m	Wed Nov 11 19:21:38 2015 -0800
+++ b/+scheme/Wave2dCurve.m	Thu Nov 12 09:47:04 2015 -0800
@@ -234,8 +234,8 @@
             % v denotes the solution in the neighbour domain
             tuning = 1.2;
             % tuning = 20.2;
-            [e_u, d_n_u, d_t_u, coeff_n_u, coeff_t_u, s_u, gamm_u, halfnorm_inv_u_n, halfnorm_inv_u_t, halfnorm_u_t, a_n_u, a_t_u] = obj.get_boundary_ops(boundary);
-            [e_v, d_n_v, d_t_v, coeff_n_v, coeff_t_v, s_v, gamm_v, halfnorm_inv_v_n, halfnorm_inv_v_t, halfnorm_v_t, a_n_v, a_t_v] = neighbour_scheme.get_boundary_ops(boundary);
+            [e_u, d_n_u, d_t_u, coeff_n_u, coeff_t_u, s_u, gamm_u, halfnorm_inv_u_n, halfnorm_inv_u_t, halfnorm_u_t] = obj.get_boundary_ops(boundary);
+            [e_v, d_n_v, d_t_v, coeff_n_v, coeff_t_v, s_v, gamm_v, halfnorm_inv_v_n, halfnorm_inv_v_t, halfnorm_v_t] = neighbour_scheme.get_boundary_ops(boundary);
 
             a_n_u = spdiag(coeff_n_u);
             a_t_u = spdiag(coeff_t_u);
@@ -314,15 +314,11 @@
                     halfnorm_inv_t = obj.Hiv;
                     halfnorm_t = obj.Hv;
                     gamm = obj.gamm_u;
-                    a_n = obj.a11;
-                    a_t = obj.a12;
                 case {'s','n'}
                     halfnorm_inv_n = obj.Hiv;
                     halfnorm_inv_t = obj.Hiu;
                     halfnorm_t = obj.Hu;
                     gamm = obj.gamm_v;
-                    a_n = obj.a22;
-                    a_t = obj.a12;
             end
         end