Mercurial > repos > public > sbplib_julia
comparison src/Grids/mapped_grid.jl @ 1560:5d32ecb98db8 feature/grids/curvilinear
Add functions for jacobian_determiant, geometric_tensor, geometric_tensor_inverse
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Thu, 25 Apr 2024 09:03:04 +0200 |
| parents | 69790e9d1652 |
| children | 64baaf29ae4e 063a2bfb03da |
comparison
equal
deleted
inserted
replaced
| 1557:9113f437431d | 1560:5d32ecb98db8 |
|---|---|
| 58 lg, | 58 lg, |
| 59 map(x,lg), | 59 map(x,lg), |
| 60 map(J,lg), | 60 map(J,lg), |
| 61 ) | 61 ) |
| 62 end | 62 end |
| 63 | |
| 64 function jacobian_determinant(g::MappedGrid) | |
| 65 return map(jacobian(g)) do ∂x∂ξ | |
| 66 det(∂x∂ξ) | |
| 67 end | |
| 68 end | |
| 69 | |
| 70 function geometric_tensor(g::MappedGrid) | |
| 71 return map(jacobian(g)) do ∂x∂ξ | |
| 72 ∂x∂ξ'*∂x∂ξ | |
| 73 end | |
| 74 end | |
| 75 | |
| 76 function geometric_tensor_inverse(g::MappedGrid) | |
| 77 return map(jacobian(g)) do ∂x∂ξ | |
| 78 inv(∂x∂ξ'*∂x∂ξ) | |
| 79 end | |
| 80 end | |
| 81 |