aktueller Status before fixing cell evaluation

added cell_values option and evaluation getInvariant and UpdateStretch in cell and added EvaluationTractionStiffness

conf/coupled/gamm-paper/coarse-biventricle.conf

+activeDeformation = 1
+WithPrestress = true
+PrestressSteps = 1
+withPressure = 1
+
+pressureED_LV = 0.001
+pressureED_RV = 0.001
+pressureES_RV = 0.001
+pressureES_RV = 0.001
+
+WithPrestress = false
+
+MechDeltaTime  = 0.001 #0.004# 0.001
+ElphyDeltaTime = 0.001 # 0.001 #0.000001
+
+#StretchFactorSheet = 2.0
+#StretchFactorNormal = 5.0
+
+linearImplicitSolverCell = 1
+cell_values = 1
+deformConcentrationStretch = 1 #1
+withPressure = 0
+timeStepExtrapolation = 0
+DampExtrapolation = 0.5;
+Mesh = KovachevaBiventricle
+#Mesh = KovachevaBiventricleInstantActivation
+ClearData = true
+
 Model = Coupled
 CoupledModel = Segregated
 
-CoupledProblem = BiventricleCoarse
-#CoupledProblem = BiventricleDirichlet
+CoupledProblem = DeformedBiventricleCoarse
+#CoupledProblem = BiventricleCoarse
 
 ElphyModelClass = MElphyModel
 ElphyModelName = TenTusscher
 TensionModelName = Rossi
-RossiViscosity = 11.0
-RossiForce = 6.0
+RossiViscosity = 7000
+RossiForce = 7
 
-ElphyLevel = 1
-MechLevel = 1
+ElphyLevel = 0 #1
+MechLevel = 0 #1
 MechPolynomialDegree = 1
 ElphyPolynomialDegree = 1
 
@@ -22,51 +50,54 @@ StartTime = 0.0
 
 #EndTime=0.00001
 EndTime = 0.8
+#EndTime = 0.8
 #EndTime = 0.02
 
-MechDeltaTime = 0.004
-ElphyDeltaTime = 0.0001
+MechDeltaTime  = 0.001 #0.004# 0.001
+ElphyDeltaTime = 0.0001 # 0.001 #0.000001
 
-DeltaTime = 0.0001
+DeltaTime = 0.001
 DeltaTimeMin = 0.00001
 DeltaTimeMax = 0.001
 
 ######################
 #      Plotting      #
 ######################
-PlottingSteps = 1
+PlottingSteps = 10
 
 #######################
 #     Coupled Data    #
 #######################
 ClearDistribution = 0
 
-
 MechDiscretization = Conforming
 MechRuntype = Default
-MechDynamics = Newmark
-#MechDynamics = RayleighNewmark
+#MechDynamics = Newmark
+MechDynamics = RayleighNewmark
+RayleighAlpha = 100
 #MechDynamics = ViscoNewmark
-
+#ViscoelasticDamping = KelvinVoigt
+#KelvinVoigtEta = 2
 #------------------------
 # Newmark Parameters
 #------------------------
 # Damped Newmark
-NewmarkBeta = 0.4
-NewmarkGamma = 0.7
+#NewmarkBeta = 0.3
+#NewmarkGamma = 0.6
 # Undamped Newmark
-#NewmarkBeta = 0.25
-#NewmarkGamma = 0.5
+NewmarkBeta = 0.25
+NewmarkGamma = 0.5
 
 
 
 # Amount of pressure steps per iteration
-PressureSteps = 1
+PressureSteps = 0
 # If Iteration with Pressuresteps fails, recursively restart PressureSolver with Steps*2 until iteration limit is reached
 PressureIterations = 10
 PressureDepth = 0
 
 WithPrestress=false
+#WithPrestress=true
 PrestressSteps=10
 
 ######################
@@ -75,16 +106,11 @@ PrestressSteps=10
 MechMeshPart = with_vessels
 MechMeshPart = vessel_base
 
-StartTime = 0.0
-EndTime = 0.4
-DeltaTime = 0.001
-DeltaTimeMin = 0.0001
-DeltaTimeMax = 0.001
-
-
-TractionK = 0.0#0.1
-TractionC = 0.0#1e-3
+#TractionK = 0
+#TractionC = 0
 
+TractionK = 0#0.1
+TractionC = 0#5e-3
 
 #------------------------
 # Stretch Parameters
@@ -96,8 +122,7 @@ StretchFactorNormal = 4.0
 ######################
 # Material Parameters#
 ######################
-Density = 0.1082
-Density = 1.082
+Density = 1082e-9
 #------------------------
 # Volumetric Penalty
 # Incompressible should be true or false
@@ -106,15 +131,13 @@ Density = 1.082
 VolumetricSplit = 0
 Incompressible = false
 QuasiCompressiblePenalty = Ciarlet #None, Ciarlet, Mixed
-VolumetricPenalty = 20
-PermeabilityPenalty = 20
+VolumetricPenalty = 1000 # ?
+PermeabilityPenalty = 1000 # ?
 
