Mercurial > repos > public > sbplib_julia

annotate test/testSbpOperators.jl @ 333:01b851161018 refactor/combine_to_one_package

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Start converting to one package by moving all the files to their correct location
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 25 Sep 2020 13:06:02 +0200
parents SbpOperators/test/runtests.jl@408c37b295c2
children f4e3e71a4ff4
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
314
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
1 using Test
219
69a6049e14d9 Create package SbpOperators
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 using SbpOperators
314
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
3 using Grids
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
4 using RegionIndices
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
5 using LazyTensors
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
6
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
7 # @testset "apply_quadrature" begin
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
8 # op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
9 # h = 0.5
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
10 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
11 # @test apply_quadrature(op, h, 1.0, 10, 100) == h
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
12 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
13 # N = 10
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
14 # qc = op.quadratureClosure
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
15 # q = h.*(qc..., ones(N-2*closuresize(op))..., reverse(qc)...)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
16 # @assert length(q) == N
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
17 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
18 # for i ∈ 1:N
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
19 # @test apply_quadrature(op, h, 1.0, i, N) == q[i]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
20 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
21 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
22 # v = [2.,3.,2.,4.,5.,4.,3.,4.,5.,4.5]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
23 # for i ∈ 1:N
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
24 # @test apply_quadrature(op, h, v[i], i, N) == q[i]*v[i]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
25 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
26 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
27
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
28 @testset "SecondDerivative" begin
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
29 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
30 L = 3.5
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
31 g = EquidistantGrid((101,), (0.0,), (L,))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
32 h_inv = inverse_spacing(g)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
33 h = 1/h_inv[1];
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
34 Dₓₓ = SecondDerivative(h_inv[1],op.innerStencil,op.closureStencils,op.parity)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
35
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
36 f0(x::Float64) = 1.
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
37 f1(x::Float64) = x
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
38 f2(x::Float64) = 1/2*x^2
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
39 f3(x::Float64) = 1/6*x^3
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
40 f4(x::Float64) = 1/24*x^4
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
41 f5(x::Float64) = sin(x)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
42 f5ₓₓ(x::Float64) = -f5(x)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
43
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
44 v0 = evalOn(g,f0)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
45 v1 = evalOn(g,f1)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
46 v2 = evalOn(g,f2)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
47 v3 = evalOn(g,f3)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
48 v4 = evalOn(g,f4)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
49 v5 = evalOn(g,f5)
219
69a6049e14d9 Create package SbpOperators
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
50
314
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
51 @test Dₓₓ isa TensorOperator{T,1} where T
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
52 @test Dₓₓ' isa TensorMapping{T,1,1} where T
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
53
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
54 # TODO: Should perhaps set tolerance level for isapporx instead?
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
55 # Are these tolerance levels resonable or should tests be constructed
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
56 # differently?
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
57 equalitytol = 0.5*1e-10
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
58 accuracytol = 0.5*1e-3
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
59 # 4th order interior stencil, 2nd order boundary stencil,
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
60 # implies that L*v should be exact for v - monomial up to order 3.
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
61 # Exact differentiation is measured point-wise. For other grid functions
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
62 # the error is measured in the l2-norm.
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
63 @test all(abs.(collect(Dₓₓ*v0)) .<= equalitytol)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
64 @test all(abs.(collect(Dₓₓ*v1)) .<= equalitytol)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
65 @test all(abs.((collect(Dₓₓ*v2) - v0)) .<= equalitytol)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
66 @test all(abs.((collect(Dₓₓ*v3) - v1)) .<= equalitytol)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
67 e4 = collect(Dₓₓ*v4) - v2
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
68 e5 = collect(Dₓₓ*v5) + v5
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
69 @test sqrt(h*sum(collect(e4.^2))) <= accuracytol
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
70 @test sqrt(h*sum(collect(e5.^2))) <= accuracytol
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
71 end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
72
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
73 @testset "Laplace2D" begin
254
4ca3794fffef Add apply_quadrature to SbpOperators
Jonatan Werpers <jonatan@werpers.com>
parents: 219
diff changeset
74 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
314
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
75 Lx = 1.5
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
76 Ly = 3.2
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
77 g = EquidistantGrid((102,131), (0.0, 0.0), (Lx,Ly))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
78
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
79 h_inv = inverse_spacing(g)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
80 h = spacing(g)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
81 D_xx = SecondDerivative(h_inv[1],op.innerStencil,op.closureStencils,op.parity)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
82 D_yy = SecondDerivative(h_inv[2],op.innerStencil,op.closureStencils,op.parity)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
83 L = Laplace((D_xx,D_yy))
254
4ca3794fffef Add apply_quadrature to SbpOperators
Jonatan Werpers <jonatan@werpers.com>
parents: 219
diff changeset
84
314
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
85 f0(x::Float64,y::Float64) = 2.
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
86 f1(x::Float64,y::Float64) = x+y
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
87 f2(x::Float64,y::Float64) = 1/2*x^2 + 1/2*y^2
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
88 f3(x::Float64,y::Float64) = 1/6*x^3 + 1/6*y^3
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
89 f4(x::Float64,y::Float64) = 1/24*x^4 + 1/24*y^4
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
90 f5(x::Float64,y::Float64) = sin(x) + cos(y)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
91 f5ₓₓ(x::Float64,y::Float64) = -f5(x,y)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
92
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
93 v0 = evalOn(g,f0)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
94 v1 = evalOn(g,f1)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
95 v2 = evalOn(g,f2)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
96 v3 = evalOn(g,f3)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
97 v4 = evalOn(g,f4)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
98 v5 = evalOn(g,f5)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
99 v5ₓₓ = evalOn(g,f5ₓₓ)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
100
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
101 @test L isa TensorOperator{T,2} where T
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
102 @test L' isa TensorMapping{T,2,2} where T
254
4ca3794fffef Add apply_quadrature to SbpOperators
Jonatan Werpers <jonatan@werpers.com>
parents: 219
diff changeset
103
314
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
104 # TODO: Should perhaps set tolerance level for isapporx instead?
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
105 # Are these tolerance levels resonable or should tests be constructed
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
106 # differently?
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
107 equalitytol = 0.5*1e-10
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
108 accuracytol = 0.5*1e-3
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
109 # 4th order interior stencil, 2nd order boundary stencil,
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
110 # implies that L*v should be exact for v - monomial up to order 3.
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
111 # Exact differentiation is measured point-wise. For other grid functions
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
112 # the error is measured in the H-norm.
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
113 @test all(abs.(collect(L*v0)) .<= equalitytol)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
114 @test all(abs.(collect(L*v1)) .<= equalitytol)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
115 @test all(collect(L*v2) .≈ v0) # Seems to be more accurate
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
116 @test all(abs.((collect(L*v3) - v1)) .<= equalitytol)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
117 e4 = collect(L*v4) - v2
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
118 e5 = collect(L*v5) - v5ₓₓ
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
119 @test sqrt(prod(h)*sum(collect(e4.^2))) <= accuracytol
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
120 @test sqrt(prod(h)*sum(collect(e5.^2))) <= accuracytol
254
4ca3794fffef Add apply_quadrature to SbpOperators
Jonatan Werpers <jonatan@werpers.com>
parents: 219
diff changeset
121 end
328
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
122
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
123 @testset "DiagonalInnerProduct" begin
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
124 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
125 L = 2.3
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
126 g = EquidistantGrid((77,), (0.0,), (L,))
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
127 h = spacing(g)
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
128 H = DiagonalInnerProduct(h[1],op.quadratureClosure)
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
129 v = ones(Float64, size(g))
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
130
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
131 @test H isa TensorOperator{T,1} where T
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
132 @test H' isa TensorMapping{T,1,1} where T
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
133 @test sum(collect(H*v)) ≈ L
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
134 @test collect(H*v) == collect(H'*v)
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
135 end
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
136
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
137 @testset "Quadrature" begin
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
138 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
139 Lx = 2.3
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
140 Ly = 5.2
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
141 g = EquidistantGrid((77,66), (0.0, 0.0), (Lx,Ly))
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
142
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
143 h = spacing(g)
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
144 Hx = DiagonalInnerProduct(h[1],op.quadratureClosure);
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
145 Hy = DiagonalInnerProduct(h[2],op.quadratureClosure);
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
146 Q = Quadrature((Hx,Hy))
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
147
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
148 v = ones(Float64, size(g))
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
149
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
150 @test Q isa TensorOperator{T,2} where T
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
151 @test Q' isa TensorMapping{T,2,2} where T
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
152 @test sum(collect(Q*v)) ≈ (Lx*Ly)
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
153 @test collect(Q*v) == collect(Q'*v)
9cc5d1498b2d Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 314
diff changeset
154 end
329
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
155
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
156 @testset "InverseDiagonalInnerProduct" begin
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
157 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
158 L = 2.3
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
159 g = EquidistantGrid((77,), (0.0,), (L,))
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
160 h = spacing(g)
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
161 H = DiagonalInnerProduct(h[1],op.quadratureClosure)
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
162
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
163 h_i = inverse_spacing(g)
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
164 Hi = InverseDiagonalInnerProduct(h_i[1],1 ./ op.quadratureClosure)
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
165 v = evalOn(g, x->sin(x))
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
166
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
167 @test Hi isa TensorOperator{T,1} where T
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
168 @test Hi' isa TensorMapping{T,1,1} where T
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
169 @test collect(Hi*H*v) ≈ v
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
170 @test collect(Hi*v) == collect(Hi'*v)
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
171 end
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
172
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
173 @testset "InverseQuadrature" begin
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
174 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
175 Lx = 7.3
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
176 Ly = 8.2
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
177 g = EquidistantGrid((77,66), (0.0, 0.0), (Lx,Ly))
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
178
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
179 h = spacing(g)
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
180 Hx = DiagonalInnerProduct(h[1], op.quadratureClosure);
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
181 Hy = DiagonalInnerProduct(h[2], op.quadratureClosure);
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
182 Q = Quadrature((Hx,Hy))
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
183
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
184 hi = inverse_spacing(g)
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
185 Hix = InverseDiagonalInnerProduct(hi[1], 1 ./ op.quadratureClosure);
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
186 Hiy = InverseDiagonalInnerProduct(hi[2], 1 ./ op.quadratureClosure);
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
187 Qinv = InverseQuadrature((Hix,Hiy))
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
188 v = evalOn(g, (x,y)-> x^2 + (y-1)^2 + x*y)
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
189
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
190 @test Qinv isa TensorOperator{T,2} where T
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
191 @test Qinv' isa TensorMapping{T,2,2} where T
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
192 @test collect(Qinv*Q*v) ≈ v
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
193 @test collect(Qinv*v) == collect(Qinv'*v)
408c37b295c2 Refactor 1D tensor mapping in inverse quadrature to separate file, InverseDiagonalNorm. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 328
diff changeset
194 end
314
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
195 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
196 # @testset "BoundaryValue" begin
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
197 # op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
198 # g = EquidistantGrid((4,5), (0.0, 0.0), (1.0,1.0))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
199 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
200 # e_w = BoundaryValue(op, g, CartesianBoundary{1,Lower}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
201 # e_e = BoundaryValue(op, g, CartesianBoundary{1,Upper}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
202 # e_s = BoundaryValue(op, g, CartesianBoundary{2,Lower}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
203 # e_n = BoundaryValue(op, g, CartesianBoundary{2,Upper}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
204 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
205 # v = zeros(Float64, 4, 5)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
206 # v[:,5] = [1, 2, 3,4]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
207 # v[:,4] = [1, 2, 3,4]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
208 # v[:,3] = [4, 5, 6, 7]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
209 # v[:,2] = [7, 8, 9, 10]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
210 # v[:,1] = [10, 11, 12, 13]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
211 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
212 # @test e_w isa TensorMapping{T,2,1} where T
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
213 # @test e_w' isa TensorMapping{T,1,2} where T
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
214 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
215 # @test domain_size(e_w, (3,2)) == (2,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
216 # @test domain_size(e_e, (3,2)) == (2,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
217 # @test domain_size(e_s, (3,2)) == (3,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
218 # @test domain_size(e_n, (3,2)) == (3,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
219 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
220 # @test size(e_w'*v) == (5,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
221 # @test size(e_e'*v) == (5,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
222 # @test size(e_s'*v) == (4,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
223 # @test size(e_n'*v) == (4,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
224 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
225 # @test collect(e_w'*v) == [10,7,4,1.0,1]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
226 # @test collect(e_e'*v) == [13,10,7,4,4.0]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
227 # @test collect(e_s'*v) == [10,11,12,13.0]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
228 # @test collect(e_n'*v) == [1,2,3,4.0]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
229 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
230 # g_x = [1,2,3,4.0]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
231 # g_y = [5,4,3,2,1.0]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
232 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
233 # G_w = zeros(Float64, (4,5))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
234 # G_w[1,:] = g_y
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
235 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
236 # G_e = zeros(Float64, (4,5))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
237 # G_e[4,:] = g_y
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
238 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
239 # G_s = zeros(Float64, (4,5))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
240 # G_s[:,1] = g_x
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
241 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
242 # G_n = zeros(Float64, (4,5))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
243 # G_n[:,5] = g_x
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
244 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
245 # @test size(e_w*g_y) == (UnknownDim,5)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
246 # @test size(e_e*g_y) == (UnknownDim,5)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
247 # @test size(e_s*g_x) == (4,UnknownDim)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
248 # @test size(e_n*g_x) == (4,UnknownDim)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
249 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
250 # # These tests should be moved to where they are possible (i.e we know what the grid should be)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
251 # @test_broken collect(e_w*g_y) == G_w
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
252 # @test_broken collect(e_e*g_y) == G_e
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
253 # @test_broken collect(e_s*g_x) == G_s
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
254 # @test_broken collect(e_n*g_x) == G_n
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
255 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
256 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
257 # @testset "NormalDerivative" begin
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
258 # op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
259 # g = EquidistantGrid((5,6), (0.0, 0.0), (4.0,5.0))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
260 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
261 # d_w = NormalDerivative(op, g, CartesianBoundary{1,Lower}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
262 # d_e = NormalDerivative(op, g, CartesianBoundary{1,Upper}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
263 # d_s = NormalDerivative(op, g, CartesianBoundary{2,Lower}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
264 # d_n = NormalDerivative(op, g, CartesianBoundary{2,Upper}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
265 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
266 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
267 # v = evalOn(g, (x,y)-> x^2 + (y-1)^2 + x*y)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
268 # v∂x = evalOn(g, (x,y)-> 2*x + y)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
269 # v∂y = evalOn(g, (x,y)-> 2*(y-1) + x)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
270 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
271 # @test d_w isa TensorMapping{T,2,1} where T
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
272 # @test d_w' isa TensorMapping{T,1,2} where T
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
273 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
274 # @test domain_size(d_w, (3,2)) == (2,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
275 # @test domain_size(d_e, (3,2)) == (2,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
276 # @test domain_size(d_s, (3,2)) == (3,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
277 # @test domain_size(d_n, (3,2)) == (3,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
278 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
279 # @test size(d_w'*v) == (6,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
280 # @test size(d_e'*v) == (6,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
281 # @test size(d_s'*v) == (5,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
282 # @test size(d_n'*v) == (5,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
283 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
284 # @test collect(d_w'*v) ≈ v∂x[1,:]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
285 # @test collect(d_e'*v) ≈ v∂x[5,:]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
286 # @test collect(d_s'*v) ≈ v∂y[:,1]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
287 # @test collect(d_n'*v) ≈ v∂y[:,6]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
288 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
289 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
290 # d_x_l = zeros(Float64, 5)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
291 # d_x_u = zeros(Float64, 5)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
292 # for i ∈ eachindex(d_x_l)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
293 # d_x_l[i] = op.dClosure[i-1]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
294 # d_x_u[i] = -op.dClosure[length(d_x_u)-i]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
295 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
296 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
297 # d_y_l = zeros(Float64, 6)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
298 # d_y_u = zeros(Float64, 6)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
299 # for i ∈ eachindex(d_y_l)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
300 # d_y_l[i] = op.dClosure[i-1]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
301 # d_y_u[i] = -op.dClosure[length(d_y_u)-i]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
302 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
303 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
304 # function prod_matrix(x,y)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
305 # G = zeros(Float64, length(x), length(y))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
306 # for I ∈ CartesianIndices(G)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
307 # G[I] = x[I[1]]*y[I[2]]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
308 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
309 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
310 # return G
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
311 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
312 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
313 # g_x = [1,2,3,4.0,5]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
314 # g_y = [5,4,3,2,1.0,11]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
315 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
316 # G_w = prod_matrix(d_x_l, g_y)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
317 # G_e = prod_matrix(d_x_u, g_y)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
318 # G_s = prod_matrix(g_x, d_y_l)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
319 # G_n = prod_matrix(g_x, d_y_u)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
320 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
321 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
322 # @test size(d_w*g_y) == (UnknownDim,6)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
323 # @test size(d_e*g_y) == (UnknownDim,6)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
324 # @test size(d_s*g_x) == (5,UnknownDim)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
325 # @test size(d_n*g_x) == (5,UnknownDim)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
326 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
327 # # These tests should be moved to where they are possible (i.e we know what the grid should be)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
328 # @test_broken collect(d_w*g_y) ≈ G_w
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
329 # @test_broken collect(d_e*g_y) ≈ G_e
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
330 # @test_broken collect(d_s*g_x) ≈ G_s
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
331 # @test_broken collect(d_n*g_x) ≈ G_n
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
332 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
333 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
334 # @testset "BoundaryQuadrature" begin
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
335 # op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
336 # g = EquidistantGrid((10,11), (0.0, 0.0), (1.0,1.0))
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
337 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
338 # H_w = BoundaryQuadrature(op, g, CartesianBoundary{1,Lower}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
339 # H_e = BoundaryQuadrature(op, g, CartesianBoundary{1,Upper}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
340 # H_s = BoundaryQuadrature(op, g, CartesianBoundary{2,Lower}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
341 # H_n = BoundaryQuadrature(op, g, CartesianBoundary{2,Upper}())
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
342 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
343 # v = evalOn(g, (x,y)-> x^2 + (y-1)^2 + x*y)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
344 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
345 # function get_quadrature(N)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
346 # qc = op.quadratureClosure
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
347 # q = (qc..., ones(N-2*closuresize(op))..., reverse(qc)...)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
348 # @assert length(q) == N
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
349 # return q
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
350 # end
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
351 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
352 # v_w = v[1,:]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
353 # v_e = v[10,:]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
354 # v_s = v[:,1]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
355 # v_n = v[:,11]
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
356 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
357 # q_x = spacing(g)[1].*get_quadrature(10)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
358 # q_y = spacing(g)[2].*get_quadrature(11)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
359 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
360 # @test H_w isa TensorOperator{T,1} where T
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
361 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
362 # @test domain_size(H_w, (3,)) == (3,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
363 # @test domain_size(H_n, (3,)) == (3,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
364 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
365 # @test range_size(H_w, (3,)) == (3,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
366 # @test range_size(H_n, (3,)) == (3,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
367 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
368 # @test size(H_w*v_w) == (11,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
369 # @test size(H_e*v_e) == (11,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
370 # @test size(H_s*v_s) == (10,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
371 # @test size(H_n*v_n) == (10,)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
372 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
373 # @test collect(H_w*v_w) ≈ q_y.*v_w
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
374 # @test collect(H_e*v_e) ≈ q_y.*v_e
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
375 # @test collect(H_s*v_s) ≈ q_x.*v_s
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
376 # @test collect(H_n*v_n) ≈ q_x.*v_n
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
377 #
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
378 # @test collect(H_w'*v_w) == collect(H_w'*v_w)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
379 # @test collect(H_e'*v_e) == collect(H_e'*v_e)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
380 # @test collect(H_s'*v_s) == collect(H_s'*v_s)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
381 # @test collect(H_n'*v_n) == collect(H_n'*v_n)
accb0876da12 Move tests for Laplace from DiffOps/test to SbpOperators/test and add test for Secondderivative
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 254
diff changeset
382 # end