[med-svn] [Git][med-team/ghmm][master] Removed unused ghmm_test.c

[Israel Komolehin (@Komolehin)]

Israel Komolehin pushed to branch master at Debian Med / ghmm


Commits:
f7306d7f by Komolehin Israel Timilehin at 2024-01-18T17:40:35+00:00
Removed unused ghmm_test.c

- - - - -


1 changed file:

- − debian/tests/ghmm_test.c


Changes:

=====================================
debian/tests/ghmm_test.c deleted
=====================================
@@ -1,132 +0,0 @@
-#include <stdio.h>
-#include <ghmm/rng.h>
-#include <ghmm/sequence.h>
-#include <ghmm/model.h>
-#include <ghmm/reestimate.h>
-
-/*
-  model with two states and three symbols
-  transition probability: 0->0: 0.9, 0->1: 0.1, 1->0: 0.1 and 1->1:0.9
-  state 0 has symbol 0 (probability 0.5) and 1 (probability 0.5)
-  state 1 has symbol 2
-*/
-int my_model()
-{
-  /* model structure, that contains states */
-  ghmm_dmodel my_model;
-  /* array of states */
-  ghmm_dstate model_states[2];
-
-  /* first state */
-  /*  probability of emmission of 0,1 or 2 */
-  double symbols_0_state[3]={0.5,0.5,0.0};
-  /* transition to which state is given in the following arrays */
-  int trans_id_0_state[2]={0,1};
-  /* transition probability from here to 0-state (self) and 1-state */
-  double trans_prob_0_state[2]={0.9,0.1};
-  /* transition probability from 0-state (self) and 1-state to this state */
-  double trans_prob_0_state_rev[2]={0.9,0.1};
-
-  /* second state , comments see above */
-  double symbols_1_state[3]={0.0,0.0,1.0};
-  int trans_id_1_state[2]={0,1};
-  double trans_prob_1_state[2]={0.1,0.9};
-  double trans_prob_1_state_rev[2]={0.1,0.9};
-  ghmm_dseq* my_output;
-
-  int pow_look[2] = {1,3};
-
-  /* flags indicating whether a state is silent */
-  /*int silent_array[2] =  {0,0};*/
-  
-  /* initialise state 0 */
-  /* start probability for this state */
-  model_states[0].pi = 0.5;
-  /* array with emission probabilities */
-  model_states[0].b=symbols_0_state;
-  /* number of fields in out_a and out_id */
-  model_states[0].out_states=2;
-  /* transition probability from this state */
-  model_states[0].out_a=trans_prob_0_state;
-  /* state ids belonging to the probability */
-  model_states[0].out_id=trans_id_0_state;
-  /* transition probability to this state */
-  /* in_states,in_id and in_a have the same function as above*/
-  model_states[0].in_states=2;
-  model_states[0].in_id=trans_id_0_state;
-  model_states[0].in_a=trans_prob_0_state_rev;
-  /* should emission probabilities be changed during reestimation? 1: no, else: yes*/
-  model_states[0].fix=0;
-  
-  /* initialise state 1 */
-  /* same meaning as above */
-  model_states[1].pi = 0.5;
-  model_states[1].b=symbols_1_state;
-  model_states[1].out_states=2;
-  model_states[1].out_a=trans_prob_1_state;
-  model_states[1].out_id=trans_id_1_state;
-  model_states[1].in_states=2;
-  model_states[1].in_id=trans_id_0_state;
-  model_states[1].in_a=trans_prob_1_state_rev;
-  model_states[1].fix=0;
-  
-    /* initialise model */
-  my_model.N=2; /* number of states, dimension of model.s */
-  my_model.M=3; /* number of symbols, dimension of states.b */
-  my_model.s=model_states; /* array of states */
-  my_model.prior=-1; /* probability of this model, used in a model array */
-
-  /*my_model.silent = silent_array;*/
-  my_model.pow_lookup = pow_look;
-  my_model.maxorder = 0;
-  my_model.model_type =0;
-  
-  
-    /* consistency check */
-  fprintf(stdout,"checking model:\n");
-  if (ghmm_dmodel_check(&my_model))
-    {
-      fprintf(stderr,"ghmm_dmodel_check failed!\n");
-      return 1;
-    }
-  fprintf(stdout,"model is ok\n");
-
-  /* print model parameters */
-  fprintf(stdout,"two_states_three_symbols model:\n");
-  ghmm_dmodel_print(stdout,&my_model);
-
-  /* generate sequences */
-  fprintf(stdout,"generating sequences:...");
-  my_output=ghmm_dmodel_generate_sequences(&my_model, /* model */
-				     0,   /* random seed */
-				     100, /* length of each sequence */
-				     100, /* no of sequences */
-		  			 100); /* maxT */ 
-  fprintf(stdout,"Done\n");
-  /*ghmm_dseq_print(stdout,my_output);*/
-
-  /* slight change of emission probabilities in state 0 */
-  symbols_0_state[0] = 0.6;
-  symbols_0_state[1] = 0.4;
-  symbols_0_state[2] = 0.0;
-
-  /* reestimation */
-  fprintf(stdout,"reestimating with Baum-Welch-algorithm...");
-  ghmm_dmodel_baum_welch(&my_model,my_output);
-
-  /* print the result */
-  fprintf(stdout,"Done\nthe result is:\n");
-  ghmm_dmodel_print(stdout,&my_model);
-
-  ghmm_dseq_free(&my_output);
-
-  return 0;
-}
-
-int main()
-{
-  /* Important! initialise rng  */
-//   ghmm_rng_init();
-
-  return my_model();
-}
\ No newline at end of file



View it on GitLab: https://salsa.debian.org/med-team/ghmm/-/commit/f7306d7fc95f7d5dba68ed001f04d055c6bd3c30

-- 
View it on GitLab: https://salsa.debian.org/med-team/ghmm/-/commit/f7306d7fc95f7d5dba68ed001f04d055c6bd3c30
You're receiving this email because of your account on salsa.debian.org.


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://alioth-lists.debian.net/pipermail/debian-med-commit/attachments/20240118/1f9a1d67/attachment-0001.htm>