#!/bin/csh
# script to generate generic neural network subroutine
# written by Alain Bellerive May 2002
#
# syntax netgen <inputfile> <outputfile> <subroutinename> 
# <inputfile>      = name of file written by JETNET3 JNDUMP subroutine
# <outputfile>     = name of outputfile for fortran/patchy program
# <subroutinename> = name of generated subroutine
#
setenv INPUTFILE $1
setenv OUTPUTFILE $2
setenv SUBNAME $3
#
echo netgen will generate subroutine $SUBNAME from data on $INPUTFILE
echo The source will be written on $OUTPUTFILE
#
# set up links for input
ln -s $INPUTFILE netgen.dat

# 
# compile and link the program
g77 -O -mcpu=i486 -mieee-fp -ffloat-store -fno-second-underscore -fno-f2c \
  -O2 -malign-double -Wall \
  netgen.f    \
  jetnet34.o  \
  /local/cern/pro/lib/libgraflib.a                 \
  /local/cern/pro/lib/libgrafX11.a                 \
  /local/cern/pro/lib/libpacklib.a                 \
  /local/cern/pro/lib/libmathlib.a                 \
  /local/cern/pro/lib/libkernlib.a                 \
 -o netgen.exe

#
# run the generation program
netgen.exe
#
rm -f netgen.dat netgen.intdat netgen.exe
#

echo Data processed...writing output
#
# write subroutine header
cat >$OUTPUTFILE <<EOF
      SUBROUTINE ${SUBNAME}(OUT,VARIN)
*  simulation of JETNET 3.4 neural net with hardwired parameters
*  routine skeleton generated by netgen - please fill in your own inputs
      implicit none
      INTEGER NVAR
      PARAMETER(NVAR=???)
      REAL OUT,VARIN(NVAR)
*
      INTEGER ILAY,IIN,IIM,IJC,IJT,I,IVAR
*
      INTEGER NLAY,NCO,NTH,NINP,MXNODE
      REAL TEMPN
*  parameters of the network size and temperature
EOF
#
# add in dimensions information on fort.11
cat netgen.dat1 >>$OUTPUTFILE
#
# dimension the arrays
cat >>$OUTPUTFILE <<EOF
      INTEGER NNODE(NLAY+1)
      REAL CO(NCO),TH(NTH),VN(MXNODE,NLAY+1),VNT
EOF
#
# add in data infomation on fort.12w
cat netgen.dat2 >>$OUTPUTFILE
#
# subroutine body
#
cat >>$OUTPUTFILE <<EOF
*
*   fill in your own input values into the first layer of the net

      DO IVAR=1,NVAR
        VN(IVAR,1)=VARIN(IVAR)
      ENDDO
*
*   loop over all the layers, calculating the neuron values for that layer
      IJC=1
      IJT=1
      DO ILAY=1,NLAY
*   loop over all the neurons in the layer being evaluated
        DO IIN=1,NNODE(ILAY+1)
*   starting value for neuron is the threshold value
          VNT=TH(IJT)
          IJT=IJT+1
*   now add in the values from the previous layer of neurons
          DO IIM=1,NNODE(ILAY)
            VNT=VNT+CO(IJC)*VN(IIM,ILAY)
*   and increment the coefficient pointer
            IJC=IJC+1
          ENDDO
*   apply the net transfer function
EOF
#
# add in the required function
cat netgen.dat3 >>$OUTPUTFILE
#
# end of the subrotuine
cat >>$OUTPUTFILE <<EOF
        ENDDO
      ENDDO
*   now we have the answer, copy to output node
      OUT=VN(1,NLAY+1)
      END
EOF
#
# all done, clear up
rm -f netgen.dat1 netgen.dat2 netgen.dat3
echo All finished
exit