Mercurial > repos > public > sbplib

changeset 326:b19e142fcae1 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed bug in setting of boundary derivative.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 27 Sep 2016 08:44:17 +0200
parents 72468bc9b63f
children d24869abc7cd
files +sbp/+implementations/d4_compatible_6.m +sbp/+implementations/d4_lonely_4_min_boundary_points.m +sbp/+implementations/d4_lonely_6_2.m +sbp/+implementations/d4_lonely_6_3.m +sbp/+implementations/d4_lonely_6_min_boundary_points.m +sbp/+implementations/d4_lonely_8_higher_boundary_order.m +sbp/+implementations/d4_lonely_8_min_boundary_points.m +sbp/+implementations/d4_variable_2.m +sbp/+implementations/d4_variable_4.m +sbp/+implementations/d4_variable_6.m
diffstat 10 files changed, 10 insertions(+), 10 deletions(-) [+]
line wrap: on
line diff
--- a/+sbp/+implementations/d4_compatible_6.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_compatible_6.m	Tue Sep 27 08:44:17 2016 +0200
@@ -25,7 +25,7 @@
 
     % Vi b?rjar med normen. Notera att alla SBP operatorer delar samma norm,
     % vilket ?r n?dv?ndigt f?r stabilitet
-    
+
     BP = 8;
     if(m<2*BP)
         error(['Operator requires at least ' num2str(2*BP) ' grid points']);
--- a/+sbp/+implementations/d4_lonely_4_min_boundary_points.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_lonely_4_min_boundary_points.m	Tue Sep 27 08:44:17 2016 +0200
@@ -32,7 +32,7 @@
 
     d1_l = sparse(m,1);
     d1_l(1:4) = 1/h*[-11/6 3 -3/2 1/3];
-    d1_r = -rot90(d1_l);
+    d1_r = -rot90(d1_l, 2);
 
     d2_l = sparse(m,1);
     d2_l(1:4) = 1/h^2*[2 -5 4 -1];
--- a/+sbp/+implementations/d4_lonely_6_2.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_lonely_6_2.m	Tue Sep 27 08:44:17 2016 +0200
@@ -41,7 +41,7 @@
 
     d1_l = sparse(m,1);
     d1_l(1:6) = [-0.137e3/0.60e2 5 -5 0.10e2/0.3e1 -0.5e1/0.4e1 0.1e1/0.5e1;]/h;
-    d1_r = -rot90(d1_l);
+    d1_r = -rot90(d1_l, 2);
 
     d2_l = sparse(m,1);
     d2_l(1:6) = [0.15e2/0.4e1 -0.77e2/0.6e1 0.107e3/0.6e1 -13 0.61e2/0.12e2 -0.5e1/0.6e1;]/h^2;
--- a/+sbp/+implementations/d4_lonely_6_3.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_lonely_6_3.m	Tue Sep 27 08:44:17 2016 +0200
@@ -36,7 +36,7 @@
 
     d1_l = sparse(m,1);
     d1_l(1:6) = [-0.137e3/0.60e2 5 -5 0.10e2/0.3e1 -0.5e1/0.4e1 0.1e1/0.5e1;]/h;
-    d1_r = -rot90(d1_l);
+    d1_r = -rot90(d1_l, 2);
 
     d2_l = sparse(m,1);
     d2_l(1:6) = [0.15e2/0.4e1 -0.77e2/0.6e1 0.107e3/0.6e1 -13 0.61e2/0.12e2 -0.5e1/0.6e1;]/h^2;
--- a/+sbp/+implementations/d4_lonely_6_min_boundary_points.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_lonely_6_min_boundary_points.m	Tue Sep 27 08:44:17 2016 +0200
@@ -39,7 +39,7 @@
 
     d1_l = sparse(m,1);
     d1_l(1:5) = [-25/12, 4, -3, 4/3, -1/4]/h;
-    d1_r = -rot90(d1_l);
+    d1_r = -rot90(d1_l, 2);
 
     d2_l = sparse(m,1);
     d2_l(1:5) = [0.35e2/0.12e2 -0.26e2/0.3e1 0.19e2/0.2e1 -0.14e2/0.3e1 0.11e2/0.12e2;]/h^2;
--- a/+sbp/+implementations/d4_lonely_8_higher_boundary_order.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_lonely_8_higher_boundary_order.m	Tue Sep 27 08:44:17 2016 +0200
@@ -38,7 +38,7 @@
 
     d1_l = sparse(m,1);
     d1_l(1:7) = [-0.49e2/0.20e2 6 -0.15e2/0.2e1 0.20e2/0.3e1 -0.15e2/0.4e1 0.6e1/0.5e1 -0.1e1/0.6e1]/h;
-    d1_r = -rot90(d1_l);
+    d1_r = -rot90(d1_l, 2);
 
     d2_l = sparse(m,1);
     d2_l(1:7) = [0.203e3/0.45e2 -0.87e2/0.5e1 0.117e3/0.4e1 -0.254e3/0.9e1 0.33e2/0.2e1 -0.27e2/0.5e1 0.137e3/0.180e3]/h^2;
--- a/+sbp/+implementations/d4_lonely_8_min_boundary_points.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_lonely_8_min_boundary_points.m	Tue Sep 27 08:44:17 2016 +0200
@@ -37,7 +37,7 @@
 
     d1_l = sparse(m,1);
     d1_l(1:6) = [-0.137e3/0.60e2 5 -5 0.10e2/0.3e1 -0.5e1/0.4e1 0.1e1/0.5e1;]/h;
-    d1_r = -rot90(d1_l);
+    d1_r = -rot90(d1_l, 2);
 
     d2_l = sparse(m,1);
     d2_l(1:6) = [0.15e2/0.4e1 -0.77e2/0.6e1 0.107e3/0.6e1 -13 0.61e2/0.12e2 -0.5e1/0.6e1;]/h^2;
--- a/+sbp/+implementations/d4_variable_2.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_variable_2.m	Tue Sep 27 08:44:17 2016 +0200
@@ -29,7 +29,7 @@
 
     d1_l = sparse(m,1);
     d1_l(1:3) = 1/h*[-3/2 2 -1/2];
-    d1_r = -rot90(d1_l);
+    d1_r = -rot90(d1_l, 2);
 
     d2_l = sparse(m,1);
     d2_l(1:3) = 1/h^2*[1 -2 1];
--- a/+sbp/+implementations/d4_variable_4.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_variable_4.m	Tue Sep 27 08:44:17 2016 +0200
@@ -29,7 +29,7 @@
 
     d1_l = sparse(m,1);
     d1_l(1:4) = 1/h*[-11/6 3 -3/2 1/3];
-    d1_r = -rot90(d1_l);
+    d1_r = -rot90(d1_l, 2);
 
     d2_l = sparse(m,1);
     d2_l(1:4) = 1/h^2*[2 -5 4 -1];
--- a/+sbp/+implementations/d4_variable_6.m	Mon Sep 26 09:55:16 2016 +0200
+++ b/+sbp/+implementations/d4_variable_6.m	Tue Sep 27 08:44:17 2016 +0200
@@ -37,7 +37,7 @@
 
     d1_l = sparse(m,1);
     d1_l(1:5) = [-25/12, 4, -3, 4/3, -1/4]/h;
-    d1_r = -rot90(d1_l);
+    d1_r = -rot90(d1_l, 2);
 
     d2_l = sparse(m,1);
     d2_l(1:5) = [0.35e2/0.12e2 -0.26e2/0.3e1 0.19e2/0.2e1 -0.14e2/0.3e1 0.11e2/0.12e2;]/h^2;