John Murray

It was my set up @joehoeller, now fixed.

Just a note that CuPy now supports float16 in user defined kernels. ``` >>> import cupy as cp >>> c = cp.random.random(1000000).astype(cp.float16) @cp.fuse() ... def squared_diff(x, y): ... return (x...

> I think ideally you could pass a `DataArray` to something that takes GPU arrays (like Numba kernels). If that doesn't make sense then perhaps passing the `DataArray.data` would be...

> Do you have to go through `__array__` (see #6845) or would accessing the underlying array using `DataArray.data` work for you? > > We could also add some properties under...

I struggled with the C++ element of the coding as I'm very rusty in it. I have it coded in Python.

> This is still worthwhile. Reopening. Happy to help @harrism @thomcom but my C++ isn't up to the standard. I have published the Python version above. I have also written...

High accuracy coordinate projection uses high precision formulae, iterative (in case of E/N to Lat/Lon), to provide surveying & civil engineering quality coordinate conversion. It is used in conjunction with...

I have implemented this already in Spatia/Rapids. It is straightforward.

I would suggest starting here. This is UK Ordnance Survey's version for UK and Ireland. https://www.ordnancesurvey.co.uk/business-government/tools-support/os-net/format The developer pack contains some very good documentation and test data with results for...

The techniques in the OS guide equally apply to NADCON (USA & Canada), Germany and other countries with Transformation requirements. Transformation grids have closely situated ground reference points and are...