Mercurial > repos > public > sbplib_julia
comparison plotDerivative2d.jl @ 141:18b3c63673b3
Merged with cell_based_test
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Thu, 21 Feb 2019 16:55:03 +0100 |
| parents | d61bfc8cf6a5 |
| children |
comparison
equal
deleted
inserted
replaced
| 134:79699dda29be | 141:18b3c63673b3 |
|---|---|
| 1 include("sbpPlot.jl") | |
| 2 | |
| 1 g = sbp.Grid.EquidistantGrid((100,75), (0.0, 0.0), (2pi, 3/2*pi)) | 3 g = sbp.Grid.EquidistantGrid((100,75), (0.0, 0.0), (2pi, 3/2*pi)) |
| 2 op = sbp.readOperator("d2_4th.txt","h_4th.txt") | 4 op = sbp.readOperator("d2_4th.txt","h_4th.txt") |
| 3 Laplace = sbp.Laplace(g, 1.0, op) | 5 Laplace = sbp.Laplace(g, 1.0, op) |
| 4 | 6 |
| 5 init(x,y) = sin(x) + cos(y) | 7 init(x,y) = sin(x) + cos(y) |
| 6 v = sbp.Grid.evalOn(g,init) | 8 v = sbp.Grid.evalOn(g,init) |
| 7 | 9 u = zero(v) |
| 8 u = zeros(length(v)) | |
| 9 | 10 |
| 10 sbp.apply!(Laplace,u,v) | 11 sbp.apply!(Laplace,u,v) |
| 11 | 12 |
| 12 #@show u | 13 #@show u |
| 13 #@show u'*u | 14 #@show u'*u |
| 14 | 15 |
| 15 sbp.Grid.plotgridfunction(g,u) | 16 plotgridfunction(g,u) |
| 16 |