[InputMesh]
  TriangleType = PlateElementMappedToShell
  ClampedBoundaryCondition = 1
  chartXAType = PlateToNonIsometricTest

  ShellType = 1
;;  ShellType = 8
  file = ../../../../data/VTK/unitsquare/unitSquareSubdiv4CrissCross.vtk
;  file = ../../../../data/VTK/lshape/Lshape_subdiv3.vtk
  tangentSpaceType = 2
  parserFileVTKPlotChart = ../../../../parameters/dkt/examplePlate/vtkPlotInterface_Plate_Undeformed.ini
  parserFileVTKPlotUndeformed = ../../../../parameters/dkt/exampleCylinder/vtkPlotInterface_Cylinder_Undeformed.ini
  parserFileVTKPlotDeformed = ../../../../parameters/dkt/exampleCylinder/vtkPlotInterface_Cylinder.ini
  parserFileVTKPlotStress = ../../../../parameters/dkt/examplePlate/vtkPlotInterface_Plate_Stress.ini

  [Force]
  numLoads = 1
  Type1 = 1
;;  Load1 = -2.5,0.0,0.0
;Load1 = 0.0,0.0,-2.0
;;Load1 = 0.0,0.0,2.5
  Load1 = 0.0,0.0,-1.0
;;  Load1 = 0.0,0.0,-4.0
  StartRange1 = -100.,-100.,-100.
;;StartRange1 = 0.5,0.5,-100.
  EndRange1 = 100.,100.,100.

;  Type1 = 1
;  Load1 = -8.0,0.0,0.0
;  StartRange1 = -100.,-100.,-100.
;  EndRange1 = 0.5,0.5,100.

;  Type2 = 1
;  Load2 = -5.0,0.0,0.0
;  StartRange2 = -100.,-100.,-100.
;  EndRange2 = 0.5,0.5,100.

  scaleForceWithAreaOfHardMaterial = 1
  scaleForceWithThickness = 0
  scaleForceWithThicknessExponent = 1.

  factor_force = 1.


;   [InputMesh]
;   TriangleType = PlateElementMappedToShell
;   ClampedBoundaryCondition = 1
;   chartXAType = PlateToHalfCone

;   ShellType = 1
;   file = ../../../../data/VTK/unitSquare/unitSquareSubdiv4.vtk
;   tangentSpaceType = 2
;   parserFileVTKPlotChart = ../../../../parameters/dkt/examplePlate/vtkPlotInterface_Plate_Undeformed.ini
;   parserFileVTKPlotUndeformed = ../../../../parameters/dkt/exampleCylinder/vtkPlotInterface_Cylinder_Undeformed.ini
;   parserFileVTKPlotDeformed = ../../../../parameters/dkt/exampleCylinder/vtkPlotInterface_Cylinder.ini
;   parserFileVTKPlotStress = ../../../../parameters/dkt/examplePlate/vtkPlotInterface_Plate_Stress.ini

;   [Force]
;    numLoads = 1
;    Type1 = 1
; ;   Load1 = -100.,0.0,0.0
;    Load1 = 0.0,0.0,-100.0
;    StartRange1 = -100,-100,-100.
;    EndRange1 = 100.,100.,100.

;    scaleForceWithAreaOfHardMaterial = 1
;    scaleForceWithThickness = 0
;    scaleForceWithThicknessExponent = 1.

;    factor_force = 1.



;  [InputMesh]
;  TriangleType = PlateElementMappedToShell
;  ClampedBoundaryCondition = 1
;  chartXAType = PlateToWobbly

;  ShellType = 1
;;  ShellType = 1111
;  file = ../../../../data/VTK/unitsquare/unitSquareSubdiv4.vtk
;;  file = ../../../../data/VTK/lshape/lshape_subdiv4.vtk
;  tangentSpaceType = 2
;  parserFileVTKPlotChart = ../../../../parameters/dkt/examplePlate/vtkPlotInterface_Plate_Undeformed.ini
;  parserFileVTKPlotUndeformed = ../../../../parameters/dkt/exampleCylinder/vtkPlotInterface_Cylinder_Undeformed.ini
;  parserFileVTKPlotDeformed = ../../../../parameters/dkt/exampleCylinder/vtkPlotInterface_Cylinder.ini
;  parserFileVTKPlotStress = ../../../../parameters/dkt/examplePlate/vtkPlotInterface_Plate_Stress.ini

;  [Force]
;  numLoads = 1
;  Type1 = 1
;  Load1 = 0.0,0.0,-2.5
;;  Load1 = 0.0,0.0,-1.0
;  StartRange1 = -100,-100,-100.
;  EndRange1 = 100.,100.,100.

;  scaleForceWithAreaOfHardMaterial = 1
;  scaleForceWithThickness = 0
;  scaleForceWithThicknessExponent = 1.

;  factor_force = 1.




[ConstraintProblemAdaptiveProlongationTypes]
 ProlongationMeshType = 1
 ProlongationPfType = 1
 ProlongationDispType = 1
 ProlongationMeshBoundaryType = 1
 constructProlongatedBoundaryByShellHandler = 0


;########################
[Material]
; #------------------------------------------------------------------------
; #####  elastic energy parameter  #####
; #------------------------------------------------------------------------
 MaterialFromDictionary = 0
# elastic materials:
 #if = 0
; ElastModulus_hard = 1.e+2
; ElastModulus_hard = 1.e+1
 ElastModulus_hard = 1.e+0
 ElastModulus_soft = 1.e+0
 PoissonRatio_hard = 0.
 PoissonRatio_soft = 0.
 thickness_hard = 1.
 thickness_soft = 1.

   factor_membraneEnergy = 0.
