I am interested in contributing to this. I am a contributor in the PINT pulsar timing package Scott mentioned. But I am new to contributing to numpy, and I am not sure where or how to begin. Are there beginner-level resources that I can look at?
Centre for Gravitation, Cosmology, and Astrophysics
University of Wisconsin Milwaukee
I think the first thing to do is have a discussion in the astropy community. Dependence on extended precision is fragile at best, it isn't available on all platforms and probably doesn't offer enough precision to be future proof. Then gather some performance data.
- What computations are needed?
- Can the computational algorithms be improved?
- Would parallel computation help?
- How much faster is extended precision?
- Where are the performance bottlenecks in the double-double time type?
- Is there hope of improving the double-double performance.
My own hope would be that a combination parallel computation and lower level code for double-double would be a viable solution.