       newArray[j].push array[i][j]
   newArray
 
- identity = [[1,0,0],[0,1,0],[0,0,1]]
+identity = [[1,0,0],[0,1,0],[0,0,1]]
 
 norm = (d)->
   # L2 norm (hypotenuse)

    e.push norm(f)

    d = (sdot(e,i) for i in axes)
    r = norm(d)

    if opts.traditionalLayout
      [y,z,x] = d

  s = axesCovariance.map Math.sqrt
  axes = transpose(axes)

  c1 = 1
  if upperHemisphere
    c1 *= -1
