Julia (normally 1-based), Polynomials package:
julia> using Polynomials
julia> p = Polynomial([1, 0, 2])
Polynomial(1 + 2*x^2)
julia> p[0]
1
As long as it is domain specific and wrapped in custom types I don't really see an issue.
Re:Again on 0-based vs. 1-based indexing
by petschge on Wednesday January 20, 2021 @11:43AM
(#25848563)
And of course there is C++ libraries to do the same. For example the blitz++ library. Negative indices to the left of element[0] are great for ghost zones in codes that do domain decomposition over different MPI ranks.
https://docs.julialang.org/en/v1/devdocs/offset-arrays/