Mercurial > repos > public > sbplib_julia

annotate test/Grids/geometry_test.jl @ 2015:5c2448d6a201 feature/grids/geometry_functions tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Structure tests a bit more
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 09 May 2025 15:57:38 +0200
parents 7895b509f9bf
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1966
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
1 using Diffinitive.Grids
2011
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
2 using Diffinitive.Grids: Line, LineSegment, linesegments, polygon_edges, Circle, TransfiniteInterpolationSurface, check_transfiniteinterpolation, arc, Arc
1966
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
3 using StaticArrays
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
4
1915
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 @testset "Line" begin
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
6 @testset "Constructors" begin
1966
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
7 @test Line([1,2],[2,3]) isa Line{SVector{2,Int}}
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
8 @test Line((1,2),(2,3)) isa Line{SVector{2,Int}}
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
9 @test Line(@SVector[1,2],[2,3]) isa Line{SVector{2,Int}}
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
10 @test Line(@SVector[1,2],@SVector[2,3]) isa Line{SVector{2,Int}}
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
11
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
12 @test Line([1,2],[2.,3]) isa Line{SVector{2,Float64}}
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
13 @test Line(@SVector[1,2.],@SVector[2,3]) isa Line{SVector{2,Float64}}
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
14 @test Line((1,2.),(2,3)) isa Line{SVector{2,Float64}}
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
15 end
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
16
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
17 @testset "Evaluation" begin
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
18 l = Line([1,2],[2,3])
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
19
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
20 @test l(0) == [1,2]
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
21 @test l(1) == [1,2] + [2,3]
478b233999c5 Add tests and better constructors for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1965
diff changeset
22 @test l(1/2) == [1,2] + [2,3]/2
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
23 end
1974
3ed7ca1f60c4 Implement Grids.jacobian for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1973
diff changeset
24
3ed7ca1f60c4 Implement Grids.jacobian for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1973
diff changeset
25 @testset "Grids.jacobian" begin
3ed7ca1f60c4 Implement Grids.jacobian for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1973
diff changeset
26 l = Line([1,2],[2,3])
3ed7ca1f60c4 Implement Grids.jacobian for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1973
diff changeset
27
3ed7ca1f60c4 Implement Grids.jacobian for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1973
diff changeset
28 @test Grids.jacobian(l,0) == [2,3]
3ed7ca1f60c4 Implement Grids.jacobian for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1973
diff changeset
29 @test Grids.jacobian(l,1) == [2,3]
3ed7ca1f60c4 Implement Grids.jacobian for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1973
diff changeset
30 @test Grids.jacobian(l,1/2) == [2,3]
3ed7ca1f60c4 Implement Grids.jacobian for Line
Jonatan Werpers <jonatan@werpers.com>
parents: 1973
diff changeset
31 end
1915
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
32 end
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
33
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
34 @testset "LineSegment" begin
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
35 @testset "Constructors" begin
1967
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
36 @test LineSegment([1,2],[2,3]) isa LineSegment{SVector{2,Int}}
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
37 @test LineSegment((1,2),(2,3)) isa LineSegment{SVector{2,Int}}
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
38 @test LineSegment(@SVector[1,2],[2,3]) isa LineSegment{SVector{2,Int}}
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
39 @test LineSegment(@SVector[1,2],@SVector[2,3]) isa LineSegment{SVector{2,Int}}
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
40
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
41 @test LineSegment([1,2],[2.,3]) isa LineSegment{SVector{2,Float64}}
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
42 @test LineSegment(@SVector[1,2.],@SVector[2,3]) isa LineSegment{SVector{2,Float64}}
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
43 @test LineSegment((1,2.),(2,3)) isa LineSegment{SVector{2,Float64}}
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
44 end
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
45
1967
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
46 @testset "Evaluation" begin
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
47 l = LineSegment([1,2],[2,3])
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
48
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
49 @test l(0) == [1,2]
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
50 @test l(1) == [2,3]
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
51 @test l(1/2) == [1,2]/2 + [2,3]/2
669361a8195a Add tests and better constructors for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1966
diff changeset
52 end
1975
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
53
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
54 @testset "Grids.jacobian" begin
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
55 a, b = [1,2], [2,3]
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
56 l = LineSegment(a,b)
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
57 d = b-a
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
58
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
59 @test Grids.jacobian(l,0) == d
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
60 @test Grids.jacobian(l,1) == d
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
61 @test Grids.jacobian(l,1/2) == d
79adad4a15d6 Implement Grids.jacobian for LineSegment
Jonatan Werpers <jonatan@werpers.com>
parents: 1974
diff changeset
62 end
1915
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
63 end
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
64
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
65 @testset "linesegments" begin
1968
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
66 a,b,c,d = [1,1],[2,2],[3,3],[4,4]
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
67 @test linesegments(a,b) == [
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
68 LineSegment(a,b),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
69 ]
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
70
1968
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
71 @test linesegments(a,b,c) == [
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
72 LineSegment(a,b),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
73 LineSegment(b,c),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
74 ]
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
75
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
76 @test linesegments(a,b,c,d) == [
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
77 LineSegment(a,b),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
78 LineSegment(b,c),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
79 LineSegment(c,d),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
80 ]
1915
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
81 end
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
82
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
83 @testset "polygon_edges" begin
1968
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
84 a,b,c,d = [1,1],[2,2],[3,3],[4,4]
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
85 @test polygon_edges(a,b) == [
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
86 LineSegment(a,b),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
87 LineSegment(b,a),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
88 ]
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
89
1968
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
90 @test polygon_edges(a,b,c) == [
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
91 LineSegment(a,b),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
92 LineSegment(b,c),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
93 LineSegment(c,a),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
94 ]
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
95
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
96 @test polygon_edges(a,b,c,d) == [
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
97 LineSegment(a,b),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
98 LineSegment(b,c),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
99 LineSegment(c,d),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
100 LineSegment(d,a),
35cb503985b6 Add tests for linesegments and polygon_edges
Jonatan Werpers <jonatan@werpers.com>
parents: 1967
diff changeset
101 ]
1915
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
102 end
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
103
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
104 @testset "Circle" begin
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
105 @testset "Constructors" begin
1969
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
106 @test Circle([1,2], 1) isa Circle{SVector{2,Int},Int}
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
107 @test Circle([1,2], 1.) isa Circle{SVector{2,Int},Float64}
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
108 @test Circle([1,2.], 1.) isa Circle{SVector{2,Float64},Float64}
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
109 @test Circle([1,2.], 1) isa Circle{SVector{2,Float64},Int}
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
110 @test Circle((1,2.), 1.) isa Circle{SVector{2,Float64},Float64}
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
111 @test Circle((1,2), 1.) isa Circle{SVector{2,Int},Float64}
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
112 @test Circle((1.,2), 1) isa Circle{SVector{2,Float64},Int}
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
113 @test Circle((1,2), 1) isa Circle{SVector{2,Int},Int}
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
114 @test Circle(@SVector[1,2], 1.) isa Circle{SVector{2,Int},Float64}
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
115 @test Circle(@SVector[1,2.], 1.) isa Circle{SVector{2,Float64},Float64}
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
116 end
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
117
1969
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
118 @testset "Evaluation" begin
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
119 c = Circle([0,0], 1)
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
120 @test c(0) ≈ [1,0]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
121 @test c(π/2) ≈ [0,1]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
122 @test c(π) ≈ [-1,0]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
123 @test c(3π/2) ≈ [0,-1]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
124 @test c(π/4) ≈ [1/√(2),1/√(2)]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
125
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
126 c = Circle([0,0], 2)
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
127 @test c(0) ≈ [2,0]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
128 @test c(π/2) ≈ [0,2]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
129 @test c(π) ≈ [-2,0]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
130 @test c(3π/2) ≈ [0,-2]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
131 @test c(π/4) ≈ [√(2),√(2)]
7f4a5146c84c Add tests and better constructor for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1968
diff changeset
132 end
1976
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
133
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
134 @testset "Grids.jacobian" begin
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
135 c = Circle([0,0], 1)
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
136 @test Grids.jacobian(c, 0) ≈ [0,1]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
137 @test Grids.jacobian(c, π/2) ≈ [-1,0]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
138 @test Grids.jacobian(c, π) ≈ [0,-1]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
139 @test Grids.jacobian(c, 3π/2) ≈ [1,0]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
140 @test Grids.jacobian(c, π/4) ≈ [-1/√(2),1/√(2)]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
141
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
142 c = Circle([0,0], 2)
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
143 @test Grids.jacobian(c, 0) ≈ [0,2]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
144 @test Grids.jacobian(c, π/2) ≈ [-2,0]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
145 @test Grids.jacobian(c, π) ≈ [0,-2]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
146 @test Grids.jacobian(c, 3π/2) ≈ [2,0]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
147 @test Grids.jacobian(c, π/4) ≈ [-√(2),√(2)]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
148
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
149 c = Circle([-1,1], 1)
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
150 @test Grids.jacobian(c, 0) ≈ [0,1]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
151 @test Grids.jacobian(c, π/2) ≈ [-1,0]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
152 @test Grids.jacobian(c, π) ≈ [0,-1]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
153 @test Grids.jacobian(c, 3π/2) ≈ [1,0]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
154 @test Grids.jacobian(c, π/4) ≈ [-1/√(2),1/√(2)]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
155
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
156 c = Circle([-1,1], 2)
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
157 @test Grids.jacobian(c, 0) ≈ [0,2]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
158 @test Grids.jacobian(c, π/2) ≈ [-2,0]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
159 @test Grids.jacobian(c, π) ≈ [0,-2]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
160 @test Grids.jacobian(c, 3π/2) ≈ [2,0]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
161 @test Grids.jacobian(c, π/4) ≈ [-√(2),√(2)]
34a7e3919e9a Implement Grids.jacobian for Circle
Jonatan Werpers <jonatan@werpers.com>
parents: 1975
diff changeset
162 end
1915
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
163 end
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
164
2011
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
165 @testset "Arc" begin
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
166 @test Arc(Circle([0,0], 1), 0, 1) isa Arc{SVector{2,Int}, Int}
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
167 @test Arc(Circle([0,0], 1.), 0, 1) isa Arc{SVector{2,Int}, Float64}
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
168 @test Arc(Circle([1., 1.], 1), 0., 1.) isa Arc{SVector{2,Float64}, Float64}
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
169 @test Arc(Circle([1., 1.], 1), 0, 1) isa Arc{SVector{2,Float64}, Int}
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
170 @test Arc(Circle([1., 1.], 1), 0, 1.) isa Arc{SVector{2,Float64}, Float64}
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
171
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
172 a = Arc(Circle([0,0], 1), 0, π/2)
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
173 @test a(0) ≈ [1,0]
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
174 @test a(1/3) ≈ [√(3)/2,1/2]
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
175 @test a(1/2) ≈ [1/√(2),1/√(2)]
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
176 @test a(2/3) ≈ [1/2, √(3)/2]
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
177 @test a(1) ≈ [0,1]
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
178
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
179 @testset "Grids.jacobian" begin
2013
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
180 c = Circle([0,0], 1)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
181
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
182 @testset "Matched to circle" begin
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
183 a = Arc(c, 0, 1)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
184 @testset for t ∈ range(0,1,8)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
185 @test jacobian(a,t) ≈ jacobian(c,t)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
186 end
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
187 end
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
188
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
189 @testset "Full circle" begin
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
190 a = Arc(c, 0, 2π)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
191 @testset for t ∈ range(0,1,8)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
192 @test jacobian(a,t) ≈ 2π*jacobian(c,t)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
193 end
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
194 end
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
195
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
196 @testset "Other" begin
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
197 a = Arc(c, π/3, 5π/4)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
198 @testset for t ∈ range(0,1,8)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
199 @test jacobian(a,t) ≈ 11π/12*jacobian(c,t)
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
200 end
7895b509f9bf Implement Grids.jacobian for Arc
Jonatan Werpers <jonatan@werpers.com>
parents: 2012
diff changeset
201 end
2011
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
202 end
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
203 end
d0b6c63c506e Add Arc type for representing partial circles
Jonatan Werpers <jonatan@werpers.com>
parents: 2010
diff changeset
204
2012
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
205 @testset "arc" begin
2015
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
206 @testset "Half circles around [0.5, 0.0]" begin
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
207 a = [0,0]
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
208 b = [1,0]
2012
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
209
2015
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
210 A = arc(a,b,1/2)
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
211 @test A(0) ≈ a atol=1e-15
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
212 @test A(1) ≈ b
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
213 @test A(0.5) ≈ [0.5, -0.5]
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
214
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
215 A = arc(a,b,-1/2)
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
216 @test A(0) ≈ a atol=1e-15
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
217 @test A(1) ≈ b
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
218 @test A(0.5) ≈ [0.5, 0.5]
5c2448d6a201 Structure tests a bit more
Jonatan Werpers <jonatan@werpers.com>
parents: 2013
diff changeset
219 end
2012
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
220
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
221 @testset "Unit arc" begin
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
222 A = arc([1,0],[0,1],1)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
223 @test A(0) ≈ [1,0]
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
224 @test A(1) ≈ [0,1]
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
225 @testset for t ∈ range(0,1,13)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
226 @test A(t) ≈ [cos(t*π/2), sin(t*π/2)]
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
227 end
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
228 end
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
229
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
230 @testset "Inverted unit arc" begin
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
231 A = arc([1,0],[0,1],-1)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
232 @test A(0) ≈ [1,0]
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
233 @test A(1) ≈ [0,1]
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
234 @testset "Inverted unit arc t=$t" for t ∈ range(0,1,13)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
235 @test A(t) ≈ [1+cos(-π/2 - t*π/2), 1+sin(-π/2 - t*π/2)]
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
236 end
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
237 end
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
238
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
239 @testset "Quarters of unit circle" begin
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
240 unitvec(θ) = [cos(θ), sin(θ)]
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
241 @testset "θ ∈ ($(i)π/4, $(i+2)π/4)" for i ∈ range(0, step=1, length=8)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
242 θ = i*π/4
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
243 @testset let θ₀ = θ, θ₁ = θ+π/2, r = 1
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
244 A = arc(unitvec(θ₀), unitvec(θ₁), r)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
245 @test A(0) ≈ unitvec(θ)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
246 @test A(1/3) ≈ unitvec(θ+π/6)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
247 @test A(1/2) ≈ unitvec(θ+π/4)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
248 @test A(2/3) ≈ unitvec(θ+π/3)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
249 @test A(1) ≈ unitvec(θ+π/2)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
250 end
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
251
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
252 @testset let θ₀ = θ+π/2, θ₁ = θ, r = -1
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
253 A = arc(unitvec(θ₀), unitvec(θ₁), r)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
254 @test A(0) ≈ unitvec(θ+π/2)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
255 @test A(1/3) ≈ unitvec(θ+π/3)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
256 @test A(1/2) ≈ unitvec(θ+π/4)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
257 @test A(2/3) ≈ unitvec(θ+π/6)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
258 @test A(1) ≈ unitvec(θ)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
259 end
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
260 end
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
261 end
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
262
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
263 @test_throws DomainError arc([-1,0], [1,0], 0.7)
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
264 end
4617e4b74b82 Add arc() for constructing circle arcs between two points
Jonatan Werpers <jonatan@werpers.com>
parents: 2011
diff changeset
265
1915
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
266 @testset "TransfiniteInterpolationSurface" begin
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
267 @testset "Constructors" begin
1973
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
268 @test TransfiniteInterpolationSurface(t->[1,2], t->[2,1], t->[0,0], t->[1,1]) isa TransfiniteInterpolationSurface
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
269
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
270 cs = polygon_edges([0,0],[1,0],[1,1],[0,1])
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
271 @test TransfiniteInterpolationSurface(cs...) isa TransfiniteInterpolationSurface
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
272 end
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
273
1973
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
274 @testset "Evaluation" begin
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
275 cs = polygon_edges([0,0],[1,0],[1,1],[0,1])
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
276 ti = TransfiniteInterpolationSurface(cs...)
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
277
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
278 @test ti(0,0) == [0,0]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
279 @test ti([0,0]) == [0,0]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
280 @test ti(1,0) == [1,0]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
281 @test ti([1,0]) == [1,0]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
282 @test ti(1,1) == [1,1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
283 @test ti([1,1]) == [1,1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
284 @test ti(0,1) == [0,1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
285 @test ti([0,1]) == [0,1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
286
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
287 @test ti(1/2, 0) == [1/2, 0]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
288 @test ti(1/2, 1) == [1/2, 1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
289 @test ti(0,1/2) == [0,1/2]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
290 @test ti(1,1/2) == [1,1/2]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
291
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
292
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
293 a, b, c, d = [1,0],[2,1/4],[2.5,1],[-1/3,1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
294 cs = polygon_edges(a,b,c,d)
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
295 ti = TransfiniteInterpolationSurface(cs...)
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
296
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
297 @test ti(0,0) == a
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
298 @test ti(1,0) == b
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
299 @test ti(1,1) == c
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
300 @test ti(0,1) == d
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
301
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
302 @test ti(1/2, 0) == (a+b)/2
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
303 @test ti(1/2, 1) == (c+d)/2
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
304 @test ti(0, 1/2) == (a+d)/2
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
305 @test ti(1, 1/2) == (b+c)/2
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
306
2010
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
307 a, b, c, d = [0,0],[1,1/2],[1,3/2],[0,1]
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
308 ti = TransfiniteInterpolationSurface(
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
309 t->@SVector[t, t^2/2],
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
310 LineSegment(b,c),
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
311 LineSegment(c,d),
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
312 LineSegment(d,a),
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
313 )
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
314
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
315 @test ti(0,0) == a
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
316 @test ti(1,0) == b
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
317 @test ti(1,1) == c
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
318 @test ti(0,1) == d
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
319
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
320 @test ti(1/2, 0) == [1/2, 1/8]
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
321 @test ti(1/2, 1) == (c+d)/2
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
322 @test ti(0, 1/2) == (a+d)/2
cf0fa2967361 Add test for transfinite interpolation with a non-linear edge
Jonatan Werpers <jonatan@werpers.com>
parents: 1982
diff changeset
323 @test ti(1, 1/2) == (b+c)/2
1973
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
324 end
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
325
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
326 @testset "check_transfiniteinterpolation" begin
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
327 cs = polygon_edges([0,0],[1,0],[1,1],[0,1])
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
328 ti = TransfiniteInterpolationSurface(cs...)
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
329
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
330 @test check_transfiniteinterpolation(ti) == nothing
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
331 @test check_transfiniteinterpolation(Bool, ti) == true
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
332
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
333 bad_sides = [
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
334 LineSegment([0,0],[1/2,0]),
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
335 LineSegment([1,0],[1,1/2]),
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
336 LineSegment([1,1],[0,1/2]),
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
337 LineSegment([0,1],[0,1/2]),
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
338 ]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
339
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
340 s1 = TransfiniteInterpolationSurface(bad_sides[1],cs[2],cs[3],cs[4])
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
341 s2 = TransfiniteInterpolationSurface(cs[1],bad_sides[2],cs[3],cs[4])
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
342 s3 = TransfiniteInterpolationSurface(cs[1],cs[2],bad_sides[3],cs[4])
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
343 s4 = TransfiniteInterpolationSurface(cs[1],cs[2],cs[3],bad_sides[4])
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
344
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
345 @test check_transfiniteinterpolation(Bool, s1) == false
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
346 @test check_transfiniteinterpolation(Bool, s2) == false
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
347 @test check_transfiniteinterpolation(Bool, s3) == false
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
348 @test check_transfiniteinterpolation(Bool, s4) == false
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
349
1982
486b3c6f919e Throw better error type for transfinite interpolation check
Jonatan Werpers <jonatan@werpers.com>
parents: 1977
diff changeset
350 @test_throws ArgumentError check_transfiniteinterpolation(s1)
486b3c6f919e Throw better error type for transfinite interpolation check
Jonatan Werpers <jonatan@werpers.com>
parents: 1977
diff changeset
351 @test_throws ArgumentError check_transfiniteinterpolation(s2)
486b3c6f919e Throw better error type for transfinite interpolation check
Jonatan Werpers <jonatan@werpers.com>
parents: 1977
diff changeset
352 @test_throws ArgumentError check_transfiniteinterpolation(s3)
486b3c6f919e Throw better error type for transfinite interpolation check
Jonatan Werpers <jonatan@werpers.com>
parents: 1977
diff changeset
353 @test_throws ArgumentError check_transfiniteinterpolation(s4)
1973
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
354 end
1977
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
355
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
356 @testset "Grids.jacobian" begin
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
357 cs = polygon_edges([0,0],[1,0],[1,1],[0,1])
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
358 ti = TransfiniteInterpolationSurface(cs...)
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
359
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
360 @test Grids.jacobian(ti, [0,0]) isa SMatrix
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
361
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
362 @test Grids.jacobian(ti, [0,0]) == [1 0; 0 1]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
363 @test Grids.jacobian(ti, [1,0]) == [1 0; 0 1]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
364 @test Grids.jacobian(ti, [1,1]) == [1 0; 0 1]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
365 @test Grids.jacobian(ti, [0,1]) == [1 0; 0 1]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
366
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
367 @test Grids.jacobian(ti, [1/2, 0]) == [1 0; 0 1]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
368 @test Grids.jacobian(ti, [1/2, 1]) == [1 0; 0 1]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
369 @test Grids.jacobian(ti, [0,1/2]) == [1 0; 0 1]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
370 @test Grids.jacobian(ti, [1,1/2]) == [1 0; 0 1]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
371
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
372
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
373 a, b, c, d = [1,0],[2,1/4],[2.5,1],[-1/3,1]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
374 cs = polygon_edges(a,b,c,d)
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
375 ti = TransfiniteInterpolationSurface(cs...)
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
376
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
377 @test Grids.jacobian(ti, [0,0]) == [b-a d-a]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
378 @test Grids.jacobian(ti, [1,0]) == [b-a c-b]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
379 @test Grids.jacobian(ti, [1,1]) == [c-d c-b]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
380 @test Grids.jacobian(ti, [0,1]) == [c-d d-a]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
381
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
382
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
383 mid(x,y) = (x+y)/2
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
384 @test Grids.jacobian(ti, [1/2, 0]) ≈ [b-a mid(c,d)-mid(a,b)]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
385 @test Grids.jacobian(ti, [1/2, 1]) ≈ [c-d mid(c,d)-mid(a,b)]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
386 @test Grids.jacobian(ti, [0, 1/2]) ≈ [mid(b,c)-mid(a,d) d-a]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
387 @test Grids.jacobian(ti, [1, 1/2]) ≈ [mid(b,c)-mid(a,d) c-b]
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
388 end
1915
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
389 end