jnshuti
/
Gen_ANN_nh3_base
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Modify the weight and set the unnormalized weight to one .i.e disabling the weighted method of David for the energy

This commit is contained in:
jean paul nshuti 2025-10-09 13:07:26 +02:00
parent 159635c494
commit c74e309497
3 changed files with 19 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
src/error.f
View File

@ -62,7 +62,7 @@
!-------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------
subroutine nnweight(wterr,pat_out) subroutine nnweight(wterr)
implicit none implicit none
! Evaluate system specific weighting for 1 pattern. ! Evaluate system specific weighting for 1 pattern.
@ -70,29 +70,30 @@
include 'nnparams.incl' include 'nnparams.incl'
include 'nncommon.incl' include 'nncommon.incl'
include 'JTmod.incl' !include 'JTmod.incl'
double precision wterr(maxpout),pat_out(maxpout) double precision wterr(maxpout)!,pat_out(maxpout)
double precision eref(nstat) !double precision eref(nstat)
double precision wten(3) !double precision wten(3)
integer j !integer j
wten=1 !wten=1
eref(1)=E0_sig !eref(1)=E0_sig
eref(2:3)=E0_pi !eref(2:3)=E0_pi
! eref(1:3)=eref(1:3)+600.d0*icm2hart ! eref(1:3)=eref(1:3)+600.d0*icm2hart
eref(1:3)=eref(1:3)+1000.d0*icm2hart !eref(1:3)=eref(1:3)+1000.d0*icm2hart
!! kill state 2 and 3 !! kill state 2 and 3
! wten(2:3)=0 ! wten(2:3)=0
do j=1,3 !do j=1,3
! weighting of energies ! weighting of energies
wterr(j)=wdamp(pat_out(j)-eref(j))*wterr(j)*wten(j) ! wterr(j)=wdamp(pat_out(j)-eref(j))*wterr(j)*wten(j)
enddo !enddo
wterr = 1.0d0
contains contains
double precision function wdamp(dE) double precision function wdamp(dE)

6
src/mkNN.f
View File

@ -408,9 +408,11 @@
write(nnunit,newline) write(nnunit,newline)
do k=1,npat do k=1,npat
call nnweight(wterr(1,k),pat_out(1,k)) !call nnweight(wterr(1,k),pat_out(1,k)) ! JP no need for the weighting schemme used for energy
! right now wter is set to 1.0 in nnweights
call nnweight(wterr(1,k))
enddo enddo
pat_out = pat_out ! to avoid unsed complain durinng compilation
total=0 total=0
do k=1,sets do k=1,sets
write(nnunit,'(A10,I6.5)') '# Scan Nr.',k write(nnunit,'(A10,I6.5)') '# Scan Nr.',k

2
src/nnparams.incl
View File

@ -48,7 +48,7 @@
!*** WARNING: maxnout should not be > maxneu, deriv-like structures !*** WARNING: maxnout should not be > maxneu, deriv-like structures
!*** assume so. !*** assume so.
parameter (maxneu=150,maxnin=15,maxnout=30) parameter (maxneu=150,maxnin=15,maxnout=30)
parameter (maxpout=15,maxset=10000) parameter (maxpout=20,maxset=10000)
parameter (maxlay=3,maxtypes=2,maxtpar=1) parameter (maxlay=3,maxtypes=2,maxtpar=1)
parameter (maxpats=50000) parameter (maxpats=50000)