title: numpy: inverting an upper triangular matrix
body: <p>In <code>numpy</code>/<code>scipy</code>, what's the canonical way to compute the inverse of an upper triangular matrix?</p>
  
  <p>The matrix is stored as 2D <code>numpy</code> array with zero sub-diagonal elements, and the result should also be stored as a 2D array.</p>
  
  <p><strong>edit</strong> The best I've found so far is <code>scipy.linalg.solve_triangular(A, np.identity(n))</code>. Is that it?</p>