// File: /workspace/github_c_makefile_star2_O0/markpeerdeman[P]mbemm/src/emath/singular.c
// Function: make_S
void make_S( const matrix S, const double beta )
{ /* IN:          Matrix S, beta
   * OUT:         nothing
   * MODIFIES:    matrix S
   * USES:        data from theElement, beta as calculated in make_D()
   * DESCRIPTION: Fills matrix S
   */

int  i;    /* counter */
double i312, /* :temporary storage: */
	 i321;

   for (i=1; i<=3; i++)
   { i312 = (n3[i]-n1[i])/theElement->z2;
	 i321 = (n3[i]-n2[i])/theElement->z1;

	/* first component: */
	 M(S,i,1) = i312 * (cosA2 + cosA3)
			  - i321 * (1.0 - cosA1);
     M(S,i,1) = M(S,i,1) / sinA3;

	/* second component: */
     M(S,i,2) = cos(theElement->alpha3-theElement->alpha1);
	 M(S,i,2) = cosA3 - M(S,i,2) + beta*sinA3*sinA3;
	 M(S,i,2) = i321 * M(S,i,2);
	 M(S,i,2) = -i312 * (1.0-cosA1) + M(S,i,2);
	 M(S,i,2) = M(S,i,2) * theElement->z2 / theElement->z1 / sinA3;

	/* third component: */
	 M(S,i,3) = cosA3 + cosA2 - beta*sinA2*sinA3;
	 M(S,i,3) = M(S,i,3) * i321;
	 M(S,i,3) = -i312 * (1.0-cosA1) + M(S,i,3);
	 M(S,i,3) = M(S,i,3) * theElement->z2 / theElement->z1 / sinA2;

   } /* for i */

}

// Assembly:
subq $0x28, %rsp
movq %rdi, 0x20(%rsp)
movsd %xmm0, 0x18(%rsp)
movl $0x1, 0x14(%rsp)
cmpl $0x3, 0x14(%rsp)
jg 0xf831
movslq 0x14(%rsp), %rcx
leaq 0xf402(%rip), %rax      # 0x1e810
movsd (%rax,%rcx,8), %xmm0
movslq 0x14(%rsp), %rcx
leaq 0xf3b1(%rip), %rax      # 0x1e7d0
subsd (%rax,%rcx,8), %xmm0
movq 0xf365(%rip), %rax      # 0x1e790
divsd 0x90(%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
movslq 0x14(%rsp), %rcx
leaq 0xf3cb(%rip), %rax      # 0x1e810
movsd (%rax,%rcx,8), %xmm0
movslq 0x14(%rsp), %rcx
leaq 0xf39a(%rip), %rax      # 0x1e7f0
subsd (%rax,%rcx,8), %xmm0
movq 0xf32e(%rip), %rax      # 0x1e790
divsd 0x88(%rax), %xmm0
movsd %xmm0, (%rsp)
movsd 0x8(%rsp), %xmm0
movsd 0xf333(%rip), %xmm2     # 0x1e7b0
addsd 0xf33b(%rip), %xmm2     # 0x1e7c0
movsd (%rsp), %xmm1
movsd 0x9b76(%rip), %xmm3     # 0x19008
subsd 0xf306(%rip), %xmm3     # 0x1e7a0
mulsd %xmm3, %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, (%rax)
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd (%rax), %xmm0
divsd 0xf2bd(%rip), %xmm0     # 0x1e7b8
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, (%rax)
movq 0xf273(%rip), %rax      # 0x1e790
movsd 0x80(%rax), %xmm0
movq 0xf264(%rip), %rax      # 0x1e790
subsd 0x70(%rax), %xmm0
callq 0x10a0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, 0x8(%rax)
movsd 0xf266(%rip), %xmm1     # 0x1e7c0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
subsd 0x8(%rax), %xmm1
movsd 0x18(%rsp), %xmm0
mulsd 0xf234(%rip), %xmm0     # 0x1e7b8
movsd 0xf22c(%rip), %xmm2     # 0x1e7b8
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, 0x8(%rax)
movsd (%rsp), %xmm0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
mulsd 0x8(%rax), %xmm0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, 0x8(%rax)
movsd 0x8(%rsp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
movsd 0x99f6(%rip), %xmm2     # 0x19008
subsd 0xf186(%rip), %xmm2     # 0x1e7a0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd 0x8(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, 0x8(%rax)
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd 0x8(%rax), %xmm0
movq 0xf113(%rip), %rax      # 0x1e790
mulsd 0x90(%rax), %xmm0
movq 0xf104(%rip), %rax      # 0x1e790
divsd 0x88(%rax), %xmm0
divsd 0xf11c(%rip), %xmm0     # 0x1e7b8
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, 0x8(%rax)
movsd 0xf100(%rip), %xmm1     # 0x1e7c0
addsd 0xf0e8(%rip), %xmm1     # 0x1e7b0
movsd 0x18(%rsp), %xmm0
mulsd 0xf0d2(%rip), %xmm0     # 0x1e7a8
movsd 0xf0da(%rip), %xmm2     # 0x1e7b8
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, 0x10(%rax)
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd 0x10(%rax), %xmm0
mulsd (%rsp), %xmm0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, 0x10(%rax)
movsd 0x8(%rsp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
movsd 0x988d(%rip), %xmm2     # 0x19008
subsd 0xf01d(%rip), %xmm2     # 0x1e7a0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd 0x10(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, 0x10(%rax)
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd 0x10(%rax), %xmm0
movq 0xefaa(%rip), %rax      # 0x1e790
mulsd 0x90(%rax), %xmm0
movq 0xef9b(%rip), %rax      # 0x1e790
divsd 0x88(%rax), %xmm0
divsd 0xefa3(%rip), %xmm0     # 0x1e7a8
movq 0x20(%rsp), %rax
movq 0x8(%rax), %rax
movl 0x14(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movq (%rax,%rcx,8), %rax
movsd %xmm0, 0x10(%rax)
movl 0x14(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x14(%rsp)
jmp 0xf3f7
addq $0x28, %rsp
retq
nopw %cs:(%rax,%rax)