!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!                                              !!!!!
!!!!!  VERIFICATION CHECK: vc-dimer-continuity-c1  !!!!!
!!!!!                                              !!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Description: Determines whether a model has a continuous energy and first
             derivative, i.e. belongs to the C^1 continuity class, for all
             possible dimers. For a model supporting N species, there are N +
             N!/(2(N-2)! distinct dimers for all possible species combinations.
             (For example if N=3, there are 3+3!/2=6 dimers.  If the species are
             A,B,C, the 6 dimers are AA, BB, CC, AB, AC, BC.) For each dimer,
             the equilibrium separation and cutoff are determined. The
             continuity across the cutoff is assessed.  Then an analysis is
             performed to detect any discontinuities from half the equilibrium
             distance to the cutoff. Although the verification check only
             requires C^1 continuity to pass, continuity up to 3rd order is
             checked and reported.

Author: Ellad Tadmor

-----------------------------------------------------------------------------------------------------
Results for KIM Model      : Morse_Shifted_GirifalcoWeizer_1959HighCutoff_Ca__MO_159753408472_003
Supported species          : Ca

-----------------------------------------------------------------------------------------------------
DIMER Ca--Ca
-----------------------------------------------------------------------------------------------------

E n e r g y  E x t r e m a
                                   distance            energy       #iter   #fn-calls    warnings
Equilibrium separation:      4.56900001e+00   -1.62283770e-01          34          69    none
Cutoff separation:           1.68661000e+01   -4.43892732e-14          28          30    none

C u t o f f  S m o o t h n e s s

                       value   continuous
Energy       -4.43892732e-14          yes
1st-Deriv     1.30704131e-05     ** NO **
2nd-Deriv    -1.05257306e-05     ** NO **
3rd-Deriv     8.47604808e-06     ** NO **

C o n t i n u i t y

Checking continuity for r = [2.28450,16.86610] at 1459 points (Delta r = 0.01000)

Local edge detection based on a normalized 5th-order local difference formula T^5
is used to determine the presence of discontinuities. The tolerance is |T^5|>1.00000.
(For details see Anne Gelb and Eitan Tadmor, J. Sci. Comp., 28:279-306, 2006.)


Checking Energy (energy vs. distance in file "dimer-energy-CaCa.dat")
... No discontinuities found.

Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-CaCa.dat")
... No discontinuities found.

Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-CaCa.dat")
... No discontinuities found.

Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-CaCa.dat")
... No discontinuities found.

=====================================================================================================

SUMMARY of Model Continuity Results Across All Dimers:

             continuous
-----------------------------
Energy              yes
1st-Deriv      ** NO **
2nd-Deriv      ** NO **
3rd-Deriv      ** NO **

=====================================================================================================
Continuity must be C^1 or higher to pass this verification check.

Grade: F

Comment: The model is C^0 continuous. This means that the model has continuous energy, but a discontinuous first derivative.