single IFO bug fix

the generic polarization updates meant that there were arrays that have to have at least two elements (h+,hx, and beyond) but those arrays historically were allocated based on the number of interferometers. That's bad when it's a single interferometer. Fixed now.

src/BayesWaveLikelihood.c

@@ -930,7 +930,7 @@ double EvaluateExtrinsicConstantLogLikelihood(UNUSED struct Network *projection,
 double EvaluateExtrinsicMarkovianLogLikelihood(struct Network *projection, double *params, double **invSnf, double *Sn, struct Wavelet **geo, double **g, struct Data *data, double fmin, double fmax)
 {
    int i, n;
-   int NI,N;
+   int NI,NP,N;
    int ifo;
    int imin,imax;
 
@@ -942,13 +942,20 @@ double EvaluateExtrinsicMarkovianLogLikelihood(struct Network *projection, doubl
 
    NI   = data->NI;
    N = data->N;
+  NP = data->Npol;
+  
+  /*
+   if Npol = 1 assume elliptical polarization, but still need room for hx
+   */
+  if(NP<2) NP = 2;
+  
    //imin = data->imin;
    //imax = data->imax;
    imin = (int)(fmin*data->Tobs);
    imax = (int)(fmax*data->Tobs);
 
    d   = data->s;
-   t   = dmatrix(0,NI-1,0,N-1); //template
+   t   = dmatrix(0,NP-1,0,N-1); //template
    h   = dmatrix(0,NI-1,0,N-1); //response
    r   = dmatrix(0,NI-1,0,N-1); //residual
 
@@ -977,7 +984,7 @@ double EvaluateExtrinsicMarkovianLogLikelihood(struct Network *projection, doubl
 
    free_dmatrix(h,0,NI-1,0,N-1);
    free_dmatrix(r,0,NI-1,0,N-1);
-   free_dmatrix(t,0,NI-1,0,N-1);
+   free_dmatrix(t,0,NP-1,0,N-1);
   
    return logL;
 }