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

Go to the documentation of this file.

/*-----------------------------------------------------------------------------
  $Id: LCDNNHFSelect.c,v 1.2 2001/05/11 21:35:52 toshi Exp $
  -----------------------------------------------------------------------------

  $Log: LCDNNHFSelect.c,v $
  Revision 1.2  2001/05/11 21:35:52  toshi
  Add CVS tags $Id: $ and $Log: $.


*/
/*********************************************************
  LCDNNHFSelect.c
  --------------------------------------------------------
  generated at Wed Feb 21 08:56:51 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[] =  {
-5.726510, 6.655230, -11.942260, -9.211320, 0.260430, 
-26.757999, 27.343969, -2.859360, 15.351010, 5.491600, 
-3.124630, -23.940670, -5.141180, 3.643160, -1.002410, 
-6.723220, 12.006890, -7.269670, -2.353060, 1.204720, 
50.095360, 4.463220, -0.085040, 1.576510, 0.387810, 
0.315170, -21.674749, -13.851600, -5.173480, -5.121190, 
10.226540, 3.393190, -4.233980, -9.292990, 3.493660, -4.339010, 

  };

  /* unit definition section (see also UnitType) */
  static UnitType Units[13] = 
  {
    { 0.0, 0.0, 0, NULL , NULL },
    { /* unit 1 (Old: 1) */
      0.0, 1.000000, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 2 (Old: 2) */
      0.0, 0.998030, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 3 (Old: 3) */
      0.0, 0.916740, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 4 (Old: 4) */
      0.0, 0.646710, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 5 (Old: 5) */
      0.0, 0.270800, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 6 (Old: 6) */
      0.0, 4.060990, 5,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 7 (Old: 7) */
      0.0, -6.909750, 5,
       &Sources[5] , 
       &Weights[5] , 
      },
    { /* unit 8 (Old: 8) */
      0.0, 11.493260, 5,
       &Sources[10] , 
       &Weights[10] , 
      },
    { /* unit 9 (Old: 9) */
      0.0, -0.767500, 5,
       &Sources[15] , 
       &Weights[15] , 
      },
    { /* unit 10 (Old: 10) */
      0.0, -21.509930, 5,
       &Sources[20] , 
       &Weights[20] , 
      },
    { /* unit 11 (Old: 11) */
      0.0, 12.785520, 5,
       &Sources[25] , 
       &Weights[25] , 
      },
    { /* unit 12 (Old: 12) */
      0.0, 1.670950, 6,
       &Sources[30] , 
       &Weights[30] , 
      }

  };



int LCDNNHFSelect(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);
}

---