// File: /workspace/github_c_makefile_star2_O1/olcaytaner[P]Isell/hmm.c
// Function: update_state_model
void update_state_model(void *model, double accum, double *accumMix, HMM_STATE_TYPE statetype)
{
 /*!Last Changed 27.03.2007*/
 Gaussianptr g;
 Gaussianmixtureptr gm;
 Multinomialptr m;
 Dirichletptr d;
 Dirichletmixtureptr dm;
 int i, j;
 double mul = (1.0 / accum);
 int numrealfeatures = current_dataset->multifeaturecount - 1;
 matrix covprior;
 switch (statetype)
  {
   case          STATE_DISCRETE:m = (Multinomialptr) model;
                                j = 0;
                                for (i = 0; i < current_dataset->featurecount; i++)
                                  if (current_dataset->features[i].type == STRING)
                                   {
                                    multiply_array(m->prob[j], mul, current_dataset->features[i].valuecount);
                                    j++;
                                   }
                                break;
   case          STATE_GAUSSIAN:g = (Gaussianptr) model;
                                multiply_array(g->mean, mul, numrealfeatures);
                                matrix_multiply_constant(&(g->covariance), mul);
                                covprior = identity(g->covariance.row);
                                matrix_multiply_constant(&(covprior), PRIOR);
                                matrix_add(&(g->covariance), covprior);
                                matrix_free(covprior);
                                break;
   case  STATE_GAUSSIAN_MIXTURE:gm = (Gaussianmixtureptr) model;
                                multiply_array(gm->weights, mul, gm->componentcount);
                                for (i = 0; i < gm->componentcount; i++)
                                 {
                                  multiply_array(gm->components[i].mean, (1.0 / accumMix[i]), numrealfeatures);
                                  matrix_multiply_constant(&(gm->components[i].covariance), (1.0 / accumMix[i]));
                                  covprior = identity(gm->components[i].covariance.row);
                                  matrix_multiply_constant(&(covprior), PRIOR);
                                  matrix_add(&(gm->components[i].covariance), covprior);
                                  matrix_free(covprior);
                                 }
                                break;
   case         STATE_DIRICHLET:d = (Dirichletptr) model;
                                multiply_array(d->alpha, mul, numrealfeatures);
                                break;
   case STATE_DIRICHLET_MIXTURE:dm = (Dirichletmixtureptr) model;
                                multiply_array(dm->weights, mul, dm->componentcount);
                                for (i = 0; i < dm->componentcount; i++)
                                  multiply_array(dm->components[i].alpha, (1.0 / accumMix[i]), numrealfeatures);
                                break;
  }
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
cmpl $0x4, %edx
ja 0x2987f
movq %rsi, %r14
movq %rdi, %rbx
movsd 0x78eb6(%rip), %xmm2    # 0xa2520
divsd %xmm0, %xmm2
leaq 0xa91fb(%rip), %r15     # 0xd2870
movq (%r15), %rax
movl 0x1c0(%rax), %ebp
decl %ebp
movl %edx, %eax
leaq 0x7bc2b(%rip), %rcx     # 0xa52b4
movslq (%rcx,%rax,4), %rax
addq %rcx, %rax
jmpq *%rax
movq (%r15), %rax
cmpl $0x0, 0x1bc(%rax)
jle 0x2987f
xorl %r14d, %r14d
movl $0x30f20, %r12d         # imm = 0x30F20
xorl %r13d, %r13d
movsd %xmm2, (%rsp)
cmpl $0x1, -0x18(%rax,%r12)
jne 0x296da
movq (%rbx), %rcx
movslq %r13d, %r13
movq (%rcx,%r13,8), %rdi
movl (%rax,%r12), %esi
movapd %xmm2, %xmm0
callq 0x378e7
movsd (%rsp), %xmm2
incl %r13d
incq %r14
movq (%r15), %rax
movslq 0x1bc(%rax), %rcx
addq $0x340, %r12            # imm = 0x340
cmpq %rcx, %r14
jl 0x296b3
jmp 0x2987f
movq 0x10(%rbx), %rdi
movl 0x8(%rbx), %esi
movapd %xmm2, %xmm0
callq 0x378e7
cmpl $0x0, 0x8(%rbx)
jle 0x2987f
xorl %r15d, %r15d
movq (%rbx), %rax
movq (%rax,%r15,8), %rdi
movsd 0x78dfc(%rip), %xmm0    # 0xa2520
divsd (%r14,%r15,8), %xmm0
movl %ebp, %esi
callq 0x378e7
incq %r15
movslq 0x8(%rbx), %rax
cmpq %rax, %r15
jl 0x29715
jmp 0x2987f
movq 0x10(%rbx), %rdi
movl 0x8(%rbx), %esi
movapd %xmm2, %xmm0
callq 0x378e7
cmpl $0x0, 0x8(%rbx)
jle 0x2987f
movl $0x8, %r12d
xorl %r13d, %r13d
leaq 0x8(%rsp), %r15
movq (%rbx), %rax
movq -0x8(%rax,%r12), %rdi
movsd 0x78da6(%rip), %xmm0    # 0xa2520
divsd (%r14,%r13,8), %xmm0
movl %ebp, %esi
callq 0x378e7
movq (%rbx), %rdi
addq %r12, %rdi
movsd 0x78d8b(%rip), %xmm0    # 0xa2520
divsd (%r14,%r13,8), %xmm0
callq 0x3e50d
movq (%rbx), %rax
movl (%rax,%r12), %edi
callq 0x3dc37
movq %rax, 0x8(%rsp)
movq %rdx, 0x10(%rsp)
movq %r15, %rdi
movsd 0x795ff(%rip), %xmm0    # 0xa2dc0
callq 0x3e50d
movq (%rbx), %rdi
addq %r12, %rdi
movq 0x8(%rsp), %rsi
movq 0x10(%rsp), %rdx
callq 0x3e64c
movq 0x8(%rsp), %rdi
movq 0x10(%rsp), %rsi
callq 0x3dbad
incq %r13
movslq 0x8(%rbx), %rax
addq $0x18, %r12
cmpq %rax, %r13
jl 0x2976a
jmp 0x2987f
movq (%rbx), %rdi
movapd %xmm2, %xmm0
movl %ebp, %esi
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x378e7
movq (%rbx), %rdi
movapd %xmm2, %xmm0
movl %ebp, %esi
movsd %xmm2, (%rsp)
callq 0x378e7
leaq 0x8(%rbx), %r14
movq %r14, %rdi
movsd (%rsp), %xmm0
callq 0x3e50d
movl 0x8(%rbx), %edi
callq 0x3dc37
leaq 0x8(%rsp), %rbx
movq %rax, (%rbx)
movq %rdx, 0x8(%rbx)
movsd 0x79564(%rip), %xmm0    # 0xa2dc0
movq %rbx, %rdi
callq 0x3e50d
movq (%rbx), %rsi
movq 0x8(%rbx), %rdx
movq %r14, %rdi
callq 0x3e64c
movq (%rbx), %rdi
movq 0x8(%rbx), %rsi
callq 0x3dbad
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq