<div dir="ltr"><div><div><div>My understanding of this is that the dtype will only hold the unit metadata. So that means units would propogate through calculations automatically, but the dtype wouldn't be able to manipulate the array data (in an in-place unit conversion for example).<br><br></div>In this world, astropy quantities and yt's YTArray would become containers around an ndarray that would make use of the dtype metadata but also implement all of the unit semantics that they already implement. Since they would become container classes and would no longer be ndarray subclasses, that avoids most of the pitfalls one encounters these days.<br><br></div>Please correct me if I'm wrong, Nathaniel.<br><br></div>-Nathan<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 26, 2017 at 5:14 PM, Marten van Kerkwijk <span dir="ltr"><<a href="mailto:m.h.vankerkwijk@gmail.com" target="_blank">m.h.vankerkwijk@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Nathaniel,<br>
<br>
Thanks for the link. The plans sounds great! You'll not be surprised<br>
to hear I'm particularly interested in the units aspect (and, no, I<br>
don't mind at all if we can stop subclassing ndarray...). Is the idea<br>
that there will be a general way for allow a dtype to define how to<br>
convert an array to one with another dtype? (Just as one now<br>
implicitly is able to convert between, say, int and float.) And, if<br>
so, is the idea that one of those conversion possibilities might<br>
involve checking units? Or were you thinking of implementing units<br>
more directly? The former would seem most sensible, if only so you can<br>
initially focus on other things than deciding how to support, say, esu<br>
vs emu units, or whether or not to treat radians as equal to<br>
dimensionless (which they formally are, but it is not always handy to<br>
do so).<br>
<br>
Anyway, do keep us posted! All the best,<br>
<br>
Marten<br>
<div class="HOEnZb"><div class="h5"><br>
On Thu, Oct 26, 2017 at 3:40 PM, Nathaniel Smith <<a href="mailto:njs@pobox.com">njs@pobox.com</a>> wrote:<br>
> On Wed, Oct 18, 2017 at 10:24 PM, Nathaniel Smith <<a href="mailto:njs@pobox.com">njs@pobox.com</a>> wrote:<br>
>> I'll also be giving a lunch talk at BIDS tomorrow to let folks locally<br>
>> know about what's going on, which I think will be recorded – I'll send<br>
>> around a link after in case others are interested.<br>
><br>
> Here's that link: <a href="https://www.youtube.com/watch?v=fowHwlpGb34" rel="noreferrer" target="_blank">https://www.youtube.com/watch?<wbr>v=fowHwlpGb34</a><br>
><br>
> -n<br>
><br>
> --<br>
> Nathaniel J. Smith -- <a href="https://vorpus.org" rel="noreferrer" target="_blank">https://vorpus.org</a><br>
> ______________________________<wbr>_________________<br>
> NumPy-Discussion mailing list<br>
> <a href="mailto:NumPy-Discussion@python.org">NumPy-Discussion@python.org</a><br>
> <a href="https://mail.python.org/mailman/listinfo/numpy-discussion" rel="noreferrer" target="_blank">https://mail.python.org/<wbr>mailman/listinfo/numpy-<wbr>discussion</a><br>
______________________________<wbr>_________________<br>
NumPy-Discussion mailing list<br>
<a href="mailto:NumPy-Discussion@python.org">NumPy-Discussion@python.org</a><br>
<a href="https://mail.python.org/mailman/listinfo/numpy-discussion" rel="noreferrer" target="_blank">https://mail.python.org/<wbr>mailman/listinfo/numpy-<wbr>discussion</a><br>
</div></div></blockquote></div><br></div>