Mercurial > repos > public > sbplib

annotate +sbp/+implementations/d4_lonely_6_3.m @ 1198:2924b3a9b921 feature/d2_compatible

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add OpSet for fully compatible D2Variable, created from regular D2Variable by replacing d1 by first row of D1. Formal reduction by one order of accuracy at the boundary point.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 16 Aug 2019 14:30:28 -0700
parents b19e142fcae1
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
324
c0cbffcf6513 Created new operator class for lonely D4 operators. Removed some output parameters of implementations.
Jonatan Werpers <jonatan@werpers.com>
parents: 322
diff changeset
1 function [H, HI, D4, e_l, e_r, M4, d2_l, d2_r, d3_l, d3_r, d1_l, d1_r] = d4_variable_6_3(m,h)
312
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
3 %%% 6:te ordn. SBP Finita differens %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
4 %%% operatorer med diagonal norm %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
5 %%% Extension to variable koeff %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
6 %%% %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
7 %%% H (Normen) %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
8 %%% D1=H^(-1)Q (approx f?rsta derivatan) %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
9 %%% D2 (approx andra derivatan) %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
10 %%% D2=HI*(R+C*D*S %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
11 %%% %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
12 %%% R=-D1'*H*C*D1-RR %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
13 %%% %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
14 %%% RR ?r dissipation) %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
15 %%% Dissipationen uppbyggd av D4: %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
16 %%% DI=D4*B*H*D4 %%%
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
17 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
246
fe26791489e0 Added a bunch of new operators. Still non-functional.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
18
312
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
19 % H?r med 7 RP ist?llet f?r 8 f?r D4 operatorn, dock samma randderivator
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
20 % Denna ?r noggrannare, och har 2a ordningens randdslutning och b?r ge 6te
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
21 % ordningens konvergens. Hade 2 fria parametrar att optimera
246
fe26791489e0 Added a bunch of new operators. Still non-functional.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
22
322
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
23 % Norm
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
24 Hv = ones(m,1);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
25 Hv(1:7) = [0.414837907e9/0.1191965760e10, 0.475278367e9/0.397321920e9, 0.13872751e8/0.12416310e8, 0.346739027e9/0.595982880e9, 0.560227469e9/0.397321920e9, 0.322971631e9/0.397321920e9, 0.616122491e9/0.595982880e9];
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
26 Hv(m-6:m) = rot90(Hv(1:7),2);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
27 Hv = h*Hv;
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
28 H = spdiag(Hv, 0);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
29 HI = spdiag(1./Hv, 0);
246
fe26791489e0 Added a bunch of new operators. Still non-functional.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
30
fe26791489e0 Added a bunch of new operators. Still non-functional.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
31
322
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
32 % Boundary operators
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
33 e_l = sparse(m,1);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
34 e_l(1) = 1;
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
35 e_r = rot90(e_l, 2);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
36
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
37 d1_l = sparse(m,1);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
38 d1_l(1:6) = [-0.137e3/0.60e2 5 -5 0.10e2/0.3e1 -0.5e1/0.4e1 0.1e1/0.5e1;]/h;
326
b19e142fcae1 Fixed bug in setting of boundary derivative.
Jonatan Werpers <jonatan@werpers.com>
parents: 325
diff changeset
39 d1_r = -rot90(d1_l, 2);
322
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
40
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
41 d2_l = sparse(m,1);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
42 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;
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
43 d2_r = rot90(d2_l, 2);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
44
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
45 d3_l = sparse(m,1);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
46 d3_l(1:6) = [-0.17e2/0.4e1 0.71e2/0.4e1 -0.59e2/0.2e1 0.49e2/0.2e1 -0.41e2/0.4e1 0.7e1/0.4e1;]/h^3;
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
47 d3_r = -rot90(d3_l, 2);
316
203afa156f59 Collected boundary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
48
203afa156f59 Collected boundary operators.
Jonatan Werpers <jonatan@werpers.com>
parents: 312
diff changeset
49
322
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
50 % Fourth derivative, 1th order accurate at first 8 boundary points
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
51 stencil = [7/240, -2/5, 169/60, -122/15, 91/8, -122/15, 169/60, -2/5, 7/240];
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
52 diags = -4:4;
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
53 M4 = stripeMatrix(stencil, diags, m);
246
fe26791489e0 Added a bunch of new operators. Still non-functional.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
54
312
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
55 M4_U = [
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
56 0.1399708478939e13/0.263487168000e12 -0.13482796013041e14/0.834376032000e12 0.344344095859e12/0.17565811200e11 -0.3166261424681e13/0.250312809600e12 0.1508605165681e13/0.333750412800e12 -0.486270829441e12/0.834376032000e12 -0.221976356359e12/0.5006256192000e13;
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
57 -0.13482796013041e14/0.834376032000e12 0.7260475818391e13/0.139062672000e12 -0.27224036353e11/0.406022400e9 0.1847477458951e13/0.41718801600e11 -0.848984558161e12/0.55625068800e11 0.247494925991e12/0.139062672000e12 0.165585445559e12/0.834376032000e12;
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
58 0.344344095859e12/0.17565811200e11 -0.27224036353e11/0.406022400e9 0.2044938640393e13/0.22250027520e11 -0.1071086785417e13/0.16687520640e11 0.502199537033e12/0.22250027520e11 -0.143589154441e12/0.55625068800e11 -0.88181965559e11/0.333750412800e12;
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
59 -0.3166261424681e13/0.250312809600e12 0.1847477458951e13/0.41718801600e11 -0.1071086785417e13/0.16687520640e11 0.628860435593e12/0.12515640480e11 -0.73736245829e11/0.3337504128e10 0.195760572271e12/0.41718801600e11 -0.81156046361e11/0.250312809600e12;
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
60 0.1508605165681e13/0.333750412800e12 -0.848984558161e12/0.55625068800e11 0.502199537033e12/0.22250027520e11 -0.73736245829e11/0.3337504128e10 0.76725285869e11/0.4450005504e10 -0.3912429433e10/0.406022400e9 0.53227370659e11/0.17565811200e11;
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
61 -0.486270829441e12/0.834376032000e12 0.247494925991e12/0.139062672000e12 -0.143589154441e12/0.55625068800e11 0.195760572271e12/0.41718801600e11 -0.3912429433e10/0.406022400e9 0.1699707221791e13/0.139062672000e12 -0.6959018412841e13/0.834376032000e12;
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
62 -0.221976356359e12/0.5006256192000e13 0.165585445559e12/0.834376032000e12 -0.88181965559e11/0.333750412800e12 -0.81156046361e11/0.250312809600e12 0.53227370659e11/0.17565811200e11 -0.6959018412841e13/0.834376032000e12 0.3012195053939e13/0.263487168000e12;
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
63 ];
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
64
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
65 M4(1:7,1:7) = M4_U;
322
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
66 M4(m-6:m,m-6:m) = rot90(M4_U, 2);
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
67 M4 = 1/h^3*M4;
246
fe26791489e0 Added a bunch of new operators. Still non-functional.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
68
322
def409c10800 Clean up of d4_variable_6_3.m
Jonatan Werpers <jonatan@werpers.com>
parents: 318
diff changeset
69 D4=HI*(M4 - e_l*d3_l'+e_r*d3_r' + d1_l*d2_l'-d1_r*d2_r');
312
9230c056a574 Fixed formatting.
Jonatan Werpers <jonatan@werpers.com>
parents: 310
diff changeset
70 end