; factor_membraneEnergy = 1.e+3
 factor_bendingEnergy = 1.e-0

 #chi(v) = ...:  -1 - harmonic, 1,2,4 - polyonmial
 pfFunctionMaterialType = 4
 #double well: 2,4 - polyonmial
 pfFunctionDoubleWellType = 4

[InitMaterial]
 #initMaterialType:
 # 1 - boundary layer
 # 2 - orthogonal layer
 # 8 - square at boundary

 designType = 1
; areaHardMaterial = 0.25
 areaHardMaterial = 1.

;########################
[saving]
 saveDirectory = ../../../../results/resultsDKT

 saveInitialization = 1
 saveSolution = 1
 plotResults = 0
 removeOldRefinementSteps = 0

 plotPointCloud = 0
 ;VTKFileType = vtk
 VTKFileType = vtp
 LATEXType = pdflatex
;  LATEXType = lualatex

 plotInterface = 1
 plotStress = 1
 thresholdInterface = 0.5
 thresholdBdr = 1.
 thresholdHardRegion = 0.5
 thresholdSoftRegion = 0.5

;########################
[DerivativeTest]
# order: 0 - no, 1 - FirstDerivativeTest, 2 - SecondDerivativeTest
 order = 0
 stepSize = 1.e-7
 tolerance = 1.e-5

 ;########################
[DerivativeTestConstraint]
# order: 0 - no, 1 - FirstDerivativeTest, 2 - SecondDerivativeTest, 12 - both
 order = 0
 stepSize = 1.e-7
 tolerance = 1.e-5


;########################
[ConstraintProblem]

; #------------------------------------------------------------------------
; #####   constraint problem (linelast )  #####
; #------------------------------------------------------------------------
 solveWithDirectSolver = 1
 #if iterative solver:
 #tolerance = sqrt( tol^2 * systemmat.col() )
 toleranceLinearSystem = 1.e-16
 #maxIters = maxItersFac * systemmat.cols()
 maxItersFacLinearSystem = 0.25

 factorIsoConstraint = 1.

 # linear:    1 - Plate Laplace, 3 - KirchhoffLove,
 # nonlinear: 11 -NonlinMembraneHinge, 21 - NonlinearMembrane linear Plate Laplace bending, 23 - NonlinearMembrane linear Plate Laplace bending, 31 -Nonlinear Membrane and Bending
 ElastEnergy = 31

#nonlinear solver: NEWTON or IPOPT or FirstTryNEWTONelseIPOPT
;  nonlinearSolver = NEWTON
 nonlinearSolver = IPOPT
;  nonlinearSolver = FirstTryNEWTONelseIPOPT
 ############### for ipopt
 # 0: MUMPS, x HSL_x x = 27,57,77,86,97
 linearSolverTypeIpopt = 97
 maxIterationsIpopt = 200
 breakConditionIpopt = 1.e-10
;   breakConditionIpopt = 1.e-12
 boundRelaxFactorIpopt = 0.
 outputLevelIpopt = 0
 ############### for newton
 # direct solver: 1 Umfpack, 2 - Cholmod, 4 - Umfpack_NEW(with Suitesparse_long)
 # iterative solver: 11 - ConjugateGradient, 24 BICGSTAB, 34 - GMRES
 linearSolver = 5
 maxIterations = 200
 breakCondition = 1.e-14

 #set initialization of displacement and multiplier to solution for current material
 updateInitialization = 1

 # 0 - no output , 1 - output if Newton fails, 2 - , 3 - Newton steps, 4 - always plot deformation when solve
 #recommended: 1 for material optimization, 3 for optimal deformation
 outputLevel = 1


 # initType : 0 - zero,
 #  Plate: 1 - simple, 2 -  twist
 #  general: 101 : h * 10^-3 (1,1,1) on interior nodes
 #           102 : h * 10^-3 f on interior nodes
 initType = 0

 # MarkingType:
#  1 - all
#  2 - all twice
# 11 - L2Gradient of material
# 12 - L2Gradient and clamped elements (ie if all nodes on boundary)
# 31 - L2Hessian of disp
# 32 - L2Hessian of disp and clamped elements (ie if all nodes on boundary)
; adaptiveMarkingType = 32
adaptiveMarkingType = 2

 #for stressvec
 adaptiveMarkingUseTotalStressVec = 0
 adaptiveMarkingUseGaussCurvatureVec = 0
 adaptiveMarkingPercentil = 0.85


 #adaptive intit material:
 #  0 - solutionMaterial
 #  1 - from file 
 #  2 - with shellHandler and designtype
 # 11 - init(x,y,z) = 0.5 ( sol(x,y,z) + sol(x,1-y,z) )
 adaptiveMaterialInitializationType = 2
 firstStepMaterialInitializationType = 0

;########################
[MaterialOptimization]
 factorComplianceCost = 1.0
 factorAreaCost       = 0.0
 factorInterfaceCost  = 0.0
 factorDoubleWell     = 1.


[Compare]
 thicknessVec = 1.0
;  factorForceVec = 1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,15.,16.,17.,18.,19.,20.,21.,22.,23.,24.,25.,26.,27.,28.,29.,30.,31.,32.
 factorForceVec = 4.,8.,16.
;  factorForceVec = 1.,2.,4.,8.,16.,32.,64.,128.

;########################
[BASH]
 numThreads = 8
 renewKerberosTicket = 0