Mercurial > repos > public > sbplib_julia

comparison test/testSbpOperators.jl @ 396:9aff09ca0052 refactor/sbp_operators_tests/collect_and_compare

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove some type annotations
author Jonatan Werpers <jonatan@werpers.com>
date Sun, 04 Oct 2020 19:39:32 +0200
parents 737e597e0e6d
children 3cecbfb3d623
comparison
equal deleted inserted replaced
395:737e597e0e6d 396:9aff09ca0052
34 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt") 34 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
35 L = 3.5 35 L = 3.5
36 g = EquidistantGrid(101, 0.0, L) 36 g = EquidistantGrid(101, 0.0, L)
37 Dₓₓ = SecondDerivative(g,op.innerStencil,op.closureStencils) 37 Dₓₓ = SecondDerivative(g,op.innerStencil,op.closureStencils)
38 38
39 f0(x::Float64) = 1. 39 f0(x) = 1.
40 f1(x::Float64) = x 40 f1(x) = x
41 f2(x::Float64) = 1/2*x^2 41 f2(x) = 1/2*x^2
42 f3(x::Float64) = 1/6*x^3 42 f3(x) = 1/6*x^3
43 f4(x::Float64) = 1/24*x^4 43 f4(x) = 1/24*x^4
44 f5(x::Float64) = sin(x) 44 f5(x) = sin(x)
45 f5ₓₓ(x::Float64) = -f5(x) 45 f5ₓₓ(x) = -f5(x)
46 46
47 v0 = evalOn(g,f0) 47 v0 = evalOn(g,f0)
48 v1 = evalOn(g,f1) 48 v1 = evalOn(g,f1)
49 v2 = evalOn(g,f2) 49 v2 = evalOn(g,f2)
50 v3 = evalOn(g,f3) 50 v3 = evalOn(g,f3)
79 Ly = 3.2 79 Ly = 3.2
80 g = EquidistantGrid((102,131), (0.0, 0.0), (Lx,Ly)) 80 g = EquidistantGrid((102,131), (0.0, 0.0), (Lx,Ly))
81 L = Laplace(g, op.innerStencil, op.closureStencils) 81 L = Laplace(g, op.innerStencil, op.closureStencils)
82 82
83 83
84 f0(x::Float64,y::Float64) = 2. 84 f0(x,y) = 2.
85 f1(x::Float64,y::Float64) = x+y 85 f1(x,y) = x+y
86 f2(x::Float64,y::Float64) = 1/2*x^2 + 1/2*y^2 86 f2(x,y) = 1/2*x^2 + 1/2*y^2
87 f3(x::Float64,y::Float64) = 1/6*x^3 + 1/6*y^3 87 f3(x,y) = 1/6*x^3 + 1/6*y^3
88 f4(x::Float64,y::Float64) = 1/24*x^4 + 1/24*y^4 88 f4(x,y) = 1/24*x^4 + 1/24*y^4
89 f5(x::Float64,y::Float64) = sin(x) + cos(y) 89 f5(x,y) = sin(x) + cos(y)
90 f5ₓₓ(x::Float64,y::Float64) = -f5(x,y) 90 f5ₓₓ(x,y) = -f5(x,y)
91 91
92 v0 = evalOn(g,f0) 92 v0 = evalOn(g,f0)
93 v1 = evalOn(g,f1) 93 v1 = evalOn(g,f1)
94 v2 = evalOn(g,f2) 94 v2 = evalOn(g,f2)
95 v3 = evalOn(g,f3) 95 v3 = evalOn(g,f3)