wrapper_cseq.i

General Hidden Markov Model Library 一个通用的隐马尔科夫模型的C代码库

%{
#include <ghmm/sequence.h>
#include <ghmm/obsolete.h>
%}

/*==========================================================================
  ===== continous sequences ================================================ */

typedef struct ghmm_cseq {
  /** sequence array. sequence[i][j] = j-th symbol of i-th seq. */
  double **seq;
  /** array of sequence length */
  int *seq_len;
#ifdef GHMM_OBSOLETE
  /**  array of sequence labels */
    long *seq_label;
#endif /* GHMM_OBSOLETE */
  /**  array of sequence IDs*/
  double *seq_id;
  /** positive! sequence weights.  default is 1 = no weight */
  double *seq_w;
  /** total number of sequences */
  long seq_number;
  /** reserved space for sequences is always >= seq_number */
  long capacity;
  /** sum of sequence weights */
  double total_w;

  /** flags (internal) */
  unsigned int flags;
} ghmm_cseq;

extern int ghmm_cseq_free(ghmm_cseq **csq);
extern ghmm_cseq* ghmm_cseq_calloc(long number);

%extend ghmm_cseq {
        ghmm_cseq(long number) { return ghmm_cseq_calloc(number); }
        ghmm_cseq(double* seq, int length) {
            ghmm_cseq* self = ghmm_cseq_calloc(1);
            self->seq[0] = seq;
            self->seq_len[0] = length;
            return self;
        }
        ghmm_cseq(double** seqs, int* lengths, int number) {
            ghmm_cseq* self = ghmm_cseq_calloc(number);
            free(self->seq);
            free(self->seq_len);
            self->seq = seqs;
            self->seq_len = lengths;
            return self;
        }
        ~ghmm_cseq() { ghmm_cseq_free(&self); }

        //ghmm_cseq** truncate(int sqd_fields, double trunc_ratio, int seed);

        ghmm_cseq* get_singlesequence(int index);

%apply SWIGTYPE * DISOWN {ghmm_cseq* self};
        int subseq_free();
%clear ghmm_cseq* self;

        int add(ghmm_cseq *source);

        void clean();

        int partition(ghmm_cseq *sqd_train, ghmm_cseq *sqd_test, double train_ratio);

        void copy_all(long t_num, ghmm_cseq *source, long s_num);

        double* getSequence(int index) { return self->seq[index]; }
        void setSequence(int seqno, double *O) { self->seq[seqno] = O; }
        double getSymbol(int seqno, int index) { return self->seq[seqno][index]; }
        void setSymbol(int seqno, int index, double value) { self->seq[seqno][index] = value; }

        int  getLength(int i) { return self->seq_len[i]; }
        void setLength(int i, int len) { self->seq_len[i] = len; }

        double getWeight(int i) { return self->seq_w[i]; }
        void   setWeight(int i, double w) { self->seq_w[i] = w; }

        void write(char* filename, int discrete = 0)
            {
                FILE* file;
                file = fopen(filename, "at");
                ghmm_cseq_print(self, file, discrete);
                fclose(file);
            }
}

%apply int *OUTPUT {int *sqd_number};
extern ghmm_cseq **ghmm_cseq_read(const char *filename, int *sqd_number);

REFERENCE_ARRAY(ghmm_cseq, cseq_ptr)