Mercurial > repos > public > sbplib_julia

annotate test/Grids/geometry_test.jl @ 1982:486b3c6f919e feature/grids/geometry_functions

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Throw better error type for transfinite interpolation check
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 26 Feb 2025 22:50:26 +0100
parents 270675bb97be
children cf0fa2967361
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
1973
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
2 using Diffinitive.Grids: Line, LineSegment, linesegments, polygon_edges, Circle, TransfiniteInterpolationSurface, check_transfiniteinterpolation
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
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
165 @testset "TransfiniteInterpolationSurface" begin
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
166 @testset "Constructors" begin
1973
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
167 @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
168
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
169 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
170 @test TransfiniteInterpolationSurface(cs...) isa TransfiniteInterpolationSurface
1965
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
171 end
eedb5d7b90b4 Add constructor testset
Jonatan Werpers <jonatan@werpers.com>
parents: 1915
diff changeset
172
1973
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
173 @testset "Evaluation" begin
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
174 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
175 ti = TransfiniteInterpolationSurface(cs...)
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
176
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
177 @test ti(0,0) == [0,0]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
178 @test ti([0,0]) == [0,0]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
179 @test ti(1,0) == [1,0]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
180 @test ti([1,0]) == [1,0]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
181 @test ti(1,1) == [1,1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
182 @test ti([1,1]) == [1,1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
183 @test ti(0,1) == [0,1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
184 @test ti([0,1]) == [0,1]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
185
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
186 @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
187 @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
188 @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
189 @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
190
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
191
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
192 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
193 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
194 ti = TransfiniteInterpolationSurface(cs...)
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
195
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
196 @test ti(0,0) == a
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
197 @test ti(1,0) == b
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
198 @test ti(1,1) == c
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
199 @test ti(0,1) == d
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
200
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
201 @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
202 @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
203 @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
204 @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
205
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
206 # TODO: Some test with curved edges?
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
207 end
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
208
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
209 @testset "check_transfiniteinterpolation" begin
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
210 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
211 ti = TransfiniteInterpolationSurface(cs...)
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
212
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
213 @test check_transfiniteinterpolation(ti) == nothing
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
214 @test check_transfiniteinterpolation(Bool, ti) == true
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
215
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
216 bad_sides = [
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
217 LineSegment([0,0],[1/2,0]),
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
218 LineSegment([1,0],[1,1/2]),
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
219 LineSegment([1,1],[0,1/2]),
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
220 LineSegment([0,1],[0,1/2]),
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
221 ]
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
222
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
223 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
224 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
225 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
226 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
227
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
228 @test check_transfiniteinterpolation(Bool, s1) == false
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
229 @test check_transfiniteinterpolation(Bool, s2) == false
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
230 @test check_transfiniteinterpolation(Bool, s3) == false
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
231 @test check_transfiniteinterpolation(Bool, s4) == false
8e9575b518a1 Add tests for TransfiniteInterpolationSurface and implemnet check_tranfiniteinterpolation
Jonatan Werpers <jonatan@werpers.com>
parents: 1969
diff changeset
232
1982
486b3c6f919e Throw better error type for transfinite interpolation check
Jonatan Werpers <jonatan@werpers.com>
parents: 1977
diff changeset
233 @test_throws ArgumentError check_transfiniteinterpolation(s1)
486b3c6f919e Throw better error type for transfinite interpolation check
Jonatan Werpers <jonatan@werpers.com>
parents: 1977
diff changeset
234 @test_throws ArgumentError check_transfiniteinterpolation(s2)
486b3c6f919e Throw better error type for transfinite interpolation check
Jonatan Werpers <jonatan@werpers.com>
parents: 1977
diff changeset
235 @test_throws ArgumentError check_transfiniteinterpolation(s3)
486b3c6f919e Throw better error type for transfinite interpolation check
Jonatan Werpers <jonatan@werpers.com>
parents: 1977
diff changeset
236 @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
237 end
1977
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
238
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
239 @testset "Grids.jacobian" begin
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
240 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
241 ti = TransfiniteInterpolationSurface(cs...)
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
242
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
243 @test Grids.jacobian(ti, [0,0]) isa SMatrix
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
244
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
245 @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
246 @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
247 @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
248 @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
249
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
250 @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
251 @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
252 @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
253 @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
254
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
255
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
256 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
257 cs = polygon_edges(a,b,c,d)
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
258 ti = TransfiniteInterpolationSurface(cs...)
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
259
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
260 @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
261 @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
262 @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
263 @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
264
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
265
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
266 mid(x,y) = (x+y)/2
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
267 @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
268 @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
269 @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
270 @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
271
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
272 # TODO: Some test with curved edges?
270675bb97be Implement Grids.jacobian for TransfiniteInterpolationSurface
Jonatan Werpers <jonatan@werpers.com>
parents: 1976
diff changeset
273 end
1915
c003685d9926 Add geometry functions removed from feature/grids/manifolds
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
274 end