Milan Klöwer

Sorry for late answer, but I'll discuss a bit the information-related questions: > Would you recommend converting 32- to 64 bit before doing the compression? That should not make a...

Do you know whether they rounded in decimal or in binary? Because while rounding in one base also quantizes in the other it doesn't round (as in make the trailing...

In fact, we already discussed that here https://github.com/nco/nco/issues/250

Just checked what Charlie means by NSD vs DSD, but the former is a relative error vs the latter is an absolute error. Great that you used nco's `ppc` option,...

At the moment we are doing https://github.com/milankl/BitInformation.jl/blob/5f3ebbd135e427c68048988fdc24f6f2d5cb71e9/src/mutual_information.jl#L103-L105 So we count the bits but if the counter is still zero afterwards an `AssertionError` is triggered. ```julia julia> using BitInformation julia> A...

> You try to calculate bitinformation along a (too) short dimension of {length(data[dim])} elements. Maybe even continuing ... bitinformation assumes that adjacents analysis elements are also physically adjacent... Feel free...

With #68 addition is almost 2x faster than it was in v0.4 ```julia julia> using SoftPosit julia> A,B = Posit16.(rand(1000,1000)),Posit16.(rand(1000,1000)); julia> using BenchmarkTools julia> @btime +($A,$B); 15.684 ms (2 allocations:...

Yeah, I mean one could define `+(::PositN,::PositN)` (and `*`,`/`) directly and without the conversion which would likely give another big speedup, but at the moment it's only important that it's...

I believe maxintfloats are ```julia Base.maxintfloat(::Type{Posit8}) = 8 # = 2^3 Base.maxintfloat(::Type{Posit16}) = 512 # = 2^9 Base.maxintfloat(::Type{Posit16_2}) = 1024 # = 2^10 Base.maxintfloat(::Type{Posit32}) = 8388608 # = 2^23 ```

`Base.maxintfloat(::AbstractPosit)` added with #68