Mercurial > repos > public > sbplib

changeset 89:53fe4b64f65e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Minor cleaning of Wave2dCurve.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 30 Nov 2015 12:20:22 +0100
parents 643bc513b8b8
children 2102af217134
files +scheme/Wave2dCurve.m
diffstat 1 files changed, 8 insertions(+), 8 deletions(-) [+]
line wrap: on
line diff
diff -r 643bc513b8b8 -r 53fe4b64f65e +scheme/Wave2dCurve.m
--- a/+scheme/Wave2dCurve.m	Mon Nov 30 11:41:34 2015 +0100
+++ b/+scheme/Wave2dCurve.m	Mon Nov 30 12:20:22 2015 +0100
@@ -294,8 +294,8 @@
                     s = -1;
 
                     I = gridMatrix(:,1);
-                    coeff_n = obj.a11(:,1);
-                    coeff_t = obj.a12(:,1);
+                    coeff_n = obj.a11(I);
+                    coeff_t = obj.a12(I);
                 case 'e'
                     e = obj.e_e;
                     d_n = obj.du_e;
@@ -303,8 +303,8 @@
                     s = 1;
 
                     I = gridMatrix(:,end);
-                    coeff_n = obj.a11(:,end);
-                    coeff_t = obj.a12(:,end);
+                    coeff_n = obj.a11(I);
+                    coeff_t = obj.a12(I);
                 case 's'
                     e = obj.e_s;
                     d_n = obj.dv_s;
@@ -312,8 +312,8 @@
                     s = -1;
 
                     I = gridMatrix(1,:)';
-                    coeff_n = obj.a22(1,:)';
-                    coeff_t = obj.a12(1,:)';
+                    coeff_n = obj.a22(I);
+                    coeff_t = obj.a12(I);
                 case 'n'
                     e = obj.e_n;
                     d_n = obj.dv_n;
@@ -321,8 +321,8 @@
                     s = 1;
 
                     I = gridMatrix(end,:)';
-                    coeff_n = obj.a22(end,:)';
-                    coeff_t = obj.a12(end,:)';
+                    coeff_n = obj.a22(I);
+                    coeff_t = obj.a12(I);
                 otherwise
                     error('No such boundary: boundary = %s',boundary);
             end