Mercurial > repos > public > sbplib_julia
annotate grid_refactor.md @ 1276:75a65db29be1 refactor/grids
Minor clean up
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Sun, 26 Feb 2023 17:09:35 +0100 |
| parents | dcd8654ca33b |
| children | 79647b60a73b |
| rev | line source |
|---|---|
| 1222 | 1 # Grids refactor |
| 2 | |
| 3 ## Mål | |
| 4 * Embedded grids | |
| 5 * Vektorvärda grid funktioner, speciellt koordinat-funktionen | |
| 6 * Olika griddar i olika riktningar? | |
| 7 * Tex periodiskt i en riktning och intervall i en. | |
| 8 * Tex ostrukturerat i ett par och strukturerat i en. | |
| 9 * Enkelt att evaluera på 0d-grid | |
| 10 | |
| 11 ## Frågor | |
| 12 | |
| 13 | |
| 14 ### Ska `Grid` vara en AbstractArray? | |
| 15 Efter som alla nät ska agera som en gridfunktion för koordinaterna måste man | |
| 16 svara på frågan hur vi hanterar generellla gridfunktioner samtidigt. | |
| 17 | |
| 18 Några saker att förhålla sig till: | |
| 19 - Multiblock nät? | |
| 20 - Unstructured? | |
| 21 - Triangular structured grids? | |
| 22 - Non-simply connected? | |
| 23 - CG/DG-nät | |
| 24 | |
| 25 | |
| 26 Om alla nät är någon slags AbstractArray så kan tillexempel ett multiblock nät vara en AbstractArray{<:Grid, 1} och motsvarande gridfunktioner AbstractArray{<:AbstractArray,1}. | |
| 27 Där man alltså först indexerar för vilket när man vill åt och sedan indexerar nätet. Tex `mg[2][32,12]`. | |
| 28 | |
| 29 Ett ostrukturerat nät med till exempel trianglar skulle vi kunna se på samma sätt som ett multiblocknät. Antagligen har de två typerna av nät olika underliggande datastruktur anpassade efter ändamålet. | |
| 30 | |
| 31 Hur fungerar tankarna ovan om man har tex tensorprodukten av ett ostrukturerat nät och en ekvidistant nät? | |
| 32 ```julia | |
| 33 m = Mesh2DTriangle(10) | |
| 34 e = EqudistantGrid(range(1:10) | |
| 35 | |
| 36 e[4] # fourth point | |
| 37 | |
| 38 m[3][5] # Fifth node in third triangle | |
| 39 m[3,5] # Fifth node in third triangle # Funkar bara om alla nät är samma, (stämmer inte i mb-fallet) | |
| 40 | |
| 41 g = TensorGrid(m, e) | |
| 42 | |
| 43 g[3,4][5] # ?? | |
| 44 g[3,4] # ?? | |
| 45 | |
| 46 g[3,5,4] # ?? | |
| 47 | |
| 48 | |
| 49 | |
| 50 ``` | |
| 51 | |
| 52 Alla grids kanske inte är AbstractArray? Måste de vara rektangulära? Det blir svårt för strukturerade trianglar och NSC-griddar. Men de ska i allafall vara indexerbara? | |
| 53 | |
| 54 Man skulle kunna utesluta MultiblockGrid i tensorgrids | |
| 55 | |
| 56 CG-nät och DG-nät blir olika. | |
| 57 På CG-nät kanske man både vill indexera noder och trianglar beroende på vad man håller på med? | |
| 58 | |
| 1235 | 59 |
| 60 Om griddarna inte ska vara AbstractArray finns det många andra ställen som blir konstiga om de är AbstractArray. TensorApplication?! LazyArrays?! Är alla saker vi jobbar med egentligen mer generella object? Finns det något sätt att uttrycka koden så att man kan välja? | |
| 61 | |
| 1254 | 62 |
| 63 Det vi är ute efter är kanske att griddarna uppfyller Iteration och Indexing interfacen. | |
| 64 | |
| 1222 | 65 #### Försök till slutsater |
| 66 * Multiblock-nät indexeras i två nivåer tex `g[3][3,4]` | |
| 67 * Vi struntar i att implementera multiblock-nät som en del av ett tensorgrid till att börja med. | |
| 68 * En grid kan inte alltid vara en AbstractArray eftersom till exempel ett NCS eller strukturerad triangel inte har rätt form. | |
| 69 * Om vi har nod-indexerade ostrukturerade nät borde de fungera med TensorGrid. | |
| 1254 | 70 * Griddar ska uppfylla Indexing och Iteration interfacen |
| 1222 | 71 |
| 72 ### Kan vi introducera 1d griddar och tensorgriddar? | |
| 73 * Vanligt intervallgrid | |
| 74 * Periodiskt grid | |
| 75 * 0-dimensionellt grid | |
| 76 | |
| 77 Det skulle kunna lösa frågan med embedded-griddar | |
| 78 och frågan om mixade grids. | |
| 79 | |
| 80 En svårighet kan vara att implemetera indexeringen för tensorgridden. Är det | |
| 81 lätt att göra den transparent för kompilatorn? | |
| 82 | |
| 83 Läs i notes.md om: Grids embedded in higher dimensions | |
| 84 | |
| 85 Periodiska gridfunktioner? Eller ska vi bara skita i det och låta användaren | |
| 86 hantera det själv? Blir det krångligt i högre dimensioner? | |
| 87 | |
| 88 | |
| 89 ### Hur ska vi hantera gridfunktioner med flera komponenter? | |
| 90 Tex koordinat-funktionen på nätet? | |
| 91 | |
| 92 Funkar det att ha StaticArray som element type? | |
| 93 Kan man köra rakt på med en operator då? | |
| 94 | |
| 95 Vi skulle också kunna använda ArraysOfArrays. Skillnaden blir att den elementtypen är Array{T,N}. Vilket betyder att den är by reference? | |
| 96 Frågan är om kompilatorn kommer att bry sig... Jag skulle luta mot StaticArrays for safety | |
| 97 | |
| 98 Läs i notes.md om: Vector valued grid functions | |
| 99 | |
|
1265
9c9ea2900250
grid_refactor.md: Add question about target dim function
Jonatan Werpers <jonatan@werpers.com>
parents:
1254
diff
changeset
|
100 ## Should Grid have function for the target manifold dimension? |
|
9c9ea2900250
grid_refactor.md: Add question about target dim function
Jonatan Werpers <jonatan@werpers.com>
parents:
1254
diff
changeset
|
101 Where would it be used? |
|
9c9ea2900250
grid_refactor.md: Add question about target dim function
Jonatan Werpers <jonatan@werpers.com>
parents:
1254
diff
changeset
|
102 In the constructor for TensorGrid |
|
1269
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
103 In eval on if we want to allow multiargument functions |
|
1265
9c9ea2900250
grid_refactor.md: Add question about target dim function
Jonatan Werpers <jonatan@werpers.com>
parents:
1254
diff
changeset
|
104 Elsewhere? |
|
9c9ea2900250
grid_refactor.md: Add question about target dim function
Jonatan Werpers <jonatan@werpers.com>
parents:
1254
diff
changeset
|
105 |
|
9c9ea2900250
grid_refactor.md: Add question about target dim function
Jonatan Werpers <jonatan@werpers.com>
parents:
1254
diff
changeset
|
106 An alternative is to analyze T in Grid{T,D} to find the answer. (See combined_coordinate_vector_type in tensor_grid.jl) |
|
9c9ea2900250
grid_refactor.md: Add question about target dim function
Jonatan Werpers <jonatan@werpers.com>
parents:
1254
diff
changeset
|
107 |
|
1269
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
108 ## Lazy version of map for our needs? |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
109 Could be used to |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
110 * evaulate functions on grids |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
111 * pick out components of grid functions |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
112 * More? |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
113 |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
114 Maybe this: |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
115 ```julia |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
116 struct LazyMappedArray <: LazyArray |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
117 f::F |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
118 v::AT |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
119 end |
|
20f42cf0800c
Add test for Grid and make them pass. Start implementing eval_on
Jonatan Werpers <jonatan@werpers.com>
parents:
1265
diff
changeset
|
120 ``` |
|
1265
9c9ea2900250
grid_refactor.md: Add question about target dim function
Jonatan Werpers <jonatan@werpers.com>
parents:
1254
diff
changeset
|
121 |
|
1275
dcd8654ca33b
Remove getcomponent stubs since we might have a cleaner way to solve it with map or lazymap
Jonatan Werpers <jonatan@werpers.com>
parents:
1269
diff
changeset
|
122 Could allow us to remove eval_on. |
|
dcd8654ca33b
Remove getcomponent stubs since we might have a cleaner way to solve it with map or lazymap
Jonatan Werpers <jonatan@werpers.com>
parents:
1269
diff
changeset
|
123 |
|
dcd8654ca33b
Remove getcomponent stubs since we might have a cleaner way to solve it with map or lazymap
Jonatan Werpers <jonatan@werpers.com>
parents:
1269
diff
changeset
|
124 ## Do we need functions like `getcomponent`? |
|
dcd8654ca33b
Remove getcomponent stubs since we might have a cleaner way to solve it with map or lazymap
Jonatan Werpers <jonatan@werpers.com>
parents:
1269
diff
changeset
|
125 Perhaps this can be more cleanly solved using map or a lazy version of map? |
|
dcd8654ca33b
Remove getcomponent stubs since we might have a cleaner way to solve it with map or lazymap
Jonatan Werpers <jonatan@werpers.com>
parents:
1269
diff
changeset
|
126 That approach would be more flexible and more general requiring few specialized functions. |
|
dcd8654ca33b
Remove getcomponent stubs since we might have a cleaner way to solve it with map or lazymap
Jonatan Werpers <jonatan@werpers.com>
parents:
1269
diff
changeset
|
127 |
|
dcd8654ca33b
Remove getcomponent stubs since we might have a cleaner way to solve it with map or lazymap
Jonatan Werpers <jonatan@werpers.com>
parents:
1269
diff
changeset
|
128 (see "Lazy version of map for our needs?" above) |
|
dcd8654ca33b
Remove getcomponent stubs since we might have a cleaner way to solve it with map or lazymap
Jonatan Werpers <jonatan@werpers.com>
parents:
1269
diff
changeset
|
129 |
| 1222 | 130 ## Notes from pluto notebook |
| 131 - Är det dåligt att använda ndims om antalet index inte matchar? | |
| 132 - Tex ostrukturerat grid | |
| 133 - Baserat på hur Array fungerar och ndims docs borde den nog returnera | |
| 134 antalet index och inte griddens dimension | |
| 135 - Å andra sidan kanske man inte behöver veta viken dimension det är? Fast det känns konstigt... | |
| 136 - Should we keep the `points` function? | |
| 137 - Maybe it's better to support `collect` on the grid? | |
| 138 - How to handle boundary identifiers and boundary grids on `TensorGrid`? | |
| 139 | |
| 140 |