// File: matthewarcus[P]mmps/transform.cpp
// Function: Perspective::ProjectInv(TransformParams const&, double, double, double&, double&) const
bool Perspective::ProjectInv(const TransformParams& params, 
			     double phi, double lambda, double& x, double& y) const
{
  (void(params));
  bool result = false;
  // Find where phi, lambda project to on the ellipsoid
  double x0,y0,z0;
  double nx,ny,nz;
  if (a == 1.0 && b == 1.0 && c == 1.0) {
    // Just a sphere
    x0 = cos(lambda) * cos(phi);
    y0 = sin(lambda) * cos(phi);
    z0 = sin(phi);
    minv.Apply(x0,y0,z0);
    nx = x0; ny = y0; nz = z0;
  } else {
    // The normal vector
    nx = cos(lambda);
    ny = sin(lambda);
    nz = tan(phi);

    // The actual vector
    double r = sqrt(1/(ia2 * sqr(nx) + ib2 * sqr(ny) + ic2 * sqr(nz)));
    x0 = nx * r * ia2;
    y0 = ny * r * ib2;
    z0 = nz * r * ic2;
    // And apply rotations
    minv.Apply(x0,y0,z0);
    minv.Apply(nx,ny,nz);
  }

  // Test for visibility - dot product of nx,ny,nz and the line
  // towards to viewpoint must be positive
  double p = (vx-x0) * nx + (vy-y0) * ny + (vz-z0) * nz;
  if (p >= 0) {
    // Project from (vx, vy, vz) through (x0,y0,z0) to (-1, y, z)
    double t = (1 + x0)/(x0 - vx);
    x = t * vy + (1-t) * y0; // Note change of axes
    y = t * vz + (1-t) * z0;
    x = (x - vy) * iscalefactor;
    y = (y - vz) * iscalefactor;
    result = true;
  }
  return result;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xc0, %rsp
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movsd %xmm0, -0x18(%rbp)
movsd %xmm1, -0x20(%rbp)
movq %rdx, -0x28(%rbp)
movq %rcx, -0x30(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0x88(%rbp)
movb $0x0, -0x31(%rbp)
movsd 0x100(%rax), %xmm0
movsd 0x92a4(%rip), %xmm1     # 0x15048
ucomisd %xmm1, %xmm0
jne 0xbeb0
jp 0xbeb0
movq -0x88(%rbp), %rax
movsd 0x108(%rax), %xmm0
movsd 0x927d(%rip), %xmm1     # 0x15048
ucomisd %xmm1, %xmm0
jne 0xbeb0
jp 0xbeb0
movq -0x88(%rbp), %rax
movsd 0x110(%rax), %xmm0
movsd 0x9256(%rip), %xmm1     # 0x15048
ucomisd %xmm1, %xmm0
jne 0xbeb0
jp 0xbeb0
movsd -0x20(%rbp), %xmm0
callq 0x5470
movsd %xmm0, -0x98(%rbp)
movsd -0x18(%rbp), %xmm0
callq 0x5470
movaps %xmm0, %xmm1
movsd -0x98(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x40(%rbp)
movsd -0x20(%rbp), %xmm0
callq 0x5080
movsd %xmm0, -0x90(%rbp)
movsd -0x18(%rbp), %xmm0
callq 0x5470
movaps %xmm0, %xmm1
movsd -0x90(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x48(%rbp)
movsd -0x18(%rbp), %xmm0
callq 0x5080
movq -0x88(%rbp), %rdi
movsd %xmm0, -0x50(%rbp)
addq $0x58, %rdi
leaq -0x40(%rbp), %rsi
leaq -0x48(%rbp), %rdx
leaq -0x50(%rbp), %rcx
callq 0x7da0
movsd -0x40(%rbp), %xmm0
movsd %xmm0, -0x58(%rbp)
movsd -0x48(%rbp), %xmm0
movsd %xmm0, -0x60(%rbp)
movsd -0x50(%rbp), %xmm0
movsd %xmm0, -0x68(%rbp)
jmp 0xc025
movsd -0x20(%rbp), %xmm0
callq 0x5470
movsd %xmm0, -0x58(%rbp)
movsd -0x20(%rbp), %xmm0
callq 0x5080
movsd %xmm0, -0x60(%rbp)
movsd -0x18(%rbp), %xmm0
callq 0x5150
movq -0x88(%rbp), %rax
movsd %xmm0, -0x68(%rbp)
movsd 0x130(%rax), %xmm0
movsd %xmm0, -0xb0(%rbp)
movsd -0x58(%rbp), %xmm0
callq 0xddd0
movq -0x88(%rbp), %rax
movsd %xmm0, -0xb8(%rbp)
movsd 0x138(%rax), %xmm0
movsd %xmm0, -0xc0(%rbp)
movsd -0x60(%rbp), %xmm0
callq 0xddd0
movsd -0xc0(%rbp), %xmm1
movsd -0xb8(%rbp), %xmm2
movq -0x88(%rbp), %rax
movaps %xmm0, %xmm3
movsd -0xb0(%rbp), %xmm0
mulsd %xmm3, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0xa0(%rbp)
movsd 0x140(%rax), %xmm0
movsd %xmm0, -0xa8(%rbp)
movsd -0x68(%rbp), %xmm0
callq 0xddd0
movsd -0xa8(%rbp), %xmm1
movaps %xmm0, %xmm2
movsd -0xa0(%rbp), %xmm0
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movsd 0x90ae(%rip), %xmm0     # 0x15048
divsd %xmm1, %xmm0
callq 0x5200
movq -0x88(%rbp), %rdi
movsd %xmm0, -0x70(%rbp)
movsd -0x58(%rbp), %xmm0
mulsd -0x70(%rbp), %xmm0
mulsd 0x130(%rdi), %xmm0
movsd %xmm0, -0x40(%rbp)
movsd -0x60(%rbp), %xmm0
mulsd -0x70(%rbp), %xmm0
mulsd 0x138(%rdi), %xmm0
movsd %xmm0, -0x48(%rbp)
movsd -0x68(%rbp), %xmm0
mulsd -0x70(%rbp), %xmm0
mulsd 0x140(%rdi), %xmm0
movsd %xmm0, -0x50(%rbp)
addq $0x58, %rdi
leaq -0x40(%rbp), %rsi
leaq -0x48(%rbp), %rdx
leaq -0x50(%rbp), %rcx
callq 0x7da0
movq -0x88(%rbp), %rdi
addq $0x58, %rdi
leaq -0x58(%rbp), %rsi
leaq -0x60(%rbp), %rdx
leaq -0x68(%rbp), %rcx
callq 0x7da0
movq -0x88(%rbp), %rax
movsd 0xd0(%rax), %xmm1
subsd -0x40(%rbp), %xmm1
movsd -0x58(%rbp), %xmm2
movsd 0xd8(%rax), %xmm0
subsd -0x48(%rbp), %xmm0
mulsd -0x60(%rbp), %xmm0
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movsd 0xe0(%rax), %xmm0
subsd -0x50(%rbp), %xmm0
movsd -0x68(%rbp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x78(%rbp)
movsd -0x78(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0xc155
movq -0x88(%rbp), %rax
movsd 0x8fb0(%rip), %xmm0     # 0x15048
addsd -0x40(%rbp), %xmm0
movsd -0x40(%rbp), %xmm1
subsd 0xd0(%rax), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0x80(%rbp)
movsd -0x80(%rbp), %xmm0
movsd 0xd8(%rax), %xmm2
movsd 0x8f80(%rip), %xmm1     # 0x15048
subsd -0x80(%rbp), %xmm1
mulsd -0x48(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq -0x28(%rbp), %rcx
movsd %xmm0, (%rcx)
movsd -0x80(%rbp), %xmm0
movsd 0xe0(%rax), %xmm2
movsd 0x8f51(%rip), %xmm1     # 0x15048
subsd -0x80(%rbp), %xmm1
mulsd -0x50(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq -0x30(%rbp), %rcx
movsd %xmm0, (%rcx)
movq -0x28(%rbp), %rcx
movsd (%rcx), %xmm0
subsd 0xd8(%rax), %xmm0
mulsd 0xc8(%rax), %xmm0
movq -0x28(%rbp), %rcx
movsd %xmm0, (%rcx)
movq -0x30(%rbp), %rcx
movsd (%rcx), %xmm0
subsd 0xe0(%rax), %xmm0
mulsd 0xc8(%rax), %xmm0
movq -0x30(%rbp), %rax
movsd %xmm0, (%rax)
movb $0x1, -0x31(%rbp)
movb -0x31(%rbp), %al
andb $0x1, %al
addq $0xc0, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)