Mercurial > repos > public > sbplib

changeset 407:e1db62d14835 feature/SBPInTimeGauss

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed bug with initial data and made gauss the default type
author Martin Almquist <martin.almquist@it.uu.se>
date Thu, 02 Feb 2017 19:30:16 +0100
parents 9fd9b1bea3d2
children ba73c9c8d1a6
files +time/SBPInTime.m
diffstat 1 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/+time/SBPInTime.m	Thu Feb 02 17:18:51 2017 +0100
+++ b/+time/SBPInTime.m	Thu Feb 02 19:30:16 2017 +0100
@@ -85,7 +85,7 @@
 
             % Pretend that the initial condition is the last level
             % of a previous step.
-            obj.v = obj.Et_r * v0;
+            obj.v =  1/(e_r'*e_r) * obj.Et_r * v0;
 
         end
 
@@ -106,7 +106,7 @@
 
     methods(Static)
         function N = smallestBlockSize(order,TYPE)
-            default_arg('TYPE','equidistant')
+            default_arg('TYPE','gauss')
 
             switch TYPE