Mercurial > repos > public > sbplib_julia
changeset 2000:889c18ad56bf feature/grids/parameter_spaces/in
Clean up of Base.in(x, s::Simplex)
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 29 Apr 2025 08:48:58 +0200 |
| parents | e6dbd4bec6cc |
| children | f45d32022e4f 607829f138a2 |
| files | src/Grids/parameter_space.jl |
| diffstat | 1 files changed, 4 insertions(+), 8 deletions(-) [+] |
line wrap: on
line diff
--- a/src/Grids/parameter_space.jl Fri Apr 25 08:16:06 2025 +0200 +++ b/src/Grids/parameter_space.jl Tue Apr 29 08:48:58 2025 +0200 @@ -170,16 +170,12 @@ v - v₁ end - A = hcat(V[2:end]...) # matrix with edge vectors as columns - b = x - v₁ + A = hcat(V[2:end]...) # Matrix with edge vectors as columns + λ = A \ (x - v₁) - # Solve Aλ = b - λ = A \ b + λ_full = (1 - sum(λ), λ...) # Full barycentric coordinates - # Compute full barycentric coordinates: first is 1 - sum(λ), then λ - λ_full = (1 - sum(λ), λ...) # Tuple of length NV - - all(λᵢ -> zero(λᵢ) ≤ λᵢ ≤ one(λᵢ), λ_full) + return all(λᵢ -> zero(λᵢ) ≤ λᵢ ≤ one(λᵢ), λ_full) end """