removed unuused constants.

This commit is contained in:
David Williams 2024-11-06 13:54:50 +01:00
parent 27d7e6915c
commit a2447f3b81
1 changed files with 3 additions and 148 deletions

View File

@ -1,89 +1,31 @@
!**** Declarations !**** Declarations
real*8 pi,infty,zero real*8 pi
real*8 scan_res
real*8 hart2eV, eV2hart real*8 hart2eV, eV2hart
real*8 hart2icm, icm2hart real*8 hart2icm, icm2hart
real*8 eV2icm, icm2eV real*8 eV2icm, icm2eV
real*8 deg2rad, rad2deg real*8 deg2rad, rad2deg
integer maxneu,maxlay,maxtypes,maxtpar integer maxnin,maxnout
integer maxpats
integer maxnin,maxnout,maxpout
integer maxwei,neucap,wbcap
integer maxxrmeta,xrcap
integer iinfty
integer iout,nnunit,perfunit,fitunit
integer ec_error,ec_read,ec_dim,ec_log
integer ec_dimrd
character*2 newline
character*8 stdfmt
character*8 nnldir
character*8 nntag
character*16 prim_tag
character*16 nnfdir,nnsdir
character*16 sline,asline,hline
character*16 mform,smform,miform
character*16 lrfmt,lifmt
character*32 nndmpfile,nnexpfile
character*32 nndatfile,nnreffile
character*32 sampfile,perfile
character*32 nnparfile,nnp10file
!********************************************************** !**********************************************************
!**** Parameters !**** Parameters
!*** maxneu: max. number of neurons per hidden layer
!*** maxnin: max. number of neurons in input layer !*** maxnin: max. number of neurons in input layer
!*** maxnout: max. number of neurons in output layer !*** maxnout: max. number of neurons in output layer
!*** maxpout: max. number of values in output pattern
!*** maxlay: max. number of layers (always >2)
!*** maxtypes: max. number of neuron types
!*** maxtpar: max. number of parameters for each neuron type
!*** maxpats: max. number of learning patterns
!*** maxxrmeta: max. number of metadata-blocks in xranges
!*** WARNING: maxnout should not be > maxneu, deriv-like structures parameter (maxnin=14,maxnout=15)
!*** assume so.
parameter (maxneu=150,maxnin=14,maxnout=15)
parameter (maxpout=15)
parameter (maxlay=3,maxtypes=2,maxtpar=1)
parameter (maxpats=50000)
parameter (maxxrmeta=3)
!**********************************************************
!**** Inferred Parameters
!*** maxwei: max. total number of weight matrix elements
!*** neucap: max. total number of neurons
!*** wbcap: max. total number of weights and biases
!*** xrcap: max. total number of used dimensions in xranges
parameter (maxwei=(maxlay-3)*maxneu**2+maxneu*(maxnin+maxnout))
parameter (neucap=(maxlay-2)*maxneu+maxnin+maxnout)
parameter (wbcap=maxwei+neucap)
parameter (xrcap=2+maxxrmeta)
!*** WARNING: maxwei may fail for 2-layered networks
!*** if maxnin*maxnout is sufficiently large!
!********************************************************** !**********************************************************
!**** Numerical Parameters !**** Numerical Parameters
!*** infty: largest possible double precision real value. !*** infty: largest possible double precision real value.
!*** iinfty: largest possible integer value. !*** iinfty: largest possible integer value.
!*** zero: sets what is considered an irrelevant difference
!*** in size. use for comarison of reals, to determine
!*** 'dangerously small' values, etc
!*** scan_res: maximum precision for geometric boundary algorithm
! 3.14159265358979323846264338327950... ! 3.14159265358979323846264338327950...
parameter (pi=3.1415926536D0) parameter (pi=3.1415926536D0)
parameter (infty=huge(1.0D0),iinfty=huge(1))
parameter (zero=1.0D-8,scan_res=1.0D-8)
!********************************************************** !**********************************************************
!**** Unit Conversion Parameters !**** Unit Conversion Parameters
!*** X2Y: convert from X to Y. !*** X2Y: convert from X to Y.
!*** !***
!**** !? currently inexact. FIX THIS.
!*** hart: hartree !*** hart: hartree
!*** eV: electron volt !*** eV: electron volt
!*** icm: inverse centimeters (h*c/cm) !*** icm: inverse centimeters (h*c/cm)
@ -99,90 +41,3 @@
parameter (icm2eV=1.0d0/eV2icm) parameter (icm2eV=1.0d0/eV2icm)
parameter (deg2rad=pi/180.0d0) parameter (deg2rad=pi/180.0d0)
parameter (rad2deg=1.0d0/deg2rad) parameter (rad2deg=1.0d0/deg2rad)
!**********************************************************
!**** I/O Parameters
!*** iout: standard output for vranf error messages
!*** nnunit: temporary UNIT for misc. output files
!*** nnuit + [0..99] are reserved for futher
!*** unspecific misc. files.
!*** perfunit: UNIT for performance logfile
!*** fitunit: UNIT added to random positive integer
!*** identifying a single core fit UNIQUELY
!***
!*** lrfmt: format for long real output
!*** lifmt: format for long integer output
!***
!*** nndatfile: filename for DATA-files
!*** (without file extension)
!*** nnreffile: filename for reference DATA-blocks
!*** (without file extension)
!*** nnparfile: filename for best fitted parameters to be
!*** written on (without file extension)
!*** nnp10file: filename for the 10th percentile parameters to
!*** be written on (without file extension)
!*** nnexpfile: filename for modified neural network parameters
!*** (without file extension)
!*** sampfile: filename for displaying sampled points in
!*** configuration space
!*** nndmpfile: filename for dumping data point pairs
!*** perfile: filename for logged fitting performances.
!*** nntag: infix for various filenames to mark their origin
!*** program should end with a trailing '_' if nonempty.
!*** prim_tag: tag added to the '***' line of primitive par-files
!*** nnfdir: directory for dumping fit files
!*** nnsdir: directory for dumping scans.
!*** nnldir: directory for dumping logfiles for each fit
parameter (nndatfile='DATA_ANN')
parameter (nnreffile='REF_ANN')
parameter (nnparfile='../nnfits/fit_pars')
parameter (nnp10file='../nnfits/fit_10p')
parameter (nnexpfile='../nnfits/exp_pars')
parameter (nndmpfile='../nnfits/fit_dump.dat')
parameter (sampfile='../scans/samples.dat')
parameter (perfile='../logs/performance.log')
parameter (nnfdir='../nnfits/',nnsdir='../scans/')
parameter (nnldir='../logs/')
parameter (nntag='',prim_tag=' Time-stamp: " "')
parameter (lrfmt='(ES20.12)',lifmt='(I)')
parameter (iout=6,perfunit=700,nnunit=800,fitunit=8000)
!**********************************************************
!**** Debugging Parameters
!*** sline: separation line
!*** asline: alternative sep. line
!*** hline: simple horizontal line
!*** newline: a single blank line
!*** mform: standard form for matrix output
!*** miform: standard form for integer matrix output
!*** smform: shortened form for matrix output
!*** stdfmt: standard format for strings
parameter (sline='(75("*"))',asline='(75("#"))')
parameter (hline='(75("-"))')
parameter (newline='()')
parameter (mform='(5ES12.4)',smform='(5ES10.2)')
parameter (miform='(5I12)')
parameter (stdfmt='(A)')
!**********************************************************
!**** Error Codes
!*** Codes should be powers of 2. Binary representation of return value
!*** should correspond to all exceptions invoked. ec_error should never
!*** be invoked with any other.
!***
!*** ec_error: generic error (catch-all, avoid!)
!*** ec_read: parsing error during les()
!*** ec_dim: dimensioning error
!*** ec_log: logic error
!***
!**** Inferred error codes
!*** ec_dimrd: ec_dim+ec_read
parameter (ec_error=1,ec_read=2,ec_dim=4,ec_log=8)
parameter (ec_dimrd=ec_dim+ec_read)