source: tag/release-8/yao8/sample/jeuvie/Biocell.h @ 1

/***************************************************************************
                 module classe Biocell.h  -  description
 ***************************************************************************/
  forward (YREAL x1, YREAL x2, YREAL x3, YREAL x4, YREAL x5, YREAL x6, YREAL x7, YREAL x8,
                 YREAL sigma, YREAL seuil, YREAL mask)
  {
    double sum=0.0;
    for(int w=0; w<YNBI_Biocell-3; ++w) sum += Yting[w];

                //cas standard discret
                //if (sum==2.0) YS1=1.0; else YS1=0.0;

                //formulation continue
    YS1 = exp(-0.5 * pow((sum-seuil),2) / pow(sigma,2));

                //mask
                YS1 = YS1 * mask;
  }

  //===========================================================================
  backward (YREAL x1, YREAL x2, YREAL x3, YREAL x4, YREAL x5, YREAL x6, YREAL x7, YREAL x8,
                  YREAL sigma, YREAL seuil, YREAL mask)
  {
    double sum=0.0, expoSum;

    for(int w=0; w<YNBI_Biocell-3; w++) sum += Yting[w];

    expoSum= exp(-0.5 * pow((sum-seuil),2) / pow(sigma,2));

    YJ1I1  = -(sum-seuil)* expoSum/pow(sigma,2)            * mask ;
    YJ1I2  = YJ1I3=YJ1I4=YJ1I5=YJ1I6=YJ1I7=YJ1I8= YJ1I1   * mask ;
    YJ1I9  = expoSum*(pow((sum-seuil),2)/pow(sigma,3))   * mask ;
    YJ1I10 = (sum-seuil)/pow(sigma,2) * expoSum         * mask ;

                YJ1I11 = expoSum;
}
//************************** FIN DU MODULE Biocell ****************************