// File: /workspace/github_c_makefile_star2_O3/markpeerdeman[P]mbemm/src/emath/post.c
// Function: write_edge_force
void write_edge_force(mesh theMesh)
{ point p;
  element e;
  double t1x, t1y, t1z, t2x, t2y, t2z, t3x, t3y, t3z;
  double fx, fy, fz;
  double fxt, fyt, fzt;
  int cnodes;
  int p1iscnode;
  int p2iscnode;
  int p3iscnode;
  /* double force1, force2; */

  printf("\n{---------------------------------------------------------------}\n");
  printf("{ POST output of the edge element forces.\n");
  printf("  format:\n");
  printf("  element_label force_x force_y force_z\n");
  printf("}\n");
  printf("edge_element_forces\n");
  fxt=0.0; fyt=0.0; fzt=0.0;

  /* the calculated tractions: */
  e = (element)element_dll(theMesh->mesh_e,START,0);
  while (e!=(element)NULL)
  { cnodes=0; p1iscnode=0; p2iscnode=0; p3iscnode=0;   
    /* get the traction forces at the three nodes of the element: */ 
	/* NODE 1: */
	p=(point)point_dll(theMesh->mesh_p,SEARCH,e->p1);
        if (p->type==CNODE) { cnodes++; p1iscnode=1; }
	t1x=p->fx; t1y=p->fy; t1z=p->fz;
	/* NODE 2: */
	p=(point)point_dll(theMesh->mesh_p,SEARCH,e->p2);
        if (p->type==CNODE) { cnodes++; p2iscnode=1; }
	t2x=p->fx; t2y=p->fy; t2z=p->fz;
	/* NODE 3: */
	p=(point)point_dll(theMesh->mesh_p,SEARCH,e->p3);
        if (p->type==CNODE) { cnodes++; p3iscnode=1; }
	t3x=p->fx; t3y=p->fy; t3z=p->fz;
  
	fx= 1.0/3.0*(t1x+t2x+t3x)*e->area*theMesh->G_SCALE;
	fy= 1.0/3.0*(t1y+t2y+t3y)*e->area*theMesh->G_SCALE;
	fz= 1.0/3.0*(t1z+t2z+t3z)*e->area*theMesh->G_SCALE;
	fxt+=fx; fyt+=fy; fzt+=fz;
        if (cnodes==2)
        { if ( (p1iscnode==1) && (p2iscnode==1) )
          { fx=fx/e->z3*theMesh->G_SCALE;
            fy=fy/e->z3*theMesh->G_SCALE;
            fz=fz/e->z3*theMesh->G_SCALE;
          }
          if ( (p1iscnode==1) && (p3iscnode==1) )
          { fx=fx/e->z2*theMesh->G_SCALE;
            fy=fy/e->z2*theMesh->G_SCALE;
            fz=fz/e->z2*theMesh->G_SCALE;
          }
          if ( (p2iscnode==1) && (p3iscnode==1) )
          { fx=fx/e->z1*theMesh->G_SCALE;
            fy=fy/e->z1*theMesh->G_SCALE;
            fz=fz/e->z1*theMesh->G_SCALE;
          }

	  printf("%4d % 8.8E % 8.8E % 8.8E (per unit length)\n", e->label,fx,fy,fz);
        }
    e=e->NEXT;
  }

  printf("end_edge_element_forces\n");

  printf("\n{---------------------------------------------------------------}\n");
  printf("{ POST output of the resulting force vector.\n");
  printf(" Force: %5.8E %5.8E %5.8E\n",fxt,fyt,fzt);
  printf(" Modulus: %5.8E\n",sqrt((double)(fxt*fxt+fyt*fyt+fzt*fzt)));

}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x48, %rsp
movq %rdi, %rbx
leaq 0xa4c4(%rip), %rdi      # 0xead4
callq 0x1060
leaq 0xa421(%rip), %rdi      # 0xea3d
callq 0x1060
leaq 0xa51c(%rip), %rdi      # 0xeb44
callq 0x1060
leaq 0xa433(%rip), %rdi      # 0xea67
callq 0x1060
leaq 0xa4d5(%rip), %rdi      # 0xeb15
callq 0x1060
leaq 0xa447(%rip), %rdi      # 0xea93
callq 0x1060
movq 0x8(%rbx), %rdi
xorl %esi, %esi
xorl %edx, %edx
callq 0xb570
testq %rax, %rax
je 0x4863
movq %rax, %r14
xorpd %xmm0, %xmm0
movapd %xmm0, 0x30(%rsp)
xorpd %xmm0, %xmm0
movsd %xmm0, 0x8(%rsp)
movq (%rbx), %rdi
movl 0x18(%r14), %edx
movl $0x2, %esi
callq 0xcc2f
movl 0x14(%rax), %r12d
xorl %r15d, %r15d
cmpl $0xf, %r12d
sete %r15b
movsd 0x48(%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movq (%rbx), %rdi
movl 0x1c(%r14), %edx
leal 0x1(%r15), %r13d
movups 0x38(%rax), %xmm0
movaps %xmm0, 0x20(%rsp)
movl $0x2, %esi
callq 0xcc2f
movl 0x14(%rax), %ebp
cmpl $0xf, %ebp
cmovnel %r15d, %r13d
movupd 0x38(%rax), %xmm0
addpd 0x20(%rsp), %xmm0
movapd %xmm0, 0x20(%rsp)
movq (%rbx), %rdi
movl 0x20(%r14), %edx
movsd 0x10(%rsp), %xmm0
addsd 0x48(%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl $0x2, %esi
callq 0xcc2f
movsd 0x10(%rsp), %xmm5
movl 0x14(%rax), %ecx
xorl %edx, %edx
cmpl $0xf, %ecx
sete %dl
addl %r13d, %edx
movupd 0x38(%rax), %xmm0
addpd 0x20(%rsp), %xmm0
movsd 0xb8(%r14), %xmm2
movsd 0x30(%rbx), %xmm1
mulpd 0x9cf7(%rip), %xmm0     # 0xe430
addsd 0x48(%rax), %xmm5
mulsd 0x9d0a(%rip), %xmm5     # 0xe450
mulsd %xmm2, %xmm5
movapd %xmm2, %xmm6
unpcklpd %xmm2, %xmm6            # xmm6 = xmm6[0],xmm2[0]
mulpd %xmm0, %xmm6
movapd %xmm1, %xmm3
unpcklpd %xmm1, %xmm3            # xmm3 = xmm3[0],xmm1[0]
mulpd %xmm3, %xmm6
mulsd %xmm1, %xmm5
cmpl $0x2, %edx
jne 0x4834
movapd %xmm5, %xmm2
movapd %xmm6, %xmm0
cmpl $0xf, %r12d
jne 0x47b3
movapd %xmm5, %xmm2
movapd %xmm6, %xmm0
cmpl $0xf, %ebp
jne 0x47b3
movsd 0x98(%r14), %xmm0
movapd %xmm5, %xmm2
divsd %xmm0, %xmm2
movapd %xmm0, %xmm4
unpcklpd %xmm0, %xmm4            # xmm4 = xmm4[0],xmm0[0]
movapd %xmm6, %xmm0
divpd %xmm4, %xmm0
mulpd %xmm3, %xmm0
mulsd %xmm1, %xmm2
cmpl $0xf, %r12d
jne 0x47db
cmpl $0xf, %ecx
jne 0x47db
movsd 0x90(%r14), %xmm4
divsd %xmm4, %xmm2
unpcklpd %xmm4, %xmm4            # xmm4 = xmm4[0,0]
divpd %xmm4, %xmm0
mulpd %xmm3, %xmm0
mulsd %xmm1, %xmm2
cmpl $0xf, %ebp
movsd %xmm5, 0x10(%rsp)
movapd %xmm6, 0x20(%rsp)
jne 0x480e
cmpl $0xf, %ecx
jne 0x480e
movsd 0x88(%r14), %xmm4
divsd %xmm4, %xmm2
unpcklpd %xmm4, %xmm4            # xmm4 = xmm4[0,0]
divpd %xmm4, %xmm0
mulpd %xmm3, %xmm0
mulsd %xmm1, %xmm2
movl 0x14(%r14), %esi
movapd %xmm0, %xmm1
unpckhpd %xmm0, %xmm1            # xmm1 = xmm1[1],xmm0[1]
leaq 0x9d28(%rip), %rdi      # 0xe549
movb $0x3, %al
callq 0x1090
movsd 0x10(%rsp), %xmm5
movapd 0x20(%rsp), %xmm6
movapd 0x30(%rsp), %xmm0
addpd %xmm6, %xmm0
movapd %xmm0, 0x30(%rsp)
movsd 0x8(%rsp), %xmm0
addsd %xmm5, %xmm0
movsd %xmm0, 0x8(%rsp)
movq 0x8(%r14), %r14
testq %r14, %r14
jne 0x467e
jmp 0x4877
xorpd %xmm0, %xmm0
movapd %xmm0, 0x30(%rsp)
xorpd %xmm0, %xmm0
movsd %xmm0, 0x8(%rsp)
leaq 0xa211(%rip), %rdi      # 0xea8f
callq 0x1060
leaq 0xa24a(%rip), %rdi      # 0xead4
callq 0x1060
leaq 0xa211(%rip), %rdi      # 0xeaa7
callq 0x1060
movaps 0x30(%rsp), %xmm0
movaps %xmm0, %xmm1
unpckhpd %xmm0, %xmm1            # xmm1 = xmm1[1],xmm0[1]
movaps %xmm1, 0x10(%rsp)
leaq 0x9c30(%rip), %rdi      # 0xe4e3
movsd 0x8(%rsp), %xmm2
movb $0x3, %al
callq 0x1090
movapd 0x10(%rsp), %xmm1
mulsd %xmm1, %xmm1
movapd 0x30(%rsp), %xmm0
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
movsd 0x8(%rsp), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm0, %xmm1
xorps %xmm0, %xmm0
sqrtsd %xmm1, %xmm0
leaq 0x9c0a(%rip), %rdi      # 0xe4fe
movb $0x1, %al
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x1090