src/warping.h

Go to the documentation of this file.

/***************************************************************************
 *   Copyright (C) 2008 by Matthias Ihrke   *
 *   mihrke@uni-goettingen.de   *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/

#ifndef WARPING_H
#define WARPING_H
#include "mathadd.h"
#include "definitions.h"
#include "time_frequency.h"
#include "regularization.h"
#include "eeg.h"

#ifdef __cplusplus
extern "C" {
#endif

  /*-----------------------------------------------------------
     -  slope constraint for DTW
     ---------------------------------------------------------*/
  typedef enum {
     SLOPE_CONSTRAINT_NONE=0,
     SLOPE_CONSTRAINT_LAX,
     SLOPE_CONSTRAINT_MEDIUM,
     SLOPE_CONSTRAINT_SEVERE
  } SlopeConstraint;

  /*-----------------------------------------------------------
     - WARPING -
     ---------------------------------------------------------*/

  typedef struct{
     int *t1; 
     int *t2; 
     int n1; 
     int n2; 
     int n; 
  } WarpPath;
    

#define warppath_CHECK( flag, w )                                                   \
  if(!( (w)->ndim==2 && (w)->dtype==UINT )){                                    \
     char *dts="";                                                                      \
     array_DTYPESTRING( dts, w->dtype );                                        \
     errprintf("not a warppath, ndim=%i, dtype=%s\n", w->ndim, dts );   \
     flag=FALSE;                                                                        \
  } else { flag=TRUE; }                                                             
  
  /* ---------------------------------------------------------------------------- 
      -- Timewarping                                                            -- 
      ---------------------------------------------------------------------------- */

  Array*  matrix_dtw_cumulate ( Array *mat, bool alloc, OptArgList *optargs );
  Array*  matrix_dtw_backtrack ( const Array *d );

  Array*  dtw_add_signals( const Array *s1, const Array *s2, const Array *path, OptArgList *opts );

  
  EEG* eeg_dtw_hierarchical( EEG *eeg_in, const double **distmatrix, 
                                      EEG *out, OptArgList *optargs );

  EEG* eeg_gibbons( EEG *eeg, int stimulus_marker, int response_marker, double k );

#ifdef EXPERIMENTAL
  Array* multiwarp( Array *in, uint window, OptArgList *opts );
  Array* multiwarp_add( Array *in, Array *times, Array *path );
#endif

  /******************************************************************/
  /******************************************************************/

  void       dtw_cumulate_matrix  ( double **d, int M, int N, OptArgList *opts );
  /* void dtw_cumulate_matrix_chiba_band( double **d, int M, int N, double restriction );*/
  WarpPath*  dtw_backtrack        ( const double **d, int M, int N, WarpPath *P );

   WarpPath* init_warppath ( WarpPath *path, int n1, int n2 );
   void      free_warppath ( WarpPath *p );
   void      reset_warppath( WarpPath *P, int n1, int n2 );
   void      print_warppath( FILE *out, WarpPath *P );
 
#ifdef __cplusplus
}
#endif

#endif