##
# c: 05.02.2008
filename_mesh = 'database/simple.mesh'

from sfepy.solvers.ts import TimeStepper
t0 = 0.0
t1 = 0.1
n_step = 11
ts = TimeStepper( t0, t1, None, n_step )

material_2 = {
    'name' : 'coef',
    'mode' : 'here',
    'region' : 'Omega',
    'val' : 0.01,
    'kind' : 'stationary', # 'stationary' or 'time-dependent'
}

field_1 = {
    'name' : 'temperature',
    'dim' : (1,1),
    'domain' : 'Omega',
    'bases' : {'Omega' : '3_4_P1'}
}

variable_1 = {
    'name' : 'T',
    'kind' : 'unknown field',
    'field' : 'temperature',
    'order' : 0,
    'history' : 'previous', # 'previous' or 'full'
}
variable_2 = {
    'name' : 's',
    'kind' : 'test field',
    'field' : 'temperature',
    'dual' : 'T',
}

region_1000 = {
    'name' : 'Omega',
    'select' : 'all',
}
region_03 = {
    'name' : 'Gamma_Left',
    'select' : 'nodes in (x < 0.00001)',
}
region_4 = {
    'name' : 'Gamma_Right',
    'select' : 'nodes in (x > 0.099999)',
}

ebc_1 = {
    'name' : 'T1',
    'region' : 'Gamma_Left',
    'dofs' : {'T.0' : 2.0},
    'value' : 2.0,
}
ebc_2 = {
    'name' : 'T2',
    'region' : 'Gamma_Right',
    'dofs' : {'T.0' : -2.0},
}

integral_1 = {
    'name' : 'i1',
    'kind' : 'v',
    'quadrature' : 'gauss_o1_d3',
}

equations = {
    'Temperature' :
    """  dw_laplace.i1.Omega( coef.val, s, T )
       + dw_mass_scalar.i1.Omega( s, dT/dt )
       = 0"""
}

solver_0 = {
    'name' : 'ls',
    'kind' : 'ls.umfpack',
}

solver_1 = {
    'name' : 'newton',
    'kind' : 'nls.newton',

    'i_max'      : 1,
    'eps_a'      : 1e-10,
    'eps_r'      : 1.0,
    'macheps'   : 1e-16,
    'lin_red'    : 1e-2, # Linear system error < (eps_a * lin_red).
    'ls_red'     : 0.1,
    'ls_red_warp' : 0.001,
    'ls_on'      : 1.1,
    'ls_min'     : 1e-5,
    'check'     : 0,
    'delta'     : 1e-6,
    'is_plot'    : False,
    'problem'   : 'linear', # 'nonlinear' or 'linear' (ignore i_max)
}

solver_2 = {
    'name' : 'ts',
    'kind' : 'ts.simple',

    't0'    : t0,
    't1'    : t1,
    'dt'    : None,
    'n_step' : n_step, # has precedence over dt!
}

options = {
    'nls' : 'newton',
    'ls' : 'ls',
    'ts' : 'ts',
    'save_steps' : -1,
}

fe = {
    'chunk_size' : 1000
}

@@ -56,119 +56,37 @@
 h'%0A%0A
-from sfepy.solvers.ts import TimeStepper%0At0 = 0.0%0At1 = 0.1%0An_step = 11%0Ats = TimeStepper( t0, t1, None, n_step )
+t0 = 0.0%0At1 = 0.1%0An_step = 11
 %0A%0Ama
@@ -650,13 +650,9 @@
 gion
-_1000
+s
  = %7B
@@ -661,105 +661,50 @@
    '
-name' : 'Omega',%0A    'select
+Omega
 ' : 
+(
 'all',
-%0A%7D%0Aregion_03 = %7B%0A    'name' : 'Gamma_Left',%0A    'selec
+ %7B%7D),%0A    'Gamma_Lef
 t' : 
+(
 'nod
@@ -728,67 +728,35 @@
 1)',
-%0A%7D%0Aregion_4 = %7B%0A    'name' : 'Gamma_Right',%0A    'selec
+ %7B%7D),%0A    'Gamma_Righ
 t' : 
+(
 'nod
@@ -781,17 +781,21 @@
 9)',
+ %7B%7D),
 %0A%7D%0A%0Aebc
-_1
+s
  = %7B
@@ -804,37 +804,14 @@
    '
-name' : 'T1',%0A    'region' 
+T1'
 : 
+(
 'Gam
@@ -819,29 +819,16 @@
 a_Left',
-%0A    'dofs' :
  %7B'T.0' 
@@ -837,75 +837,22 @@
 2.0%7D
-,%0A    'value' : 2.0,%0A%7D%0Aebc_2 = %7B%0A    'name' : 'T2',%0A    'region' 
+),%0A    'T2'
 : 
+(
 'Gam
@@ -865,21 +865,8 @@
 ht',
-%0A    'dofs' :
  %7B'T
@@ -876,16 +876,17 @@
  : -2.0%7D
+)
 ,%0A%7D%0A%0Aint

shortened input