Mercurial > repos > public > sbplib_julia
comparison src/Grids/parameter_space.jl @ 2001:f45d32022e4f feature/grids/manifolds
Merge feature/grids/parameter_spaces/in
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 29 Apr 2025 08:59:29 +0200 |
| parents | 889c18ad56bf |
| children |
comparison
equal
deleted
inserted
replaced
| 1999:a1b2453c02c9 | 2001:f45d32022e4f |
|---|---|
| 168 v₁ = s.verticies[1] | 168 v₁ = s.verticies[1] |
| 169 V = map(s.verticies) do v | 169 V = map(s.verticies) do v |
| 170 v - v₁ | 170 v - v₁ |
| 171 end | 171 end |
| 172 | 172 |
| 173 A = hcat(V[2:end]...) # matrix with edge vectors as columns | 173 A = hcat(V[2:end]...) # Matrix with edge vectors as columns |
| 174 b = x - v₁ | 174 λ = A \ (x - v₁) |
| 175 | 175 |
| 176 # Solve Aλ = b | 176 λ_full = (1 - sum(λ), λ...) # Full barycentric coordinates |
| 177 λ = A \ b | 177 |
| 178 | 178 return all(λᵢ -> zero(λᵢ) ≤ λᵢ ≤ one(λᵢ), λ_full) |
| 179 # Compute full barycentric coordinates: first is 1 - sum(λ), then λ | |
| 180 λ_full = (1 - sum(λ), λ...) # Tuple of length NV | |
| 181 | |
| 182 all(λᵢ -> zero(λᵢ) ≤ λᵢ ≤ one(λᵢ), λ_full) | |
| 183 end | 179 end |
| 184 | 180 |
| 185 """ | 181 """ |
| 186 verticies(s::Simplex) | 182 verticies(s::Simplex) |
| 187 | 183 |