src/common/filter_net.c

/*********************************************************
  test_1/filter_net.c
  --------------------------------------------------------
  generated at Tue Jul 19 14:46:44 2005
  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[10];
  /* Sources definition section */
  static pUnit Sources[] =  {
Units + 1, Units + 2, Units + 3, 
Units + 1, Units + 2, Units + 3, 
Units + 1, Units + 2, Units + 3, 
Units + 1, Units + 2, Units + 3, 
Units + 1, Units + 2, Units + 3, 
Units + 4, Units + 5, Units + 6, Units + 7, Units + 8, 

  };

  /* Weigths definition section */
  static float Weights[] =  {
2.075310, 3.843340, 0.061740, 
-0.027020, 0.112820, 0.606700, 
-0.672260, 0.391710, 0.799930, 
0.888670, 3.404130, 0.143620, 
-0.815620, -0.186080, -0.807600, 
-4.194760, -0.058950, -1.101900, -3.117450, 1.956400, 

  };

  /* unit definition section (see also UnitType) */
  static UnitType Units[10] = 
  {
    { 0.0, 0.0, 0, NULL , NULL },
    { /* unit 1 (unit) */
      0.0, 1.000000, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 2 (unit) */
      0.0, 0.998030, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 3 (unit) */
      0.0, 0.916740, 0,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 4 (unit) */
      0.0, 0.406360, 3,
       &Sources[0] , 
       &Weights[0] , 
      },
    { /* unit 5 (unit) */
      0.0, -0.057010, 3,
       &Sources[3] , 
       &Weights[3] , 
      },
    { /* unit 6 (unit) */
      0.0, -0.434500, 3,
       &Sources[6] , 
       &Weights[6] , 
      },
    { /* unit 7 (unit) */
      0.0, 0.663390, 3,
       &Sources[9] , 
       &Weights[9] , 
      },
    { /* unit 8 (unit) */
      0.0, 0.311790, 3,
       &Sources[12] , 
       &Weights[12] , 
      },
    { /* unit 9 (unit) */
      0.0, 0.392790, 5,
       &Sources[15] , 
       &Weights[15] , 
      }

  };



int filter_net(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[3] = {Units + 1, Units + 2, Units + 3}; /* members */

  static pUnit Hidden1[5] = {Units + 4, Units + 5, Units + 6, Units + 7, Units + 8}; /* members */

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

  static int Output[1] = {9};

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

  for (member = 0; member < 5; 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);
}

---