This repository was archived by the owner on Jul 19, 2023. It is now read-only.
File tree Expand file tree Collapse file tree 2 files changed +4
-4
lines changed Expand file tree Collapse file tree 2 files changed +4
-4
lines changed Original file line number Diff line number Diff line change @@ -19,8 +19,8 @@ function fourth_deriv_approx_stencil(N)
1919 A[1 ,1 : 8 ] = [3.5 - 56 / 3 42.5 - 54.0 251 / 6 - 20.0 5.5 - 2 / 3 ]
2020 A[2 ,1 : 8 ] = [2 / 3 - 11 / 6 0.0 31 / 6 - 22 / 3 4.5 - 4 / 3 1 / 6 ]
2121
22- A[N- 1 ,N- 5 : end ] = reverse ([3.5 - 56 / 3 42.5 - 54.0 251 / 6 - 20.0 5 .5 - 2 / 3 ], dims= 2 )
23- A[N,N- 5 : end ] = reverse ([2 / 3 - 11 / 6 0.0 31 / 6 - 22 / 3 4 .5 - 4 / 3 1 / 6 ], dims= 2 )
22+ A[N- 1 ,N- 5 : end ] = reverse ([2 / 3 - 11 / 6 0.0 31 / 6 - 22 / 3 4 .5 - 4 / 3 1 / 6 ], dims= 2 )
23+ A[N,N- 5 : end ] = reverse ([3.5 - 56 / 3 42.5 - 54.0 251 / 6 - 20.0 5 .5 - 2 / 3 ], dims= 2 )
2424
2525 for i in 3 : N- 2
2626 A[i,i- 2 : i+ 4 ] = [- 1 / 6 2.0 - 13 / 2 28 / 3 - 13 / 2 2.0 - 1 / 6 ]
Original file line number Diff line number Diff line change @@ -4,8 +4,8 @@ function fourth_deriv_approx_stencil(N)
44 A = zeros (N,N+ 2 )
55 A[1 ,1 : 8 ] = [3.5 - 56 / 3 42.5 - 54.0 251 / 6 - 20.0 5.5 - 2 / 3 ]
66 A[2 ,1 : 8 ] = [2 / 3 - 11 / 6 0.0 31 / 6 - 22 / 3 4.5 - 4 / 3 1 / 6 ]
7- A[N- 1 ,N- 5 : end ] = reverse ([3.5 - 56 / 3 42.5 - 54.0 251 / 6 - 20.0 5 .5 - 2 / 3 ], dims= 2 )
8- A[N,N- 5 : end ] = reverse ([2 / 3 - 11 / 6 0.0 31 / 6 - 22 / 3 4 .5 - 4 / 3 1 / 6 ], dims= 2 )
7+ A[N- 1 ,N- 5 : end ] = reverse ([2 / 3 - 11 / 6 0.0 31 / 6 - 22 / 3 4 .5 - 4 / 3 1 / 6 ], dims= 2 )
8+ A[N,N- 5 : end ] = reverse ([3.5 - 56 / 3 42.5 - 54.0 251 / 6 - 20.0 5 .5 - 2 / 3 ], dims= 2 )
99 for i in 3 : N- 2
1010 A[i,i- 2 : i+ 4 ] = [- 1 / 6 2.0 - 13 / 2 28 / 3 - 13 / 2 2.0 - 1 / 6 ]
1111 end
You can’t perform that action at this time.
0 commit comments