Mercurial > repos > public > sbplib

comparison +time/SBPInTime.m @ 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 38173ea263ed b5e5b195da1e
comparison
equal deleted inserted replaced
406:9fd9b1bea3d2 407:e1db62d14835
83 % LU factorization 83 % LU factorization
84 [obj.L,obj.U,obj.P,obj.Q] = lu(obj.M); 84 [obj.L,obj.U,obj.P,obj.Q] = lu(obj.M);
85 85
86 % Pretend that the initial condition is the last level 86 % Pretend that the initial condition is the last level
87 % of a previous step. 87 % of a previous step.
88 obj.v = obj.Et_r * v0; 88 obj.v = 1/(e_r'*e_r) * obj.Et_r * v0;
89 89
90 end 90 end
91 91
92 function [v,t] = getV(obj) 92 function [v,t] = getV(obj)
93 v = obj.Et_r' * obj.v; 93 v = obj.Et_r' * obj.v;
104 end 104 end
105 end 105 end
106 106
107 methods(Static) 107 methods(Static)
108 function N = smallestBlockSize(order,TYPE) 108 function N = smallestBlockSize(order,TYPE)
109 default_arg('TYPE','equidistant') 109 default_arg('TYPE','gauss')
110 110
111 switch TYPE 111 switch TYPE
112 112
113 case 'equidistant' 113 case 'equidistant'
114 switch order 114 switch order