heartcorr.c

Go to the documentation of this file.

/*****************************************************************************/
#include "libtpcidi.h"
/*****************************************************************************/
 
/*****************************************************************************/

int heartcorr(
  int verbose
) {
  if(verbose>0) printf("heartcorr()\n");
 
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/
int heartRecoverySpilloverCorrectionFactors(
  double R,
  double d,
  double s,
  double Vb,
  double *FMM,
  double *FMB,
  double *FBM,
  double *FBB
) {
  double hp;
 
  if(R<=0.0 || d<=0.0 || s<=0.0) return(1);
  if(FMM==NULL || FMB==NULL || FBM==NULL || FBB==NULL) return(2);
  if(Vb<0.0 || Vb>=1.0) return(3);
  hp = exp(-R*R/(2.0*s*s));
  *FMB = hp - exp(-(R+d)*(R+d)/(2.0*s*s));
  *FMM = erf( (d/2.0) / (M_SQRT2*s) );
  *FBB = 1.0 - hp;
  *FBM = Vb + 0.5*(1.0 - *FMM);  //0.6 - 0.5 * *FMM;
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/