mrl.c File Reference

Calculation of maximum run length. More...

#include "libtpcmodext.h"

Go to the source code of this file.

Functions
int	mrl_between_tacs (double *y1, double *y2, int n)

Detailed Description

Calculation of maximum run length.

Author

Vesa Oikonen

Definition in file mrl.c.

Function Documentation

mrl_between_tacs()

int mrl_between_tacs	(	double *	y1,
		double *	y2,
		int	n )

Return the maximum run length between given n length arrays of data

Parameters

y1	Array of data1; may contain NaNs
y2	Array of data2; may contain NaNs
n	Nr of samples in array 1 and 2

Definition at line 13 of file mrl.c.

  {
  int i, mrl=0, rl=0;
  char last_sign=0, sign;
 
  if(n<1 || y1==NULL || y2==NULL) return(0);
  for(i=0; i<n; i++) {
    if(isnan(y1[i]) || isnan(y2[i])) continue;
    if(y1[i]>y2[i]) sign=1; else if(y1[i]<y2[i]) sign=-1; else sign=0;
    if(sign!=last_sign) {
      rl=0; last_sign=sign;
    } else {
      if(sign!=0) {rl++; if(rl>mrl) mrl=rl;}
    }
  }
  return(mrl);
}