----------------------------------------------------------------------------
   
   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**
   
   KIM Model name: EAM_Dynamo_Angelo_Moody_NiAlH__MO_418978237058_004
   
   Maximum cut-off distance:      5.6500000000

 **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       = ACTIVATED

 **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.8918999667458973
   
   Updated lattice vector after uniform optimization of lattice parameter:
        3.5675998669835893     0.0000000000000000     0.0000000000000000
        0.0000000000000000     3.5675998669835893     0.0000000000000000
        0.0000000000000000     0.0000000000000000     3.5675998669835893
   Updated basis atom positions after uniform optimization of lattice parameter:
        0.0000000000000000     0.0000000000000000     0.0000000000000000
        0.0000000000000000     1.7837999334917947     1.7837999334917947
        1.7837999334917947     0.0000000000000000     1.7837999334917947
        1.7837999334917947     1.7837999334917947     0.0000000000000000
   
   dwde for CheckEquilibrium
        0.0000000000000000     0.0000000000000001    -0.0000000000000001
        0.0000000000000000    -0.0000000000000001     0.0000000000000001
        0.0000000000000000    -0.0000000000000001    -0.0000000000000001
        0.0000000000000001    -0.0000000000000001    -0.0000000000000001
   * CB-KIM calculation successfully completed