-DGSign=-1.0
-DGPenalty=20.0
 
 ElphyModel = SemiImplicit
-OrderStaggeredScheme=wcV
-ReassembleRHSonce=1
+OrderStaggeredScheme = wcV
+ReassembleRHSonce = 0
 
 CrankNicolsonTheta = 0.5
 ElphySplittingMethod =Godunov  #Strang, Stein,Godunov
@@ -122,7 +145,7 @@ ElphySplittingMethod =Godunov  #Strang, Stein,Godunov
 #Iext splitted in PDE =1, Iext in ODE =0
 IextInPDE=1
 #Glaettung von Anregung  in Zeit: Stepfunction,Arctan,Tanh
-ExternalCurrentSmoothingInTime = Arctan
+ExternalCurrentSmoothingInTime = Stepfunction#Arctan #Stepfunction
 #Glättung Anregung im Raum: discrete, linear
 ExternalCurrentSmoothingInSpace = linear
 
@@ -133,14 +156,14 @@ ElphyMeshPart = ventriclesonly
 SurfaceToVolumeRatio= 140#1/mm
 MembraneCapacitance=0.01#muF/mm^2
 
-IntraLongitudinal = 0.170
+IntraLongitudinal = 0.170 #0
 ExtraLongitudinal = 0.620
 IntraTransversal = 0.019
-ExtraTransversal = 0.240
+ExtraTransversal = 0.240  #0
 
 ElphyIntegrator=ExplicitEuler
 CellModelScheme = ExponentialIntegrator
-IonScheme = RungeKutta2 # ExplicitEuler,RungeKutta2,RungeKutta4
+IonScheme = Quadratic # ExplicitEuler,RungeKutta2,RungeKutta4
 
 ThresholdVForActivation=-40.0
 ScaleExcitationAmplitude=0.04101
@@ -155,6 +178,7 @@ LinearEpsilon = 1e-8
 LinearReduction = 1e-8
 LinearReduction = 1e-11
 LinearSteps = 10000
+#LinearPrintSteps = 1000
 BasePreconditioner = GaussSeidel
 Preconditioner = GaussSeidel
 Preconditioner = LIB_PS
@@ -169,22 +193,22 @@ GatingReduction = 1e-12
 # === Elphy Solver ================================================
 ElphySolver = gmres
 ElphyPreconditioner=Jacobi
-ElphyEpsilon = 1e-8
-ElphyReduction = 1e-12
+ElphyEpsilon = 1e-6
+ElphyReduction = 1e-10
 
 # === Mechanics Solver ============================================
 MechSolver = gmres
 MechPreconditioner = GaussSeidel
-MechEpsilon = 1e-8
-MechReduction = 1e-12
+MechEpsilon = 1e-10
+MechReduction = 1e-10 #1e-12
 MechSteps = 20000
 
 # === Newton Method === #
-NewtonEpsilon = 1e-6
-NewtonReduction = 1e-12
+NewtonEpsilon = 1e-8 #1e-8
+NewtonReduction = 1e-12 #1e-5
 NewtonSteps = 300
 NewtonLineSearchSteps = 10
-NewtonLinearizationReduction = 1e-12
+NewtonLinearizationReduction = 1e-4 #1e-12
 
 #Smoother = GaussSeidel
 #Smoother = PointBlockGaussSeidel
@@ -192,35 +216,35 @@ NewtonLinearizationReduction = 1e-12
 presmoothing = 3
 postsmoothing = 3
 
-ConfigVerbose = 1
+ConfigVerbose = 10
 AssembleVerbose = 1
 
 # = Solver Verbose =
 MultigridVerbose = 0
 BaseSolverVerbose = -1
-LinearVerbose = -5
-NewtonVerbose = 0
+LinearVerbose = -1
+NewtonVerbose = 1
 EulerVerbose = 0
-LinearImplicitVerbose = -5
-ElphyLinearSolverVerbose =-1
+LinearImplicitVerbose = 1
+ElphyLinearSolverVerbose =1
 MonodomainVerbose =-1
 NonLinearVerbose = -5
-ElphySolverVerbose = 0
+ElphySolverVerbose = 1
 CoupledSolverVerbose = 4
 CellModelVerbose = 0
 
 # = Main Verbose =
 ElphyVerbose = 1
 MechVerbose = 1
-
+PrintVSteps = 10
 # = Plot Verbose =
-ElphyVTK = 0 # 0 Disables all plots from Electrophysiology
+ElphyVTK = 1 # 0 Disables all plots from Electrophysiology
 MechVTK = 1 # 0 Disables all plots from Elasticity
 
 # = Plot Verbose =
 PressureSolverVTK = 0
-DynamicSolverVTK = 1
-ElphySolverVTK = 0
-CoupledSolverVTK = 1
-
-PlotVTK = 1
+DynamicSolverVTK = 10
+ElphySolverVTK = 10
+CoupledSolverVTK = 10
+PlotCalcium=1
+PlotVTK = 1
\ No newline at end of file