----------------------------------------------------------------------------
   
   c a u c h y - b o r n   c o n s t i t u t i v e   m o d e l
   
   This program computes the finite deformation constitutive response
   of an arbitrary multilattice single crystal.  The code is based on
   the local quasicontinuum (QC) method (http://qcmethod.org).
   
   The program uses interatomic potentials compatible with the
   Knowledgebase of Interatomic Models (KIM) project (https://openkim.org).
   
   Copyright (C) 2015 E. B. Tadmor
   
   ----------------------------------------------------------------------------
   
   Test calculation: Stillinger-Weber silicon ground state and compression

 **setting factor**
   * Default globaltol set to  0.10000E-15
   CG global tolerance = 0.10000E-15

 **reading material data**
   
   Material name 'Al-Ni(L12)'
   
   material information
   
   Bravais lattice vectors:
   a1 =    4.00000   0.00000   0.00000
   a2 =    0.00000   4.00000   0.00000
   a3 =    0.00000   0.00000   4.00000
   Basis atoms:
   Number of basis atoms =  4
   spec    mass        reference shift vector
     13    1.00       0.00000   0.00000   0.00000
     28    1.00       0.00000   2.00000   2.00000
     28    1.00       2.00000   0.00000   2.00000
     28    1.00       2.00000   2.00000   0.00000

 **reading constitutive information**
* Error: at line 169 in mod_poten_kim.F90
	Message: kim_api_string_init
	KIM_STATUS_MSG: unsuccessful completion
   
   KIM Model name: IMD_EAM_Schopf_AlNiCo_A__MO_122703700223_002
   
   Match failed; removing process_d2Edr2 and trying again.
   
   Moduli calculation disabled.
   
   Preconditioning disabled in lattice parameter optimization.
   
   Maximum cut-off distance:      7.2000000000

 **global parameter settings**
   
   
   Number of spatial dimensions =       3
   Number of dofs per node      =       3
   Number of nodes per element  =       4
   Shuffles optimization                    = ACTIVATED
   Prefer multilattice stress correction    = OFF
   Optimize F and shuffles concurrently     = ACTIVATED
   Use preconditioned CG optimization       = OFF

 **check equilibrium of multilattice**
   
   Input lattice vectors for CheckEquilibrium
        4.0000000000000000     0.0000000000000000     0.0000000000000000
        0.0000000000000000     4.0000000000000000     0.0000000000000000
        0.0000000000000000     0.0000000000000000     4.0000000000000000
   Input basis atom positions for CheckEquilibrium
        0.0000000000000000     0.0000000000000000     0.0000000000000000
        0.0000000000000000     2.0000000000000000     2.0000000000000000
        2.0000000000000000     0.0000000000000000     2.0000000000000000
        2.0000000000000000     2.0000000000000000     0.0000000000000000
   
   scaling parameter:      0.9316432643371826
   
   Updated lattice vector after uniform optimization of lattice parameter:
        3.7265730573487303     0.0000000000000000     0.0000000000000000
        0.0000000000000000     3.7265730573487303     0.0000000000000000
        0.0000000000000000     0.0000000000000000     3.7265730573487303
   Updated basis atom positions after uniform optimization of lattice parameter:
        0.0000000000000000     0.0000000000000000     0.0000000000000000
        0.0000000000051714     1.8632865287242122     1.8632865286875797
        1.8632865286631897     0.0000000000051711     1.8632865286224738
        1.8632865287242175     1.8632865286631950    -0.0000000000051714
   
   dwde for CheckEquilibrium
        0.0000000000000000    -0.0000000000688929    -0.0000000007094164
        0.0000000000000000     0.0000000009742882    -0.0000000000688937
        0.0000000000000000     0.0000000007656732    -0.0000000010306955
       -0.0000000000688929    -0.0000000007094164     0.0000000009744395
   * CB-KIM calculation successfully completed