$NON_LINEAR_SOLVER
------------------------------



--> OLD FORMATS
; method -- local tolerance -- global tolerance -- max_iterations -- relaxation
  NEWTON    1.e-3                1.e-6             20                0.0
  
; method -- tolerance -- max_iterations -- relaxation
  PICARD    1.e-3        20                0.0
  
!Problems with this:
 1) No selection of error method.
 2) Tolerance only allowed for 1 degree of freedom (this is the vital change)
 
 --> NEW FORMATS
; method -- error_method -- max_iterations -- relaxation -- local_tolerance -- tolerance[1:dof]
  NEWTON    1               20                0.0           1.e-3              1.e-6  --> number of dof's
  
; method -- error_method -- max_iterations -- relaxation -- tolerance[1:dof]
  PICARD    1               20                0.0           1.e-6  --> number of dof's

  
FOR EXAMPLE, A PROCESS WITH 3 PRIMARY VARIABLES (Pw, Snw, and T, for instance) WOULD BE:

; method -- error_method -- max_iterations -- relaxation -- local_tolerance -- tolerance_Pw -- tolerance_Snw -- tolerance_T
  NEWTON    1               20                0.0           1.e-3              1.0e-2          1.0e-7           1.0e-5
  
; method -- error_method -- max_iterations -- relaxation -- tolerance_Pw -- tolerance_Snw -- tolerance_T
  PICARD    1               20                0.0           1.0e-2          1.0e-7           1.0e-5
  

  
---------------------------------------
CHANGES TO COUPLING SCHEME
---------------------------------------
  
here is the input file for a 3 process inner and outer coupling. Ask me if you have questions, it is easiest to explain by providing the input file.


GeoSys-NUM: Numerical Parameter ----------------------------------------
$OVERALL_COUPLING
; min_iter -- max_iter
  1           2

#NUMERICS
 $PCS_TYPE
  MULTI_PHASE_FLOW
 $LINEAR_SOLVER 
; method error_tolerance max_iterations theta precond storage
  2      2 1.e-10       2000           1.0   100       4
 $NON_LINEAR_SOLVER 
; method - error method - max iter - relax - tol_Pc - tol_Pnw
  PICARD   1              30         0.0     1.0e-4   1.0e-4
 $COUPLING_CONTROL
; error_method -- tol_Pc -- tol_Pnw
  1               1.0e-3    1.0e-3
 $COUPLED_PROCESS
; process    --    min_iter -- max_iter
  HEAT_TRANSPORT   1           3
 $ELE_MASS_LUMPING
  0
 $ELE_GAUSS_POINTS
  2  
  
#NUMERICS
 $PCS_TYPE
  HEAT_TRANSPORT
 $LINEAR_SOLVER 
; method error_tolerance max_iterations theta precond storage
  2      1 1.e-012       1000           0.5   100       4
 $NON_LINEAR_SOLVER 
; method - error method - max iter - relax - tol_T
  PICARD   1              100        0.0     1.0e-6
 $COUPLING_CONTROL
; error_method -- tol_T
  1               1.0e-6 
   
#NUMERICS
 $PCS_TYPE
  DEFORMATION
 $NON_LINEAR_SOLVER 
; method - error method - max iter - relax - local_tol -- tol_u
  NEWTON   1              100        0.0     1.0e-2       1.0e-10
 $LINEAR_SOLVER 
; method error_tolerance max_iterations theta precond storage
  2      0 1.e-011       2000         1.0   100     4
#STOP