title: Is there a way to make a complex number in NumPy with the accuracy of a Decimal type?
body: <p>I'm working in Python with NumPy arrays of complex numbers that extend well past the normal floating point limits of NumPy’s default <em>Complex</em> type (numbers greater than 10^500). I wanted to know if there was some way I could extend NumPy so that it will be able to handle complex numbers with this sort of magnitude. For example, is there a way to make a NumPy complex type that uses functionality from the Decimal module?</p>
  <p>I know there are resources available such as <a href="https://code.google.com/p/mpmath/" rel="nofollow noreferrer">mpmath</a> that could probably do what I need. However, it is a requirement of my project that I use NumPy.</p>
  <p>For anyone that's interested in why I need these enormous numbers, it's because I'm working on a <a href="https://en.wikipedia.org/wiki/Numerical_relativity" rel="nofollow noreferrer">numerical relativity</a> simulation of the early universe.</p>