# 
ChampGeoLin    "ChampGeo" ""

booleen Reassemble f
booleen VerboseInstationnaire  v
booleen VerbosePointFixe v
booleen E_InfVieillissant f
booleen TempsImpose f

scalaire TempsInitial 0
scalaire LongueurPasDeTemps .1
scalaire NombreDePasDeTemps 420
scalaire LongueurMaximalePasDeTemps 0.1

scalaire TempsFinal   41
scalaire "IndicateurEtat_Inf" 0
#scalaire "Integrale_bInf" 1

v3dquad  "U" [0,0,0]
v3dquad  "UPrec" [0,0,0]

to4elasticite "ERef_Infijkl"  [1e5,1e5,1e5,0.3,0.3,0.3,38462,38462,38462]
to4elasticite "E_Infijkl"  [1e5,1e5,1e5,0.3,0.3,0.3,38462,38462,38462]

scalaire "NombreBrancheViscoElastique" 1
{ Branche_00
  to4elasticite "ERef_ijkl"     [4e5,4e5,4e5,0.3,0.3,0.3,153846,153846,153846]
  empilementto4 "FonctElasViscoRef_ijkl" [Un,Zero,Un, Un,Zero,Un, Zero,Zero,Zero,Un,Zero,Zero,Zero,Zero,Un,Un,Zero,Un,Zero,Zero,Un]
  #1 0 1 1 0 1 0 0 0 1 0 0 0 0 1 1 0 1 0 0 1 
  booleen MateriauVieillissant f
  #booleen MateriauVieillissant v
  scalaire "Multiplicateur_b" 1
  scalaire "Multiplicateur_bPrecedent" 1
  scalaire "DeriveeMultiplicateur_b" 0
  scalaire "DeriveeMultiplicateur_bPrecedent" 0
  scalaire "MultiplicateurLambda_l" 1
  scalaire "MultiplicateurLambda_lPrecedent" 1
  #scalaire "IndicateurEtat" f(M,MP)=(M-MP)
  #scalaire "Multiplicateur_b" f(One,Alpha,M,M12)=One+Alpha*(M-M12)
  #scalaire "Multiplicateur_bPrecedent" f(One,Alpha,MP,M12)=One+Alpha*(MP-M12)
}

# ceci donne des pas de temps 
scalaire NbPasDeTemps    500
scalaire freqCalcul      1
scalaire freqImpression  1

#
# valeur specifique au solveurInstNlinPETSc et solveurLinPETSc
#
# Tolerance accrue pour le solveur
scalaire ToleranceLin 1e-12
scalaire ToleranceLinMec 1e-12
scalaire ToleranceLinMecR 1e-40
scalaire Divergence 1e+10
scalaire TolerancePointFixe 1e-8
scalaire NombreMaximalIterationsPointFixe 400

# nombre maximum d'iterations pour les solveur
scalaire MaxItLinMec 6000

#
# on veut la trace de la solution en certains points
# NbEval nombre de points d'evaluations
# Coord_i le point i pour i=0 a NbEval-1
#
# en mettant NbEval à zéro on ne fait plus les evaluations


scalaire NbEval 1
#ajouter les autres coins
vectoriel3d Coord_0 [20, 0, 0]       
vectoriel3d Coord_1 [20, .5, 0]    
vectoriel3d Coord_2 [20.,-.5, 0.0] 
vectoriel3d Coord_3 [1., 1., 0] 

# On choisit le type de solveur qu'on utilise

#solveurlinpetsc      typesolveur      "SolveurLineaire" gradient_conjugue
solveurlinpetsc      typesolveur      "SolveurLineaire" direct
#solveurlinpetsc      typeprecond      "SolveurLineaire" sor
#solveurlinpetsc      configprecondsor "SolveurLineaire" [v,0.7]
solveurlinpetsc      suiviconvergence "SolveurLineaire" defaut
solveurlinpetsc      tolerance        "SolveurLineaire" [ToleranceLin,ToleranceLin,Divergence,MaxItLinMec]

scalaire "H1" f()=-1*((0<=t)&(t<=10))

scalaire "nul" f()=0


scalaire "Theta" 1
scalaire "Coef" 1