/Home/aooliver/hepsoft/lcdroot/LCDRootApps/PhUtil/src/LCDNNTrackAttachGhostZpole.c

Go to the documentation of this file.

/*********************************************************
  LCDNNTrackAttachGhostZpole.c
  --------------------------------------------------------
  generated at Mon Jul 30 18:39:57 2001
  by snns2c ( Bernward Kett 1995 ) 
*********************************************************/

#include <math.h>

#define Act_Logistic(sum, bias)  ( (sum+bias<10000.0) ? ( 1.0/(1.0 + exp(-sum-bias) ) ) : 0.0 )
#ifndef NULL
#define NULL (void *)0
#endif

typedef struct UT {
          float act;         /* Activation       */
          float Bias;        /* Bias of the Unit */
          int   NoOfSources; /* Number of predecessor units */
   struct UT   **sources; /* predecessor units */
          float *weights; /* weights from predecessor units */
        } UnitType, *pUnit;

  /* Forward Declaration for all unit types */
  static UnitType Units[13];
  /* Sources definition section */
  static pUnit Sources[] =  {
Units + 1, Units + 2, Units + 3, Units + 4, Units + 5, 
Units + 1, Units + 2, Units + 3, Units + 4, Units + 5, 
Units + 1, Units + 2, Units + 3, Units + 4, Units + 5, 
Units + 1, Units + 2, Units + 3, Units + 4, Units + 5, 
Units + 1, Units + 2, Units + 3, Units + 4, Units + 5, 
Units + 1, Units + 2, Units + 3, Units + 4, Units + 5, 
Units + 6, Units + 7, Units + 8, Units + 9, Units + 10, Units + 11, 

  };

  /* Weigths definition section */
  static float Weights[] =  {
-25.288549, -4.046800, 5.016060, 0.513620, 0.626160, 
2.070870, -0.100070, 2.819970, 25.874180, -4.198960, 
7.204980, 6.118720, -0.542130, 0.093370, 5.098960, 
-4.418000, 2.607140, 7.771570, -2.675650, 18.379610, 
-9.906980, 6.152710, 2.138680, -0.293120, 36.910759, 
8.016760, -2.452140, -6.188510, -11.696310, -6.958790, 
-4.582510, -3.758210, -8.969600, 2.015460, 7.999760, -1.914720, 

  };

  /* unit definition section (see also UnitType) */
  static UnitType Units[13] = 
  {
    { 0.0, 0.0, 0, NULL , NULL },
    { /* unit 1 (Old: 1) */
      0.0, 0.858420, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 2 (Old: 2) */
      0.0, -0.234780, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 3 (Old: 3) */
      0.0, 0.597270, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 4 (Old: 4) */
      0.0, 0.585270, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 5 (Old: 5) */
      0.0, 0.605490, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 6 (Old: 6) */
      0.0, -5.710260, 5,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 7 (Old: 7) */
      0.0, -0.678110, 5,
       &Sources[5] , 
       &Weights[5] , 
      },
    { /* unit 8 (Old: 8) */
      0.0, -1.901810, 5,
       &Sources[10] , 
       &Weights[10] , 
      },
    { /* unit 9 (Old: 9) */
      0.0, -2.729140, 5,
       &Sources[15] , 
       &Weights[15] , 
      },
    { /* unit 10 (Old: 10) */
      0.0, -6.891270, 5,
       &Sources[20] , 
       &Weights[20] , 
      },
    { /* unit 11 (Old: 11) */
      0.0, 4.918040, 5,
       &Sources[25] , 
       &Weights[25] , 
      },
    { /* unit 12 (Old: 12) */
      0.0, 2.315280, 6,
       &Sources[30] , 
       &Weights[30] , 
      }

  };



int LCDNNTrackAttachGhostZpole(float *in, float *out, int init)
{
  int member, source;
  float sum;
  enum{OK, Error, Not_Valid};
  pUnit unit;


  /* layer definition section (names & member units) */

  static pUnit Input[5] = {Units + 1, Units + 2, Units + 3, Units + 4, Units + 5}; /* members */

  static pUnit Hidden1[6] = {Units + 6, Units + 7, Units + 8, Units + 9, Units + 10, Units + 11}; /* members */

  static pUnit Output1[1] = {Units + 12}; /* members */

  static int Output[1] = {12};

  for(member = 0; member < 5; member++) {
    Input[member]->act = in[member];
  }

  for (member = 0; member < 6; member++) {
    unit = Hidden1[member];
    sum = 0.0;
    for (source = 0; source < unit->NoOfSources; source++) {
      sum += unit->sources[source]->act
             * unit->weights[source];
    }
    unit->act = Act_Logistic(sum, unit->Bias);
  };

  for (member = 0; member < 1; member++) {
    unit = Output1[member];
    sum = 0.0;
    for (source = 0; source < unit->NoOfSources; source++) {
      sum += unit->sources[source]->act
             * unit->weights[source];
    }
    unit->act = Act_Logistic(sum, unit->Bias);
  };

  for(member = 0; member < 1; member++) {
    out[member] = Units[Output[member]].act;
  }

  return(OK);
}

---