Some bugs have been worked out, but problems still remain. The correct eigenvalues are not found and the code crashes with an error (libc++abi.dylib: terminating with uncaught exception of type std::runtime_error: Unable to diagonalize matrix Abort trap: 6) about 60% of the time.