// File: /workspace/github_c_makefile_star2_O2/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 0x26733
movl %edx, %r15d
movq %rsi, %r14
movq %rdi, %rbx
movsd 0x6efe0(%rip), %xmm1    # 0x95520
divsd %xmm0, %xmm1
leaq 0xa7325(%rip), %r12     # 0xcd870
movq (%r12), %rax
movl 0x1c0(%rax), %ebp
decl %ebp
movl %edx, %ecx
leaq 0x71d2c(%rip), %rdx     # 0x9828c
movslq (%rdx,%rcx,4), %rcx
addq %rdx, %rcx
jmpq *%rcx
movl $0x30f20, %r14d         # imm = 0x30F20
xorl %r13d, %r13d
movsd %xmm1, (%rsp)
movslq 0x1bc(%rax), %rcx
cmpq %rcx, %r13
jge 0x26733
cmpl $0x1, -0x18(%rax,%r14)
jne 0x265b2
movq (%rbx), %rcx
movslq %r15d, %r15
movq (%rcx,%r15,8), %rdi
movl (%rax,%r14), %esi
movapd %xmm1, %xmm0
callq 0x33a4f
movsd (%rsp), %xmm1
incl %r15d
movq (%r12), %rax
incq %r13
addq $0x340, %r14            # imm = 0x340
jmp 0x26577
movq 0x10(%rbx), %rdi
movl 0x8(%rbx), %esi
movapd %xmm1, %xmm0
callq 0x33a4f
xorl %r15d, %r15d
movslq 0x8(%rbx), %rax
cmpq %rax, %r15
movsd 0x6ef40(%rip), %xmm0    # 0x95520
jge 0x26733
movq (%rbx), %rax
movq (%rax,%r15,8), %rdi
divsd (%r14,%r15,8), %xmm0
movl %ebp, %esi
callq 0x33a4f
incq %r15
jmp 0x265d1
movq 0x10(%rbx), %rdi
movl 0x8(%rbx), %esi
movapd %xmm1, %xmm0
callq 0x33a4f
pushq $0x8
popq %r12
xorl %r13d, %r13d
leaq 0x8(%rsp), %r15
movslq 0x8(%rbx), %rax
cmpq %rax, %r13
movsd 0x6eef6(%rip), %xmm0    # 0x95520
jge 0x26733
movq (%rbx), %rax
movq -0x8(%rax,%r12), %rdi
divsd (%r14,%r13,8), %xmm0
movl %ebp, %esi
callq 0x33a4f
movq (%rbx), %rdi
addq %r12, %rdi
movsd 0x6eecd(%rip), %xmm0    # 0x95520
divsd (%r14,%r13,8), %xmm0
callq 0x39a42
movq (%rbx), %rax
movl (%rax,%r12), %edi
callq 0x3930e
movq %rax, 0x8(%rsp)
movq %rdx, 0x10(%rsp)
movq %r15, %rdi
movsd 0x6f749(%rip), %xmm0    # 0x95dc8
callq 0x39a42
movq (%rbx), %rdi
addq %r12, %rdi
movq 0x8(%rsp), %rsi
movq 0x10(%rsp), %rdx
callq 0x39b81
movq 0x8(%rsp), %rdi
movq 0x10(%rsp), %rsi
callq 0x3929f
incq %r13
addq $0x18, %r12
jmp 0x2661b
movq (%rbx), %rdi
movapd %xmm1, %xmm0
movl %ebp, %esi
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x33a4f
movq (%rbx), %rdi
movapd %xmm1, %xmm0
movl %ebp, %esi
movsd %xmm1, (%rsp)
callq 0x33a4f
leaq 0x8(%rbx), %r14
movq %r14, %rdi
movsd (%rsp), %xmm0
callq 0x39a42
movl 0x8(%rbx), %edi
callq 0x3930e
leaq 0x8(%rsp), %rbx
movq %rax, (%rbx)
movq %rdx, 0x8(%rbx)
movsd 0x6f6b8(%rip), %xmm0    # 0x95dc8
movq %rbx, %rdi
callq 0x39a42
movq (%rbx), %rsi
movq 0x8(%rbx), %rdx
movq %r14, %rdi
callq 0x39b81
movq (%rbx), %rdi
movq 0x8(%rbx), %rsi
callq 0x3929f
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq