user: What would the C source look like for the assembly code provided here:

```asm
<GPS_Math_Bonne_EN_To_LatLon>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x1c0,%rsp
  movsd  %xmm0,-0x178(%rbp)
  movsd  %xmm1,-0x180(%rbp)
  mov    %rdi,-0x188(%rbp)
  mov    %rsi,-0x190(%rbp)
  movsd  %xmm2,-0x198(%rbp)
  movsd  %xmm3,-0x1a0(%rbp)
  movsd  %xmm4,-0x1a8(%rbp)
  movsd  %xmm5,-0x1b0(%rbp)
  movsd  %xmm6,-0x1b8(%rbp)
  movsd  %xmm7,-0x1c0(%rbp)
  mov    -0x198(%rbp),%rax
  movq   %rax,%xmm0
  call   <GPS_Math_Deg_To_Rad@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x198(%rbp)
  mov    -0x1a0(%rbp),%rax
  movq   %rax,%xmm0
  call   <GPS_Math_Deg_To_Rad@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1a0(%rbp)
  mov    -0x198(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x148(%rbp)
  movsd  D169(%rip),%xmm0
  movsd  %xmm0,-0x140(%rbp)
  movsd  D490(%rip),%xmm0
  movsd  %xmm0,-0x138(%rbp)
  movsd  -0x1a0(%rbp),%xmm0
  comisd D288(%rip),%xmm0
  jbe    L33
  movsd  -0x1a0(%rbp),%xmm0
  subsd  -0x140(%rbp),%xmm0
  movsd  %xmm0,-0x1a0(%rbp)
L33:
  movsd  -0x1b8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x130(%rbp)
  movsd  -0x1c0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x128(%rbp)
  movsd  -0x130(%rbp),%xmm0
  subsd  -0x128(%rbp),%xmm0
  divsd  -0x130(%rbp),%xmm0
  movsd  %xmm0,-0x120(%rbp)
  movsd  -0x120(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x118(%rbp)
  movsd  -0x120(%rbp),%xmm0
  mulsd  -0x118(%rbp),%xmm0
  movsd  %xmm0,-0x110(%rbp)
  movsd  -0x110(%rbp),%xmm1
  movsd  D22(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D255(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x108(%rbp)
  movsd  -0x118(%rbp),%xmm1
  movsd  D167(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x100(%rbp)
  movsd  -0x120(%rbp),%xmm0
  movsd  D226(%rip),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  D157(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  -0x100(%rbp),%xmm1
  movsd  D40(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0x110(%rbp),%xmm2
  movsd  D383(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  movsd  D44(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf8(%rbp)
  movsd  -0x120(%rbp),%xmm1
  movsd  D167(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D261(%rip),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  -0x100(%rbp),%xmm0
  movsd  D113(%rip),%xmm2
  divsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x108(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  movsd  -0x118(%rbp),%xmm1
  movsd  D295(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D44(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0x108(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  movsd  -0x110(%rbp),%xmm1
  movsd  D520(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D123(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  mov    -0x198(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,-0xd8(%rbp)
  movsd  -0x120(%rbp),%xmm0
  mulsd  -0x148(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x148(%rbp),%xmm1
  movsd  D157(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movsd  D456(%rip),%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movapd %xmm0,%xmm1
  movsd  -0xd8(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  movsd  -0xf8(%rbp),%xmm0
  mulsd  -0x198(%rbp),%xmm0
  movsd  %xmm0,-0xc8(%rbp)
  movsd  -0x198(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0xf0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc0(%rbp)
  movsd  -0x198(%rbp),%xmm1
  movsd  D226(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0xe8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xb8(%rbp)
  movsd  -0x198(%rbp),%xmm1
  movsd  D150(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0xe0(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  movsd  -0xc8(%rbp),%xmm0
  subsd  -0xc0(%rbp),%xmm0
  addsd  -0xb8(%rbp),%xmm0
  subsd  -0xb0(%rbp),%xmm0
  movsd  -0x1b8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  D157(%rip),%xmm0
  subsd  -0x120(%rbp),%xmm0
  movq   %xmm0,%rax
  movsd  D456(%rip),%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,-0xa0(%rbp)
  movsd  D157(%rip),%xmm0
  subsd  -0xa0(%rbp),%xmm0
  movsd  -0xa0(%rbp),%xmm2
  movsd  D157(%rip),%xmm1
  addsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0x98(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0x90(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0x88(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x98(%rbp),%xmm1
  movsd  D167(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D476(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0x88(%rbp),%xmm2
  movsd  D364(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  movsd  D113(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0x90(%rbp),%xmm1
  movsd  D558(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D27(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0x80(%rbp),%xmm2
  movsd  D21(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  movsd  D113(%rip),%xmm2
  divsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0x88(%rbp),%xmm1
  movsd  D310(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D178(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x80(%rbp),%xmm1
  movsd  D409(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D565(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  pxor   %xmm0,%xmm0
  ucomisd -0x148(%rbp),%xmm0
  jp     L167
  pxor   %xmm0,%xmm0
  ucomisd -0x148(%rbp),%xmm0
  jne    L167
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x168(%rbp)
  jmp    L11
L167:
  movsd  -0x1b8(%rbp),%xmm0
  mulsd  -0xd0(%rbp),%xmm0
  divsd  -0x148(%rbp),%xmm0
  movsd  %xmm0,-0x168(%rbp)
L11:
  movsd  -0x178(%rbp),%xmm0
  subsd  -0x1a8(%rbp),%xmm0
  movsd  %xmm0,-0x160(%rbp)
  movsd  -0x180(%rbp),%xmm0
  subsd  -0x1b0(%rbp),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x168(%rbp),%xmm0
  subsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,-0x158(%rbp)
  movsd  -0x160(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0x158(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movsd  D456(%rip),%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x150(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x198(%rbp),%xmm0
  jbe    L138
  movsd  -0x150(%rbp),%xmm0
  movq   D472(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x150(%rbp)
L138:
  movsd  -0x168(%rbp),%xmm0
  addsd  -0xa8(%rbp),%xmm0
  subsd  -0x150(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x1b8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xf8(%rbp),%xmm1
  movsd  -0x50(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x48(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x78(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x48(%rbp),%xmm1
  movsd  D226(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x70(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x48(%rbp),%xmm1
  movsd  D150(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x68(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x48(%rbp),%xmm1
  movsd  D261(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x60(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x48(%rbp),%xmm0
  addsd  -0x40(%rbp),%xmm0
  addsd  -0x38(%rbp),%xmm0
  addsd  -0x30(%rbp),%xmm0
  addsd  -0x28(%rbp),%xmm0
  mov    -0x188(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  D7(%rip),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x138(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x20(%rbp),%xmm1
  mov    -0x188(%rbp),%rax
  movsd  (%rax),%xmm0
  movq   D575(%rip),%xmm2
  andpd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  jbe    L8
  mov    -0x188(%rbp),%rax
  movsd  (%rax),%xmm0
  movq   D575(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  -0x138(%rbp),%xmm0
  addsd  -0x20(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jbe    L8
  mov    -0x190(%rbp),%rax
  movsd  -0x1a0(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L34
L8:
  mov    -0x188(%rbp),%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x188(%rbp),%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x10(%rbp)
  movsd  -0x120(%rbp),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D157(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movsd  D456(%rip),%xmm0
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movapd %xmm0,%xmm1
  movsd  -0x18(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x198(%rbp),%xmm0
  jbe    L49
  movsd  -0x160(%rbp),%xmm0
  movq   D472(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x160(%rbp)
  movsd  -0x158(%rbp),%xmm0
  movq   D472(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x158(%rbp)
L49:
  movsd  -0x158(%rbp),%xmm0
  mov    -0x160(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <atan2@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  mulsd  -0x150(%rbp),%xmm0
  movsd  -0x1b8(%rbp),%xmm1
  mulsd  -0x8(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  addsd  -0x1a0(%rbp),%xmm0
  mov    -0x190(%rbp),%rax
  movsd  %xmm0,(%rax)
L34:
  mov    -0x188(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd -0x138(%rbp),%xmm0
  jbe    L17
  mov    -0x188(%rbp),%rax
  movsd  -0x138(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L141
L17:
  mov    -0x188(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  -0x138(%rbp),%xmm0
  movq   D472(%rip),%xmm2
  xorpd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  jbe    L141
  movsd  -0x138(%rbp),%xmm0
  movq   D472(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  mov    -0x188(%rbp),%rax
  movsd  %xmm0,(%rax)
L141:
  mov    -0x190(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd D288(%rip),%xmm0
  jbe    L105
  mov    -0x190(%rbp),%rax
  movsd  (%rax),%xmm0
  subsd  -0x140(%rbp),%xmm0
  mov    -0x190(%rbp),%rax
  movsd  %xmm0,(%rax)
L105:
  mov    -0x190(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D475(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L143
  mov    -0x190(%rbp),%rax
  movsd  (%rax),%xmm0
  addsd  -0x140(%rbp),%xmm0
  mov    -0x190(%rbp),%rax
  movsd  %xmm0,(%rax)
L143:
  mov    -0x190(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd D288(%rip),%xmm0
  jbe    L0
  mov    -0x190(%rbp),%rax
  movsd  D288(%rip),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L12
L0:
  mov    -0x190(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D475(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L12
  mov    -0x190(%rbp),%rax
  movsd  D475(%rip),%xmm0
  movsd  %xmm0,(%rax)
L12:
  mov    -0x190(%rbp),%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movq   %xmm0,%rax
  mov    -0x190(%rbp),%rdx
  mov    %rax,(%rdx)
  mov    -0x188(%rbp),%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movq   %xmm0,%rax
  mov    -0x188(%rbp),%rdx
  mov    %rax,(%rdx)
  nop
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: tag: D169
type_infered: byte[8]
name: retrieve_data
arguments: tag: D490
type_infered: byte[8]
name: retrieve_data
arguments: tag: D288
type_infered: f64[2]
name: retrieve_data
arguments: tag: D22
type_infered: f64
name: retrieve_data
arguments: tag: D255
type_infered: f64
name: retrieve_data
arguments: tag: D167
type_infered: f64
name: retrieve_data
arguments: tag: D226
type_infered: f64
name: retrieve_data
arguments: tag: D157
type_infered: f64
name: retrieve_data
arguments: tag: D40
type_infered: f64
name: retrieve_data
arguments: tag: D383
type_infered: f64
name: retrieve_data
arguments: tag: D44
type_infered: f64
name: retrieve_data
arguments: tag: D261
type_infered: f64
name: retrieve_data
arguments: tag: D113
type_infered: f64
name: retrieve_data
arguments: tag: D295
type_infered: f64
name: retrieve_data
arguments: tag: D520
type_infered: f64
name: retrieve_data
arguments: tag: D123
type_infered: f64
name: retrieve_data
arguments: tag: D456
type_infered: f64
name: retrieve_data
arguments: tag: D150
type_infered: f64
name: retrieve_data
arguments: tag: D476
type_infered: f64
name: retrieve_data
arguments: tag: D364
type_infered: f64
name: retrieve_data
arguments: tag: D558
type_infered: f64
name: retrieve_data
arguments: tag: D27
type_infered: f64
name: retrieve_data
arguments: tag: D21
type_infered: f64
name: retrieve_data
arguments: tag: D310
type_infered: f64
name: retrieve_data
arguments: tag: D178
type_infered: f64
name: retrieve_data
arguments: tag: D409
type_infered: f64
name: retrieve_data
arguments: tag: D565
type_infered: f64
name: retrieve_data
arguments: tag: D472
type_infered: f32[2]
name: retrieve_data
arguments: tag: D7
type_infered: f64
name: retrieve_data
arguments: tag: D575
type_infered: qword
name: retrieve_data
arguments: tag: D475
type_infered: byte[8]

tool: {"D169": [0x1e,0x26,0x44,0x54,0xfb,0x21,0x19,0x40]}

tool: {"D490": [0x1e,0x26,0x44,0x54,0xfb,0x21,0xf9,0x3f]}

tool: {"D288": [3.141592653589,45.0]}

tool: D22: 45

tool: D255: 1024

tool: D167: 3

tool: D226: 4

tool: D157: 1

tool: D40: 64

tool: D383: 5

tool: D44: 256

tool: D261: 8

tool: D113: 32

tool: D295: 15

tool: D520: 35

tool: D123: 3072

tool: D456: 0.5

tool: D150: 6

tool: D476: 2

tool: D364: 27

tool: D558: 21

tool: D27: 16

tool: D21: 55

tool: D310: 151

tool: D178: 96

tool: D409: 1097

tool: D565: 512

tool: {"D472": [0.0,-0.0]}

tool: D7: 1e-05

tool: D575: 9223372036854775807

tool: {"D475": [0x1e,0x26,0x44,0x54,0xfb,0x21,0x09,0xc0]}

assistant: ```c
void GPS_Math_Bonne_EN_To_LatLon(double E, double N, double *phi, double *lambda, double phi0, double M0, double E0, double N0, double a, double b) {
  double p2;
  double po2;
  double a2;
  double b2;
  double e2;
  double e4;
  double e6;
  double M1;
  double m1;
  double c0;
  double c1;
  double c2;
  double c3;
  double A0;
  double A1;
  double A2;
  double A3;
  double j;
  double te4;
  double E1;
  double E2;
  double E3;
  double E4;
  double x;
  double phi0s;
  double lat;
  double phi0c;
  double phi0s2;
  double phi0s4;
  double phi0s6;
  double as;
  double phis;
  double phic;
  double dx;
  double dy;
  double mu;
  double mm;
  double MM;
  double asdy;
  double rho;
  double smu2;
  double smu4;
  double smu6;
  double smu8;
  double tol;
  phi0 = GPS_Math_Deg_To_Rad(phi0);
  M0 = GPS_Math_Deg_To_Rad(M0);
  phi0s = sin(phi0);
  p2 = (double)3.141592653589 * (double)2.0;
  po2 = (double)3.141592653589 / (double)2.0;
  if (M0 > 3.141592653589)
    M0 -= p2;
  a2 = a * a;
  b2 = b * b;
  e2 = (a2 - b2) / a2;
  e4 = e2 * e2;
  e6 = e2 * e4;
  j = (double)45.0 * e6 / (double)1024.0;
  te4 = (double)3.0 * e4;
  c0 = (double)1.0 - e2 / (double)4.0 - te4 / (double)64.0 - (double)5.0 * e6 / (double)256.0;
  c1 = (double)3.0 * e2 / (double)8.0 + te4 / (double)32.0 + j;
  c2 = (double)15.0 * e4 / (double)256.0 + j;
  c3 = (double)35.0 * e6 / (double)3072.0;
  phi0c = cos(phi0);
  m1 = phi0c / pow(((double)1.0 - e2 * phi0s * phi0s), (double)0.5);
  lat = c0 * phi0;
  phi0s2 = c1 * sin((double)2.0 * phi0);
  phi0s4 = c2 * sin((double)4.0 * phi0);
  phi0s6 = c3 * sin((double)6.0 * phi0);
  M1 = a * (lat - phi0s2 + phi0s4 - phi0s6);
  x = pow((double)1.0 - e2, (double)0.5);
  E1 = ((double)1.0 - x) / ((double)1.0 + x);
  E2 = E1 * E1;
  E3 = E2 * E1;
  E4 = E3 * E1;
  A0 = (double)3.0 * E1 / (double)2.0 - (double)27.0 * E3 / (double)32.0;
  A1 = (double)21.0 * E2 / (double)16.0 - (double)55.0 * E4 / (double)32.0;
  A2 = (double)151.0 * E3 / (double)96.0;
  A3 = (double)1097.0 * E4 / (double)512.0;
  if (!phi0s)
    as = (double)0.0;
  else
    as = a * m1 / phi0s;
  dx = E - E0;
  dy = N - N0;
  asdy = as - dy;
  rho = pow(dx * dx + asdy * asdy, (double)0.5);
  if (phi0 < (double)0.0)
    rho = -rho;
  MM = as + M1 - rho;
  mu = MM / (a * c0);
  smu2 = A0 * sin((double)2.0 * mu);
  smu4 = A1 * sin((double)4.0 * mu);
  smu6 = A2 * sin((double)6.0 * mu);
  smu8 = A3 * sin((double)8.0 * mu);
  *phi = mu + smu2 + smu4 + smu6 + smu8;
  tol = (double)0.00001;
  if (((po2 - tol) < fabs(*phi)) && (fabs(*phi) < po2 + tol))
    *lambda = M0;
  else {
    phic = cos(*phi);
    phis = sin(*phi);
    mm = phic / pow(((double)1.0 - e2 * phis * phis), (double)0.5);
    if (phi0 < (double)0.0) {
      dx = -dx;
      asdy = -asdy;
    }
    *lambda = M0 + rho * (atan2(dx, asdy)) / (a * mm);
  }
  if (*phi > po2)
    *phi = po2;
  else if (*phi < -po2)
    *phi = -po2;
  if (*lambda > 3.141592653589)
    *lambda -= p2;
  if (*lambda < -3.141592653589)
    *lambda += p2;
  if (*lambda > 3.141592653589)
    *lambda = 3.141592653589;
  else if (*lambda < -3.141592653589)
    *lambda = -3.141592653589;
  *lambda = GPS_Math_Rad_To_Deg(*lambda);
  *phi = GPS_Math_Rad_To_Deg(*phi);
  return;
}
```