void filter_boundary_MEX(double *out,
			 double *in,  int lin,
			 double *num, int lnum, 
			 int den_length, 
			 double *d1, int ld1, 
			 double *d2, int ld2, 
			 int condition);