/* 
  This subroutine accomplishes reweighting of posterior simulator output
  and/or construction of new functions of interest.  The code given here
  leaves the prior and functions of interest unchanged, and therefore
  accomplishes nothing.  This code must be replaced by code specific to
  the client's interests.
  Inputs:
  npar    Number of functions of interest in the investigator's posterior
          simulation file.
  thetai  Vector of dimension npar, containing the parameter values read from
          the investigator's simulation file in the current record.
  pri     Logarithm of the investigator's prior evaluated at thetai
  Outputs:
  nc      Number of functions of interest in the client's posterior simulation
          file.
  gc      Vector of dimension nc, containing the values of the client's
          function of interest
  priorc  Logarithm of the client's prior evaluated at thetai
*/
#include <math.h>
void client(int npar, double *thetai, double pri, 
            int *nc, double *gc, double *priorc)
{
   int i=1;
   *nc=npar;
   *priorc=pri;
   dcopy_(nc,thetai,&i,gc,&i);
}